Fw: [請問] 請問我這個迴圈的CODE，寫得對嗎 @@??

看板Programming作者Ecampus時間6年前 (2018/05/23 17:23)推噓0(0推 0噓 0→)

留言0則, 0人參與討論串1/1

※ [本文轉錄自 MATLAB 看板 #1R1J6OY3 ] 作者: Ecampus () 看板: MATLAB 標題: [請問] 請問我這個迴圈的CODE，寫得對嗎 @@?? 時間: Wed May 23 17:17:08 2018 https://imgur.com/a/xAtl4kf 這是一個迴圈矩陣，從 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之後，結果跟想像的有落差，所以好奇想發問一下，是否這個迴圈矩陣寫錯了?? -- ※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 125.224.101.239 ※ 文章網址: https://www.ptt.cc/bbs/MATLAB/M.1527067032.A.883.html ※ 發信站: 批踢踢實業坊(ptt.cc) ※ 轉錄者: Ecampus (125.224.101.239), 05/23/2018 17:23:42 ※ 編輯: Ecampus (125.224.101.239), 05/23/2018 17:27:15