| 查看: 412 | 回复: 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 ] |
回复此楼» 猜你喜欢
面上的水很深
已经有25人回复
-
国社科通讯评审
已经有7人回复
-
5个%2F，非固定段22个字母
已经有13人回复
-
小木虫没落了，除了祈祷帖子，几乎看不到有价值的帖子
已经有16人回复
-
生命口C13面上会评时间
已经有11人回复
-
生命口面上会评结束了吗大家有知道结果吗？现在还没有结果正常不？
已经有9人回复
-
投英文期刊如何查重
已经有6人回复
-
伯胺上甲基求助
已经有7人回复
-
生生命学部会评
已经有4人回复
-
H口，面上，时间戳7月7日。
已经有3人回复
sunyuanxin
禁虫 (著名写手)
|
本帖内容被屏蔽 |
2楼2009-03-10 11:28:32