文档介绍：基于通信故障的多传感器系统融合估计算法研究报告
基于通信故障的多传感器系统融合估计算法研究报告
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
差值
均值

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=;
q=;
h=;
r=;
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); %观测方程
end
end
end
%————————————Kalman滤波估计————————————————————
for k=1:l
for i=1:N
if k==1
x_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