改进的变步长维纳系统盲源分离

1、 改进的变步长维纳系统盲源分离方法研究摘要：在基于非线性盲源分离的维纳系统算法中，采用固定步长导致算法的收敛速度和稳态误差之间存在矛盾，直接影响分离算法的性能。为了解决该问题，本文提出了基于非线性函数的变步长维纳系统盲源分离方法，该方法将更新的步长以非线性函数的形式引入到分离算法中，使得稳态时参数更新的步长尽可能小，以避免发生震荡。变步长算法在分离过程中的每次更新都会使步长自动做出合理的调整，提高维纳系统盲源分离算法的性能。通过实验仿真表明，本文提出的维纳系统盲源分离方法可以更好的分离出信源信号，误差较低。关键词：盲源分离；维纳系统；非线性盲源分离；变步长算法中图分类号：TN911.7 文献标

2、识码：A引言盲源分离（Blind Source Separation, BSS ）是指在信源信号和信道参数都未知的条件下，从观测到的混合信号中估计出信源信号，被广泛用于多种信号处理和分析领域。目前的研究仍然主要集中于线性瞬时混合信号的盲源分离问题，但在许多的实际系统中，非线性混合模型更为常见。为此，近年来许多学者提出了非线性盲源分离问题。非线性盲源分离是一种针对非线性混合信号的盲源分离方法，已经被广泛应用于医学信号处理、通信信号处理、图像处理及故障诊断等方面1-4。目前具有代表性的非线性盲源分离方法主要有以下三类：基于互信息最小化的非线性盲源分离方法5-7，该方法采用互信息作为衡量相互独立性的

3、标准，互信息越小，分离效果越好；基于贝叶斯的非线性盲源分离方法8-9，成功利用贝叶斯网络理论处理非线性混合模型中各个变量和参数间的关系；基于参考信号的非线性盲源分离方法10-11，它是一种运用信源信号的先验信息作为参考信号的分析方法。维纳系统被应用于信号处理、生物、金融、社会以及心理分析等多方面，针对盲源分离问题，研究学者提出了基于非线性盲源分离的维纳系统BSS方法12-13。为了克服固定步长非线性盲源分离算法收敛性能差的问题，本文提出变步长和基于马尔可夫原理的后置非线性盲源分离算法，提高维纳系统BSS的性能。1 非线性BSS数学模型假设n个相互独立的未知源信号，首先经过未知的线性混合矩阵A(

4、nn维)，得到线性混合信号,再将分别通过一个非线性混合系统，得到观测信号： (1)盲源分离中的解混和混合是一个互逆的过程。此非线性混合系统的解混由两部分组成：第一部分是对非线性混合函数的求逆，即它是一个非线性反变换函数，用来补偿混合过程中的非线性失真；第二部分为线性解混矩阵，系统的输出信号可以定义为： (2)式中：，为对角矩阵，为置换矩阵，则非线性解混系统为： (3)后置非线性BSS混合-分离如图1所示。 图1 后置非线性BSS混合-分离结构框图Figure 1 Mixing and separating structure diagram of post nonlinear blind so

5、urce separation非线性盲源分离算法的关键是根据分离信号的相互统计独立性来对非线性函数和矩阵A求逆。2 基于非线性BSS的维纳系统基于后置非线性BSS的维纳系统将后置非线性BSS的线性混合矩阵A用一个线性滤波器来代替，其信号混合和分离数学模型如图2所示12-14，是信源信号，是未知的可逆滤波器，是未知的可逆无记忆非线性函数，是观测信号，是解混非线性函数，是解混矩阵，是对的估计。图2 维纳系统BSS混合模型和分离模型Figure 2 Mixing and separating model structure diagram of wiener system blind source

6、separation则维纳系统的输出为14： (4) 马尔可夫过程是一典型的随机过程,设是一个随机过程，当过程在时刻所处的状态为已知时，时刻所处的状态与过程在时刻之前的状态无关，这个无后效性的随机过程称为马尔可夫过程。由于后置非线性混合过程中观测信号是一个瞬时的混合过程，因此也满足马尔可夫过程。本文以最小化互信息作为衡量相互独立的标准，它被定义为： (5)其中，代表了分离信号各个分量的信息熵。对于q阶马尔可夫模型，条件互信息可以表示为： (6) 其中：表示均值。因为概率密度满足： (7) 式(6)中：是解混系统中非线性反变换函数对的求导。所以推得：(8)式中，不依赖于矩阵B及非线性函数g的参数

7、，所以被省略。目标函数变为：(9)盲源分离算法的相互独立性判据是互信息量传输最小化原则，所以通过对参数B和的调整，使输出的互信息尽可能小，已达到最佳的分离效果。B、的梯度计算式如下： (10) 其中：,。 (11)式中：因此可得： (12) (13)其中：为步长。3 基于非线性函数变步长方法一般来说，大的步长收敛速度比较快，但最小均方误差较大；相反，小的步长收敛速度比较慢，但最小均方误差较小。本文提出的基于非线性函数变步长算法，使步长随着分离算法的实际情况做动态变化。该方法通过误差的平方值去调节步长，表达式如下： (14)式中：是稳态误差，和为算法中的参数，这两个参数值的选取对分离算法的性能会

8、产生一定的影响。对参数B的调整，算法收敛后应该满足： (15)假设： (16) 则： (17)因此，可得到误差的一种表达式为： (18)由(13)和（17）可以得到变步长算法为： （19）如果选取的太大，对分离算法的稳定性会造成很大的影响。因此，给出的上界： （20）式（19）中：为的上界。的取值为： （21）同理，可以求出对参数调整时的步长。4 MATLAB仿真实验为了验证本文算法的有效性，选取随机信号作为源信号，其波形图如图3所示。图3源信号的波形图Figure 3 Waveform of source signal源信号先后经过一个滤波器（）和一个非线性系统（）后得到观测信号如图4。图4

9、观测信号的波形图Figure 4 Waveform of mixing signal本文通过变步长算法分离出的信号波形图如图5所示。对比图5与图3的波形图可以看出：两信号波形基本一致，说明了该算法很好地分离出了源信号。 图5 分离出的信号波形图Figure 5 Waveform of separating signal为验证所提出算法的性能，用最小均方误差评价分离效果，其值越接近于零，说明算法的分离性能越好。图6是基于非线性函数的变步长算法与固定步长算法的MMSE比较图。图6 最小均方误差性能对比图Figure 6 Performance comparison of minimum Mean

10、Squared Error由图6可知，固定步长算法分离出源信号需要更新199次左右，而变步长算法只需130次左右，收敛速度明显加快。采用固定步长算法收敛时的最小均方误差为1.3313，而变步长算法的最小均方误差为0.7335，可见误差性能有很大的改善。变步长算法与固定步长算法相比，分离效果有了明显的改善。5 结论本文对维纳系统盲源分离中的收敛速度和误差问题进行了研究，提出了非线性函数变步长的维纳系统盲源分离方法，该方法相对传统方法可以更好的分离出信源信号，而且收敛速度比较快，系统整体性能较好。参考文献1 任东晓, 叶茂，殷英. 基于互信息最小的非线性混合胎儿心电信号提取方法J.电子测量与仪器学

11、报，2010,24(7):680-685R Dongxiao, Ye Mao, Yin Ying. FECG extraction from nonlinear mixture based on minimization of mutual informationJ. Journal of electronic measurement and instrument, 2010,24(7):680-6852 K Diamantaras, T Papadimitriou, G Vranou. “Blind separationofmultiple binary sourcesfromone non

12、linear mixture ”C. 2011 IEEE International conference on acoustics, speech and signal processing.2011,22-173 M Almeida, L Almeida. Nonlinear separation of show through image mixtures using a physical model trained with ICAJ. Signal Processing,2012,29(4):872-8844 王晓伟, 石林锁. 自适应非线性BSS及其在齿轮故障诊断中的应用J.振动与

13、冲击,2012,31(10):45-48W Xiaowei, S linsuo. Adaptive nonlinear blind source separation and its application in gear faults diagnosisJ. Journal of vibration and shock,2012,31(10):45-485 LB Almeida. Linearandnonlinear ICA basedonmutual informationtheMISEPmethodJ.Signal processing, 2004,84(2):231-2456 Z. L

14、. Sun An extension of MISEP for post-nonlinear-linear mixture separationJ, IEEE Trans. Circuits Syst.2009,56(8):654 -658.7 L. T. Duarte and C. Jutten A mutual information minimization approach for a class of nonlinear recurrent separating systems, Proc. IEEE MLSP 2007:122 -127.8 C Wei, L C Khor, W L

15、 Woo. Post-nonlinear underdetermined ICA by Bayesian statisticsC. ICA 2006.IEEE Press, 2006: 773-780.9 Duarte, L. T.; Jutten, C.; Moussaoui, S. A Bayesian Nonlinear Source Separation Method for Smart Ion-Selective Electrode Arrays, Sensors Journal, IEEE, On page(s): 1763 1771, Volume: 9, Issue: 12,

16、Dec. 2009.10 R Dongxiao, Y Mao. Extracting post-nonlinear signal with referenceJ. Conputers and electrical engineering,2011:1171-118111 R Dongxiao, Y Mao. Extracting post-nonlinear signal with specific kurtosis rangeJ. Applied mathematics and computation.2012:5726-573812 S Zhangli, H Deshuang, Z Chu

17、nhou.Blind inversion of wiener system for single source using nonlinear blind source separationC.Conference on neural networks.2005:1235-123813 A Table, J Sole. Quasi-nonprarmetric blind inversion of wiener systemsJ. IEEE Trans on signal processing,2001，59(5):917-92414 Sole-Casals, J. and Jutten, C.

18、 “Post-Nonlinear Mixtures and Beyond”, vol. 16, pp. 67-74, 2004.Research on improved Wiener system Blind Source Separation methodAbstract：In the nonlinear blind source separation algorithm of Wiener system,the fixed step size lead to the contradiction between convergence speed and steady-state error

19、 of the algorithm,which directly affect the performance of the separation algorithm. To solve the problem, an algorithm based on variable step size blind source separation method of Wiener system is proposed in this paper.This method will update the step length in the form of a nonlinear function is introduced into the separation algorithm,makes the steady-state pa