| 查看: 427 | 回复: 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 ] |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有8人回复
-
面上意见出来了
已经有12人回复
-
为什么到现在没收到通知？
已经有4人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有10人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有6人回复
-
有没有仍没收到信息的
已经有5人回复
-
导师吐槽：我怎么摊上了这么个极品研究生！
已经有9人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有5人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有4人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有10人回复
sunyuanxin
禁虫 (著名写手)
|
本帖内容被屏蔽 |
2楼2009-03-10 11:28:32
投票:
30