※ [本文轉錄自 MATLAB 看板 #1R1J6OY3 ]
作者: Ecampus () 看板: MATLAB
標題: [請問] 請問我這個迴圈的CODE，寫得對嗎 @@??
時間: Wed May 23 17:17:08 2018

這是一個迴圈矩陣，從 0分鐘 迭代到 160分鐘
I＝eye(360,360);
[Xt]是360*1的矩陣、
[A]是360*360的矩陣、
u＝0
[Ｐt+1]跟[Ｐt]都是360*360的矩陣、
[Kg]是360*1的矩陣、
[H]是1*360的矩陣、
yt+1是純量，是實測值的data、
下面是CODE：
K=zeros(360,1);
for f=2:N % f，代表從第2次迭代～到第N次迭代
x(:,1,f) = A*x(:,1,f-1);
P(:,:,f) = A*P(:,:,f-1)*A';
K(:,1,f) = P(:,:,f)*H' / (H*P(:,:,f)*H');
x(:,1,f)=x(:,1,f)+K(:,1,f)*(y(f)-H*x(:,1,f-1));
x0(1,f)=x(230,1,f);
% x0的物理意義是把卡爾曼濾波算法得到的每一k時刻的位移值放在x0中,
於是x0就是每一時刻的濃度值
P(:,:,f)=(I-K(:,1,f)*H)*P(:,:,f);;
end
y
x0
t=1:N;
plot(t,y,t,x0);
我把參數輸入自己的CODE之後，結果跟想像的有落差，
所以好奇想發問一下，是否這個迴圈矩陣寫錯了??