基于通信故障的多传感器系统融合估计算法研究报告.doc

基于通信故障的多传感器系统融合估计算法研究报告.doc

ID:55296319

大小:197.00 KB

页数:13页

时间:2020-05-09

基于通信故障的多传感器系统融合估计算法研究报告.doc_第1页
基于通信故障的多传感器系统融合估计算法研究报告.doc_第2页
基于通信故障的多传感器系统融合估计算法研究报告.doc_第3页
基于通信故障的多传感器系统融合估计算法研究报告.doc_第4页
基于通信故障的多传感器系统融合估计算法研究报告.doc_第5页
资源描述:

《基于通信故障的多传感器系统融合估计算法研究报告.doc》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库

1、基于通信故障的多传感器系统融合估计算法研究报告基于通信故障的多传感器系统融合估计算法研究报告1.原理简介在实际的多传感器系统进行数据融合时,传输线路由于各种各样的原因总会一些故障,而导致某些传感器的数据无法及时传送到中信处理器,为了检测哪些传感器传输出现故障,首先给出经过改进后的无反馈层次结构的融合模型。该模型就是在正常的最优联合融合算法的基础上增加一个标志寄存器Bi,该操作根据硬件电路的设计来实现。寄存器Bi初始状态为低电平,当中心处理器接收到传感器i的信息时,立刻给相应的寄存器一个高电平信号,则相应的寄存器Bi被置为

2、1,否则Bi保持0。利用Kalman滤波的方法对数据进行预测估计,得到结果信息融合。利用公式(1)(2)2.仿真结果从图1中可以看出,故障情况下最优融合估计比正常情况下的融合估计较不稳定。在某些时间点估计偏差较大。可以看出当传输出现故障时,IOCFA算法的精确度低于正常情况。但他具有实时检测以及故障处理的功能。图1NOCFA与IOCFA算法融合估计效果图根据附件1传输故障统计表可以看出在第55、66次采样时传输故障为0个由图可以看出,此时图形IOCFA与NOCFA的值相等,且误差为0,如图2、3所示。图2局部放大图图3N

3、OCFA与IOCFA算法融合估计误差图图4IOCFA算法跟踪到(部分)出现通信故障的传感器示意图求取两种算法的曲线均值,得到两种的误差,如表1所示。表1两种算法的据对误差均值算法NOCFAIOCFA差值均值2.90461.97440.93033.存在疑问老师给的论文中,给出了IOCFA算法跟踪到出现故障的传感器示意图。这个图在程序中怎么表示?附录1仿真程序代码%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%通信故障时

4、Kalman最优滤波融合算法%%公式步骤:%%(1)在k+1时刻产生1*N的随机向量(u=randn(1,N));%%(2)对每个随机变量中的元素进行判断,如果u(i)>=E(E是一个任意的实数,在本节仿%%真中取E=0),则将相应的标志寄存器Bi置1;否则为0;%%(3)如果Bi全为0,则转到(1)重新产生随机向量;否则继续;%%(4)继续融合算法%%要求:%%能够比较滤波后估计与预测值的差别,能够比较出它的优越性。%%作者:时间:2013年8月7日%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

5、%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%clearcloseallclcx0=5;x_gj0=5;p0=10;fai=0.98;q=0.5;h=0.98;r=0.2;N=5;l=100;E=0;s=0;s1=0;%————————————产生仿真信号数据————————————————————fork=1:la=0;while(1)u=randn(1,N);fori=1:Nifu(i)

6、tinueelsebreakendendfori=1:Nw(k,i)=sqrt(q)*u(i);v(k,i)=sqrt(r)*u(i);ifk==1x(k,i)=fai*x0+w(k,i);%%表示第i个传感器的k时刻的状态z(k,i)=h*x0+v(k,i);elsex(k,i)=fai*x(k-1,i)+w(k,i);%状态方程z(k,i)=h*x(k,i)+v(k,i);%观测方程endendend%————————————Kalman滤波估计————————————————————fork=1:lfori=1:N

7、ifk==1x_yc(k,i)=fai*x_gj0;p_yc(k,i)=fai*p0*fai'+q;K(k,i)=p_yc(k,i)*h'*inv(h*p_yc(k,i)*h'+r);x_gj(k,i)=x_yc(k,i)+K(k,i)*(z(k,i)-h*x_yc(k,i));p_gj(k,i)=(1-K(k,i)*h)*p_yc(k,i);elsex_yc(k,i)=fai*x_gj(k-1,i);p_yc(k,i)=fai*p_gj(k-1,i)*fai'+q;K(k,i)=p_yc(k,i)*h'*inv(h*p

8、_yc(k,i)*h'+r);x_gj(k,i)=x_yc(k,i)+K(k,i)*(z(k,i)-h*x_yc(k,i));p_gj(k,i)=(eye(1)-K(k,i)*h)*p_yc(k,i);endendend%————————————————信息融合————————————————————fork=1:l

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。