24小时热门版块排行榜    

查看: 376  |  回复: 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=[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)方差')
  
  
   什么原因啊,本人是菜鸟,希望大侠多多帮忙,只有两个金币,全部贡献了。。
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

11jxli1

禁虫 (著名写手)

感谢参与,应助指数 +1
本帖内容被屏蔽

3楼2012-12-22 23:57:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 3 个回答

csgt0

荣誉版主 (著名写手)

彩色挂图

【答案】应助回帖

感谢参与,应助指数 +1
Z你的给出啊,
在命令窗口中先载入Z,然后用
[X1 X2 K]=my_kalman(Z)运行
showmethemoney
2楼2012-12-21 14:27:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +7 布布和一二 2026-08-10 9/450 2026-08-11 22:12 by 医学老男孩
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 F4zcb7xRCEJA 2026-08-11 3/150 2026-08-11 21:51 by pMlf9MbdMHyc
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 yPv6ndDtiRxQ 2026-08-11 4/200 2026-08-11 21:15 by thgx1vqMqrLw
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +17 Lanmanbaby 2026-08-09 27/1350 2026-08-11 21:12 by klxdzh
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 yPv6ndDtiRxQ 2026-08-10 7/350 2026-08-11 20:39 by thgx1vqMqrLw
[公派出国] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 yPv6ndDtiRxQ 2026-08-10 8/400 2026-08-11 19:56 by thgx1vqMqrLw
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 GpvwGNVshCOp 2026-08-10 7/350 2026-08-11 19:39 by thgx1vqMqrLw
[硕博家园] 读博的好处 +3 lnee 2026-08-11 3/150 2026-08-11 18:10 by 希望我好好的
[基金申请] 什么时候出结果,有咨询渠道??? +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 分享一下我之前已中青C的计划书的filecode +3 布布和一二 2026-08-11 4/200 2026-08-11 15:20 by aasahr
[论文投稿] 职称评审,求友友推荐见刊最快的期刊 +3 工厂打螺丝 2026-08-08 3/150 2026-08-11 13:42 by 会议编辑
[基金申请] 听说今天filecode变了 +26 布布和一二 2026-08-06 49/2450 2026-08-11 13:18 by WH3796
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 关于代码变化问题,想知道的进来 +17 且听虎啸 2026-08-07 24/1200 2026-08-10 18:35 by zhangduo2008
[基金申请] 关于Filecode分析方法 +3 majunge000 2026-08-10 3/150 2026-08-10 15:46 by lch2012
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
信息提示
请填处理意见