| 查看: 374 | 回复: 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)方差') 什么原因啊，本人是菜鸟，希望大侠多多帮忙，只有两个金币，全部贡献了。。 |
回复此楼» 猜你喜欢
如何从铁电相到顺电相。
已经有1人回复
-
各位大佬，求一份最新Wien2k版本作为学习用途
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有295人回复
-
南科大活性流体和软物质课题组诚招2027级博士、硕士研究生和博后
已经有42人回复
-
瑞典林雪平大学博士后招聘｜PEC Water Splitting 方向
已经有24人回复
-
PRB投稿，Acknowledgment sent to author状态十天了
已经有3人回复
-
PVA溶解
已经有2人回复
-
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助 matlab交点问题
已经有19人回复
- 求助matlab问题
已经有5人回复
- 求助有关Matlab有约束非线性最优化问题
已经有9人回复
- 求助matlab-simulink运行时报错问题
已经有3人回复
- 【求助】matlab 中，几个m文件调用
已经有14人回复
- 求助matlab的安装问题
已经有3人回复
- 【求助】大家看看，matlab安装过程的截图，找出错
已经有14人回复
- 【求助】matlab 积分的问题
已经有5人回复
- 【求助】matlab中保存计算数据的问题
已经有8人回复
- 【求助】修改MATLAB支持向量机核函数程序遇到的疑问？求助
已经有5人回复
- 【求助】使用Matlab预估动力学方程问题
已经有13人回复
2楼2012-12-21 14:27:43
感谢参与，应助指数 +1
|
本帖内容被屏蔽 |
3楼2012-12-22 23:57:31
100