多传感器大数据融合算法

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

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

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

4、方法基本上可概括为随机和人工智能两大类，随机类方法有加权平均法、卡尔曼滤波法、 多贝叶斯估计法、产生式规则等；而人工智能类则有模糊逻辑理论、神经网络、粗集理论、专家系统等。可以预见，神经网络和人工智能等新概念、新技术在多传感器数据融合中将起 到越来越重要的作用。数据融合存在的问题(1) 尚未建立统一的融合理论和有效广义融合模型及算法；(2) 对数据融合的具体方法的研究尚处于初步阶段；(3) 还没有很好解决融合系统中的容错性或鲁棒性问题；(4) 关联的二义性是数据融合中的主要障碍；(5) 数据融合系统的设计还存在许多实际问题。二、算法介绍：2.1多传感器数据自适应加权融合估计算法：设有n个传感器

5、对某一对象进行测量，如图 1所示，对于不同的传感器都有各自不同的加 权因子，我们的思想是在总均方误差最小这一最优条件下，根据各个传感器所得到的测量值以自适应的方式寻找各个传感器所对应的最优加权因子，使融合后的X值达到最优。图】多传感器数据自适应加权融合估汁模型最优加权因子及所对应的均方误差：21 , b22，b2n ； 所要估计的真（多传感器方法的理论依据：设n个传感器的方差分别为b值为X，各传感器的测量值分别为X1,X2，Xn，它们彼此互相独立，并且是 X的无偏估计；各传感器的加权因子分别为W1,W 2，W n ,则融合后的X值和加权因子满足以nXWpXp,nWp1p 1总均方误差为nWp2

6、 Xp 1XpnWpWqx1,q 1XpX Xq因为X 1 , X 2-,X n彼此独立，并且为X的无偏估计,所以E (X-Xp)(X-Xq) =0,(p FP =12，n;q =1, 2，n），故b2可写成n少n2 2 2 2EWp X XpWp 2p 1p 1从式可以看出,总均方误差b 2是关于各加权因子的多元二次函数,因此b2必然存在最小值。该最小值的求取是加权因子 W1 , W2，Wn满足式约束条件的多元函数极值求取。根据多元函数求极值理论，可求出总均方误差最小时所对应的加权因子:Wp 1/2丄Psw 1,2,L , ni 1 i此时对应的最小均方误差为：min 1/p 1 pX为常量

7、时，则以上是根据各个传感器在某一时刻的测量值而进行的估计，当估计真值 可根据各个传感器历史数据的均值来进行估计。设1 k nXp k Xp ip 1,2,L ,n 此时估计值为 XWpXp kk i 1p 1总均方误差为-2 E X )? 2nEWp2 XXp kp 1n2WpWqX Xp k X Xq kp 1,q 1同理，因为X1 , X2，X n为X的无偏估计,所以X 1(k) , X 2(k)，,X n(k)也一定p qn2 2是X的无偏估计，故一 EWp X Xpp 1自适应加权融合估计算法的线性无偏最小方差性1）线性估计由式可以看出，融合后的估计是各传感器测量值或测量值样本均值的线

8、性函数。2）无偏估计因为Xp（p =1 , 2，n）为X的无偏估计，即 EX-Xp =0（p =1, 2，n），所以可得nEWpp 10，X为无偏估计。nX XpWpE X Xp 1同理,由于Xp(p =12，n）为X的无偏估计，所以Xp（k）也一定是X的无偏估计。WpX XpnWpE X Xp k最小均方误差估计 在推导过程中，是以均方误差最小做为最优条件，因而该估计算法的均方误差一定是最的。与用单个传感器均值做估计和用多传为了进一步说明这一点，我们用所得的均方误差b 2Lmin感器均值平均做估计的均方误差相比较。我们用n个传感器中方差最小的传感器L做均值估计,设传感器L的方差b 2Lmin

9、为测量数据的个数为k，则一2h/k 一初山2n 1所以今12 minI 1minp 1 pp LF面我们讨论与用多个传感器均值平均做估计均方误差相比较的情况。所谓用多个传感器均值平均做估计是用n个传感器测量数据的样本平均再做均值处理而得1 n到的估计，即 求Xp k此时均方误差为n p i?同理,E2Xp(k) 定为2Xp k的无偏估计，可得?2丄n n2 p iXp k E XXq kiq i p q2X Xp ki nn2k p i2min若我们事先已经将各个传感器的方差进行排序，且不妨设222? i n 2 i0 i 2 L n，则根据契比雪夫不等式得p 2 imin n p ip各传感

10、器方差b p2的求取从以上分析可以看出，最佳加权因子Wp*决定各个传感器的方差bP2。般不是已知的，我们可根据各个传感器所提供的测量值，依据相应的算法，将它们求出。设有任意两个不同的传感器p、q ,其测量值分别为 X p、X q ,所对应观测误差分别为V、V q，即Xp X Vp； Xq X Vq，其中，Vp、V q为零均值平稳噪声，则传感器p方差E Vp2，因为Vp、V q互不相关，与X也不相关，所以X p、X q的互协方差函数 Rpq2满足Rpq E XpXqE X , X p的自互协方差函数 Rpp满足E XpXpE X2 E Vp2 作差得 2 E VpRppRpq对于Rpp、Rpq的

11、求取，可由其时间域估计值得出。设传感器测量数据的个数为k, Rpp的时间域估计值为Rpp (k) , Rpq的时间域估计值为Rpq (k)，则Rpp ki Xp i Xp ik i ik iiRpp k i -Xp k Xp kkkjRpqk iiXp k Xq k如用传感器q(q工p ;q =i , 2，n)与传感器p做相关运算，则可以得到Rpq (k)(qzp ;q =i , 2，n)值。因而对于 R pq可进一步用R pq (k)的均值Rp(k)来做为它的估计，即实用标准文案T7q1Rpqq p由此，我们依靠各个传感器的测量值求出了Rpp与Rpq的时间域的估计值，从而可估计出各 个传感器

12、的方差。22基于最小二乘原理的多传感器加权融合算法以存在随机扰动环境中的不同参数多传感器为研究对象，基于最小二乘原理， 提出了一种加权融合算法，推导出各传感器的权系数与测量方差的关系。并且根据测量信息，提出了一种方差估计学习算法，实现对各传感器测量方差的估计，从而对各传感器的权值进行合理的分配。该算法简单，能快速、准确的估计出待测物理量的状态信息。同种类型不同参数的多个传感器对存在随机扰动环境中的某一状态进行测量时，如何使状态的估计值在统计意义上更加接近于状态的真实值，针对这一问题进行了研究。 依据最小二乘原理，推导出了多传感器的加权融合公式，并且在最优原则下， 得出测量过程中各传感器的测量方

13、差与其权系数的关系。针对以上不足，充分利用多传感器测量这一特点，将传感器内部噪声与环境干扰综合考虑，提出了一种对各传感器测量方差及待测物理量状态进行实时估计的算法。设n个传感器对某系统状态参数的观测方程为：Y Hx e，式中，x为一维状态量；Y为n维测量向量，设 丫 y2L yn 丁，e为n维测量噪声向量，包含传感器的内部噪声及环T境干扰噪声，设e e 62L en ，h为已知n维常向量。采用加权最小二乘法从测量向量Y中估计出状态量x的估计量。加权最小二乘法估计的准则是使加权误差平方和Jw ? Y H)?TW Y H)?取最小值。其中 W是一个正定对角加权阵，设Jw XttW diag w,

14、L wn，对之求偏导，令一HWWYHxO 得n到加权最小一乘估计：1W| yi? H TWHH TWYnWi 1 状态的估计时，其估计方差 ：Wii 1对测量噪声作如下假设：(1)各传感器的测量噪声为相互独立的白噪声；(2)由于测量噪声是传感器内部噪声和环境干扰等多种相互独立因素引起的，利用概率知识可以证明：多个相互独立的随机变量相加的和接近正态分布。因而可以假设测量噪声的分布规律也是正态的。2 2以E e2Ex yiR I 1,2,L ,n写作矩阵形式:E e 0, E ee R diag R R2LRn其中，Ri为第i个传感器的测量方差，R为测文档量方差矩阵。可得估计方差:2Ex)?nWi

15、yix 口nWii 1Wixyi由于i不等于E x ?j 时 ei、nWinWii 1ejyi相互独立，故WjWjnWii 12 xyix yjWnWii 1WinWii 1R令偏导数为零得1,2,L , n得估计方差为1n丄i 1 Ri不难看出，采用加权融合的估计方差比任何一个传感器的测量方差都小。当以算术平均作为2 R，可以证明n i 11Ri说明加权融合的效果要i 1 R|T1 tT1 T优于算术平均估计。可得 HtWH HTWe, E% HtWH H tWE e 0 可知基于最小二乘原理的加权融合算法是一种无偏估计算法。通过以上的推导，公式）即为基于最小二乘原理的加权融合算法的计算公式

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

17、象，测量方差是传感器内部噪声与环境干扰的一种综合属性，这一属性始终存在于测量的全过程中，因此要将单个传感器历次采样时的方差分配与当前方差分配的算术平均值作为当前测量方差的实时估算。 亦即在此提出了方差估计学习算法。基于以上分析，方差估计学习算 法如下：设ymi表示第i个传感器第m次采样的结果，则第m次采样时各传感器测量算术平均值为:ym1ymi。第i个传感器第m次采样时测量方差的估计分配Rmi为:n i i_ 2Rmi求算术平均值ymiym对各传感器测量方差在历次采样时的估计分配值Rmi1Rji此式即为第m次采样时第i个传感器测量方差的估计值，写成递推公式形m j i为:R0iSi m 1,2

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

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

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

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

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

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

24、xi Xji, j 1,2, , n其中，0 f 1，为连续下降函数。一般给出融合上限1XiXjmijjmj 0，令bj若b ij=0，认为第i个传感器与第j0xi xj mij个相互不信任，若bij=l ,则认为二者间信任。若一个传感器不被其他传感器信任，或只被少数传感器信任，则该传感器的读数在进行数据融合时即被删掉。这样处理不利于对实际情况做出客观判别，进而使融合结果受主观因素的影响过大。e lx xj x x m改进方法将bij设为指数函数，bji J ，即设定当二者差值大0人 Xj M于上限值M，二者不再信任bij=0。将bij定义成满足模糊性的指数函数形式这样既充分利用了模糊理论中隶

25、属度函数范围确定的优点，又避免了数据之间相互信任程度的绝对化，更加符合实际问题的真实性，同 时便于具体实施，可以使融合的结果更加精确和稳定。设有n个传感器测量同一参数，根据测得数据间的信任度函数bij，建立信任度矩阵Bbiib12binBb21b22b2nbn1bn2bnnn对于B中第i行元素来说，若bj较大，表明第i个传感器的测得数据被多数传感j i器信任；反之，第i个传感器的测得数据为真实数据的可能性较小。数据融合过程用Wi表示第i个传感器测得数据 Xi在融合过程中所占权重。 由于Wi值的大小反映了其它传感器测的数据对第i个传感器测得数据xi的综合信任程度，可以利用wi对xin进行加权求和

26、，得到数据融合的表达式?Wi Xii 1,2, ,ni 1n其中，权系数满足Wi10 Wi 1i 1在信任度矩阵B中，信任度函数bj仅仅表示测得数据Xj对Xi的信任程度，并不能反映系统中所有传感器的测得数据对Xi的信任程度，而Xi的真实程度实际上应该由bi1 , bi2，bin综合来体现。Wi应综合一个关于Xi的信任度系统中，各子系统bi1, bi2 , . , bin的全部信息，所以需要求出一组非负数a1,a2，an，使得Wi aa a2bi2an*i 1,2, ,n 改写为矩阵形式 W BA式中，W W1,W2, ,WnT , A aa2， a. T因为bij0，所以信任度矩阵B是一个非负

27、矩阵，并且该对称矩阵存在最大模特征值入0，使得 A BA。求出X及对应特征向量 A，满足ai0，则WWiaiWjajWiaiaa?可以作为对可以作为各传感器测得数据间综合信任程度的度量，即i, j 1,2, ,n对Wi进行归一化处理，得到naxi得到对所有传感器测得数据融合估计的最终结果为)?丄ai a2an2.5提高测量可靠性的多传感器数据融合有偏估计方法为了提高测量数据可靠性，多传感器数据融合在过程控制领域得到了广泛应用。本文基于有偏估计能够减小最小二乘无偏估计方差的思想，提出采用多传感器有偏估计数据融合改善测量数据可靠性的方法。首先，基于岭估计提出了有偏测量过程，并给出了测量数据可靠性定

28、量表示方法，同时证明了有偏测量可靠度优于无偏测量可靠度。其次，提出了多传感器有偏估计数据融合方法，证明了现有集中式与分布式无偏估计数据融合之间的等价性。最后，证明了多传感器有偏估计数据融合收敛于无偏估计数据融合。证明了方法的有效性。目前单传感器测量数据的处理方法主要有三种 ：平均值法1、加权平均法2和递推滤 波算法3.通过理论推导， 发现这些方法都是特殊形式的最小二乘估计(Least squareestimation ， LS)。基于模型y Hx w，x的最小二乘估计为:?LS H TH H T ymyiy为观测矢量，H为观测矩阵，未知矢量 x，当H=l，可化简为x_s 丄m可知，平均值法与其

29、具有相同的表达形式.采用类似的分析过程，可得另外两种方法与最小二乘估计是等价的.由于最小二乘估计是一种无偏估计，所以这种等价关系也说明上述三种数据处理方法具有无偏性，本文称之为无偏测量过程。无偏测量过程可以采用方差直接衡量测量可靠性，即方差越小测量可靠性越高。为了提高测量可靠性，国内外学者提出了多传感器数据融合的方法，旨在减小测量方差。目前多传感器数据融合常用的理论方法为线性无偏估计理论(简称多传感器无偏估计数据融合)，其中又以最小二乘估计应用最为 广泛.但是现有多传感器无偏估计数据融合方法存在两方面问题：1)融合结果可靠性均为定 性说明而无法量化表示， 即只能通过比较不同融合结果的方差定性地

30、判断融合结果可靠性 的优劣 ;2) 虽然多传感器无偏估计数据融合具有无偏性的优良性质， 但是并不能由此认为它的测量 结果一定是高可靠的因为根据高斯-马尔科夫定理可知，最小二乘估计方差有下界，所以此时无偏估计数据融合具有最小的方差， 但是当这个最小方差本身却很大时， 那么无偏 估计数据融合将不能保证测量数据的可靠性一定是可接受的 .但值得一提的是，无偏测量过程与最小二乘估计之间的等价关系为线性有偏估计算法用于提高测量可靠性成为可能. 如James-Stein 估计、 压缩最小二乘估计、 岭估计 (Ridge estimation ， RE)18 - 19 等 . 其 中岭估计是应用最为广泛的改进

31、最小二乘估计方法 . 本文以岭估计为基础提出多传感器有 偏估计数据融合方法， 岭估计长期以来一直是广泛用于改善最小二乘估计方差的有偏估计 方法 . 由于无偏测量与最小二乘估计之间是等价的，所以本文借鉴岭估计的思想通过引入较小的偏差改善无偏测量数据的方差， 并称之为有偏测量过程 . 在此基础上解决有偏测量 与无偏测量的可靠性定量表示问题.这种方法引入的偏差是可知的固定性偏差，且可以在一定程度上减小估计值的方差，其余并没有创新，不详细介绍了。2.6 基于小波去噪和数据融合的多传感器数据重建算法 为了从被噪声干扰的各个传感器测量值中获得更准确的测量结果，提出了一种基于小波去噪和多传感器数据融合的传感

32、器数据重建算法。 仿真和实验的结果都表明 : 由该算法重建得 到的各个传感器的重建数据的方差低于传感器测量值的方差。 可以认为多传感器数据重建算 法给出了对每一个传感器的更为准确的测量结果。一个传感器组， 利用每一个传感器的测量值对其加权， 进而对这组传感器的测量结果 进行数据融合以达到提高测量精度的目的.具体方法是在方差基本定义的基础上提出递归的估计方差的算法， 利用估计的方差估计出每个数据的权值，进而对电磁流量计的流量进行 递归估计 ，从而达到提高精度的目的。为了从受到不同噪声干扰的各个传感器测量值中获得更准确的各个传感器数据 ，本文 提出了一种基于小波去噪和多传感器数据融合的传感器数据重

33、建算法.该方法首先将每个传感器的测量值用小波阈值的方法去噪， 减小噪声对传感器测量值的影响 . 为了更好的重建传 感器信号， 先将各个传感器测量值进行归一化处理 ，再将归一化后的各个传感器测量值做 基于最小均方的数据融合 . 多传感器数据融合目的在于用较大的数据量， 充分利用对被测 目标的在时间与空间的信息 ，获得对被测量的描述。来自多传感器的信号所提供的信息具 有相关性、互补性和冗余性 ，将同源数据进行组合，可得到统计上的优势。基于小波去噪及多传感器数据融合的传感器数据重建算法：假设 N 个传感器在不同位置对同一测量值 Y测量，每个传感器测量值记为 Xj (j=1 , 2，N )由于测量中，

34、存在内部 外部噪声影响，测量值表示为Xj n Sn 勺n j 1,2, , N。 S(n)为被测量，ej(n)为第 j 个传感器在时刻加性噪声， Xj(n) 为第 j 个传感器在 n 时刻观测值。信号小波消噪方法主要通过设置阈值 .通过信号的离散小波变换， 计算所有小波系数 ， 然后剔除被认为跟噪声有关的小波系数。 例如通常的方法是设置阈值， 将小于阈值的小波系 数去掉 .最后，然后通过小波变换的逆变换来得到信号.数据融合数据重建算法：首先对 每 一 个 传 感 器 获得的一 组测 量 值 用这组数据中 的 最大测 量 值 归一 ：Yj n Yj n /MaxYj n j 1,2, ,N;n

35、1,2 ,I其中，MaxYj n是在估计长度I内第j个传感器的最大测量值。 Yj n为第j个传感器在n 时刻归一化后的测量值， 由于每个传感器收到噪声干扰程度不同， 所以偏离真实被测量程度不同， 对每个传感器根据一定原则确定权值，可从 N 个传感器得到估计值 Y。NNYWjYj ,Wj 1j 1j 1由于各传感器之间受到噪声干扰的程度不同，所以各传感器测量值的方差并不一致 ， 即各传感器测量值的可信度是不同的 若将较大的权值赋予可信度高的传感器 ，将较小的权 值赋予可信度小的传感器 ，就可以使估计值更精确地描述原信号。1,2, ,N，归一化权值为 Wj,N1Wjjj 器测量系统的数据融合问题进

36、行研究就具有更加广泛的应用价值。对Y反归一化,得到各传感器重建数据:Yj n Y nMax Yj n j 1,2, N; n1,2,1算法步骤：1置估计长度I； 2对各传感器测量值作小波阈值去噪处理；3采用MA模型用递归方法估计方差；4计算Wj ； 5计算Y,计算各传感器重建数据。2.7测量噪声相关情况下的多传感器数据融合对于测量噪声相关的多传感器测量模型,利用Cholesky 分解和单位下三角阵的求逆方法，将其转化为测量噪声互不相关的等价的多传感器伪测量模型,然后基于M arkov 估计，提出了一种测量噪声相关情况下多传感器数据融合的新方法。与直接利用原始传感器测量值的Markov 估计数据

37、融合方法相比，两者的计算精度相同,但新方法的计算复杂度却大大降低。 数值仿真实验进一步验证了新方法的有效性。所谓多传感器数据融合，就是将来自多个同类或异类传感器的数据（信息）进行综合处理，以获得比单一传感器更为准确可靠的结果。已有的多传感器数据融合方法,般利用含有加性噪声的线性测量方程来估计未知常值参数，大多假设各传感器的测量噪声之间互不相关。但是在实际应用中，由于各传感器通常处于同一测量环境，所以传感器的测量结果中除由于传感器自身精度限制而引入的测量误差外，共同的环境噪声的影响也不容忽略，而这往往会导致各传感器的测量噪声之间相关，所以对测量噪声相关情况下多传感为了解决测量噪声相关情况下的多传

38、感器测量数据融合问题，文献在最小二乘准则下，利用 Lagrange 乘子条件极值方法 ， 给出了一种最佳的线性数据融合方法， 但是仅适用 于被测参数为标量的情况， 无法直接扩展到参数为矢量的情况， 另外， 由于需要对累积观 测矢量的自相关阵直接求逆， 所以计算复杂度非常大 ; 文献则利用实对称矩阵的正交相似 变换实现了多传感器测量噪声互协方差阵的对角化 ，从而实现了各传感器测量噪声之间的 去相关 ，法在实际应用时比较困难。本文首先利用但是一般来说 ，这种对角化不能在有限步中完成， 只能通过迭代步骤求近似值， 所以该方Cholesky 分解和单位下三角阵的求逆方法将多然后基于与直接利用原，但新方

39、法的计算传感器的测量模型转化成各传感器的测量噪声互不相关的等价的伪测量模型，Markov 估计提出了一种测量噪声相关情况下的多传感器数据融合的新方法。始传感器测量值的 Markov 估计数据融合方法相比， 两者的结果相同复杂度大大降低。数值仿真实验进一步验证了本文方法的有效性。采用 N 个传感器对同一常值参数进行线性测量模型一般表示成ziH i x vi 测量噪声vi 服从均值为 0，方差为 Rii 的高斯分布z Hx v ，z z1 z2 zN ，HH1T H2T HTN假定各传感器的测量噪声相关，即R cov v,vR11 R12R21 R22R1NR2NRN1 RN 2的非对角块不全为

40、0，且 R 正定。噪声相关的多传感器数据融合对于上式所描述的多传感器线性测量系统偏估计 ，且对应的均方误差阵最小，由于 Markov 估计 9 为被测常值参数的无 所以根据 Markov 估计可得这 N 个传感器对于常值参数的最优融合估计为 ？FT 11 T 1H R H H R z，相应的融合均方误差Pf传感器噪声去相关：由于R = rij是一正定的实对称阵，根据矩阵的 Cholesky 分解可知，R可以唯一地分解成 Rldlt。其中，L =叩 为单位下三角阵D =diag( di 12，dNm)且正定，diriidkhki 1,2, , Nmk 1j 1rijdkl ikl jk口j 1,2, ,i 1djlij1j i0j i 1,i2, Nm对于单位下三角阵 l ,其逆阵存在，且仍为单位下三角阵，记m L 1 mj令 j=1，2，.，Nm，则0 i1,2,ji1j1mj1“ji 1ij 1l ik mkjk ji j2,j3, Nm将M分块阵表示M11O12O1NMM 21M 22O2N为单位下三角阵；MiiOij为零矩阵。在原格式两边左乘MM N1M N2M NN则