多传感器数据融合算法

1、、背景介绍:多传感器数据融合是一种信号处理、辨识方法，可以与神经网络、小波变换、kalman滤波技术结合进一步得到研究需要的更纯净的有用信号。如信号处理、估计理论、不确定性理论、 多传感器数据融合比较确切的定义可概括 采用计算机技术对按时间序列获得的多多传感器数据融合涉及到多方面的理论和技术, 最优化理论、模式识别、神经网络和人工智能等。 为：充分利用不同时间与空间的多传感器数据资源，传感器观测数据，在一定准则下进行分析、综合、支配和使用，获得对被测对象的一致性解释与描述，进而实现相应的决策和估计，使系统获得比它的各组成部分更充分的信息。多传感器信息融合技术通过对多个传感器获得的信息进行协调、

2、组合、互补来克服单个传感器的不确定和局限性，并提高系统的有效性能，进而得出比单一传感器测量值更为精确 的结果。数据融合就是将来自多个传感器或多源的信息在一定准则下加以自动分析、综合以完成所需的决策和估计任务而进行的信息处理过程。当系统中单个传感器不能提供足够的准确度和可靠性时就采用多传感器数据融合。数据融合技术扩展了时空覆盖范围，改善了系统的可靠性，对目标或事件的确认增加了可信度，减少了信息的模糊性， 这是任何单个传感器做不到的。实践证明：与单传感器系统相比， 运用多传感器数据融合技术在解决探测、跟踪和目标识别等问题方面，能够增强系统生存能力，提高整个系统的可靠性和鲁棒性，增强数据的可信度，并

3、提高精度，扩展整个系统的时间、 空间覆盖率，增加系统的实时性和信息利用率等。 信号级融合方法最简单、最直观方法是加权平均法，该方法将一组传感器提供的冗余信息进 行加权平均，结果作为融合值，该方法是一种直接对数据源进行操作的方法。卡尔曼滤波主要用于融合低层次实时动态多传感器冗余数据。该方法用测量模型的统计特性递推，决定统计意义下的最优融合和数据估计。多传感器数据融合虽然未形成完整的理论体系和有效的融合算法，但在不少应用领域根据各自的具体应用背景， 已经提出了许多成熟并且有效的融合方法。多传感器数据融合的常用方法基本上可概括为随机和人工智能两大类，随机类方法有加权平均法、卡尔曼滤波法、 多贝叶斯估

4、计法、产生式规则等 ；而人工智能类则有模糊逻辑理论、神经网络、粗集理论、 专家系统等。可以预见，神经网络和人工智能等新概念、新技术在多传感器数据融合中将起 到越来越重要的作用。数据融合存在的问题(1) 尚未建立统一的融合理论和有效广义融合模型及算法；(2) 对数据融合的具体方法的研究尚处于初步阶段；(3) 还没有很好解决融合系统中的容错性或鲁棒性问题；(4) 关联的二义性是数据融合中的主要障碍；(5) 数据融合系统的设计还存在许多实际问题。二、算法介绍：2.1多传感器数据自适应加权融合估计算法：设有n个传感器对某一对象进行测量，如图1所示，对于不同的传感器都有各自不同的加权因子，我们的思想是在

5、总均方误差最小这一最优条件下，根据各个传感器所得到的测量值以自适应的方式寻找各个传感器所对应的最优加权因子，使融合后的X值达到最优。图I多件感器数据白适应加权融合估计模电最优加权因子及所对应的均方误差：(多传感器方法的理论依据：设n个传感器的方差分别为,022,，C2n；所要估计的真值为X，各传感器的测量值分别为 X1, X2,，Xn,它们彼此互相独立，并且是X的无偏估计；各传感器的加权因子分别为W, W2,，W,则融合后的X值和加权因子满足以下两式：nnX 八 WpXp, Wp =1p Ap 42 n 2 2 n总均方误差为 b2 =E Wp2(X _Xp ) +2 送 WpWq(X _Xp

6、 X X _Xq )pTqV_因为X 1 , X 2 ,，X n彼此独立，并且为 X的无偏估计，所以 E (X-Xp)(X-Xq) =0(p 为;P =1 , 2 ,，n;q =1 , 2 ,，n)，故 0 可写成-n 2 n/ = E 匡 W；(X -Xp ) “ W芯-p#从式可以看出，总均方误差0是关于各加权因子的多元二次函数，因此0必然存在最小值。该最小值的求取是加权因子W1 , W2 ,，Wn满足式约束条件的多元函数极值求取。根据多元函数求极值理论，可求出总均方误差最小时所对应的加权因子：n 1、W； =1/ 确瓦(Psw=1,2,川,n)2 n 1此时对应的最小均方误差为：二min

7、 =1/2p仝口 p以上是根据各个传感器在某一时刻的测量值而进行的估计，当估计真值X为常量时，则1 kXpg=k2Xp(i)可根据各个传感器历史数据的均值来进行估计。设np =1,2,HI,n此时估计值为 刃八 WpXp kpm总均方误差为-E十2W；(XXp(k)2严一n2、WpWqX-Xp k X-Xq kph,q 巧p刊同理，因为X1 , X2 ,，X n为X的无偏估计，所以X 1(k) , X 2(k),X n(k)也一定2 一 2H 1 n 2 2是X的无偏估计，故T2 =E巨 Wp2(X Xp(k)卜一 W冷2 用一 k pa自适应加权融合估计算法的线性无偏最小方差性1) 线性估计

8、由式可以看出，融合后的估计是各传感器测量值或测量值样本均值的线性函数。2) 无偏估计因为Xp(p =1 , 2,，n)为X的无偏估计，即 EX-Xp =0(p =1 , 2，n)，所以可得-nnE X 灭=E 区 Wp(X Xp )=送 WpE X Xp| =0 , X 为无偏估计。.PA_ P=1同理，由于Xp(p =1 , 2，n)为X的无偏估计，所以Xp(k)也一定是X的无偏估计。nnE 石 Wp(X Xp (k = Z WpE X Xp(k )=0 j 二-p4最小均方误差估计在推导过程中，是以均方误差最小做为最优条件，因而该估计算法的均方误差一定是最的。为了进一步说明这一点，我们用所

9、得的均方误差*Lmin与用单个传感器均值做估计和用多传感器均值平均做估计的均方误差相比较。我们用n个传感器中方差最小的传感器L做均值估计，设传感器L的方差o2Lmin为测量数据2(J.所以甘- min的个数为k,则千乂 h/k, 盂=1/n 11 VmQ - 1P八pp=L下面我们讨论与用多个传感器均值平均做估计均方误差相比较的情况。所谓用多个传感器均值平均做估计是用n个传感器测量数据的样本平均再做均值处理而得1 n到的估计，即 乂二丄a Xp k此时均方误差为1门2kJ；则n pz1cP2 = E2(X-刃)=集E2(X-Xp(k)+ $ EXn pmnP=1,qz1p=q-Xp k E X

10、 - Xq k1 n:?2二 2minN n=丄Q2 丄 p ln P=1同理，Xp(k)一定为X的无偏估计，可得 炉=2送E(xXp(k) n p=11我们事先已经将各个传感器的方差进行排序，且不妨设 八n pcp222cP21 n 2 10 ：匚1汀2 -川-；n，则根据契比雪夫不等式得2P飞勺amin ng ap各传感器方差cp2的求取从以上分析可以看出，最佳加权因子Wp*决定各个传感器的方差c2。一般不是已知的，我们可根据各个传感器所提供的测量值，依据相应的算法，将它们求出。设有任意两个不同的传感器p、q，其测量值分别为 X p、X q，所对应观测误差分别为Vp、V q，即Xp二XVp

11、； Xq二XVq，其中，Vp、V q为零均值平稳噪声，则传感器P方差二E V；，因为Vp、V q互不相关，与X也不相关，所以X p、Xq的互协方差函数 Rpq满足 Rpq 二 E X pXq二E X2，X p的自互协方差函数 Rpp满足-Rpp - RpqRpp =EXpXp = E X2|+E Vp2 作差得 p=E Vk, Rpp的时对于Rpp、Rpq的求取，可由其时间域估计值得出。设传感器测量数据的个数为 间域估计值为Rpp(k)，Rpq的时间域估计值为Rpq(k)，则12k-11Rppk 編；XpiXpi -pRppk1 匚XpkXpkk -11Rpq =Rpq k-1?Xp k Xq

12、 kkk如用传感器q(q Mp ；q =1 , 2,，n)与传感器p做相关运算，则可以得到 Rpq(k)(q和；q =1 , 2,，n)值。因而对于 R pq可进一步用R pq(k)的均值Rp(k)来做为它的估计，即1 n Rpq rRp k = Rpq kn 1 q 吕q=p由此，我们依靠各个传感器的测量值求出了Rpp与Rpq的时间域的估计值，从而可估计出各个传感器的方差。基于最小二乘原理， 提出了一 并且根据测量信息， 提出了 从而对各传感器的权值进行合理2.2基于最小二乘原理的多传感器加权融合算法以存在随机扰动环境中的不同参数多传感器为研究对象， 种加权融合算法，推导出各传感器的权系数与

13、测量方差的关系。 一种方差估计学习算法，实现对各传感器测量方差的估计， 的分配。该算法简单，能快速、准确的估计出待测物理量的状态信息。同种类型不同参数的多个传感器对存在随机扰动环境中的某一状态进行测量时，如何使状态的估计值在统计意义上更加接近于状态的真实值，针对这一问题进行了研究。依据最小二乘原理，推导出了多传感器的加权融合公式，并且在最优原则下， 得出测量过程中各传感器的测量方差与其权系数的关系。针对以上不足，充分利用多传感器测量这一特点，将传感器内部噪声与环境干扰综合考虑，提出了一种对各传感器测量方差及待测物理量状态进行实时估计的算法。设n个传感器对某系统状态参数的观测方程为：Y = Hx

14、 e，式中，x为一维状态量；Y为n维测量向量，设 Y - “ 丫2丨1|丫丁，e为n维测量噪声向量，包含传感器的内部噪声及环境干扰噪声，设e =1,H为已知n维常向量。采用加权最小二乘法从测量向量Y中估计出状态量x的估计量。加权最小二乘法估计的准则是使加权误差平方和Jw (刃)=(丫 一 HXWWH)?)取最小值。其中 W是一个正定对角加权阵，设W =diag w1 wJHwn，对之求偏导，令：?=-Ht W WT 丫 _ Hx = 0得二Z Wiyi到加权最小二乘估计：9= H TWH 一 H TWY二毘一Z Wii 4对测量噪声作如下假设：（1）各传感器的测量噪声为相互独立的白噪声；（2）

15、由于测量噪声是传感器内部噪声和环境干扰等多种相互独立因素引起的，利用概率知识可以证明：多个相互独立的随机变量相加的和接近正态分布。因而可以假设测量噪声的分布规律也是正态的。所以 Ee2 = E （xyi ）2= R（i =1,2,川,n）写作矩阵形式：E le I - 0, E eeT yRdiag RR2|Rn其中，Ri为第i个传感器的测量方差，R为测量方差矩阵。可得估计方差:E(x)2=E= eZi4n Wiyi i吕nWi 4Win工Wi由于i不等于j时ei、ejr 7nE i4x yi相互独立，故、2Wn wii Wx yiWi x-yiwii 4WjWj72I - Wi(x yi )

16、(X yj )winwii壬R令偏导数为零得W i =1,2,H),n得估计方差为E x-x2不难看出，采用加权融合的估计方差比任何一个传感器的测量方差都小。当以算术平均作为1 n 1 1 n状态的估计时，其估计方差七 R，可以证明R说明加权融合的效果要 n i m丁1 n i mi =1 R优于算术平均估计。可得 x=HTWH HTWe, Exi- H TWH 1 H TWE l.e - 0可知基于最小二乘原理的加权融合算法是一种无偏估计算法。通过以上的推导，公式）即为基于最小二乘原理的加权融合算法的计算公式。测量方差阵R的计算方法：进行测量方差的估计时， 把传感器的内部噪声与环境干扰综合考

17、虑，将得出一个随不确定因素而变化的测量方差阵 R的估计方法。在对测量方差进行估计之前，先作如下分析：（1）横向 分析（针对多个传感器一次采样结果的分析）:多个传感器单次采样结果的算术平均值是该采 样时刻状态的无偏估计。基于这个原理，各传感器测量方差的估计可先基于算术平均值作一 个粗略的分配估算；以每个传感器的测量值与该次采样时各传感器测量算术平均值的偏差平 方作为各传感器该次采样的方差分配。横向分析中利用了多传感器在某一采样时刻的测量信息。（2）纵向分析（针对一个传感器多次采样结果的分析）:以单个传感器为研究对象，测量方差是传感器内部噪声与环境干扰的一种综合属性，这一属性始终存在于测量的全过程

18、中，因此要将单个传感器历次采样时的方差分配与当前方差分配的算术平均值作为当前 测量方差的实时估算。亦即在此提出了方差估计学习算法。基于以上分析，方差估计学习算法如下：设ymi表示第i个传感器第m次采样的结果，则第m次采样时各传感器测量算术平1 n均值为：ym =丄V ymi。第i个传感器第m次采样时测量方差的估计分配Rmi为：n i 4Rmi二lymyj2对各传感器测量方差在历次采样时的估计分配值代求算术平均值1 mFU二丄Rji此式即为第m次采样时第i个传感器测量方差的估计值，写成递推公式形m y式为：RmiR(m J.)iRoi=01-Rmi m=1,2|l|m将结果代入，便得测量过程中各

19、传感器的权系数。由测量方差估计的计算过程可以看出，每次新的测量数据都对各传感器的测量方差有调节作用，但这种调节作用将越来越小。这是因为把传感器与测量环境综合起来考虑，测量向量从统计意义上说， 它的概率分布是确定的。方差估计学习算法实际上是随着采样时刻的推移，对测量向量分布特性的学习过程，而在学习过程中，最初的几个采样时刻是对测量向量分布特性从无到有的认识，因而学习速度较快，体现在对测量方差的估计中是相邻采样点间各传感器测量方差估计值的变化率较大。而随着采样的进行，这种学习过程将趋于稳定，体现在对测量方差的估计中是每次新的测量数据对 各传感器测量方差的估计只起微小的调节作用，相邻采样点间各传感器

20、测量方差估计值的变化率较小。2.3同类多传感器自适应加权估计的数据级融合算法研究针对同类多传感器测量中含有的噪声，提出了多传感器数据自适应加权融合估计算法，该算法不要求知道传感器测量数据的任何先验知识，依据估计的各传感器的方差的变化，及时调整参与融合的各传感器的权系数，使融合系统的均方误差始终最小，并在理论上证明了该估计算法的线性无偏最小方差性仿真结果表明了本算法的有效性，其融合结果在精度、容错性方面均优于传统的平均值估计算法。同类多传感器数据的测量可以看作是从含有噪声的大量测量数据中估计一个非随机量， 由于测量数据中存在着噪声，那么根据这些测量数据所得到的估计值也存在估计误差，然而这种估计误

21、差是随机量，一般用均方误差来评价测量方法的优劣，而影响估计值均方误差的主要因素是传感器自身的均方误差。在单一传感器测量时，为了减少估计值的均方误差就必须增大测量数据的数量，这必然降低实时性。为了提高测量的实时性和精度，就需要用同种类的多个传感器同时测量一个物理量。数据一致性检验设有m个传感器对某一对象进行测量，首先对 X i（i =1 ，2，m）进行数据检验，检验准则是Xi , X 2 ,，Xm的相邻两值之差不应超过给定门限& 根据传感器精度确定。即X2-:,X3 -X2兰名，Xm _Xmd兰名自适应加权融合估计算法理论：与2.2完全相同算法流程：1）根据递推式算出采样时刻k的Rppk与Rpq

22、k ; 2）计算k时刻Rpk ;3）计算21k 11k时刻二2 ； 4）求出各传感器k时刻均值Xp kXp iXpk-1 Xpk ;k ykk5）求出此时刻各传感器最优加权因子W； ; 6）得出此时刻估价式 X。从以上运算流程可以看出，对于每个传感器所对应的最优加权因子，只是根据各个传感器的测量数据以自适应的 方式将它们求取出来，因而，称该算法为多传感器数据自适应加权融合估计算法。2.4基于信任度的多传感器数据融合及其应用针对多传感器信息采集系统中的数据不确定性问题，提出了一种基于信任度的多传感器数据融合方法。该方法首先定义一个模糊型指数信任度函数，对两传感器测得数据 间的信任程度进行量化处理

23、，并通过信任度矩阵度量各传感器测得数据的综合信任程 度，以合理地分配测得数据在融合过程中所占权重，得到数据融合估计的最终表达式， 从而实现了多传感器数据的融合。在多传感器信息采集系统中，由于不可避免会受到传感器精度、传输误差、环境噪声和人为干扰等因素的影响，将使得它们的测得数据产生不确定性因此在数据融合过程中，必须首先确定被融合数据的可信程度：若某些数据表现异常， 就不能作为被融合的数据；若某些数据相互接近，则可以把它们融合在一起，从而提高融合结果的精确度和稳定性。针对上述问题，本文充分利用模糊集合理论中隶属度函数范围确定的优点，定义了一种模糊型指数信任度函数，对传感器测得数据间的信任程度进行

24、量化处理，并通过信任度矩阵度量各传感器测得数据的综合信任程度，合理地分配测得数据在融合过程中所占权重，得出数据融合估计的最终表达式，进而得到一种对多个传感器测得数据进行融合处理的简便有效的方法。设多个传感器测量同一参数，第 i个传感器和第j个传感器测得的数据分别为 xi和xj。 如果xi的真实性越高，xi被其余数据所信任的程度就越高。所谓 xi被xj信任程度，即从xj 看xi为真实数据的可能程度，多传感器测得数据间的这种信任程度被称为信任度。为了对测得数据间的信任度进行进一步地统一量化处理，定义一个信任度函数bij，表示xi被xj信任程度。q = f氐xj ） i, j =1,2,，n其中，0

25、兰f兰1，为连续下降函数。般给出融合上限 口耳=0，令勺=Xj Xj Emu若bj=0，认为第i个传感器与第jXi Xj mij个相互不信任，若 bj=1，则认为二者间信任。若一个传感器不被其他传感器信任，或只被少数传感器信任，则该传感器的读数在进行数据融合时即被删掉。这样处理不利于对 实际情况做出客观判别，进而使融合结果受主观因素的影响过大。erx 亠 x _ x I 0所以信任度矩阵 B是一个非负矩阵，并且该对称矩阵存在最大模特征 值 40,使得KA=BA。求出 入及对应特征向量 A，满足aiQ则W =A可以作为对可以作为各传感器测得数据间综合信任程度的度量，即Wi色Wj ajWiaia1

26、 aani, j =1,2, ,n对Wi进行归一化处理，得到n a/得到对所有传感器测得数据融合估计的最终结果为空印 + a? +an2.5提高测量可靠性的多传感器数据融合有偏估计方法为了提高测量数据可靠性，多传感器数据融合在过程控制领域得到了广泛应用。本文基于有偏估计能够减小最小二乘无偏估计方差的思想，提出采用多传感器有偏估计数据融合改善测量数据可靠性的方法。首先，基于岭估计提出了有偏测量过程，并给出了测量数据可靠性定量表示方法，同时证明了有偏测量可靠度优于无偏测量可靠度。其次，提出了多传感器有偏估计数据融合方法，证明了现有集中式与分布式无偏估计数据融合之间的等价性。最后，证明了多传感器有偏

27、估计数据融合收敛于无偏估计数据融合。证明了方法的有效性。目前单传感器测量数据的处理方法主要有三种：平均值法1、加权平均法2和递推滤波算法3.通过理论推导，发现这些方法都是特殊形式的最小二乘估计(Least squareestimation， LS)。基于模型y = Hx w，x的最小二乘估计为 ？_s = H TH H T ym瓦yiy为观测矢量，H为观测矩阵，未知矢量 x，当H=l，可化简为5?ls =m可知，平均值法与其具有相同的表达形式.采用类似的分析过程，可得另外两种方法与最小二乘估计是等价的.由于最小二乘估计是一种无偏估计，所以这种等价关系也说明上述三种数据处理方法具有无偏性，本文称

28、之为无偏测量过程。无偏测量过程可以采用方差直接衡量测量可靠性，即方差越小测量可靠性越高。为了提高测量可靠性，国内外学者提出了多传感器数据融合的方法，旨在减小测量方差。目前多传感器数据融合常用的理论方法为线性无偏估计理论(简称多传感器无偏估计数据融合)，其中又以最小二乘估计应用最为 广泛.但是现有多传感器无偏估计数据融合方法存在两方面问题：1)融合结果可靠性均为定性说明而无法量化表示，即只能通过比较不同融合结果的方差定性地判断融合结果可靠性的优劣；2)虽然多传感器无偏估计数据融合具有无偏性的优良性质，但是并不能由此认为它的测量结果一定是高可靠的.因为根据高斯-马尔科夫定理可知，最小二乘估计方差有

29、下界，所以此时无偏估计数据融合具有最小的方差，但是当这个最小方差本身却很大时，那么无偏估计数据融合将不能保证测量数据的可靠性一定是可接受的.但值得一提的是，无偏测量过程与最小二乘估计之间的等价关系为线性有偏估计算法用于提高测量可靠性成为可能.如James-Stein估计、压缩最小二乘估计、岭估计(Ridge estimation , RE)18-19 等.其中岭估计是应用最为广泛的改进最小二乘估计方法.本文以岭估计为基础提出多传感器有偏估计数据融合方法，岭估计长期以来一直是广泛用于改善最小二乘估计方差的有偏估计方法.由于无偏测量与最小二乘估计之间是等价的，所以本文借鉴岭估计的思想通过引入较小的

30、偏差改善无偏测量数据的方差，并称之为有偏测量过程.在此基础上解决有偏测量与无偏测量的可靠性定量表示问题.这种方法引入的偏差是可知的固定性偏差，且可以在一定 程度上减小估计值的方差，其余并没有创新，不详细介绍了。2.6基于小波去噪和数据融合的多传感器数据重建算法为了从被噪声干扰的各个传感器测量值中获得更准确的测量结果，提出了一种基于小波去噪和多传感器数据融合的传感器数据重建算法。仿真和实验的结果都表明：由该算法重建得到的各个传感器的重建数据的方差低于传感器测量值的方差。可以认为多传感器数据重建算法给出了对每一个传感器的更为准确的测量结果。一个传感器组，利用每一个传感器的测量值对其加权，进而对这组

31、传感器的测量结果进行数据融合以达到提高测量精度的目的具体方法是在方差基本定义的基础上提出递归的估计方差的算法，利用估计的方差估计出每个数据的权值，进而对电磁流量计的流量进行递归估计，从而达到提高精度的目的。为了从受到不同噪声干扰的各个传感器测量值中获得更准确的各个传感器数据，本文提出了一种基于小波去噪和多传感器数据融合的传感器数据重建算法该方法首先将每个传感器的测量值用小波阈值的方法去噪，减小噪声对传感器测量值的影响为了更好的重建传感器信号，先将各个传感器测量值进行归一化处理，再将归一化后的各个传感器测量值做基于最小均方的数据融合多传感器数据融合目的在于用较大的数据量，充分利用对被测目标的在时

32、间与空间的信息，获得对被测量的描述。来自多传感器的信号所提供的信息具有相关性、互补性和冗余性，将同源数据进行组合，可得到统计上的优势。基于小波去噪及多传感器数据融合的传感器数据重建算法：假设N个传感器在不同位置对同一测量值 Y测量，每个传感器测量值记为Xj(j=1，2, .N )由于测量中，存在内部外部噪声影响，测量值表示为Xj n =Sn 0 n j =1,2,N。S(n)为被测量，ej(n)为第j个传感器在时刻加性噪声，Xj(n)为第j个传感器在n时刻观测值。信号小波消噪方法主要通过设置阈值通过信号的离散小波变换，计算所有小波系数 ，然后剔除被认为跟噪声有关的小波系数。例如通常的方法是设置

33、阈值，将小于阈值的小波系数去掉最后，然后通过小波变换的逆变换来得到信号数据融合数据重建算法： 首先对每一个传感器获得的一组测量值用这组数据中的最大测量值归一：Yj n 二Y n /MaxYj n j =1,2, ,N;n =1,2 ,I其中，MaxYj n是在估计长度I内第j个传感器的最大测量值。Yj n为第j个传感器在n时刻归一化后的测量值， 由于每个传感器收到噪声干扰程度不同，所以偏离真实被测量程度不同，对每个传感器根据一定原则确定权值，可从N个传感器得到估计值 丫。NNY 7 WjYj, W 胡j #由于各传感器之间受到噪声干扰的程度不同，所以各传感器测量值的方差并不一致，即各传感器测量

34、值的可信度是不同的 若将较大的权值赋予可信度高的传感器，将较小的权值赋予可信度小的传感器 ，就可以使估计值更精确地描述原信号。1 1 1Wjj =1,2/ ,N，归一化权值为 WjN i =1,2/ ,N6 j 1对Y反归一化，得到各传感器重建数据：YjnAYn MaxYj n j -1,2/ , N;n =1,2, ,1算法步骤：1置估计长度I; 2对各传感器测量值作小波阈值去噪处理；3采用MA模型用递归方法估计方差；4计算Wj ; 5计算Y，计算各传感器重建数据。2.7测量噪声相关情况下的多传感器数据融合对于测量噪声相关的多传感器测量模型，利用Cholesky分解和单位下三角阵的求逆方法，

35、将其转化为测量噪声互不相关的等价的多传感器伪测量模型，然后基于M arkov估计，提出了一种测量噪声相关情况下多传感器数据融合的新方法。与直接利用原始传感器测量值的Markov估计数据融合方法相比，两者的计算精度相同，但新方法的计算复杂度却大大降低。数值仿真实验进一步验证了新方法的有效性。所谓多传感器数据融合，就是将来自多个同类或异类传感器的数据（信息）进行综合处理，以获得比单一传感器更为准确可靠的结果。已有的多传感器数据融合方法，一般利用含有加性噪声的线性测量方程来估计未知常值参数，大多假设各传感器的测量噪声之间互不相关。但是在实际应用中，由于各传感器通常处于同一测量环境，所以传感器的测量结

36、果中除由于传感器自身精度限制而引入的测量误差外，共同的环境噪声的影响也不容忽略，而这往往会导致各传感器的测量噪声之间相关，所以对测量噪声相关情况下多传感器测量系统的数据融合问题进行研究就具有更加广泛的应用价值。为了解决测量噪声相关情况下的多传感器测量数据融合问题，文献在最小二乘准则下， 利用Lagrange乘子条件极值方法，给出了一种最佳的线性数据融合方法，但是仅适用于被测参数为标量的情况，无法直接扩展到参数为矢量的情况，另外，由于需要对累积观测矢量的自相关阵直接求逆，所以计算复杂度非常大；文献则利用实对称矩阵的正交相似变换实现了多传感器测量噪声互协方差阵的对角化，从而实现了各传感器测量噪声之

37、间的去相关， 但是一般来说，这种对角化不能在有限步中完成，只能通过迭代步骤求近似值，所以该方法在实际应用时比较困难。本文首先利用 Cholesky分解和单位下三角阵的求逆方法将多传 感器的测量模型转化成各传感器的测量噪声互不相关的等价的伪测量模型，然后基于Markov估计提出了一种测量噪声相关情况下的多传感器数据融合的新方法。与直接利用原 始传感器测量值的 Markov估计数据融合方法相比，两者的结果相同，但新方法的计算复杂度大大降低。数值仿真实验进一步验证了本文方法的有效性。Zi = HiX vi测量噪声vi服从均值为0，方差为Rii的高斯分布zN 】，h 二 HIh： hN 1采用N个传感

38、器对同一常值参数进行线性测量模型一般表示成假定各传感器的测量噪声相关，即-RiiDRi2RiNR = cov v,v 】=R219DR22aR2N-_RniDRN2Rnn -的非对角块不全为0,且R正定。噪声相关的多传感器数据融合对于上式所描述的多传感器线性测量系统，由于Markov估计9为被测常值参数的无偏估计，且对应的均方误差阵最小，所以根据Markov估计可得这N个传感器对于常值参数的最优融合估计为XF = HTRH iHTRz，相应的融合均方误差PF = HTRH匸传感器噪声去相关：由于R = rij是一正定的实对称阵，根据矩阵的 Cholesky分解可知，R可以唯一地分解成 R= L

39、DLt。其中，L =叩 为单位下三角阵 ，D =diag（ d 1 ，i J2k=L，dNm)且正定，di dJik i =1,2, Nmj JLj ,2, ,i -Li 二 iij _ 送 dJikljk k)djlj = * L0 j =i +L,i +2,，Nm对于单位下三角阵 L ,其逆阵存在，且仍为单位下三角阵，记M =L，= Rj】令 j=L , 2, ., Nm，则mij-liji A八 likmkjkdi “,2,j -Li = ji = j Lj 2,j3, Nm将M分块阵表示M LLIm2L0L2M 220L N02NMii为单位下三角阵；0ij为零矩阵。在原格式两边左乘M NLM N2M NN则 y 二 Gx ， y 二 Mz 二 j： y； yNG 二 MH 二 GTg； GN T，= Mv = T N