| 查看: 387 | 回复: 1 | |||
| 当前主题已经存档。 | |||
[交流]
【求助】精通kalman滤波的大侠请进.
|
|||
|
我编写了关于kalman滤波的matlab程序,程序不长,可结果就是不对,请大侠帮忙看看程序有什么问题.多谢. 跟踪一个纯复数信号y(t)=exp(i*w*t); clear all; t=0:0.02/64:0.02*20; s=2*(cos(2*pi*50*t)+sqrt(-1)*sin(2*pi*50*t));%被分析信号 y=s; tx1(1)=1;tx2(1)=1;pk=[1 0;0 1];kk=[0.3 1]';%设初值 h=[0 1]; for k=2:1:100 px1(k)=tx1(k-1); px2(k)=tx1(k-1)*tx2(k-1); pkx=[px1(k) px2(k)]'; %求x(k|k-1) txk=pkx+kk*(y(k)-h*pkx); %求x(k|k) fk=[1 0;tx2(k-1) tx1(k-1)]; mk=fk*pk*(conj(fk))'; kk=mk*(conj(h))'*inv(h*mk*(conj(h))'+1);%求增益k tx1(k)=txk(1); tx2(k)=txk(2); pk=mk-kk*h*mk; end; [ Last edited by 2007骑猪逛街 on 2009-1-28 at 00:51 ] |
回复此楼» 猜你喜欢
323求调剂
已经有6人回复
-
一志愿北京化工大学 070300 学硕 336分 求调剂
已经有4人回复
-
352求调剂
已经有3人回复
-
一志愿东华大学化学070300，求调剂
已经有8人回复
-
277材料科学与工程080500求调剂
已经有7人回复
-
317求调剂
已经有18人回复
-
293求调剂
已经有5人回复
-
280分求调剂 一志愿085802
已经有7人回复
-
0854电子信息求调剂
已经有3人回复
-
263求调剂
已经有4人回复
sunyuanxin
禁虫 (著名写手)
|
本帖内容被屏蔽 |
2楼2009-03-10 11:28:32
