I write this code for calculating Kalman Gain and show them for 50 iterations on a plot.
I want to get a plot like here in the attachment. I need to show (0,0) and show all the points.
Please, help me 🙂
clearclcP = 1;R = 0.1;XK1 =0;hold onZ=random('Normal',0,2,1,50);randn('seed',0)uzunluk=length(Z);XK_sakla=[1:uzunluk]; KG_sakla=[1:1:uzunluk];a=0;for y = [1:uzunluk] a=a+1; KG = P/(P+ R);XK = XK1 + (KG*(Z(a)-XK1));P = (1-KG)*P;XK1 = XK;XK_sakla(a) = XK;KG_sakla(a) = KG; endplot (1:uzunluk , KG_sakla,'blue','Linewidth',1 );grid on;grid minor;xlabel('Number of Observations');ylabel('Kalman Gain');
Best Answer