| 3 | 1/1 | 返回列表 |
| 查看: 352 | 回復(fù): 2 | |||
[求助]
求助,matlab運行錯誤
|
|
錯誤信息如下: Error using my_kalman (line 15) Not enough input arguments. 代碼如下: function [X1 X2 K]=my_kalman(Z) x1=0;x2=0; %初始條件 F=[1 1;0 1]; %狀態(tài)轉(zhuǎn)移矩陣 T=[1 1]'; %激勵轉(zhuǎn)移矩陣 H=[1 0;0 1]; %輸出轉(zhuǎn)移矩陣 P0=[10 0;0 10]; X=zeros(2,11); X(1,1)=x1; X(2,1)=x2; R=[1 0;0 2]; for k=1:10 P1=F*P0*F'+T*T'; %P1為單步預(yù)測誤差協(xié)方差矩陣即P(k+1|k) K=P1*H'*inv(H*P1*H'+R); P0=(eye(2,2)-K*H)*P1;%P0為誤差協(xié)方差矩陣即P(k+1|k+1) X(:,k+1)=F*X(:,k)+K*(Z(:,k)-H*F*X(:,k)); fangcha=zeros(2,10); fangcha(1,k)=P0(1,1); fangcha(2,k)=P0(2,2); end K X1=X(1, ![]() X2=X(2, %X1,X2為狀態(tài)矢量估計值X11=[0,1.65428714,3.50300702,5.997852924,9.15040740,12.50873910,... 16.92192594,21.34483352,25.89335144,31.54135330,36.93605670]; X22=[0,1.65428714,1.84871988,2.47552222,3.17187816,3.35833170,... 4.41318684,4.42290758,4.54851792,5.64800186,5.39447034]; %X11,X22為狀態(tài)矢量真實值 wucha1=X11-X1; wucha2=X22-X2; %wucha1,wucha2分別為x1(k),x2(k)誤差值 k=0:10; subplot(3,1,1),plot(k,X1,'r-',k,X11,'g-',k,X2,'b-',k,X22,'y-') title('x1(k),x2(k)估計值和真實值'),xlabel('k') legend('x1(k)估計值','x1(k)真實值','x2(k)估計值','x2(k)真實值') subplot(3,1,2),plot(k,wucha1,'r-',k,wucha2,'b-') title('x1(k),x2(k)誤差'),xlabel('k') legend('x1(k)誤差','x2(k)誤差') k=1:10; subplot(3,1,3),plot(k,fangcha(1, ,'r-',k,fangcha(2, ,'b-')title('x1(k),x2(k)方差'),xlabel('k') legend('x1(k)方差','x2(k)方差') 什么原因啊,本人是菜鳥,希望大俠多多幫忙,只有兩個金幣,全部貢獻了。。 |
榮譽版主 (著名寫手)
彩色掛圖
![]() |
專家經(jīng)驗: +2 |

|
本帖內(nèi)容被屏蔽 |
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 352求調(diào)劑 +4 | 大米飯! 2026-03-22 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 292求調(diào)劑 +9 | 鵝鵝鵝額額額額?/a> 2026-03-25 | 10/500 |
|
|
[考研] 0856調(diào)劑 +4 | 求求讓我有書讀?/a> 2026-03-26 | 5/250 |
|
|
[考研] 0856求調(diào)劑 +8 | zhn03 2026-03-25 | 9/450 |
|
|
[考研] 303求調(diào)劑 +6 | 藍(lán)山月 2026-03-25 | 6/300 |
|
|
[考研] 材料求調(diào)劑 +4 | .m.. 2026-03-25 | 4/200 |
|
|
[考研] 機械學(xué)碩總分317求調(diào)劑。。! +4 | Acaciad 2026-03-25 | 4/200 |
|
|
[考研] 290分調(diào)劑求助 +3 | 吉祥止止陳 2026-03-25 | 3/150 |
|
|
[考研] 333求調(diào)劑 +6 | 87639 2026-03-21 | 11/550 |
|
|
[考研] 347求調(diào)劑 +4 | L when 2026-03-25 | 4/200 |
|
|
[考研] 求調(diào)劑323材料與化工 +4 | 1124361 2026-03-24 | 4/200 |
|
|
[考研] 考研化學(xué)308分求調(diào)劑 +10 | 你好明天你好 2026-03-23 | 11/550 |
|
|
[考研] 一志愿北化315 求調(diào)劑 +3 | akrrain 2026-03-24 | 3/150 |
|
|
[考研] 305分求調(diào)劑(食品工程) +5 | Sxy112 2026-03-21 | 7/350 |
|
|
[考研] 材料與化工085600,總分304,本科有兩篇sci參與,求調(diào)劑 +4 | 幸運的醬醬 2026-03-22 | 5/250 |
|
|
[考研] 一志愿華中農(nóng)業(yè)071010,總分320求調(diào)劑 +5 | 困困困困坤坤 2026-03-20 | 6/300 |
|
|
[考研]
|
Grand777 2026-03-21 | 3/150 |
|
|
[考研] 求調(diào)劑 +3 | .m.. 2026-03-21 | 4/200 |
|
|
[考研] 0805材料320求調(diào)劑 +3 | 深海物語 2026-03-20 | 3/150 |
|
|
[考研] 求調(diào)劑 +3 | @taotao 2026-03-20 | 3/150 |
|