24小时热门版块排行榜

>论坛更新日志 (4617)
>虫友互识 (537)
>文献求助 (439)
>休闲灌水 (401)
>导师招生 (399)
>考博 (136)
>硕博家园 (133)
>博后之家 (104)
>基金申请 (101)
>招聘信息布告栏 (100)
>论文投稿 (95)
>考研 (88)
>教师之家 (65)
>找工作 (50)
>公派出国 (48)
>绿色求助(高悬赏) (46)

返回列表

后知后觉的夜

新虫 (初入文坛)

应助: 0 (幼儿园)
帖子: 1
在线:
虫号: 1997005
注册: 2012-09-13

[求助] 求助，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=[1 1]'; %激励转移矩阵
H=[1 0;0 1]; %输出转移矩阵
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为单步预测误差协方差矩阵即P（k+1|k）
   K=P1*H'*inv(H*P1*H'+R);
   P0=(eye(2,2)-K*H)*P1;%P0为误差协方差矩阵即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为状态矢量估计值
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为状态矢量真实值
  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)方差')


什么原因啊，本人是菜鸟，希望大侠多多帮忙，只有两个金币，全部贡献了。。

回复此楼