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

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

  1. 1、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
  2. 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  3. 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
  4. 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
  5. 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们
  6. 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
  7. 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
基于通信故障多传感器系统融合估计算法研究报告

基于通信故障的多传感器系统融合估计算法研究报告 基于通信故障的多传感器系统融合估计算法研究报告 1. 原理简介 在实际的多传感器系统进行数据融合时,传输线路由于各种各样的原因总会一些故障,而导致某些传感器的数据无法及时传送到中信处理器,为了检测哪些传感器传输出现故障,首先给出经过改进后的无反馈层次结构的融合模型。 该模型就是在正常的最优联合融合算法的基础上增加一个标志寄存器Bi,该操作根据硬件电路的设计来实现。寄存器Bi初始状态为低电平,当中心处理器接收到传感器i的信息时,立刻给相应的寄存器一个高电平信号,则相应的寄存器Bi被置为1,否则Bi保持0。利用Kalman滤波的方法对数据进行预测估计,得到结果信息融合。利用公式 (1) (2) 2. 仿真结果 从图1中可以看出,故障情况下最优融合估计比正常情况下的融合估计较不稳定。在某些时间点估计偏差较大。可以看出当传输出现故障时,IOCFA算法的精确度低于正常情况。但他具有实时检测以及故障处理的功能。 图1 NOCFA与IOCFA算法融合估计效果图 根据附件1传输故障统计表可以看出在第55、66次采样时传输故障为0个由图可以看出,此时图形IOCFA与NOCFA的值相等,且误差为0,如图2、3所示。 图2 局部放大图 图3 NOCFA与IOCFA算法融合估计误差图 图4 IOCFA算法跟踪到(部分)出现通信故障的传感器示意图 求取两种算法的曲线均值,得到两种的误差,如表1所示。 表1 两种算法的据对误差均值 算法 NOCFA IOCFA 差值 均值 2.9046 1.9744 0.9303 3. 存在疑问 老师给的论文中,给出了IOCFA算法跟踪到出现故障的传感器示意图。这个图在程序中怎么表示? 附录1 仿真程序代码 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %% 通信故障时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日 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% clear close all clc x0=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; %————————————产生仿真信号数据———————————————————— for k=1:l a=0; while(1) u=randn(1,N); for i=1:N if u(i)E B(k,i)=0; a=a+1; else B(k,i)=1; end end if a==N a=0; continue else break end end for i=1:N w(k,i)=sqrt(q)*u(i); v(k,i)=sqrt(r)*u(i); if k==1 x(k,i)=fai*x0+w(k,i); %%表示第i个传感器的k时刻的状态 z(k,i)=h*x0+v(k,i); else x(k,i)=fai*x(k-1,i)+w(k,i); %状态方程 z(k,i)=h*x(k,i)+v(k,i);

文档评论(0)

phl805 + 关注
实名认证
文档贡献者

建筑从业资格证持证人

该用户很懒,什么也没介绍

领域认证 该用户于2023年05月12日上传了建筑从业资格证

1亿VIP精品文档

相关文档