3.5.22 · D1 · Physics › Guidance, Navigation & Control (GNC) › Kalman gain — minimizes trace of covariance
Tumhare paas do blurry guesses hain — ek model se, ek sensor se — aur tumhe unhe mila ke ek sharper guess banana hai. Kalman gain bas ek recipe hai ki sensor ko model ke comparison mein kitna weight do, aisa choose karo ki final guess jitni ho sake utni un-blurry ho.
Is page par kuch bhi assume nahi kiya gaya. Parent note Kalman gain padhne se pehle, tumhe uske har symbol ko padhna aana chahiye: x , x ^ , e , P , E [ ⋅ ] , tr ( ⋅ ) , H , z , v , R , K , aur chhote superscripts − , + , ⊤ . Hum inhe ek-ek karke build karenge, har ek ko ek picture se anchor karke.
x
State x unhi numbers ki list hai jo is waqt jo cheez tum track kar rahe ho usse poori tarah describe karti hai. Ek chalte drone ke liye yeh position aur velocity ho sakti hai: x = [ p v ] .
Picture: ek number line par (ya plane mein agar do numbers hain) true location par baitha ek akela dot.
Kyun zaroori hai: baaki sab kuch is true dot se kitna door hai, usi se define hota hai. Lekin yahan ek pakad hai — tum x ko actually kabhi dekh nahi sakte. Yeh chhipi hui sachhai hai.
x ^ (padho "x-hat")
Hat ^ ka matlab hamesha hota hai "humare hisaab se best guess." Toh x ^ humara guessed state hai — hum sochte hain x yeh hai.
Picture: ek doosra dot, true dot ke paas lekin uske upar nahi.
Kyun zaroori hai: poora filter ek machine hai jo acche hats banata hai. Dot aur uske hat ke beech ka gap woh dushman hai jisse hum lad rahe hain.
Figure s01 (neeche): ek horizontal number line jisme ek green dot dayi taraf true state x mark kar raha hai aur ek blue dot uske baayi taraf humara estimate x ^ mark kar raha hai; ek red arrow blue dot se green dot tak jaata hai, labelled e = x − x ^ — woh error jo hum abhi define karne wale hain.
Hat = guess. Koi bhi letter jo chhota hat pehne hue hai woh "our estimate of" hai, kabhi truth nahi.
Intuition "Average mein" kyun?
Kisi ek run mein, e kismet se bada ya chhota ho sakta hai. Jo cheez hum pin down kar sakte hain woh hai kai runs mein uska typical size . "Typical" ki baat karne ke liye humein expectation symbol chahiye.
Worked example Chhoti si expectation
Agar ek error + 2 half the time aur − 2 half the time ho, toh E [ e ] = 2 1 ( + 2 ) + 2 1 ( − 2 ) = 0 — average error zero. Lekin error clearly hamesha zero nahi hai! Isliye hum e 2 dekhte hain, e nahi: E [ e 2 ] = 2 1 ( 4 ) + 2 1 ( 4 ) = 4 . Squaring plus aur minus ko cancel hone se rokta hai.
Abhi tak hamare paas sirf estimation error e hai; sensor aur uska noise baad mein §7 mein aayega. Toh filhaal zero-mean idea sirf e ke liye bataya ja raha hai (sensor noise define hone ke baad hum wahan bhi isko repeat karenge).
Definition Unbiased (zero-mean) error
Ek accha filter unbiased hota hai: average mein uski guess exactly truth par padti hai, isliye
E [ e ] = 0.
Error dots truth ke aaspaas scatter hote hain bina kisi consistent lean ke.
Picture: error cloud ka centre-of-mass exactly true dot par baitha hai — na baayi taraf, na dayi taraf.
Kyun zaroori hai: iske bina, "spread" aur "covariance" (agle §4 mein define) galat centre ke aaspaas measure ki jaati.
E [ e 2 ] sirf tab jab mean zero ho
Variance ki general definition pehle mean subtract karti hai: Var ( e ) = E [ ( e − E [ e ] ) 2 ] . Sirf jab E [ e ] = 0 ho tabhi yeh woh clean E [ e 2 ] mein collapse hoti hai jo humne upar use ki. Jab hum §4 mein matrix covariance build karenge, yahi zero-mean assumption hi hai jo ise true covariance (mean ke aaspaas spread) banati hai, na ki raw second moment. Isliye is assumption ko apni jeb mein rakh lo — parent note ka har formula quietly isi par depend karta hai.
Definition Variance (ek number)
Upar wale zero-mean assumption ke saath, variance = E [ e 2 ] = average squared error = error cloud ki spread . Badi variance = fuzzy guess; chhoti variance = sharp guess.
Jab state mein kai components hain (position aur velocity), toh ek component per variance kaafi nahi hoti — alag-alag components ki errors linked ho sakti hain. Hum inhe sab ek matrix mein collect karte hain. (⊤ "transpose" symbol jo abhi neeche use hoga, woh poori tarah §6 mein define hai; filhaal e e ⊤ ko "errors ka column times wahi errors row mein laid out" padho, jo ek square table banata hai.)
P hamesha symmetric aur positive semidefinite hoti hai
P = E [ e e ⊤ ] se do properties seedhi nikalti hain aur yahi "ellipse" picture ko valid banati hain:
Symmetric (P = P ⊤ ): kyunki e e ⊤ har e ke liye symmetric hai (entry i , j equals entry j , i , dono e i e j hain), isliye uska average bhi symmetric hai. Ek symmetric matrix ke perpendicular axes hote hain — ellipse ke axes well-defined hote hain.
Positive semidefinite (P ⪰ 0 ): kisi bhi direction u ke liye, u ki direction mein error ki variance u ⊤ P u = E [( u ⊤ e ) 2 ] ≥ 0 hai — ek squared cheez negative nahi ho sakti. Iska matlab ellipse ki "negative width" kabhi nahi hogi, jo nonsense hota.
Ek matrix jo symmetric aur positive semidefinite hai woh exactly wohi hai jo ek sensible ellipse draw kar sakti hai — isliye P ko is tarah picture kiya ja sakta hai.
Figure s02 (neeche): ek wide blue ellipse aur ek small green ellipse, dono ek white dot x ^ par centred, "state component 1" aur "state component 2" axes ke upar drawn. Wide ellipse labelled "big blur = big P (unsure)"; tight wali "small blur (sure)" — dikhata hai ki ellipse ka size hi uncertainty P ki matra hai.
Intuition Chhoti ellipse = better
Poori Kalman story yeh hai: badi ellipse se shuru karo, ek measurement blend karo, chhoti ellipse pe khatam karo. Gain K aise choose kiya jaata hai ki final ellipse jitni ho sake utni chhoti ho. Geometry ke liye dekho Covariance Matrices and Uncertainty .
Trace se pehle, notation ka ek chhota piece jo parent note mein kaam aata hai:
Definition Euclidean norm
∥ ⋅ ∥
Ek vector w = w 1 w 2 ⋮ ke liye, Euclidean norm uski ordinary straight-line length hai:
∥ w ∥ = w 1 2 + w 2 2 + ⋯ , ∥ w ∥ 2 = w ⊤ w = w 1 2 + w 2 2 + ⋯
Picture: Pythagoras — origin se w ki tip tak arrow ki length.
Kyun zaroori hai: ∥ x − x ^ ∥ 2 hai guess aur truth ke beech squared distance , humara ek "kitna galat" number. Squared component-errors ko sum karna hi ek pure error vector ko ek honest scalar mein badalta hai.
Common mistake "Determinant kyun nahi use karte?"
Kyun sahi lagta hai: determinant ellipse ka area hai, ek aur size measure.
Fix: trace = total mean-squared error — ek physically meaningful "expected wrongness." Woh clean matrix-calculus jo K deta hai, woh bhi sirf trace ke liye kaam karta hai. Isliye trace dono meaningful aur tractable hai.
A ⊤
A ⊤ matrix ko uski diagonal par flip karta hai — rows columns ban jaati hain. Ek column vector e ek row vector e ⊤ ban jaata hai.
Picture: ledger ko side pe tilt karo.
Kyun zaroori hai: e e ⊤ (column times row) poora P table build karta hai; e ⊤ e (row times column) single number ∥ e ∥ 2 mein collapse ho jaata hai. Order matter karta hai, aur transpose hi humein choose karne deta hai. (Yahi woh symbol hai jo humne §4 mein explain karne ka promise kiya tha.)
Definition Measurement matrix
H
H state-space ko sensor-space mein translate karta hai. State x mein position aur velocity dono ho sakti hain, lekin ek ruler sirf position padh sakta hai. H woh parts pick out / mix karta hai jo sensor actually dekhta hai.
Picture: state world se measurement world ki taraf point karta ek arrow labelled "convert."
Kyun zaroori hai: sensor reading ko x se directly compare nahi kar sakte — alag worlds hain. H x ^ hai "agar humari guess sahi ho toh sensor kya padhega, humari expectation."
Figure s03 (neeche): bayi taraf ek blue box labelled "state-space" jisme x = ( p , v ) (position aur velocity) hai, ek yellow arrow labelled "H (convert)" dayi taraf pointing karta hai, aur dayi taraf ek green box labelled "measurement-space" jisme "H x = position only" hai — dikhata hai ki H state components ko drop/mix karke waisi cheez mein le jaata hai jo sensor actually padhta hai.
z
z = H x + v
Woh number jo sensor actually nikalta hai: true state H se dekha, plus ek random smudge v .
Picture: true reading H x , phir upar se ek random wobble add.
v aur uski covariance R
v woh random error hai jo har reading ke andar hota hai. Bilkul §3 ke estimation error ki tarah, hum assume karte hain ki yeh zero-mean hai, E [ v ] = 0 (sensor true reading ke aaspaas wobble karta hai bina kisi built-in offset ke). Uska typical size tab true covariance R = E [ v v ⊤ ] hai — sensor ki khud ki uncertainty ellipse (symmetric aur positive semidefinite, usi argument se jo §4 mein P ko woh properties deta hai).
Picture: true reading ke aaspaas ek fuzzy halo; wide halo = noisy sensor (bada R ), tight halo = precise sensor (chhota R ).
Kyun zaroori hai: sensor par kitna trust karna hai woh completely uski noisiness par depend karta hai. R wahi noisiness hai, aur woh gain formula ke andar hi baitha hai. Dekho Innovation and Residuals ki z − H x ^ kaisa "surprise" ban ke filter react karta hai.
− aur posterior +
x ^ − , P − (minus / "prior"): naya measurement dekhne se pehle humari guess aur uska blur — seedha model ki prediction se. Dekho Kalman Filter — Predict Step .
x ^ + , P + (plus / "posterior"): measurement fold in karne ke baad humari guess aur uska blur. Dekho Kalman Filter — Update Step .
Picture: ek wide "before" ellipse aur ek narrower "after" ellipse. Measurement minus ko plus mein squeeze karti hai.
Minus = measure karne se pehle, Plus = measure karne ke baad.
Intuition Kya measuring hamesha ellipse shrink karti hai?
Linear-Gaussian assumptions ke under jo pura topic use karta hai (linear H , zero-mean uncorrelated Gaussian noise, aur optimal gain), plus-ellipse kabhi minus-ellipse se badi nahi hoti — nayi information sirf help kar sakti hai. Yeh un assumptions ka consequence hai, koi universal law nahi: nonlinear model ke saath (dekho Extended Kalman Filter (EKF) in GNC ), bura chosen gain, ya non-Gaussian noise ke saath, posterior actually worse bhi aa sakta hai. "Measuring ellipse shrink karti hai" wali picture rakh lo, lekin yaad raho ki yeh assumptions se milti hai, free nahi hai.
K
K woh weight hai jo hum sensor ki surprise par rakhte hain jab apni guess correct karte hain:
x ^ + = x ^ − + K ( z − H x ^ − ) .
Picture: 0 se 1 tak ek slider. K = 0 par hum sensor ignore karte hain aur prediction rakhte hain; K = 1 ke paas hum almost poori tarah measurement ki taraf jump karte hain.
Kyun zaroori hai: yahi ek cheez hai jo hum choose kar sakte hain. Baaki sab (P − , H , R ) given hai; K free dial hai. Parent note ka poora kaam woh K dhundna hai jo tr ( P + ) ko sabse zyada shrink kare.
Intuition Woh tug-of-war jo tum abhi milne wale ho
K badhao → sensor par zyada rely karo (accha agar model bura hai, bura agar sensor noisy hai). K ghataao → model par zyada rely karo. Kahin beech mein ek sweet spot hai — ek U-shaped valley ka bottom. Wahi valley hai jo parent note ka derivative hunt karta hai. Yeh do sources ka blend karna exactly Sensor Fusion in Navigation hai, aur mathematically yeh Least Squares Estimation ki ek recursive form hai.
Zero-mean assumption - E of e and v are 0
Covariance P - the blur ellipse
Expectation - average over luck
Euclidean norm - straight-line length
Trace - total blur as one number
Transpose - flip rows and columns
Noise v inside each reading
Gain K - the trust slider
Kalman gain - minimizes trace of covariance
Test karo apne aap ko — tum parent note ke liye ready ho sirf tab jab bina dekhe har sawaal ka jawab de sako.
x ^ mein hat, matlab kya hota hai?"Humare hisaab se best estimate / guess of" — kabhi true value nahi.
Estimation error likho aur batao yeh kahan point karta hai. e = x − x ^ ; humari guess se wapas truth ki taraf gap/arrow.
E [ ⋅ ] kya compute karta hai aur humein isko kyun chahiye?Infinitely many random repetitions mein average; noise random hota hai isliye hum strategies ko average behaviour se compare karte hain.
Error ko average karne se pehle square kyun karte hain? Taaki positive aur negative errors cancel na hon — E [ e ] 0 ho sakta hai jabki guess phir bhi galat ho; E [ e 2 ] true spread capture karta hai.
Filter kaun si zero-mean assumptions karta hai, aur kyun? E [ e ] = 0 aur E [ v ] = 0 (unbiased); sirf tab P = E [ e e ⊤ ] aur R = E [ v v ⊤ ] true covariances hain na ki raw second moments.
∥ w ∥ ka matlab kya hai?Euclidean (straight-line) length
w 1 2 + w 2 2 + ⋯ ; aur
∥ w ∥ 2 = w ⊤ w .
P words mein kya hai, aur iske saath kaun si picture jaati hai?Error covariance P = E [ e e ⊤ ] ; guess ke aaspaas ek uncertainty ellipse — uska size humara blur hai.
P mein hamesha kaun si do properties hoti hain, aur yeh kyun matter karti hain?Symmetric (P = P ⊤ ) aur positive semidefinite (u ⊤ P u ≥ 0 ); yahi exactly woh hain jo P ko ek sensible uncertainty ellipse draw karne deti hain.
Dikhao kyun tr ( P ) = E [ ∥ x − x ^ ∥ 2 ] . tr ( P ) = tr ( E [ e e ⊤ ]) = E [ tr ( e e ⊤ )] = E [ e ⊤ e ] = E [ ∥ e ∥ 2 ] .
Matrix H kya karta hai? State-space ko measurement-space mein map karta hai; H x ^ hai "agar humari guess sahi ho toh sensor kya padhega."
v aur R kya hain?v har reading ke andar random sensor noise hai; R = E [ v v ⊤ ] uska size hai (sensor ki khud ki uncertainty).
Superscripts − aur + mein farq batao. − = prior, measurement use karne se pehle ; + = posterior, measurement use karne ke baad .
Kaun si assumptions mein posterior ellipse prior se kabhi badi nahi hoti? Linear-Gaussian setting mein optimal gain ke saath; nonlinear models, bure gains, ya non-Gaussian noise ke saath fail ho sakta hai.
Ek sentence mein, Kalman gain K kya hai? Woh 0-to-1 slider setting jo decide karta hai ki sensor ki surprise ka kitna part humari guess mein fold karein, choose kiya jaata hai after-blur minimize karne ke liye.