一种基于稀疏重构的网格外目标波达方向估计方法 |
|||||||
申请号 | CN201410169416.6 | 申请日 | 2014-04-25 | 公开(公告)号 | CN103941220A | 公开(公告)日 | 2014-07-23 |
申请人 | 电子科技大学; | 发明人 | 段惠萍; 王艳艳; | ||||
摘要 | 一种基于稀疏重构的网格外目标 波达方向 估计方法,包括:1)使用一种基于字典平滑的新网格不匹配模型:通过对空间区域的 角 度进行网格划分,形成两组角度成平移等差关系的子字典,用网格上角度对应的导向矢量以线性组合的方式表示非网格角度上的导向矢量,得到新的基于字典平滑的网格不匹配模型。2)基于该模型在稀疏重构 框架 下,用联合稀疏重构方法求解稀疏 信号 ,用最小二乘方法求解方向偏差,求解 稀疏信号 和求解方向偏差的过程交替 迭代 直至收敛,最后将稀疏信号中的非零元素 位置 映射到角度空间并根据方向偏差进行补偿得到信号的波达方向。本 发明 比传统的泰勒模型更高的波达方向估计 精度 ,平滑后的字典维数降低从而使稀疏重构 算法 的运算量减少。 | ||||||
权利要求 | 1.一种基于稀疏重构的网格外目标波达方向估计方法,包括: |
||||||
说明书全文 | 一种基于稀疏重构的网格外目标波达方向估计方法技术领域[0001] 本发明属于阵列信号处理领域,具体涉及稀疏重构框架下的信号波达方向估计技术,针对字典不匹配引起的波达方向估计性能下降问题提出一种新网格不匹配模型及稀疏重构方法对网格外目标的波达方向进行估计。 背景技术[0002] 波达方向(Direction of Arrival,DOA)估计是阵列信号处理领域研究的一个重要课题,在雷达、声呐、地震探测、通信、医学等领域都有广泛的应用。经典的波达方向估计方法包括参数估计方法和非参数估计方法;参数估计方法主要基于最大似然估计算法(Maximum Likelihood,ML)得到,其中又包括确定性最大似然估计和随机性最大似然估计。最大似然估计算法的局限性在于初始值的选择直接影响到能否获得全局最优解,而且计算复杂度高。非参数估计法之一是Bartlett波束形成法,该方法的分辨率取决于阵列孔径。 1979年Schmidt(施密特)等人提出多重信号分类(Multiple Signal Classification,MUSIC)算法,这是一类子空间类的算法,可以达到超分辨,但是子空间类算法对相干或相关信源波达方向无法获得理想的估计结果。 [0003] 近年来随着压缩感知理论的发展,出现了一类基于稀疏重构框架的信号波达方向估计方法,这类方法依据压缩感知理论,对空间区域进行网格划分,由此形成基于网格的字典。在恢复信号的过程中,使用l1范数约束来保证信号的稀疏性。最经典的一个方法是Malioutov提出的l1-SVD方法,该方法为了降低计算的复杂性,用奇异值分解(Singular Value Decomposition,SVD)来降低信号的维数。l1-SVD算法可以分辨出间隔很近的两个相关信号,在低信噪比或者快拍数比较少的情况下,也可以达到很好的效果。一般稀疏重构框架下的信号波达方向估计方法均假设信号方向恰好落在网格上,但是实际上这个假设往往不能得到满足,当信号角度不在网格上时,l1-SVD算法就不能准确地估计出信号波达方向。此后,Aris Gretsistas等人提出了泰勒估计模型,将网格外导向矢量在最近网格处作一阶泰勒展开,通过估计方向偏差,对网格外目标的波达方向进行补偿,弥补了l1-SVD算法的不足。但是,泰勒估计模型的估计精度还较低,且受字典维数影响稀疏重构算法的运算量大,这些问题都有待进一步改进。 发明内容[0004] 本发明的目的在于改善稀疏重构框架下字典不匹配带来的信号波达方向估计性能降低问题,通过设计新的网格不匹配模型及相应的稀疏重构算法提高信号波达方向估计的精度。 [0005] 为了实现上述目的,本发明设计了如下技术方案: [0006] 一种基于稀疏重构的网格外目标波达方向估计方法,包括: [0007] 准备步骤:通过对空间区域的角度进行网格划分,形成两组角度成平移等差关系的子字典,依据字典平滑的思想,用网格上角度对应的导向矢量以线性组合的方式表示非网格角度上的导向矢量,建立一个基于字典平滑的网格不匹配模型; [0008] 估计步骤:基于网格不匹配模型,通过交替迭代的过程实现信号波达方向的估计和方向偏差的补偿:在求解信号的波达方向时,用联合稀疏重构的方法求解稀疏信号,用最小二乘的方法求解方向偏差,求解稀疏信号和求解方向偏差的过程交替迭代直至收敛,最后将稀疏信号的非零行位置映射到角度空间并根据方向偏差进行补偿得到估计的信号波达方向。 [0009] 该方法的具体步骤是: [0010] 1)、获取阵列接收数据Y,初始化方向偏差矢量Λ=diag(δ); [0011] 2)、根据输入的子字典构造平滑字典A=[I-Λk/r].A1+Λk/r.A2; [0012] 3)、根据接收数据维数用稀疏重构或联合稀疏重构算法求解稀疏信号X,优化准则为: [0013] 在使用联合稀疏重构算法时可以选择使用l1-SVD法降低数据维数; [0014] 4)、用最小二乘法求解δ,优化准则为: 求解公式为: k + k k [0015] δ =Cvec(R)=C+vec(Y-A1X),其中C的第n列是 [0016] 5)、判断是否收敛或迭代次数是否达到最大,如果不是重复执行步骤2)-4); [0017] 6)、在结束步骤2)-4)的循环后输出估计的稀疏信号 和方向偏差 将非零行的 [0018] 位置映射到角度空间并根据方向偏差矢量 进行补偿得到估计的信号波达方向。 [0019] 所述准备步骤的详细步骤是: [0020] 设信号波达方向为θ,对空间区域进行网格划分,离θ最近的两个网格是θ1和θ2,对a(θ1),a(θ2)分别在θ处进行一阶泰勒展开,得: [0021] a(θ1)=a(θ)+a′(θ)(θ1-θ),a(θ2)=a(θ)+a′(θ)(θ2-θ)[0022] 如果角度θ处对应的导向矢量a(θ)可以用a(θ1),a(θ2)线性表示: [0023] a(θ)=k1a(θ1)+k2a(θ2) [0024] 把a(θ1),a(θ2)的表达式代入上式,得方程: [0025] a(θ)=k1(a(θ)+a′(θ)(θ1-θ))+k2(a(θ)+a′(θ)(θ2-θ)[0026] =(k1+k2)a(θ)+[k1(θ1-θ)+k2(θ2-θ)]a′(θ) [0027] 上面方程对任意角度θ成立的条件是: [0028] k1+k2=1 [0029] k1(θ1-θ)+k2(θ2-θ)=0 [0030] 设网格间距为r,θ-θ1=δ,θ2-θ=r-δ,求解出k1,k2取值分别为: [0031] [0032] 把k1,k2代入a(θ)的表达式得到: [0033] [0034] δ表示方向偏差,是信号波达方向和最近网格之间的偏差;如果每个潜在目标,或者说字典中的每 [0035] 一个导向矢量都采用上面的线性组合表达方式,那么网格不匹配模型可以表达为如下形式: [0036] Y={[I-Λ/r].A1+Λ/r.A2}X+E [0037] 其中Y∈CM×T中包含了M个阵元测量得到的T次快拍,X表示T次快拍对应的稀疏信号矩阵,X [0038] 的每一列具有相同的稀疏结构,将X的非零行位置映射到角度空间就可以找到目标信号的波达方向; [0039] 其中A1,A2∈CM×N-1表示两个网格子集对应的子字典: [0040] A1=[a(θ1)a(θ2)...a(θN-1)] [0041] A2=[a(θ2)a(θ3)...a(θN)], [0042] N为网格总数,每个网格子集中包含N-1个网格; [0043] Λ=diag(δ)为对角矩阵包含了所有可能的方向偏差; [0044] δ为方向偏差矢量:δ=[δ1,δ2?,δN-1]T;E∈CM×T表示加性噪声。 [0045] 所述估计步骤的详细步骤是: [0046] 设方向偏差矢量δ为某一个固定值,采用联合稀疏重构的方法基于如下优化准则对最优稀疏信号进行求解: [0047] [0048] 其中: [0049] X2,1表示对X各行先求l2范数再对得到的列向量求l1范数; [0050] λ为正则化参数,λ值的选取与信噪比有关; [0051] 采用任何一种联合稀疏重构的方法都可以获得稀疏信号的估计 [0052] 在获得稀疏信号的估计 之后,对方向偏差矢量δ进行估计: [0053] 基于字典平滑的新网格不匹配模型经过整理可以重写为: [0054] [0055] 因此可以采用如下优化准则对方向偏差矢量δ进行估计: [0056] [0057] 假设 那么 上面的优化准则可以重新表达为: [0058] [0059] 其中bn表示矩阵B的第n列, 表示 的第n行,vec(R)表示将R各列堆叠实现矩阵矢量化,矩阵C的第n列是 该优化准则意味着δ可以用最小二乘法进行估计: 以上稀疏信号和方向偏差矢量估计的过程采用如下步骤交替迭代直至收敛: [0060] 步骤1:输入子字典A1,A2,观测矩阵Y,来波信号的个数K,最大迭代次数Max_iters,最大误差ε;k [0061] 步骤2:初始化循环次数和方向偏差矢量:k=1,δ =0;k k k k [0062] 步骤3:构造平滑字典A=[I-Λ/r].A1+Λ/r.A2,其中Λ =diag(δ),根据接收数据维数用稀疏重构或联合稀疏重构算法求解稀疏信号X,优化准则为: [0063] 步骤4:用最小二乘法求解δ,优化准则为: 求解k + k + k 公式为:δ =Cvec(R)=Cvec(Y-A1X),其中C的第n列是 k k-1 [0064] 步骤5:k=k+1,判断:如果k>Max_iters或δ-δ >ε则退出循环到步骤6,否则回到步骤3; [0065] 步骤6:在结束步骤3-5的循环后输出稀疏信号 和方向偏差 将非零行的位置映射到角度空间并根据方向偏差矢量 进行补偿得到估计的信号波达方向。 [0066] 本发明与现有技术相比具有如下优点:对于网格外目标,本发明设计的基于字典平滑的新网格不匹配模型可以比传统的泰勒模型提供更高的估计精度,平滑后的字典维数降低从而使稀疏重构算法的运算量减少。附图说明 [0067] 图1:本发明的流程图; [0068] 图2:本发明与传统一阶泰勒模型的方向估计误差随目标方向变化的比较图; [0069] 图3:本发明与传统一阶泰勒模型的方向估计均方根误差随信噪比变化的比较图。 具体实施方式[0070] 下面结合图1、图2、图3和仿真结果详细说明本发明的内容和效果: [0071] 一.基于字典平滑的新网格不匹配模型: [0072] 1)假设信号波达方向为θ,对空间区域进行网格划分,离θ最近的两个网格是θ1和θ2,对a(θ1),a(θ2)分别在θ处进行一阶泰勒展开,可以得到: [0073] a(θ1)=a(θ)+a′(θ)(θ1-θ),a(θ2)=a(θ)+a′(θ)(θ2-θ); [0074] 如果角度θ处对应的导向矢量a(θ)可以用a(θ1),a(θ2)线性表示: [0075] a(θ)=k1a(θ1)+k2a(θ2) [0076] 把a(θ1),a(θ2)的表达式代入上式: [0077] a(θ)=k1(a(θ)+a′(θ)(θ1-θ))+k2(a(θ)+a′(θ)(θ2-θ)[0078] =(k1+k2)a(θ)+[k1(θ1-θ)+k2(θ2-θ)]a′(θ) [0079] 上面方程对任意角度θ成立的条件是: [0080] k1+k2=1 [0081] k1(θ1-θ)+k2(θ2-θ)=0 [0082] 假设网格间距为r,θ-θ1=δ,θ2-θ=r-δ,可以求解出k1,k2取值分别为: [0083] [0084] 把k1,k2代入a(θ)的表达式得到: [0085] [0086] 注意δ表示方向偏差,是信号波达方向和最近网格之间的偏差。如果每个潜在目标,或者说字典中的每一个导向矢量都采用上面的线性组合表达方式,那么新的网格不匹配模型可以表达为如下形式:Y={[I-Λ/r].A1+Λ/r.A2}X+EM×T [0087] 其中Y∈C 中包含了M个阵元测量得到的T次快拍,X表示T次快拍对应的稀疏信号矩阵,X的每一列具有相同的稀疏结构,将X的非零行位置映射到角度空间就可以找到目标信号的波达方向。 [0088] 由此式可以看出,如果每个潜在目标方向对应的导向矢量用两个最近网格对应的导向矢量的线性组合来表示,最终稀疏重构的字典就是对两个子字典的平滑,因此申请人提出的新模型叫做字典平滑的网格不匹配模型。需要注意的是,在字典平滑的网格不匹配模型中,网格间距不能太大。M×N-1 [0089] 另外,A1,A2∈C 表示两个网格子集对应的子字典: [0090] A1=[a(θ1)a(θ2)?a(θN-1)] [0091] A2=[a(θ2)a(θ3)?a(θN)], [0092] N为网格总数,每个网格子集中包含N-1个网格。Λ=diag(δ)为对角矩阵包含T M×T了所有可能的方向偏差,δ为方向偏差矢量:δ=[δ1,δ2?,δN-1]。E∈C 表示加性噪声。 [0093] 可以看出在字典平滑的网格不匹配模型中,字典的列数为N-1,方向偏差向量的长度为N-1,比泰勒模型中减少一维,所以基于新提出的网格不匹配模型,稀疏重构的运算量减少。 [0094] 二.基于新网格不匹配模型的稀疏重构算法 [0095] 基于本发明设计的新网格不匹配模型,稀疏信号X和方向偏差矢量δ采用交替迭代的过程进行求解。 [0096] 首先,假设方向偏差矢量δ为某一个固定值,采用联合稀疏重构的方法基于如下优化准则对最优稀疏信号进行求解: [0097] [0098] 其中X2,1表示对X各行先求l2范数再对得到的列向量求l1范数。λ为正则化参数,λ值的选取与信噪比有关。以上优化问题为多测量矢量(Multiple Measurement Vectors,MMV)的联合稀疏重构问题,通过任何一种联合稀疏重构的方法都可以获得稀疏信号的估计 假设已知目标信号个数K,先通过l1-SVD的方法对多快拍的数据先进行降维,然后采用联合稀疏重构的方法,比如同时正交匹配追踪法(Simultaneous Orthogonal Matching Pursuit,SOMP)求解稀疏信号。 [0099] 其次,在获得稀疏信号的估计 之后,对方向偏差矢量δ进行估计: [0100] 基于字典平滑的新网格不匹配模型经过整理可以重写为: [0101] [0102] 因此可以采用如下优化准则对方向偏差矢量δ进行估计: [0103] [0104] 假设 那么 上面的优化准则可以重新表达为: [0105] [0106] 其中bn表示矩阵B的第n列, 表示 的第n行,vec(R)表示将R各列堆叠实现矩阵矢量化,矩阵C的第n列是 该优化准则意味着δ可以用最小二乘法进行估计: 以上稀疏信号和方向偏差矢量估计的过程采用如下步骤交替迭代直至收敛: [0107] 步骤1:输入子字典A1,A2,观测矩阵Y,来波信号的个数K,最大迭代次数Max_iters,最大误差ε; [0108] 步骤2:初始化循环次数和方向偏差矢量:k=1,δk=0; [0109] 步骤3:构造平滑字典A=[I-Λk/r]?A1+Λk/r.A2,其中Λk=diag(δk),根据接收数据维数用稀疏重构或联合稀疏重构算法求解稀疏信号X,优化准则为: [0110] 步骤4:用最小二乘法求解δ,优化准则为: 求解k + k k 公式为:δ =Cvec(R)=C+vec(Y-A1X),其中C的第n列是 [0111] 步骤5:k=k+1,判断:如果k>Max_iters或δk-δk-1>ε则退出循环到步骤6,否则回到步骤3; [0112] 步骤6:在结束步骤3-5的循环后输出稀疏信号 和方向偏差 将非零行的位置映射到角度空间并根据方向偏差矢量 进行补偿得到估计的信号波达方向。 [0113] 本发明的效果通过以下仿真进行说明: [0114] 仿真条件:阵元个数为M=8,阵元间距为半波长,网格划分的区域为[0°:180°],划分间隔为r=2°,假设有一个目标:K=1,快拍数T=200。 [0115] 仿真1:与传统一阶泰勒模型的方向估计误差随目标方向的变化进行比较。假设目标方向在[81°:171°]之间变化,以2°为间隔从而保证仿真的目标方向都在网格外,图2所示为方向估计误差随信源方向变化的曲线,可以看出对于网格外目标,本发明设计的基于字典平滑的网格不匹配模型能够提供优于传统泰勒模型的方向估计精度。 [0116] 仿真2:假设目标方向固定为121°,信噪比在[-10:15]dB范围内以1dB为间隔变化,图3所示为目标方向估计的均方根误差RMSE随信噪比SNR的变化曲线。可以看出与传统的泰勒模型相比,在不同的信噪比环境下本发明中的新网格不匹配模型以及基于该模型的稀疏重构算法可以提供具有更低均方根误差的波达方向估计结果。 [0117] 对于网格外目标,本发明提出的稀疏重构算法由于采用了一种基于字典平滑的新网格不匹配模型,可以比传统泰勒模型下的稀疏重构算法提供更高的波达方向估计精度,平滑后的字典维数降低从而使稀疏重构算法的运算量减少。 |