一种雷达均匀线性阵列幅度和相位的估计方法技术

本发明专利技术公开了一种雷达均匀线性阵列幅度和相位的估计方法，涉及阵列信号处理领域。包括：步骤1，设定雷达天线阵列具有M个阵元，其中前g个阵元的幅度和相位不存在误差；步骤2，将幅度、相位的预估计值作为均匀线性阵列的幅度、相位扰动初值；通过幅相扰动初值建立均匀线性阵列的阵列流型；步骤3，构建观测矩阵的子空间；步骤4，利用阵列流型及观测矩阵的子空间求解目标稀疏矩阵；步骤5，利用阵列流型及观测矩阵的子空间以及目标稀疏矩阵对阵列流型进行优化；步骤6，由优化后阵列流型矩阵对角线元素求出均匀线性阵列的幅度估计值和相位估计值。本发明专利技术在低信噪比且快拍数量较少的情况下，有效的对均匀线性阵列的幅度、相位进行估计。

【技术实现步骤摘要】
一种雷达均匀线性阵列幅度和相位的估计方法
本专利技术属于雷达阵列信号处理领域，涉及一种雷达均匀线性阵列幅度和相位的估计方法。
技术介绍
阵列信号处理是现代信号处理的一个重要研究分支。但在实际工程应用中，由于各种不可避免的误差，实际阵列流型往往会出现一定程度的偏差或扰动，这样，高分辨率的谱估计计算性能便会严重恶化，因而阵列的误差估计是阵列信号处理中很重要的研究方向之一。由于各种扰动、误差的影响，当实际模型与假设模型不符合时，许多基于理想模型基础上的谱估计性能将严重下降，因而必须采取补偿或者阵列误差校准的措施。实际工程应用时，通常会出现快拍数量有限且信噪比低的问题。当快拍数少时或信噪比低时，信号与噪声没有足够的时间解相关，噪声协方差矩阵也没有收敛，因此会对后期各种参数估计性能造成影响。目前，大多数阵列误差校准方法是通过对阵列扰动进行建模，将阵列误差校准逐渐转化为一个参数估计的问题的思想实现的。参数类的阵列校准方法可分为有源校准和自校准两类；有源校准通过在空间设置方位已知的辅助信源对阵列扰动参数进行离线估计，但这类校准算法对辅助信号源有较高的方位信息的要求，所以当辅助信号的方位信息有偏差时，这类校准算法会带来较大的偏差；而自校准通常根据某种优化函数对空间信源的方位与阵列的扰动参数进行联合估计，可以不需要方位已知的辅助信源，可以在线完成实际参数估计，所以校准的度较高，但由于误差参数与方位参数之间的耦合和某些病态的阵列结构，参数估计的唯一辨识往往无法保证，且其参数联合估计对应的高维、多模非线性优化问题带来了庞大的运算量。在期刊IEEETrans中，FriedlanderB和WeissAJ基于子空间原理，提出一种信源方位、阵列互耦、阵元增益及相位扰动交替迭代估计的阵列自校准技术，但该技术需要求解高维非线性的优化问题，运算量庞大、收敛速度慢，全局收敛性无法保证，而且对于均匀线阵，阵列扰动的参数的估计存在模糊问题。在期刊IEEETrans中，SeeCMS和PohBK通过对阵列观测空间内一系列离散方位上阵列导向矢量的测量，基于最大似然准则，提出了一种信源方位、校准矩阵以及阵元位置同时联合优化估计的方法。虽然该方法通过参数联合估计和高斯—牛顿梯度算法与遗传算法相混合的优化算法来提高参数估计的收敛性和全局搜索性能，但其对应的非线性优化问题有很高的维数导致了更庞大的运算量。《空间谱估计理论与算法》一书中提出了用辅助阵元法(ISM)来进行幅度相位误差的估计，通过引入少量校准的辅助阵元对阵元幅度相位误差进行无模糊联合估计，该方法只需要参数的一维搜索，不存在参数联合估计的局部收敛问题。但当快拍数据有限、信噪比低时，该估计结果会出现一些偏差。
技术实现思路
本专利技术的目的在于针对上述已有技术的不足，提出了一种雷达均匀线性阵列幅度和相位的估计方法，实现了在低信噪比且快拍数量较少的情况下，有效的对均匀线性阵列的幅度值、相位值进行估计，具有良好的稳健性。为达到上述目的，本专利技术采用以下技术方案予以实现。一种雷达均匀线性阵列幅度和相位的估计方法，其特征在于，包括以下步骤：步骤1，设定雷达天线阵列为均匀线性阵列，均匀线性阵列具有M个阵元，其中前g个阵元的幅度和相位不存在误差；步骤2，确定均匀线性阵列的导向矢量L和均匀线性阵列的幅度相位扰动矢量Γ；均匀线性阵列的导向矢量L表达式为：L＝[l1l2,...,lh,...,lM]，其中，lh表示均匀线性阵列第h个阵元的导向矢量，1≤h≤M，均匀线性阵列的幅度相位扰动矢量Γ表达式Γ＝[111,...,ΓM-g-1ΓM-g,...,ΓM]T，其中，扰动矢量Γ中前g个1代表幅度和相位不存在误差的阵元，后M-g个元素代表幅度和相位存在误差的阵元；M表示均匀线性阵列阵元数目，g表示幅度和相位不存在误差的阵元数目；对均匀线性阵列的幅度相位扰动矢量Γ进行预估计，得到均匀线性阵列的幅度相位扰动预估计矢量Γ′；均匀线性阵列的幅度相位扰动预估计矢量Γ′表达式为：Γ′＝[111,...,Γ′M-g-1Γ′M-g,...,Γ′M]T；(·)T表示转置操作，扰动矢量Γ′中前g个1代表幅度和相位不存在误差的阵元，后M-g个元素代表幅度和相位存在误差的阵元；利用预估计矢量Γ′和均匀线性阵列的导向矢量L构建均匀线性阵列的阵列流型D，即D＝diag(Γ′)·L，diag(.)表示矢量对角化，·表示矩阵点乘；步骤3，雷达天线接收T个快拍的目标回波数据，利用目标回波数据以及阵列流型D建立观测矩阵Y*；通过观测矩阵Y*求解观测矩阵Y*的子空间X；步骤4，利用阵列流型D构造矩阵Φ＝DA(θ)，A(θ)是目标的导向矢量；再利用观测矩阵Y*的子空间X求取矩阵Φ的支撑集Ω；构建支撑集Ω对应矩阵Φ中不为零的列向量组成的矩阵ΦΩ，并利用观测矩阵Y*的子空间X构建矩阵SΩ；矩阵SΩ组成目标稀疏矩阵S的非零列，目标的稀疏矩阵S的其余列为零，得到目标的稀疏矩阵S；步骤5，利用观测矩阵Y*的子空间X和目标的稀疏矩阵S对阵列流型D进行优化，得到优化后的阵列流型D*；步骤6，根据优化后的阵列流型D*求取均匀线性阵列的幅度估计值和相位估计值，即对优化后的阵列流型D*对角线上的数据求绝对值，得到均匀线性阵列的幅度估计值，对优化后的阵列流型D*对角线上的数据求角度，即求取对角线上的数据的实部与虚部的反正切值，得到均匀线性阵列的相位估计值。上述技术方案的特点和进一步改进在于：(1)步骤3包括以下子步骤：3a)接收T个快拍的目标回波数据，利用目标回波数据和阵列流型D建立观测矩阵Y*：Y*＝DA(θ)S+E其中，观测矩阵Y*∈CM×T,噪声矩阵E∈CM×T，稀疏矩阵S∈CN×T，A(θ)是目标的导向矢量，其中D表示阵列流型；M为该均匀线性阵列的阵元数，T为快拍数，N为采样点数，C为基矩阵；在[-90°,90°]间每隔1°进行采样，则采样点数N＝181；3b)设定X＝USW1/2为观测矩阵Y*的子空间；其中，W＝(ΛS-σe2IK)2ΛS-1为渐进最佳权值，US表示观测矩阵Y*的K个目标的奇异值对应的奇异向量组成的矩阵，ΛS为K个目标的奇异值组成的对角矩阵，且1≤K≤M；M表示均匀线性阵列阵元数目；λP表示观测矩阵Y*的第P个噪声的奇异值，且1≤P≤M；σe2表示噪声功率；IK为K×K的单位矩阵，(·)-1表示矩阵取逆操作。(2)步骤4包括以下子步骤：4a)设定矩阵Φ＝DA(θ)，A(θ)是目标的导向矢量，Ω为矩阵Φ的支撑集，支撑集Ω中包含了矩阵Φ中不为零的列的标号，元素β是支撑集Ω中的一个元素，通过求解以下公式得到第i次迭代中的元素βi：其中，表示在第i-1次迭代中由支撑集Ωi-1对应于矩阵Φ中不为零的列向量组成的矩阵，设定Ri-1表示在第i-1次迭代时观测矩阵Y*的子空间X的矩阵，设定R0为观测矩阵Y*的子空间X；(·)H表示共轭转置操作，p表示目标的稀疏矩阵S的行的序号，0<p≤N，N为采样点数；q表示目标的稀疏矩阵S的列的序号，0<q≤T，T为快拍数；a本文档来自技高网...

【技术保护点】
一种雷达均匀线性阵列幅度和相位的估计方法，其特征在于，包括以下步骤：步骤1，设定雷达天线阵列为均匀线性阵列，均匀线性阵列具有M个阵元，其中前g个阵元的幅度和相位不存在误差；步骤2，确定均匀线性阵列的导向矢量L和均匀线性阵列的幅度相位扰动矢量Γ；均匀线性阵列的导向矢量L表达式为：L＝[l1l2,...,lh,...,lM]，其中，lh表示均匀线性阵列第h个阵元的导向矢量，1≤h≤M，均匀线性阵列的幅度相位扰动矢量Γ表达式Γ＝[1 1 1,...,ΓM‑g‑1 ΓM‑g,...,ΓM]T，其中，扰动矢量Γ中前g个1代表幅度和相位不存在误差的阵元，后M‑g个元素代表幅度和相位存在误差的阵元；M表示均匀线性阵列阵元数目，g表示幅度和相位不存在误差的阵元数目；对均匀线性阵列的幅度相位扰动矢量Γ进行预估计，得到均匀线性阵列的幅度相位扰动预估计矢量Γ′；均匀线性阵列的幅度相位扰动预估计矢量Γ′表达式为：Γ′＝[1 1 1,...,Γ′M‑g‑1 Γ′M‑g,...,Γ′M]T；(·)T表示转置操作，扰动矢量Γ′中前g个1代表幅度和相位不存在误差的阵元，后M‑g个元素代表幅度和相位存在误差的阵元；利用预估计矢量Γ′和均匀线性阵列的导向矢量L构建均匀线性阵列的阵列流型D，即D＝diag(Γ′)·L，diag(.)表示矢量对角化，·表示矩阵点乘；步骤3，雷达天线接收T个快拍的目标回波数据，利用目标回波数据以及阵列流型D建立观测矩阵Y*；通过观测矩阵Y*求解观测矩阵Y*的子空间X；步骤4，利用阵列流型D构造矩阵Φ＝DA(θ)，A(θ)是目标的导向矢量；再利用观测矩阵Y*的子空间X求取矩阵Φ的支撑集Ω；构建支撑集Ω对应于矩阵Φ中不为零的列向量的矩阵ΦΩ，并利用观测矩阵Y*的子空间X构建矩阵SΩ；矩阵SΩ组成目标稀疏矩阵S的非零列，目标的稀疏矩阵S的其余列为零，得到目标的稀疏矩阵S；步骤5，利用观测矩阵Y*的子空间X和目标的稀疏矩阵S对阵列流型D进行优化，得到优化后的阵列流型D*；步骤6，根据优化后的阵列流型D*求取均匀线性阵列的幅度估计值和相位估计值，即对优化后的阵列流型D*对角线上的数据求绝对值，得到均匀线性阵列的幅度估计值，对优化后的阵列流型D*对角线上的数据求角度，即求取对角线上的数据的实部与虚部的反正切值，得到均匀线性阵列的相位估计值。...

【技术特征摘要】
1.一种雷达均匀线性阵列幅度和相位的估计方法，其特征在于，包括以下步骤：步骤1，设定雷达天线阵列为均匀线性阵列，均匀线性阵列具有M个阵元，其中前g个阵元的幅度和相位不存在误差；步骤2，确定均匀线性阵列的导向矢量L和均匀线性阵列的幅度相位扰动矢量Γ；均匀线性阵列的导向矢量L表达式为：L＝[l1l2，...，lh，...，lM]，其中，lh表示均匀线性阵列第h个阵元的导向矢量，1≤h≤M，均匀线性阵列的幅度相位扰动矢量Γ表达式Γ＝[111，...，ΓMg-1ΓM-g，...，ΓM]T，其中，扰动矢量Γ中前g个1代表幅度和相位不存在误差的阵元，后M-g个元素代表幅度和相位存在误差的阵元；M表示均匀线性阵列阵元数目，g表示幅度和相位不存在误差的阵元数目；对均匀线性阵列的幅度相位扰动矢量Γ进行预估计，得到均匀线性阵列的幅度相位扰动预估计矢量Γ′；均匀线性阵列的幅度相位扰动预估计矢量Γ′表达式为：Γ′＝[111，...，Γ′Mg-1Γ′M-g，...，Γ′M]T；(·)T表示转置操作，扰动矢量Γ′中前g个1代表幅度和相位不存在误差的阵元，后M-g个元素代表幅度和相位存在误差的阵元；利用预估计矢量Γ′和均匀线性阵列的导向矢量L构建均匀线性阵列的阵列流型D，即D＝diag(Γ′)·L，diag(.)表示矢量对角化，·表示矩阵点乘；步骤3，雷达天线接收T个快拍的目标回波数据，利用目标回波数据以及阵列流型D建立观测矩阵Y*；通过观测矩阵Y*求解观测矩阵Y*的子空间X；步骤3包括以下子步骤：3a)接收T个快拍的目标回波数据，利用目标回波数据和阵列流型D建立观测矩阵Y*：Y*＝DA(θ)S+E其中，观测矩阵Y*∈CM×T，噪声矩阵E∈CM×T，稀疏矩阵S∈CN×T，A(θ)是目标的导向矢量，其中D表示阵列流型；M为该均匀线性阵列的阵元数，T为快拍数，N为采样点数，C为基矩阵；在[-90°，90°]间每隔1°进行采样，则采样点数N＝181；3b)设定X＝USW1/2为观测矩阵Y*的子空间；其中，W＝(ΛS-σe2IK)2ΛS-1为渐进最佳权值，US表示观测矩阵Y*的K个目标的奇异值对应的奇异向量组成的矩阵，ΛS为K个目标的奇异值组成的对角矩阵，且1≤K≤M；M表示均匀线性阵列阵元数目；λP表示观测矩阵Y*的第P个噪声的奇异值，且1≤P≤M；σe2表示噪声功率；IK为K×K的单位矩阵，(·)-1表示矩阵取逆操作；步骤4，利用阵列流型D构造矩阵Ф＝DA(θ)，A(θ)是目标的导向矢量；再利用观测矩阵Y*的子空间X求取矩阵Ф的支撑集Ω；构建支撑集Ω对应于矩阵Ф中不为零的列向量的矩阵ФΩ，并利用观测矩阵Y*的子空间X构建矩阵SΩ；矩阵SΩ组成目标稀疏矩阵S的非零列，目标的稀疏矩阵S的其余列为零，得到目标的稀疏矩阵S；步骤5，利用观测矩阵Y*的子空间X和目标的稀疏矩阵S对阵列流型D进行优化，得到优化后的阵列流型D*；步骤5包括以下子步骤：5a)利用观测矩阵Y*的子空间X，将均匀线性阵列的幅度相位估计转换为优化问题，也就是公式(1)：其中，||·||F为Frobenius范数运算符，X为观测矩阵Y*的子空间，A(θ)是目标的导向矢量，S为目标的稀疏矩阵，D为均匀线性阵列的阵列流型，||S||∞，0定义为||·||∞，0表示混合范数，p表示目标的稀疏矩阵S的行的序号，0＜p≤N，N为采样点数；q表示目标的稀疏矩阵S的列的序号，0＜q≤T，T为快拍数，λ＞0为正则化参数；5b)将目标的稀疏矩阵S代入优化公式(1)得到优化求解阵列流型D的公式(2)；s.t.||S||∞，0＝常数(2)5c)将优化求解阵列流型D的公式(2)进行矢量化，得到阵列流型D的矢量化形式d，对阵列流型D的矢量化形式d进行迭代，得到阵列流型D的矢量化形式...

【专利技术属性】
技术研发人员：陈伯孝，杨明磊，朱芳芳，雷文英，
申请(专利权)人：西安电子科技大学，
类型：发明
国别省市：陕西;61