书签 分享 收藏 举报 版权申诉 / 13
立即下载 加入VIP,免费下载
当前位置： 首页 > 物理 > 计算；推算；计数 > 基于张量正则分解的欠定盲源分离中的混合矩阵识别方法.pdf

基于张量正则分解的欠定盲源分离中的混合矩阵识别方法.pdf

  • 上传人：a3
  • 文档编号：4704260
  • 上传时间：2018-10-28
  • 格式：PDF
  • 页数：13
  • 大小：3.01MB
    • 摘要
    申请专利号：

    CN201410613325.7

    		 申请日：

    2014.11.04
    公开号：

    CN104375976A

    		 公开日：

    2015.02.25
    当前法律状态：

    授权

    		 有效性：

    有权
    法律详情：

    授权|||实质审查的生效IPC(主分类):G06F 17/16申请日:20141104|||公开
    IPC分类号：

    G06F17/16; G10L21/0272(2013.01)I

    		 主分类号：

    G06F17/16
    申请人：

    西安电子科技大学
    发明人：

    罗勇江; 艾小凡; 汤建龙; 赵国庆; 杨松涛
    地址：

    710071陕西省西安市太白南路2号
    优先权：

    专利代理机构：

    陕西电子工业专利中心61205

    		 代理人：

    王品华; 王喜媛
    PDF完整版下载： PDF下载
    内容摘要

    本发明公开了一种基于张量正则分解的欠定盲源分离中的混合矩阵识别方法，主要解决现有技术在估计混合矩阵时受特定条件限制的问题。其实现步骤是：(1)对源信号进行采样得到观测数据；(2)利用观测数据的四阶累积量计算在不同时延下的四阶协方差矩阵；(3)将不同时延下的四阶协方差矩阵扩展成三阶张量的形式；(4)对三阶张量进行张量正则分解得到待识别混合矩阵的Khatri-Rao乘积矩阵；(5)利用特征值分解的方法对该乘积矩阵进行处理，得到混合矩阵的估计值。本发明具有识别精度高的优点，可用于语音、通信、雷达及生物医学领域源信号在时频混叠条件下的欠定盲源分离。

    权利要求书

    权利要求书
    1.  一种基于张量正则分解的欠定盲源分离中的混合矩阵识别方法，包括如下步骤：
    (1)在接收端对源信号进行采样，得到观测信号；
    (2)计算观测信号的四阶协方差矩阵其中，τ1＝0，τ2＝0，τ3为整数且τ3∈[0,R-1]，M为观测信号个数，R是大于P的正整数，取值为2*P，P为源信号的个数；
    (3)将四阶协方差矩阵Q(0,0,0)，Q(0,0,1)，…，Q(0,0,R-1)扩展成三阶张量Ti,j,k＝[Q(0,0,k)]i,j，1≤i,j≤M2，0≤k<R；
    (4)对三阶张量T进行张量正则分解，得到待识别混合矩阵的Khatri-Rao乘积源信号的四阶统计特性矩阵D∈CR×P及AQ的共轭矩阵
    (5)将待识别混合矩阵的Khatri-Rao乘积AQ的第e列元素be表示为矩阵的形式Be，
    其中，Be的每个元素为：Be[i,j]＝be((i-1)M+j),1≤i,j≤M，1≤e≤P，然后对Be进行特征值分解，其中最大的特征值对应的特征向量即为识别出的混合矩阵的第e列。

    2.  根据权利要求1所述的基于张量正则分解的欠定盲源分离中的混合矩阵识别方法，其中步骤(2)所述的计算观测信号的四阶协方差矩阵按如下步骤进行：
    (2.1)计算观测信号的四阶累积量：
    Ci,j,k,l(τ1,τ2,τ3)=m^i,j,k,l(τ1,τ2,τ3)-Rx(τ)Rx(τ3-τ2)-Rx(τ2)Rx(τ3-τ1)-Rx(τ3)Rx(τ2-τ1),]]>
    其中，τ1,τ2,τ3为时延，m^i,j,k,l(τ1,τ2,τ3)=1TΣt=1Txi(t)xj*(t+τ1)xk*(t+τ2)xl(t+τ3),]]>*表示复数共轭，Rx(τ1)=1TΣtTxi(t)xj*(t+τ1),]]>Rx(τ3-τ2)=1TΣtTxj*(t+τ2)xl(t+τ3),]]>Rx(τ2)=1TΣtTxi(t)xk*(t+τ2),]]>Rx(τ3-τ1)=1TΣtTxj*(t+τ1)xl(t+τ3),]]>Rx(τ3)=1TΣtTxi(t)xl(t+τ3),]]>Rx(τ2-τ1)==1TΣtTxj*(t+τ1)xk*(t+τ2).]]> 1≤i,j,k,l≤M，xi(t)表示第i路观测信号，当t≤0或者t≥T时，均取xi(t)＝0；
    (2.2)通过四阶累积量Ci,j,k,l(τ1,τ2,τ3)计算观测信号的四阶协方差矩阵Q(τ1,τ2,τ3)：
    Q(τ1,τ2,τ3)[M(i-1)+j，M(k-1)+l]＝Ci,j,k,l(τ1,τ2,τ3)，
    (2.3)在四阶协方差矩阵Q(τ1,τ2,τ3)中取τ1＝0，τ2＝0，τ3∈[0,R-1]，获得不同时延下的四阶协方差矩阵：Q(0,0,0)，Q(0,0,1)，…,Q(0,0,R-1)，其中R是大于P的正整数，取值为2*P，P为源信号的个数。

    3.  根据权利要求1所述的基于张量正则分解的欠定盲源分离中的混合矩阵识别方法，其中步骤(3)所述的将观测信号的四阶协方差矩阵扩展成张量T，按如下步骤进行：
    (3.1)根据源信号是统计独立的特性，将四阶协方差矩阵Q(τ1,τ2,τ3)表示为：
    Q(τ1,τ2,τ3)=Σp=1Pcp(τ1,τ2,τ3)[ap&CircleTimes;ap*][ap&CircleTimes;ap*]H=AQCs(τ1,τ2,τ3)AQH,]]>
    其中，ap表示混合矩阵的第p列，AQ=Δ[a1&CircleTimes;a1*,...,ap&CircleTimes;ap*]&Element;CM2×P,]]>为Kronecker乘积，Cs(τ1,τ2,τ3)是P×P的对角矩阵，Cs(τ1,τ2,τ3)＝diag[c1(τ1,τ2,τ3),…,cp(τ1,τ2,τ3)]，cp的定义为：
    cp(τ1,τ2,τ3)=ΔCum(sp(t),sp(t-τ1)*,sp(t-τ2)*,sp(t-τ3)),]]>
    其中，sp表示第p路未知源信号；
    (3.2)将四阶协方差矩阵Q(0,0,0)，Q(0,0,1)，…，Q(0,0，R-1)表示为：
    Q(0,0,0)=AQCs(0,0,0)AQH,Q(0,0,1)=AQCs(0,0,1)AQH,...Q(0,0,R-1)=AQCs(0,0,R-1)AQH,]]>
    其中，表示AQ的酉矩阵；
    (3.3)将四阶协方差矩阵Q(0,0,0)，Q(0,0,1)，Q(0,0,2)，…，Q(0,0,R-1)扩展成三阶张量T，其中，T的第(i,j,k)个元素为Ti,j,k＝[Q(0,0,k)]i,j，1≤i,j≤M2，0≤k<R-1，
    式中Q(0,0,0)为T的第一维切片，Q(0,0,1)为T的第二维切片，Q(0,0,2)为T的第三维切片，以此类推，Q(0,0,R-1)为T的第R维切片。

    4.  根据权利要求1所述的基于张量正则分解的欠定盲源分离中的混合矩阵识别方法，其中步骤(4)所述的对三阶张量进行张量正则分解，按如下步骤进行：
    (4.1)根据三阶张量的定义，求解目标矩阵1≤r≤3，使得代价函数f(U(1),U(2),U(3))最小，其中代价函数f(U(1),U(2),U(3))的定义由下式给出：

    其中，1≤i,j≤M2，0≤k<R-1，和分别为U(1)、U(2)和U(3)的第r个列向量，和分别为和的第i、j、k个元素；
    (4.2)采用线性搜索迭代最小二乘LS_ALS算法，对上述代价函数f(U(1),U(2),U(3))进行优化，得到目标矩阵U(r)，1≤r≤3，其中U(1)为待识别混合矩阵的Khatri-Rao乘积U(2)为AQ的共轭矩阵U(3)为源信号的四阶统计特性矩阵D∈CR×P。

    说明书

    说明书基于张量正则分解的欠定盲源分离中的混合矩阵识别方法
    技术领域
    本发明属于通信技术领域，特别涉及一种混合矩阵的识别方法，可用于语音、通信、雷达及生物医学领域源信号在时频混叠条件下的欠定盲源分离中。
    背景技术
    盲源分离BSS是指在未知的传输通道和源信号的条件下，仅通过传感器收到的观测信号来达到分离源信号的目的，该方法已经广泛地应用于语音信号处理、图像处理、雷达、通信及生物医学等各个领域。作为盲源分离的经典算法，独立分量分析ICA及其扩展算法大多用于解决观测信号数量等于或大于源信号数量条件下的问题，这种盲源分离称为正定或超定的盲源分离，但在实际过程中，往往需要解决源信号数量小于观测信号数量的问题，即欠定盲源分离UBSS。欠定盲源分离系统的线性瞬时模型为X(t)＝AS(t)+W(t)，其中X(t)∈CM表示观测信号，M为观测信号个数；S(t)∈CP为未知的源信号，P为源信号个数；W(t)∈CM表示加性噪声；未知混合矩阵A＝[a1,a2,…,aP]∈CM×P。在欠定盲分离系统中，观测信号个数小于源信号个数，即M<P。欠定条件下的混合矩阵盲识别，就是在未知的混合矩阵A和源信号S(t)的条件下从观测信号X(t)中识别出混合矩阵，是盲源分离问题的一个难点。
    目前，稀疏成分分析SCA是解决欠定盲源分离的问题的主要方法，大多数算法都是通过“两步法”来完成盲源分离，第一步则是估计出未知的传输通道，即混合矩阵模型，然后利用识别出的混合矩阵和稀疏分解的方法完成源信号的恢复，因此混合矩阵的识别在盲源分离问题中非常关键，其精度影响到后续源信号的恢复。部分学者利用信号的稀疏性，采用聚类的方法进行混合矩阵的识别，当源信号在时域上不满足稀疏性时，则利用傅立叶变换或小波变换等工具将信号变换到稀疏的频域上，然后利用聚类或势函数的方法识别混合矩阵，例如，NgutyenLin-Trung,ABelouchrani,KarimA-M.Separatingmoresourcesthansensorsusing time-frequencydistributions.EURASIPJournalonAppliedSignalProcessing,2005,17,pp.2828-2847，当源信号在时频域均混叠的情况下，该方法的性能则不理想。部分学者利用时频的方法，例如，陆凤波,黄知涛,彭耿,等,“基于时频分布的欠定混叠盲分离”，电子学报,2011,39(9),pp.2067-2072，该方法对观测信号进行时频处理，然后提取信号的自源时频点，利用自 源时频点构造张量模型并对该模型进行张量正则分解，从而完成混合矩阵的识别，但是在频域交叠比较严重的情况下自源时频点的提取并不理想，因此会影响混合矩阵的识别性能；还有一部分学者利用信号的统计特性，例如DeLathauwerL,CastaingJ,CardosoJ,“Fourth-order cumulant-basedblindidentificationofunderdeterminedmixtures”,IEEETransactionsonSignal Processing,2007,55(6),pp.2965-2973，该方法无需源信号满足稀疏的特性，只需要源信号是统计独立的非高斯信号，在实际过程中，这个条件往往是容易满足的，但该算法在求解的过程中，需要假设源信号具有相同符号的峭度，即反映振动信号分布特性的数值统计量，是归一化的四阶中心矩，而在源信号先验知识不足的情况下该条件往往很难满足，从而影响混合矩阵的识别精度。
    发明内容
    本发明的目的在于克服上述已有技术的不足，提出了一种基于张量正则分解的欠定盲源分离中的混合矩阵识别方法，以在无需特定条件的前提下，提高识别精度。
    本发明的技术方案是：在接收传感器处对经过未知通道的非高斯独立统计源信号进行采样，得到观测信号；利用观测信号的四阶累积量构建四阶协方差矩阵，并将四阶协方差矩阵表示成三阶张量模型；采用张量正则分解对三阶张量模型进行求解，对求解后得到的矩阵进行特征值分解，完成混合矩阵的识别。其实现步骤如下:
    (1)在接收端对源信号进行采样，得到观测信号；
    (2)计算观测信号的四阶协方差矩阵其中，τ1＝0，τ2＝0，τ3为整数且τ3∈[0,R-1]，M为观测信号个数，R是大于P的正整数，取值为2*P，P为源信号的个数；
    (3)将四阶协方差矩阵Q(0,0,0)，Q(0,0,1)，…,Q(0,0,R-1)扩展成三阶张量Ti,j,k＝[Q(0,0,k)]i,j，1≤i,j≤M2，0≤k<R；
    (4)对三阶张量T进行张量正则分解，得到待识别混合矩阵的Khatri-Rao乘积源信号的四阶统计特性矩阵D∈CR×P及AQ的共轭矩阵
    (5)将待识别混合矩阵的Khatri-Rao乘积AQ的第e列元素be表示为矩阵的形式Be，其中，Be的每个元素为：Be[i,j]＝be((i-1)M+j),1≤i,j≤M，1≤e≤P，然后对Be进行特征值分解，其中最大的特征值对应的特征向量即为识别出的混合矩阵的第e列。
    与现有技术相比，本发明具有以下优点：
    第一，本发明利用观测信号的四阶统计特性解决欠定盲源分离中的混合矩阵识别的问题，克服现有技术在混合矩阵识别时要求源信号满足稀疏性的缺点，使得本发明可以解决源信号在时频均混叠的条件下的欠定盲源分离中混合矩阵的识别问题。
    第二，本发明引入张量正则分解的方法对观测信号的四阶协方差矩阵扩展的三阶张量模型进行求解，克服现有技术中对观测信号的自源时频点提取困难的缺点，使得本发明提高了欠定盲源分离中混合矩阵的识别精度。
    第三，本发明将不同时延下的四阶协方差矩阵扩展成三阶张量模型，然后再进行求解，克服现有技术中需要源信号具有相同符号的峭度的缺点，从而可以解决源信号先验知识不足的条件下欠定盲源分离中混合矩阵的识别。
    附图说明
    图1为本发明的实现流程图；
    图2为本发明仿真实验中设置的4路源信号时域波形图；
    图3为本发明仿真试验中设置的4路源信号时频图；
    图4为用3个接收传感器对4路源信号的线性混合信号采样得到的3路观测信号；
    图5为本发明和现有方法对图4仿真得到的混合矩阵估计性能随信噪比变化的曲线图。
    具体实施方式
    以下参照附图对本发明作进一步详细的描述。
    参照图1，本发明实现步骤如下：
    步骤1：在接收端对源信号进行采样得到观测信号。
    M个传感器在t时刻对源信号进行等间隔采样，得到观测信号xi(t)，其中，1≤i≤M，t∈[1,2,…,N]，N为采样数据长度。
    步骤2：计算观测信号的四阶协方差矩阵。
    (2.1)计算观测信号的四阶矩：
    m^i,j,k,l(τ1,τ2,τ3)=1TΣt=1Nxi(t)xj*(t+τ1)xk*(t+τ2)xl(t+τ3),]]>
    其中，1≤i,j,k,l≤M，τ1,τ2,τ3分别为第j路，第k路以及第l路观测信号的时延；
    (2.2)计算观测信号的互相关：
    计算第i路观测信号xi(t)与第j路观测信号xj(t)在时延τ1下的互相关为：
    Bx(τ1)=1TΣtNxi(t)xj*(t+τ1),]]>
    计算第j路观测信号xj(t)与第l路观测信号xl(t)在时延τ3-τ2下的互相关为：
    Bx(τ3-τ2)=1TΣtNxj*(t+τ2)xl(t+τ3),]]>
    计算第i路观测信号xi(t)与第k路观测信号xk(t)在时延τ2下的互相关为：
    Rx(τ2)=1TΣtNxi(t)xk*(t+τ2),]]>
    第j路观测信号xj(t)与第l路观测信号xl(t)在时延τ3-τ1下的互相关为：
    Rx(τ3-τ2)=1TΣtNxj*(t+τ1)xl(t+τ3),]]>
    第i路观测信号xi(t)与第l路观测信号xl(t)在时延τ3下的互相关为：
    Bx(τ3)=1TΣtNxi(t)xl(t+τ3),]]>
    第j路观测信号xj(t)与第k路观测信号xk(t)在时延τ2-τ1下的互相关为：
    Rx(τ2-τ1)=1TΣtNxj*(t+τ1)xk*(t+τ2),]]>
    其中，*表示复数共轭；
    (2.3)计算观测信号的四阶累积量：
    Ci,j,k,l(τ1,τ2,τ3)=m^i,j,k,l(τ1,τ2,τ3)-Rx(τ1)Rx(τ3-τ2)-Rx(τ2)Rx(τ3-τ1)-Rx(τ3)Rx(τ2-τ1);]]>
    (2.4)通过四阶累积量Ci,j,k,l(τ1,τ2,τ3)计算观测信号的四阶协方差矩阵Q(τ1,τ2,τ3)：
    Q(τ1,τ2,τ3)[M(i-1)+j，M(k-1)+l]＝Ci,j,k,l(τ1,τ2,τ3)，
    (2.5)四阶协方差矩阵Q(τ1,τ2,τ3)中取τ1＝0，τ2＝0，τ3∈[0,R-1]，获得不同时延下的四阶协方差矩阵：Q(0,0,0)，Q(0,0,1)，…,Q(0,0,R-1)，其中R是大于P的正整数，取值为2*P，P为源信号的个数。
    步骤3：将观测信号的四阶协方差矩阵Q(0,0,0)，Q(0,0,1)，…,Q(0,0,R-1)扩展成三阶张量Ti,j,k＝[Q(0,0,k)]i,j(1≤i,j≤M2,0≤k<R)；
    (3.1)根据源信号是统计独立的特性，将四阶协方差矩阵Q(τ1,τ2,τ3)表示为：
    Q(τ1,τ2,τ3)=AQCS(τ1,τ2,τ3)AQH,]]>
    其中，ap表示混合矩阵的第p列，为Kronecker乘积，CS(τ1,τ2,τ3)是P×P的对角矩阵；
    (3.2)将四阶协方差矩阵Q(0,0,0)，Q(0,0,1)，…,Q(0,0,R-1)表示为：
    Q(0,0,0)=AQCS(0,0,0)AQH,Q(0,0,1)=AQCS(0,0,1)QQH,...Q(0,0,R-1)=AQCS(0,0,R-1)AQH,]]>
    其中，表示AQ的酉矩阵；
    (3.3)将四阶协方差矩阵Q(0,0,0)，Q(0,0,1)，Q(0,0,2)，…,Q(0,0,R-1)扩展成三阶张量T，其中，T的第(i,j,k)个元素为Ti,j,k＝[Q(0,0,k)]i,j，1≤i,j≤M2，0≤k<R-1，Q(0,0,0)为T的第一维切片，Q(0,0,1)为T的第二维切片，Q(0,0,2)为T的第三维切片，以此类推，Q(0,0,R-1)为T的第R维切片。
    步骤4：对三阶张量T进行张量正则分解：
    (4.1)根据三阶张量的定义，求解目标矩阵1≤r≤3，使得代价函数f(U(1),U(2),U(3))最小，其中代价函数f(U(1),U(2),U(3))的定义如下式：

    其中，1≤i,j≤M2，0≤k<R-1，和分别为U(1)、U(2)和U(3)的第r个列向量，和分别为和的第i、j、k个元素；
    4.2)采用线性搜索迭代最小二乘LS_ALS算法，对上述代价函数f(U(1),U(2),U(3))进行优 化，得到目标矩阵U(r)，1≤r≤3，其中U(1)为待识别混合矩阵的Khatri-Rao乘积U(2)为AQ的共轭矩阵U(3)为源信号的四阶统计特性矩阵D∈CR×P。
    步骤5：对待识别混合矩阵的Khatri-Rao乘积AQ进行特征值分解，得到混合矩阵估计值：
    (5.1)将待识别混合矩阵的Khatri-Rao乘积AQ的每一列be表示为矩阵的形式Be，其中，Be的每个元素表示为：Be[i,j]＝be((i-1)M+j),1≤i,j≤M，1≤e≤P；
    (5.2)对Be进行特征值分解，其中最大的特征值对应的特征向量即为识别的混合矩阵的第e列，1≤e≤P。
    本发明的效果可通过以下仿真进一步说明：
    为验证本发明的有效性和正确性，采用三种已有的混合矩阵识别方法与本发明方法同时对某混合矩阵进行识别。所有仿真实验均在Windows8.1操作系统下采用Matlab2012b软件实现。
    1)仿真参数
    采样率200MHz，采样点数为1024。4个源信号是时域上和频域上混叠的LFM信号，各个源信号归一化频率范围为[0.5,0]、[0,0.4]、[0.5,0.24]和[0.35,0.15]，其时域波形图和时频分布如图2和图3所示。
    根据空间自由传输模型，混合矩阵A中每一个元素的定义为：
    amp＝exp(2πj(αmcos(θp)cos(φp)+βmcos(θp)sin(φp)))，
    其中，αm＝(Ra/λ)cos(2π(m-1)/M)，βm＝(Ra/λ)sin(2π(m-1)/M)，Ra/λ＝0.55。信号入射方位角为θ1＝3π/10、θ2＝3π/10、θ3＝2π/5和θ4＝0；俯仰角φ1＝7π/10、φ2＝9π/10、φ3＝3π/5和φ4＝4π/5。
    为了评价算法对混合矩阵的识别性能，采用平均相对误差EA作为评价因子：
    EA=E{||A-A^||F||A||F},]]>
    式中，表示识别的混合矩阵，||·||F这里表示的是F范数。
    2)仿真内容
    在信噪比变化为0-35dB的范围，采用基于聚类的混合矩阵识别算法Cluster_based、基于 时频分布的混合矩阵识别算法TFDs_based、基于联合对角化的四阶累积量混合矩阵估计算法FOOBI以及本发明方法对图4中的3路观测信号做100次MonteCarlo仿真实验，得到混合矩阵的性能随信噪比的变化曲线，如图5所示。
    从图5可见，本发明方法在混合矩阵识别精度上要优于其他三种算法，从而验证了本发明方法对欠定盲源分离中的混合矩阵识别的有效性和正确性。

    举报
    举报     版权申诉
    版权申诉     PDF格式文档无特别注明外均可编辑修改；预览文档经过压缩，下载后原文更清晰！ 立即下载
    关 键  词：
    基于 张量 正则 分解 欠定盲源 分离 中的 混合 矩阵 识别 方法
      专利查询网所有文档均是用户自行上传分享，仅供网友学习交流，未经上传用户书面授权，请勿作他用。
    0条评论

    还可以输入200字符

    暂无评论，赶快抢占沙发吧。

    关于本文
    本文标题：基于张量正则分解的欠定盲源分离中的混合矩阵识别方法.pdf
    链接地址：https://www.zhuanlichaxun.net/p-4704260.html
    关于我们 - 网站声明 - 网站地图 - 资源地图 - 友情链接 - 网站客服 - 联系我们

    copyright@ 2017-2018 zhuanlichaxun.net网站版权所有
    经营许可证编号:粤ICP备2021068784号-1