联邦滤波matlab程序,联邦滤波器仿真
不知道为什么,滤波器仿真效果很好,但是三个组合到一起效果就很差,是我联邦滤波算法的问题吗?请教各位大神,有做这个的麻烦帮我看看代码问题所在,或者有一个简单的联邦滤波程序能麻烦发给我一下吗?小弟不胜感激。
for j=2:step
%%
%预测
B=chol(3*P_,'lower');
Chi=B*Xi+repmat(x_,1,2*n);
Chi_=feval(ffun,Chi);
lamda=feval(hfun1,Chi);
x__=1/(2*n)*sum(Chi_,2);
y_=1/(2*n)*sum(lamda,2);
P__=1/(2*n)*Chi_*Chi_'-x__*(x__)'+3*Q;
Pyy=1/(2*n)*lamda*lamda'-y_*(y_)'+R1;%行向量
Pxy=1/(2*n)*Chi_*lamda'-x__*y_'+SC1;%Pxzi,i-1|i-2%矩阵按列3*3
M=Pxy*inv(Pyy);%矩阵按列3*3
x1=x__+M*(Y(1,j-1)-y_);
P1=P__-M*Pyy*(M)';
XX1(:,j)=x1;
PP1(:,j)=trace(P1);
B=chol(3*P_,'lower');
Chi=B*Xi+repmat(x_,1,2*n);
Chi_=feval(ffun,Chi);
lamda=feval(hfun2,Chi);
x__=1/(2*n)*sum(Chi_,2);
y_=1/(2*n)*sum(lamda,2);
P__=1/(2*n)*Chi_*Chi_'-x__*(x__)'+3*Q;
Pyy=1/(2*n)*lamda*lamda'-y_*(y_)'+R2;%行向量
Pxy=1/(2*n)*Chi_*lamda'-x__*y_'+SC2;%Pxzi,i-1|i-2%矩阵按列3*3
M=Pxy*inv(Pyy);%矩阵按列3*3
x2=x__+M*(Y(2,j-1)-y_);
P2=P__-M*Pyy*(M)';
XX2(:,j)=x2;
PP2(:,j)=trace(P2);
B=chol(3*P_,'lower');
Chi=B*Xi+repmat(x_,1,2*n);
Chi_=feval(ffun,Chi);
lamda=feval(hfun3,Chi);
x__=1/(2*n)*sum(Chi_,2);
y_=1/(2*n)*sum(lamda,2);
P__=1/(2*n)*Chi_*Chi_'-x__*(x__)'+3*Q;
Pyy=1/(2*n)*lamda*lamda'-y_*(y_)'+R3;
Pxy=1/(2*n)*Chi_*lamda'-x__*y_'+SC3;%Pxzi,i-1|i-2
M=Pxy*inv(Pyy);
x3=x__+M*(Y(3,j-1)-y_);
P3=P__-M*Pyy*(M)';
XX3(:,j)=x3;
PP3(:,j)=trace(P3);
P_=inv(inv(P1)+inv(P2)+inv(P3));%三个传感器融合
x_=P_*(inv(P1)*x1+inv(P2)*x2+inv(P3)*x3);
xx(:,j)=x_;
PP(:,j)=trace(P_);
end
联邦滤波matlab程序,联邦滤波器仿真相关推荐
- matlab 多 带阻,matlab程序之——滤波器(带通-带阻
<matlab程序之--滤波器(带通-带阻>由会员分享,可在线阅读,更多相关<matlab程序之--滤波器(带通-带阻(8页珍藏版)>请在人人文库网上搜索. 1.matlab程 ...
- MATLAB实现卡尔曼滤波器仿真
MATLAB实现卡尔曼滤波器仿真 本仿真是在学习B站DR_CAN的视频之后的练习.联系针对的是第五个视频中的二维示例. 视频连接为 [[卡尔曼滤波器]1_递归算法_Recursive Processi ...
- matlab程序 地震 相干噪声_SAR滤波MATLAB程序
SAR滤波MATLAB程序 所属分类:图形图像处理 开发工具:matlab 文件大小:1613KB 下载次数:28 上传日期:2018-03-20 17:21:52 上 传 者:pumpking 说明 ...
- 算术平均值滤波matlab程序,基于S7-1200 AD采样的高效数字滤波算法的设计与实践...
在工业控制中,常常会伴随大量的模拟量数据采样.模拟量采样时,由于待采量本身.传感器和传输过程中的外界干扰,特别是非稳态干扰信号等因素的影响,采样值中通常含有各种周期性和非周期性的噪声和干扰.为得到稳定 ...
- 做微波有源电路的设计仿真 一般都用什么软件 matlab可以吗,MATLAB在有源滤波器仿真设计中的应用...
MA TLAB 在有源滤波器仿真设计中的应用Ξ王金星 王庆平 贾长朱 杨 刚 陈超英 (天津大学 天津 300072) APPL ICATION OFMATLAB IN DESIGN AND SIM ...
- matlab实现中值滤波程序,中值滤波matlab程序
中值滤波3*3滤波窗口 程序如下: >> i=imread('Winter.jpg'); %读入图像 subplot(2,2,1); imshow(i); %显示原始图像 title('o ...
- 分享matlab程序之——滤波器篇(高通,低通)
快毕业了,把自己写的现成的matlab函数分享给有需要的人,由于个人水平有限,写的不好请见谅,愿意拍砖的尽管拍好了.目前还不考虑读博,所以写的程序仍了可惜,所以就拿出来分享.好了不废话了,开始正题. ...
- 指数高通滤波器代码_分享matlab程序之——滤波器篇(高通,低通)
快毕业了,把自己写的现成的matlab函数分享给有需要的人,由于个人水平有限,写的不好请见谅,愿意拍砖的尽管拍好了.目前还不考虑读博,所以写的程序仍了可惜,所以就拿出来分享.好了不废话了,开始正题. ...
- 谐波均值滤波matlab,谐波均值滤波器.PPT
谐波均值滤波器 数字图像处理Digital Image Processing 第6章 图像复原(Image Restoration) 第6章 图像复原(Image Restoration) 6.1 F ...
最新文章
- python编程例子 输入 输出-推荐 :手把手教你用Python创建简单的神经网络(附代码)...
- Date 和 SimpleDateFormat 类表示时间
- Xen server三 为xen添加存储及创建虚拟机
- UOJ#80 二分图最大权匹配 [模板题]
- OpenCV形态学变换函数morphologyEx()开运算的使用
- mac下安装配置mongodb
- SVT-AV1:开源编解码最新进展
- 为什么越来越多的开发者选择使用Spring Boot
- [XSY] 简单的博弈题(博弈+dp+组合数+容斥)
- [Flink] Flink运行报错The number of requested virtual cores for application master
- Linux性能分析工具汇总
- 源码nginx+php
- Linux基金会呼吁厂商放弃微软FAT文件系统
- FullScreen
- 一般家用路由器买多大的合适_家用路由器选多少M(兆)的合适
- linux 网卡绑定解绑,网卡绑定开机绑定,关机解绑
- 微信小程序 语音录制功能和文件(ppt、word、excel、pdf、txt格式)上传
- 微信服务号和订阅号的五大区别,如何选择申请微信公众号
- 使众人行:如何带领一群人把事做成?
- win7/win10 + VS2015 + cocos2dx-3.13的win32环境配置 与 项目编译时间优化和项目精简