24小时热门版块排行榜    

查看: 377  |  回复: 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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +4 Tide man 2026-08-11 4/200 2026-08-12 00:39 by liuqidi
[基金申请] 是这周出结果还是下周出结果? +3 yuleib84 2026-08-11 3/150 2026-08-11 21:53 by jnhyjjm
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 什么时候出结果,有咨询渠道??? +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[论文投稿] 职称评审,求友友推荐见刊最快的期刊 +3 工厂打螺丝 2026-08-08 3/150 2026-08-11 13:42 by 会议编辑
[基金申请] FileCode能看出啥? +7 要乐观耀哥 2026-08-10 17/850 2026-08-11 10:50 by 要乐观耀哥
[基金申请] filecode +9 documentary 2026-08-10 9/450 2026-08-11 10:08 by wk7465
[基金申请] 静等基金结果 +5 gjjjzhong 2026-08-10 16/800 2026-08-10 17:19 by Tide man
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 好奇怪的filecode +4 布布和一二 2026-08-08 5/250 2026-08-10 14:20 by 冰心玉壶晴
[基金申请] filecode与中标关系的预测 +5 布布和一二 2026-08-07 5/250 2026-08-09 16:15 by 袁向阳007
[基金申请] 关于豆爷回答的JTJC与%2F数量 +5 yang182083 2026-08-06 7/350 2026-08-08 18:29 by zhanghaozhu
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变,有变的大神么? +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] 固定端突然变了,今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] 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 医学老男孩
[基金申请] filecode +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
[基金申请] 影响面上的因素 +8 布布和一二 2026-08-05 11/550 2026-08-06 10:41 by 宝贝虫子
信息提示
请填处理意见