一种雷达均匀线性阵列幅度和相位的估计方法
【专利摘要】本发明公开了一种雷达均匀线性阵列幅度和相位的估计方法，涉及阵列信号处理领域。包括：步骤1，设定雷达天线阵列具有M个阵元，其中前g个阵元的幅度和相位不存在误差；步骤2，将幅度、相位的预估计值作为均匀线性阵列的幅度、相位扰动初值；通过幅相扰动初值建立均匀线性阵列的阵列流型；步骤3，构建观测矩阵的子空间；步骤4，利用阵列流型及观测矩阵的子空间求解目标稀疏矩阵；步骤5，利用阵列流型及观测矩阵的子空间以及目标稀疏矩阵对阵列流型进行优化；步骤6，由优化后阵列流型矩阵对角线元素求出均匀线性阵列的幅度估计值和相位估计值。本发明在低信噪比且快拍数量较少的情况下，有效的对均匀线性阵列的幅度、相位进行估计。
【专利说明】 一种雷达均匀线性阵列幅度和相位的估计方法
【技术领域】
[0001]本发明属于雷达阵列信号处理领域，涉及一种雷达均匀线性阵列幅度和相位的估计方法。
【背景技术】
[0002]阵列信号处理是现代信号处理的一个重要研究分支。但在实际工程应用中，由于各种不可避免的误差，实际阵列流型往往会出现一定程度的偏差或扰动，这样，高分辨率的谱估计计算性能便会严重恶化，因而阵列的误差估计是阵列信号处理中很重要的研究方向之一 O
[0003]由于各种扰动、误差的影响，当实际模型与假设模型不符合时，许多基于理想模型基础上的谱估计性能将严重下降，因而必须采取补偿或者阵列误差校准的措施。实际工程应用时，通常会出现快拍数量有限且信噪比低的问题。当快拍数少时或信噪比低时，信号与噪声没有足够的时间解相关，噪声协方差矩阵也没有收敛，因此会对后期各种参数估计性能造成影响。
[0004]目前，大多数阵列误差校准方法是通过对阵列扰动进行建模，将阵列误差校准逐渐转化为一个参数估计的问题的思想实现的。参数类的阵列校准方法可分为有源校准和自校准两类；有源校准通过在空间设置方位已知的辅助信源对阵列扰动参数进行离线估计，但这类校准算法对辅助信号源有较高的方位信息的要求，所以当辅助信号的方位信息有偏差时，这类校准算法会带来较大的偏差；而自校准通常根据某种优化函数对空间信源的方位与阵列的扰动参数进行联合估计，可以不需要方位已知的辅助信源，可以在线完成实际参数估计，所以校准的度较高，但由于误差参数与方位参数之间的耦合和某些病态的阵列结构，参数估计的唯一辨识往往无法保证，且其参数联合估计对应的高维、多模非线性优化问题带来了庞大的运算量。
[0005]在期刊IEEE Trans中，Friedlander B和Weiss A J基于子空间原理，提出一种信源方位、阵列互耦、阵元增益及相位扰动交替迭代估计的阵列自校准技术，但该技术需要求解高维非线性的优化问题，运算量庞大、收敛速度慢，全局收敛性无法保证，而且对于均匀线阵，阵列扰动的参数的估计存在模糊问题。在期刊IEEE Trans中，See C M S和Poh BK通过对阵列观测空间内一系列离散方位上阵列导向矢量的测量，基于最大似然准则，提出了一种信源方位、校准矩阵以及阵元位置同时联合优化估计的方法。虽然该方法通过参数联合估计和高斯一牛顿梯度算法与遗传算法相混合的优化算法来提高参数估计的收敛性和全局搜索性能，但其对应的非线性优化问题有很高的维数导致了更庞大的运算量。
[0006]《空间谱估计理论与算法》一书中提出了用辅助阵元法(ISM)来进行幅度相位误差的估计，通过引入少量校准的辅助阵元对阵元幅度相位误差进行无模糊联合估计，该方法只需要参数的一维搜索，不存在参数联合估计的局部收敛问题。但当快拍数据有限、信噪比低时，该估计结果会出现一些偏差。
【发明内容】
[0007]本发明的目的在于针对上述已有技术的不足，提出了一种雷达均匀线性阵列幅度和相位的估计方法，实现了在低信噪比且快拍数量较少的情况下，有效的对均匀线性阵列的幅度值、相位值进行估计，具有良好的稳健性。
[0008]为达到上述目的，本发明采用以下技术方案予以实现。
[0009]一种雷达均匀线性阵列幅度和相位的估计方法，其特征在于，包括以下步骤:
[0010]步骤1，设定雷达天线阵列为均匀线性阵列，均匀线性阵列具有M个阵元，其中前g个阵元的幅度和相位不存在误差；
[0011]步骤2，确定均匀线性阵列的导向矢量L和均匀线性阵列的幅度相位扰动矢量 ；均匀线性阵列的导向矢量L表达式为:L = [I1I2,...，lh，...，1M]，其中，Ih表示均匀线性阵列第h个阵元的导向矢量，l<h<M，均匀线性阵列的幅度相位扰动矢量表达式 =[I I ι，...， Μτ1 rM_g，...， Μ]Τ，其中，扰动矢量中前g个I代表幅度和相位不存在误差的阵元，后M-g个元素代表幅度和相位存在误差的阵元;M表示均匀线性阵列阵元数目，g表示幅度和相位不存在误差的阵元数目；
[0012]对均匀线性阵列的幅度相位扰动矢量进行预估计，得到均匀线性阵列的幅度相位扰动预估计矢量 ，;均匀线性阵列的幅度相位扰动预估计矢量 ，表达式为: ，=[I I i，...，r' Mf1 r' M_g，...，r' Μ]τ;(.)τ表示转置操作，扰动矢量 r'中前g个I代表幅度和相位不存在误差的阵元，后M-g个元素代表幅度和相位存在误差的阵元；
[0013]利用预估计矢量 ，和均匀线性阵列的导向矢量L构建均匀线性阵列的阵列流型D,即D = diag(r ' ).L, diag(.)表示矢量对角化，.表示矩阵点乘；
[0014]步骤3，雷达天线接收T个快拍的目标回波数据，利用目标回波数据以及阵列流型D建立观测矩阵Y* ;通过观测矩阵Y*求解观测矩阵Y*的子空间X ；
[0015]步骤4，利用阵列流型D构造矩阵Φ = DA( Θ )，Α( Θ )是目标的导向矢量；再利用观测矩阵Y*的子空间X求取矩阵Φ的支撑集Ω ;构建支撑集Ω对应矩阵Φ中不为零的列向量组成的矩阵φΩ，并利用观测矩阵Y*的子空间X构建矩阵Si2 ;矩阵Si2组成目标稀疏矩阵S的非零列，目标的稀疏矩阵S的其余列为零，得到目标的稀疏矩阵S ；
[0016]步骤5，利用观测矩阵Y*的子空间X和目标的稀疏矩阵S对阵列流型D进行优化，得到优化后的阵列流型D* ；
[0017]步骤6，根据优化后的阵列流型D*求取均匀线性阵列的幅度估计值和相位估计值，即对优化后的阵列流型D*对角线上的数据求绝对值，得到均匀线性阵列的幅度估计值，对优化后的阵列流型D*对角线上的数据求角度，即求取对角线上的数据的实部与虚部的反正切值，得到均匀线性阵列的相位估计值。
[0018]上述技术方案的特点和进一步改进在于:
[0019](I)步骤3包括以下子步骤:
[0020]3a)接收T个快拍的目标回波数据，利用目标回波数据和阵列流型D建立观测矩阵Y*:
[0021]Y* = DA ( Θ ) S+E
[0022]其中，观测矩阵Y* e Cmxt,噪声矩阵E e Cmxt,稀疏矩阵S e CnxtjA ( Θ )是目标的导向矢量，其中D表示阵列流型…为该均匀线性阵列的阵元数，T为快拍数，N为采样点数，C为基矩阵；在[-90°，90° ]间每隔1°进行采样，则采样点数N= 181 ；
[0023]3b)设定X = UsW172为观测矩阵Y*的子空间；
[0024]其中，W = (As-O62Ik)2As-1为渐进最佳权值，Us表示观测矩阵Y*的K个目标的奇异值对应的奇异向量组成的矩阵，As为K个目标的奇异值组成的对角矩阵，且I < K < M ;
M表示均匀线性阵列阵元数目
【权利要求】
1.一种雷达均匀线性阵列幅度和相位的估计方法，其特征在于，包括以下步骤: 步骤1，设定雷达天线阵列为均匀线性阵列，均匀线性阵列具有M个阵元，其中前g个阵元的幅度和相位不存在误差； 步骤2，确定均匀线性阵列的导向矢量L和均匀线性阵列的幅度相位扰动矢量I' ;均匀线性阵列的导向矢量L表达式为:L = [I1I2,...，lh，...，1M]，其中，Ih表示均匀线性阵列第h个阵元的导向矢量，I ≤h<M，均匀线性阵列的幅度相位扰动矢量I'表达式I' = [I II,..., I'Μτ1 rM_g，...，I'Μ]Τ，其中，扰动矢量I'中前g个I代表幅度和相位不存在误差的阵元，后M-g个元素代表幅度和相位存在误差的阵元…表示均匀线性阵列阵元数目，g表示幅度和相位不存在误差的阵元数目； 对均匀线性阵列的幅度相位扰动矢量I'进行预估计，得到均匀线性阵列的幅度相位扰动预估计矢量I'，;均匀线性阵列的幅度相位扰动预估计矢量I'，表达式为:I' ' = [II i，...，r' Mf1 r' M_g，...，r' Μ]τ;(.)τ表示转置操作，扰动矢量I' '中前g个I代表幅度和相位不存在误差的阵元，后M-g个元素代表幅度和相位存在误差的阵元； 利用预估计矢量I''和均匀线性阵列的导向矢量L构建均匀线性阵列的阵列流型D，即D = diag(r ' ).L, diag(.)表示矢量对角化，.表示矩阵点乘； 步骤3，雷达天线接收T个快拍的目标回波数据，利用目标回波数据以及阵列流型D建立观测矩阵r;通过观测矩阵Y*求解观测矩阵Y*的子空间X ； 步骤4，利用阵列流型D构造矩阵Φ = DA( Θ )，Α( Θ )是目标的导向矢量；再利用观测矩阵Y*的子空间X求取矩阵Φ的支撑集Ω ;构建支撑集Ω对应于矩阵Φ中不为零的列向量的矩阵ΦΩ，并利用观测矩阵Y*的子空间X构建矩阵Si2 ;矩阵Si2组成目标稀疏矩阵S的非零列，目标的稀疏矩阵S的其余列为零，得到目标的稀疏矩阵S ； 步骤5，利用观测矩阵Y*的子空间X和目标的稀疏矩阵S对阵列流型D进行优化，得到优化后的阵列流型D* ； 步骤6，根据优化后的阵列流型D*求取均匀线性阵列的幅度估计值和相位估计值，即对优化后的阵列流型D*对角线上的数据求绝对值，得到均匀线性阵列的幅度估计值，对优化后的阵列流型D*对角线上的数据求角度，即求取对角线上的数据的实部与虚部的反正切值，得到均匀线性阵列的相位估计值。
2.根据权利要求1所述的一种雷达均匀线性阵列幅度和相位的估计方法，其特征在于，步骤3包括以下子步骤: 3a)接收T个快拍的目标回波数据，利用目标回波数据和阵列流型D建立观测矩阵Y*: Y* = DA ( Θ ) S+E 其中，观测矩阵Y* e Cmxt,噪声矩阵E e Cmxt，稀疏矩阵S ∈ CNXT, A( Θ )是目标的导向矢量，其中D表示阵列流型；M为该均匀线性阵列的阵元数，T为快拍数，N为采样点数，C为基矩阵；在[-90°，90° ]间每隔1°进行采样，则采样点数N = 181 ； 3b)设定X = UsWv2为观测矩阵Y*的子空间； 其中，W = (As-O62Ik)2As-1为渐进最佳权值，Us表示观测矩阵Y*的K个目标的奇异值对应的奇异向量组成的矩阵，As为K个目标的奇异值组成的对角矩阵，且I≤K≤M ;M表示均匀线性阵列阵元数目
,λ P表示观测矩阵Y*的第P个噪声的奇异值，且I≤P≤M ; σ e2表示噪声功率；IK为KXK的单位矩阵，(.)-1表示矩阵取逆操作。
3.根据权利要求1所述的一种雷达均匀线性阵列幅度和相位的估计方法，其特征在于，步骤4包括以下子步骤: 4a)设定矩阵Φ =?Α(θ)，Α(θ)是目标的导向矢量，Ω为矩阵Φ的支撑集，支撑集Ω中包含了矩阵Φ中不为零的列的标号，元素β是支撑集Ω中的一个元素，通过求解以下公式得到第i次迭代中的元素β i:
其中，％,4表示在第1-Ι次迭代中由支撑集Ω H对应于矩阵Φ中不为零的列向量组成的矩阵，设记-Q0 = 0; Rh表示在第i_l次迭代时观测矩阵Y*的子空间X的矩阵，设定Rtl为观测矩阵Y*的子空间X ； (.)Η表示共轭转置操作，P表示目标的稀疏矩阵S的行的序号，0〈ρ ( N, N为采样点数；q表示目标的稀疏矩阵S的列的序号，0〈q ^ Τ,Τ为快拍数；arg(.)表示求解最优化，1 = 1,2, 3,…，K, K为目标数，I I.I I =O表示无穷范数运算符； 4b)利用第i次迭代中的元素β i求解第i次迭代时支撑集Qi=QHUiii; 4c)利用第i次迭代时支撑集Ω i和观测矩阵Y*的子空间X得到第i次迭代时矩阵Ri表达式为下式:
其中φ犮示在第i次迭代中由支撑集Qi对应于矩阵Φ中不为零的列向量组成的矩阵；χ为观测矩阵Y*的子空间；(.)H表示共轭转置操作，(.)_1表示矩阵取逆操作； 4d)令迭代次数i增加I，重复迭代以上几步4a)-4c)，直到i等于K，得到Ωκ;设定支撑集Ω = Ωκ，求解出支撑集Ω ； 4e)构建支撑集Ω对应于集合Φ中不为零的列向量的矩阵ΦΩ，再利用观测矩阵Y*的子空间X构建矩阵Si2，矩阵Si2公式为
其中，ΦΩ表示支撑集Ω对应于矩阵Φ中的列向量集合，(.)Η表示共轭转置操作，X为观测矩阵Y*的子空间，(.)-1表示矩阵取逆操作； 4f)矩阵Si2组成目标稀疏矩阵S的非零列，目标的稀疏矩阵S的其余列为零，得到目标的稀疏矩阵S。
4.根据权利要求1所述的一种雷达均匀线性阵列幅度和相位的估计方法，其特征在于，步骤5包括以下子步骤: 5a)利用观测矩阵Y*的子空间X，将均匀线性阵列的幅度相位估计转换为优化问题，也就是公式⑴:
其中，Il.||F*FrobeniUS范数运算符，X为观测矩阵Y*的子空间，Α(θ)是目标的导向矢量，S为目标的稀疏矩阵，D为均匀线性阵列的阵列流型，I |S| I μ定义为
，M.11 =0,0表示混合范数，P表示目标的稀疏矩阵s的行的序号，ο〈ρ ≤N,N为米样点数；q表不目标的稀疏矩阵S的列的序号，0〈q ≤ T, T为快拍数，λ >0为正则化参数； 5b)将目标的稀疏矩阵S代入优化公式(I)得到优化求解阵列流型D的公式(2);
5c)将优化求解阵列流型D的公式(2)进行矢量化，得到阵列流型D的矢量化形式d，对阵列流型D的矢量化形式d进行迭代，得到阵列流型D的矢量化形式d的第j次迭代值Clj的表达式为下式(3):
其中，vec(.)表示矢量化，Cltl为阵列流型D对角线元素，j表示迭代次数，Il.||2表示2范数运算符，X为观测矩阵Y*的子空间，Α(θ)是目标的导向矢量，S表示目标的稀疏矩阵，arg(.)表示求解最优化； 5d)将第j次迭代值4进行I范数加权得到第j次迭代值范数加权值(，d*j的表达
其中，炉=Im 0W, W表示观测矩阵Y*中行的互协方差矩阵的非零对角矩阵，其中Im表示MXM的单位矩阵，M表示均匀线性阵列阵元数目，vec(.)表示矢量化，| |.| |工表示I范数运算符，A ( Θ )是目标的导向矢量，S表示目标的稀疏矩阵，d*0由Cltl加权I范数得至Ij, arg(.)表示求解最优化； 5e)求解公式⑷，若达到终止条件||d*j-d*j-1|| 2〈 ε，取ε = 10_4，设定第j次迭代值范数加权值为阵列流型D矢量化并加权I范数矢量(T ;若没有达到终止条件，则令j增加I，继续求解公式⑷； 5f)通过阵列流型D矢量化并加权I范数矢量(T得到优化后的阵列流型D%即D* =diag(cf) ；diag(.)表示矢量对角化。
【文档编号】G01S7/02GK104076332SQ201410290130
【公开日】2014年10月1日 申请日期:2014年6月25日 优先权日:2014年6月25日
【发明者】陈伯孝, 杨明磊, 朱芳芳, 雷文英 申请人:西安电子科技大学