互耦条件下基于四阶累积量稀疏表示的mimo雷达波达方向估计方法

互耦条件下基于四阶累积量稀疏表示的mimo雷达波达方向估计方法
【专利摘要】本发明属于单基地MIMO雷达系统技术领域，具体涉及互耦条件下基于四阶累积量稀疏表示的MIMO雷达波达方向估计方法。本发明包括：发射阵列发射相互正交的相位编码信号，接收端进行匹配滤波处理后获得接收数据，并利用发射和接收阵列都具有的互耦矩阵带状对称Toeplitz结构特点，通过线性变换消除未知互耦的影响；构造降维转换矩阵，对消除互耦后的数据进行降维处理，进而基于新的数据矩阵构造具有特殊形式的四阶累积量矩阵。本发明由于四阶累积量技术和加权稀疏表示框架的应用，成功地抑制了色噪声，在高斯色噪声情况下，本发明提供精确的波达方向估计，具有更高的角度分辨率和更好的角度估计性能。
【专利说明】
互耦条件下基于四阶累积量稀疏表示的ΜI MO雷达波达方向估 计方法
技术领域
[0001]本发明属于单基地ΜΜ0雷达系统技术领域，具体涉及互耦条件下基于四阶累积量 稀疏表示的ΜΙΜΟ雷达波达方向估计方法。
【背景技术】
[0002] 多输入多输出（ΜΙΜΟ)雷达由新型的阵列结构构成，发射相互正交的波形，由于比 传统相控阵雷达具有更高分辨率和更好的参数可识别性等很多优势，而吸引越来越多的关 注。ΜΙΜΟ雷达可以分为两类，统计ΜΙΜΟ雷达和相干ΜΙΜΟ雷达（IEEESignalProcessingMagazi ne，2007,24(5) :106-114)。统计]\ΠΜ0雷达通过大天线距获得空间增益，而相干ΜΜ0雷达通 过紧密的天线距形成虚拟阵列大孔径，获得更高的空间分辨率和更多的自由度。相干ΜΜ0 雷达包括双基地和单基地ΜΜ0雷达，前者发射阵列和接收阵列是互相分开的，后者是紧密 共置的。本发明中，我们研究的是单基地ΜΙΜ0雷达中的波达方向估计问题。
[0003] 角度估计在信号处理和雷达应用中具有重要作用，利用ΜΜ0雷达和传统阵列信号 模型的相似之处，大量关于ΜΙΜ0雷达D0A估计的文献涌现出来，其中最具有代表性的是基于 角度搜索或旋转不变特性的子空间角度估计方法，具体包括MUSIC，Capon和ESPRIT方法等。 MUSIC和Capon具有很大的计算复杂度而不利于信号实时处理。ESPRIT算法具有高运算效率 和更好的角度估计性能。派生的方法包括降维Capon，降维ESPRIT (RD-ESPRIT， ElectronicsLetters，2011，47(4) :283-284)和酉ESPRIT，RD-ESPRIT和酉ESPRIT比ESPRIT 方法具有更低的计算量。发射波束域能量集中技术(IEEETransactionsonSignalProcessin g，2011，59(6) :2669-2682)，使信噪比（SNR)增益得到最大化。然而，当考虑互耦时，以上子 空间方法性能下降甚至失效。为了消除互耦影响，一些方法被提出来，相似MUSIC (SignalProcessing，2012,92(2) :517-522)通过角度搜索对角度进行估计，求根MUSIC具有 更低的计算复杂度，相似ESPRIT方法（3丨811 &1?仰〇688丨即，2012,92(12):3039-3048)具有良 好的角度估计性能，同时避免了空间峰值搜索。
[0004] 在阵列信号处理中，经常遇到空间高斯色噪声的情况而不是白噪声，以上提到方 法的波达方向估计性能进一步下降。高阶累积量在D0A估计上的应用能够解决这个问题，因 为四阶累积量(F0C)能够抑制任意方差的高斯噪声。基于这个特性，F0C-MUSIC，F0C-ESPRIT 等方法被提出来。互親条件下，一种改进的F0C-MUSIC方法（SignalProcessing，2009,89: 1839-1843)在高斯白噪声和色噪声条件下都具有良好的估计性能。
[0005] 新兴的稀疏表示在信号分析领域吸引越来越多的关注，通过找到数据最稀疏的表 示能够将其应用到D0A估计中。所有仿真结果证实，与传统子空间方法相比，稀疏表示方法 具有显著优点，能够更好地适应特殊情况，此外提供更高角度分辨率并且更少地依赖入射 信号数目的先验信息。为了估计波达方向一些方法被提出来，例如1 i-SVD，F0CUSS，1 rSRACV 和W-li-SRACV等。W-li-SRACV(IEEESignalProcessingLetters，2012,19(3) :155-158)通过 设计阵列协方差矢量的加权1:范数最小化框架而具有很好的分辨率。实值稀疏表示算法 (IEEEAntennasandWirelessPropagationLetters，2013，12:376-379)降低了 计算量。互親 条件下，改进的 li_SVD( IEEEAntennasandWirelessPropagationLetters，2〇12，11:1 21〇-1213)角度估计方法被提出来，未知互耦影响得以消除。然而，当考虑互耦或者噪声为高斯 色噪声时，以上方法的估计性能下降。

【发明内容】

[0006] 本发明的目的在于提供一种互耦条件下基于四阶累积量稀疏表示的ΜΜ0雷达波 达方向估计方法。
[0007] 本发明的目的是这样实现的：
[0008] 互耦条件下基于四阶累积量稀疏表示的Μ頂0雷达波达方向估计方法，包括如下步 骤：
[0009] (1)发射阵列发射相互正交的相位编码信号，接收端进行匹配滤波处理后获得接 收数据，并利用发射和接收阵列都具有的互耦矩阵带状对称Toeplitz结构特点，通过线性 变换消除未知互耦的影响；
[0010] (2)构造降维转换矩阵，对消除互耦后的数据进行降维处理，进而基于新的数据矩 阵构造具有特殊形式的四阶累积量矩阵；
[0011] (3)对四阶累积量观测矩阵进行降维处理，获得稀疏表示框架下的相应模型，并利 用导向矢量和相应噪声子空间的正交性，设计权值矩阵以加强稀疏解；
[0012] (4)设计加权h范数约束最小化的稀疏表示框架，利用编程软件包S0C计算方法， 获得恢复矩阵，寻找恢复矩阵中的非零行，实现对高斯白噪声和色噪声条件下ΜΜ0雷达系 统中目标D0A的精确估计。
[0013] 步骤（1)中利用单基地ΜΜ0雷达互耦矩阵带状对称Toeplitz结构特点，通过线性 变换消除未知互耦的影响按如下步骤：
[0014] (1.1)在快拍下，接收数据为
[0015] x(t) =CAs(t)+n(t)
[0016] 其中<^ = €；0<：,.，(^和(^分别是发射和接收阵列的互耦矩阵，0:是&〇11扣1?^积操 作。 A=[a(01)，· · ·，a(0p)]且 0($):=七(|)?~(&)，ρ = 1，2, ..·Ρ，Ρ 是远场目标总数。 a,(^) = [l,e^，…，，("'1>8&'外和~(6^ = [1,，8^，…，分别是发射和接收导 向矢量，ΘΡ是第p个目标的D0A，M和N分别是发射和接收天线数;s(t) = [S1(t)，...，sP(t)]T 是非高斯信号，n(t)是零均值高斯白噪声或色噪声；
[0017] (1.2)消除未知互耦影响，构造选择矩阵J，新的数据矢量为
[0018]
[0019] 其中 J = <=[0^ 心 0財]，Λ=[0μ 心 表示 QXQ 维单位 矩阵，丑=鉍-2^，反=#- 2&，1^+1为非零互耦系数个数。吨)是新噪声，^ =[砂1)，...,封4)]， _P) =以⑴，先e'、nQ-]i"]τ，足的)=[1，0^^ 包含未知互耦信息的对角矩阵。
[0020] 所述的步骤(2)中构造降维转换矩阵，对消除互耦后的数据进行降维处理，进而基 于新的数据矩阵构造具有特殊形式的四阶累积量矩阵按如下步骤：
[0021] (2.1)对消除互耦后的数据矩阵进行降维处理，构造降维转换矩阵J3如下
[0022]
[0023]
[0024] (2.2)将降维转换矩阵J3与数据矢量y(t)相乘，新的数据矢量为
[0025]
[0026] 其中
(ΘΡ)]且^(。^，，^，…^^^^":^豆⑴是去除互耦和冗杂行后的新噪声向量； [0027] (2.3)构造具有特殊形式的四阶累积量矩阵：
[0028]
[0029] 其中心氏為)表示CV中第k2行且第kdU的元素，通过收集iKO的J个快拍，获得夕 的四阶累积量矩阵CV的样本估计A中每一项
[0030]
[0031] 其中.?是J的第h个兀尜，丨丨，/ = 1,2,3,4，在收集】个快拍条件下， J
^ /=i
[0032] 所述的步骤(3)中对四阶累积量观测矩阵进行降维处理，获得稀疏表示框架下的 相应模型，并利用导向矢量和相应噪声子空间的正交性，设计权值矩阵以加强稀疏解按如 下步骤：
[0033] (3.1)对四阶累积量观测矩阵进行降维处理：
[0034]
[0035] 其中Vs由对应P个大奇异值的右奇异向量组成，U和V分别由奇异值分解的左和右 奇异向量组成，:(斗. Dmc和Cs分别为包含奇异值，未知互親 和信号四阶累积量的对角矩阵，?ι: > ?2. 2…> Mm ;.Sfb中第i行第j列兀素为SFB(i，j)= [卩13(0」）]1|2[卩13(0」）]1，[卩13(0」）] 1是卩13(0」）第i个值。
[0036] (3.2)利用导向矢量和噪声子空间的正交性，设计权值矩阵加强稀疏解：
[0037]
[0038]其中r =(#；α/,.,，鳥么>]是稀疏表示框架下的完备字典，攝t 为所有感兴趣D0A的离散样本网格，L>>p，if包含完备字典&中可能对应目标真实DOAs 的P个导向矢量。％ 由馬中剩下导向矢量组成;Un为噪声子空间，由U中第 (p+i)到(双+#-1)列组成;r(w =(巧?/")(/:)潜在地包含阶p和; :.r(w的第i个值是w中 第i行2-范数，[默似1 =?? ιτ0.，:)ιμ
[0039] 所述步骤(4)中设计加权h范数约束最小化的稀疏表示框架，利用编程软件包SOC 计算方法，获得恢复矩阵，寻找恢复矩阵中的非零行，实现对高斯白噪声和色噪声条件下 Μ頂0雷达系统中目标D0A的精确估计按如下步骤：
[0040] 设计加权h范数约束最小化的稀疏表示框架，
[0041]
[0042] 其中4是匀真实值的样本估计;7Γ是LX1维向量，[0 =|丨^仏為，稀疏矩阵 巧满足g 正则化参数选择为在四阶累积量估计误差的方差均一化基础上，具有 (i? +反-1)Ρ自由度高概率l-ε置信区间的卡方分布的上限值，ε =0.001足够；利用公式 4/2/?ν(1-5,(及+及-1)乃计算ri;最后，通过S0C程序计算，测绘寻找Ρ个峰值获得目标 的DOAs。
[0043] 本发明的有益效果在于：
[0044] 1、本发明由于四阶累积量技术和加权稀疏表示框架的应用，成功地抑制了色噪 声，在高斯色噪声情况下，本发明提供精确的波达方向估计，比相似ESPRIT，F0C-MUSIC和 h-SVD方法具有更高的角度分辨率和更好的角度估计性能。
[0045] 2、本发明由于加权稀疏表示技术的应用，权值矩阵强化了稀疏解，无论噪声是高 斯白噪声还是色噪声，本发明估计性能均优于相似ESPRIT，FOC-MUSIC和h-SVD方法，具有 更低的SNR阈值，提供更高的角度分辨率；
[0046] 3、本发明由于降维转换技术的应用和所设计的四阶累积量矩阵的特殊形式，计算 量大大降低，本发明的计算复杂度比F0C-MUSIC方法更合理。
【附图说明】
[0047]图1本发明的整体框架图；
[0048] 图2本发明互耦K = 2和K = 3时，SNR = -10dB和SNR = 0dB不同情况下三个目标角度 估计的空间谱；
[0049]图3不同方法互耦白噪声时三个目标角度估计均方根误差和信噪比关系；
[0050]图4不同方法互耦色噪声时三个目标角度估计均方根误差和信噪比关系；
[0051]图5不同方法互耦白噪声时两个目标角度估计均方根误差和角度间隔关系；
[0052]图6不同方法互耦色噪声时两个目标角度估计均方根误差和角度间隔关系；
[0053]图7不同方法互耦白噪声时三个目标角度估计均方根误差和快拍数关系；
[0054]图8不同方法互耦色噪声时三个目标角度估计均方根误差和快拍数关系；
[0055]图9不同方法互耦白噪声时三个目标角度估计分辨概率和信噪比关系；
[0056]图10不同方法互耦色噪声时三个目标角度估计分辨概率和信噪比关系。
【具体实施方式】
[0057]下面结合附图对本发明做进一步描述。
[0058]本发明提供一种未知互耦误差条件下基于四阶累积量稀疏表示的单基地多输入 多输出（multiple-input multiple_output，MIMO)雷达波达方向（direction of arrival， 简称DOA)估计方法，主要是为了解决目前在稀疏表示框架下，单基地MMO雷达系统中的角 度估计方法存在不适用于高斯色噪声和互耦误差情况，估计精度不理想等缺点。首先通过 线性转换消除互耦影响，构造降维转换矩阵减少计算复杂度。然后利用四阶累积量矩阵的 有利形式，设计加权h范数约束最小化稀疏表示框架获得目标的D0A。其过程为:获得单基 地ΜΙΜΟ雷达系统的接收信号，利用发射和接收阵列都具有的互耦矩阵(MCM)的特殊结构，通 过线性变换消除未知互耦的影响;然后构造降维转换矩阵，对消除互耦后的数据进行降维， 进而基于新的数据矩阵构造四阶累积量观测矩阵，并进一步进行降维处理;获得稀疏表示 模型，利用导向矢量和噪声子空间的正交性，设计权值矩阵;最后设计加权h范数约束最小 化稀疏表示框架，获得恢复矩阵，寻找恢复矩阵中的非零行，实现对高斯白噪声或色噪声情 况下MM0雷达系统中目标D0A的精确估计。与互耦误差条件下单基地MM0雷达波达方向估 计传统的相似ESPRIT，FOC-MUSIC和h-SVD方法相比，本发明在高斯白噪声和色噪声条件下 都提供精确的D0A估计，同时具有更高的角度分辨率和更好的角度估计性能。
[0059]本发明的目的在于克服上述方法的缺陷，提供一种新的未知互耦误差条件下基于 四阶累积量稀疏表示的单基地ΜΜ0雷达波达方向估计方法。本发明利用发射和接收阵列都 具有的互耦矩阵(MCM)带状对称Toeplitz结构，消除未知互耦的影响;然后基于新的接收数 据，构造降维转换矩阵，进而构造具有特殊形式的四阶累积量矩阵以降低稀疏信号恢复的 计算复杂度;最后设计加权h范数约束最小化的稀疏表示框架，通过找到恢复矩阵中的非 零行获得目标的D0A。与未知互耦误差条件下单基地ΜΜ0雷达角度估计的传统方法相比，由 于利用了四阶累积量技术和加权稀疏表示框架，本发明在高斯白噪声和色噪声条件下都提 供精确的D0A估计，同时具有更高的角度分辨率和更好的角度估计性能。本发明波达方向估 计主要包括以下几个方面：
[0060] 1、根据互耦条件下单基地ΜΜ0雷达接收数据的结构和发射、接收阵列都具有的互 親矩阵带状对称Toeplitz结构特点，通过线性变换消除未知互親的影响。
[0061] 考虑一个窄带单基地ΜΜ0雷达系统，发射阵列和接收阵列分别由Μ和N个天线组 成，阵列均为半波长空间距的均匀线性阵列(ULAs)。在单基地ΜΙΜ0雷达中发射和接收阵列 紧密共置，因此对于一个远场目标可以看作它们具有相同的角度（即波达方向（D0A))。发射 阵列利用Μ个天线发射Μ个具有相同带宽和中心频率的不同窄带正交波，假设目标数目为Ρ， ΘΡ表示第ρ个目标的D0A。在有K = k+1个非零互耦系数和min{M，N}>2k的情况下，发射阵列 和接收阵列都考虑互耦的影响。于是，在某个快拍下接收端匹配滤波器的输出表示为
[0062] Λ：(?) = As(t) + n{t)(\)
[0063] 其中x(i) e €λ/λ_χ1，s(i)=卜丨⑴4(?),…七⑴f e CPxl是非高斯信号，《(〇 e C·1是 具有零均值的高斯白噪声或色噪声向量。和 = 分别是发射和接收导向矢量，p = 1，2，· · · P。 ⑷）? C>,.⑷…,CX吼）? C>,.(仏)]e C1Λ x"'，其中Ct和匕分别是发射和接收阵列 的互耦矩阵，互耦矩阵为带状对称Toeplitz矩阵并且表示为
[0064]
[0065] 其中(^(1 = ^;」=0，1，...，1〇是互耦系数，与两个天线间的距离有关，并且满足 0〈 | Cik|〈. . .〈 | Cil I〈 I Ci〇| =1。基于J_的结构特点和Kronecker积操作的特性，未知互親条件 下Μ頂0雷达的接收数据为
[0066]
[0067]其中 C = <C, ?C,.)€C·咖，J 是导向矩阵，并且 (义关于信号和噪声引入如下假设：
[0068] 1)信号是零均值平稳非高斯信号，彼此相互独立；
[0069] 2)噪声是零均值高斯噪声，可能为白噪声或空间相关色噪声；
[0070] 3)噪声和信号互相独立。
[0071] 稀疏表示方法由于未知互耦矩阵的存在而失效。为构造有效的四阶累积量稀疏表 示框架，首先需要消除发射和接收阵列的互耦影响。利用公式（2)中互耦矩阵带状对称 Τ 〇 e ρ 1 i t ζ结构特点，对于C t a t ( Θ p)，p = 1，2，. . .，Ρ，定义Μ X Μ维选择矩阵 心U以选择Ct的中间Μ -2Η亍。令z = Zsin(9，有
[0072]
[0073] 其中辦,）=作+ 1>"_(ζ? + ζ--3]是标量，并且承(4) = (1之:2…定义另 一个选择矩阵/2=[0w 4 #:=JV-类似公式(4)，则
[0074]
[0075] 其中
，因此进一步推导出
[0076]
[0077] X(t)的左边，得到
[0078]
[0079]其中y(t)是新接收数据矢量，对角矩阵D = diag[ I c%) I，· · ·，I c(0p) I ]包含未知 互耦的信息。^ =[如及.):] e 和岭)=i?(〇分别是新导向矩阵和噪声。由于互耦 矩阵C转变成为了对角矩阵，因此互耦不再对导向矩阵存在影响，即公式(7)中互耦误差影 响得到补偿。
[0080] 2、构造降维转换矩阵，对消除互耦后的数据进行降维处理，进而基于新的数据矩 阵构造具有特殊形式的四阶累积量矩阵。
[0081] 构造降维转换矩阵，对消除互耦后的数据进行降维处理，如下所示
[0082] 在消除未知互耦影响后，为了降低下面四阶累积量计算的复杂度，基于单基地 μ頂〇雷达的特殊结构构造降维转换。J的第ρ列4(1)?)元(4) e 为
[0083]
[0084] 可以发现中包含很多重复项，并且满足
[0085] α(θ.:) = 0){θρ) (9)
[0086] 其中G和b(0p)的具体表达式如下所示
[0087]
[0088] 其中 6巧），G e C_^-V，且 m = 1。重复项在下面的四阶累积量计算中将导致大量的冗杂信息，进而在稀疏信 号恢复中导致很大的计算量。为了解决这个问题，构造降维转换如下
[0089]
[0090] 因此，将降维转换矩阵J3与数据矢量y(t)相乘，有
[0091]
[0092] 其中5 = ⑷),和恥)分别是消除互耦影响和去除 冗杂行后的新数据向量和噪声。同时，F计算为
[0093]
[0094] 其中A = j。降维转换不仅将接收数据矩阵的维数从涵?X/显著减少 到(反+友-1)X J，而且避免了附加空间色噪声，J为快拍数。
[0095] 基于新的数据矩阵构造具有特殊形式的四阶累积量矩阵，如下所示
[0096] 在补偿了未知互耦误差影响和降低了接收数据维数之后，通过收集的J个快 拍，数据f的四阶累积量為，夂4)定义为
[0097]
[0098] 其中爲是f的第1^个元素，十」?-1，f=l,2,3,4。在J个快拍条件下， i', (m為)的估计中的每一项按如下公式获得
[0099] ' " J ?=ι
-
[0100] 根据公式(14)，设计四阶累积量矩阵Cy如下所示
[0101]
[0102] 其中CV(AA)表示?-,中第k2行且第1^列的元素。基于四阶累积量的特性和信号Sl (i = l，2,. . .，P)之间相互独立等假设，推导出
[0103]
p-i
[0104] 其中气(4)和, (6,)分别是Μθρ)中第k2和第ki个元素，鳥=0_：(1^,<，>^， >^):是信 号8[)的四阶累积量。构造四阶累积量矩阵$后，(#+友-1)2个值被合理地排列在?-中。由于 中包含1^和1?的所有可能的排列置换，推导出?.如下所示
[0105]
[0106] 其中Cf g C(A/+A-?μ〇υ+ΛΓ-11 〇Dmc = diag( | c(0〇 |4, · · ·，| c(0P) |4)和Cs = diag(Pi，· · ·， fr)为对角矩阵。同时，SFB中的第i行且第j列元素为
[0107] SFB(i,j)=|[Fb(0j)]i|2[Fb(0j)]i(19)
[0108] 其中丨_</外4"?|夕-丨），1彡」彡卩，[卩13(0」）1是卩13(00中的第1个值。
[0109] 3、对四阶累积量观测矩阵进行降维处理，获得稀疏表示框架下的相应模型，并利 用导向矢量和相应噪声子空间的正交性，设计权值矩阵以加强稀疏解。
[0110] 对四阶累积量观测矩阵进行降维处理，如下所示
[0111] 为了进一步降低D0A估计稀疏表示方法的计算复杂度，四阶累积量观测矩阵CV能 够从(# + # -1) X (及+万-1)维转换成为(# +友-1>χΡ维，如下所示
[0112]
[0113]其中包含绝大部分的信号能量，可以用来代替对DOAs进行精确的估 计。1由对应着大奇异值的右奇异向量组成，此外，，u和V分别由左奇异向量和右 奇异向量组成，4 = diagiq，,.. ·, ?λ7+λμ )是对角矩阵，且碑2…是奇异值。令 Γ = ，Dm。和Cs均为对角矩阵。当不考虑互耦时原始的恢复矩阵模型中的某些行是 零向量时，T中的对应行也必然为零，即使Dmc和*的组成是完全未知的。这表明0^， 的引入不会改变稀疏解，因此通过公式(7)的线性转换消除互耦影响后，在公式(20)现有的 形式中，稀疏表示方法能够成功地对波达方向进行估计。
[0114] 获得稀疏表示框架下的相应模型，如下所示
[0115] 为了应用稀疏表示方法估计D0A，令旧丨^为所有感兴趣DOAs的离散样本网格。可能 的D0A数目远远大于目标数目，即L>>P，进而为可能DOAs构造完备字典
>%满足A = ，当空域某点有真实目标存在 时，％中对应行为非零向量而其他为零向量。因此，％与原始矩阵T具有相同的行支持，我们 能够将D0A估计问题转换成为找到完备字典中对应％非零行的位置。为估计非零行的最小 数目，直接的稀疏度量是1〇_范数处罚，但1〇_范数最小化是非凸优化且NP困难而无法解决。 1:-范数处罚能够用来合理地解决这个问题。为获得稀疏矩阵巧，构造如下1:-范数约束最小 化框架
[0116]
[0117] 其中[d%(/,:)112，且?f) eCixl是正则化参数，通过收集J个快拍ξ是真实 值<^的样本估计，从公式(14) (15) (16)和(20)获得。
[0118] 利用导向矢量和相应噪声子空间的正交性，设计权值矩阵以加强稀疏解，如下所 示
[0119] 为了加强h-范数处罚以更好地接近1〇_范数处罚，我们提出如下加权h-范数最小 化稀疏表示方法。将完备字典&按列分为两部分，一部分包含可能对应目标真实DOAs的P个 导向矢量，另一部分由剩下的对应着恢复矩阵％中零行的导向矢量组成。因此，^能够表示 为為^啤"肩 21]，于是选择U中从(P+1)到(M + #-1)的共(凝+及-1-P).列作为噪声子空 间1^，^￡/"表示为
[0120]
[0121] 兵甲白J - 叮，匕泣骰址明Wi甲的母一项趋于零，因此基于此特性，设计一个加 权矩阵如下所示
[0122]
[0123]其中坏f2)(/)/max(?f:))比灰^)(y)/max(Wf:丨）小，且 0 < )(、/)/max()fd< 1， 同时当①时，七')/ 4 0，>(/)和％@(/)分别是iff >和iff1的第i项和第 j项潜在地包含两部分，故〇^?/,,)(/，替在地包含fff jPfff'Wr能够实现大权值处罚稀 疏向量中可能为零的元素，小权值保留值更大的项。换言之，加权矩阵的引入能够加强 稀疏解，提高DOA估计的精确度。
[0124] 4、设计加权h范数约束最小化的稀疏表示框架，利用编程软件包S0C(二阶锥)计 算方法，获得恢复矩阵，寻找恢复矩阵中的非零行，实现对高斯白噪声和色噪声条件下Μ頂0 雷达系统中目标D0A的精确估计。
[0125] 构造加权稀疏表示框架如下所示
[0126]
[0127 ]其中η是设置误差量的正则化参数，η数值的选择在最后的doa估计中起到很重要 的作用。令見7；为四阶累积量估计误差，为了获得vec^A)的分布特性，首先推 导vec(ACV)的分布，并且M：f. =4 -g，^是真实值G的样本估计，因此= 。基于 4的元素结构，可以认为vec(^)是一系列#的四阶累积量估计值。根据文献（Signals， SystemsandComputers，1993,2:1186-1190)，对于统计过程,，四阶累积量估计vecf,..)中所 有项.％(m先)具有(A:,，H λ:) H & )]满足渐近正态分布的特性。此外， 氏，&，&，&)] = c;%名名名），因此多快拍情况下，Vec(ACf)是具有零均值的渐近正
态分布:ΗΠ
[0128]
[0129] 其中AsN(y，V)表示均值为μ，协方差矩阵为的渐近正态分布。 并且根据文献可知四阶累积量估计误差的协方差如下所示
[0130]
[0131] 其中1 =让1，1?，1?，1^4，0 = 1^1'，1^'2，1?'，1^'4，缩写表示〇44 = 〇44(丁，0)，1112(入）=1112(1^1， 00. λ)，且m2(λ，）=m2 (h，，λ，）。根据=丨,吐、/COY. .j"V Η Σ cov 认(')， ξ.-^co if {I ,^?/i (〇 = J'i, (〇·· ,V/,; 4, ./；.(/) = V·^ (/)>', (/)···yK (/),n < 4 0kn = k2,ki2 = k3，kl3 = k4 ; k21 = k3，k22 = k2，k23 = k4 ; k31 = k4，k32 = k3，k33 = k2，同理k ' 具有相同定义。同时， 利用公式(15)获得k阶样本矩
[0132] 当让2 =让3 =让4且让'2 =让3'=让'4，￥由【':(.4/)=.。〇4以左1，是2為，灸 2)為(灸1為，是2，女2)}构 成，A来自于4。推导出i和h，匕的关系如下所示
[0133]
[0134] 其中[· ]rem表示值?/(?? + /?-1)的余数，[· ]ceii表示不小于值/7〇7 + V-1)的最 小的整数。MPkAk's的关系与公式(27)相同。由于·^二從#，利用矢量化操作的特性， 可以推导出
[0135] \
.(/.Z
[0136] 其中五v e C(nilPx(i?+〃-1):。基于公式(28)和vec(Aq冲渐近正态分布的线性转换 不变特性，推导出vec(Afy)中关于均值及和协方差矩阵f的结论如下所示
[0137]
[0138] 其中elC(m)i3x(A7*--μ。即vec(A(? f)满足均值为零且协方差矩阵为EvV(Ev)h的渐 近正态分布。因此，根据文献（IEEETransactionsonSignalProcessing，2005，53(8): 3010-3022)，本发明中正则化参数τι选择为在四阶累积量估计误差的方差均一化基础上，具有 (W + ,? -.丨自由度高概率l-ε置信区间的卡方分布的上限值，ε = 0.001足够。利用Mat lab 软件，公式?'Λ/_2/·/?ν，(! - 6%(Μ + /V -丨)/3)可以用来计算η。另一方面，公式(24)可以通过S0C(二 阶锥)程序软件包计算，例如SeDuMi和CVX。最后，通过解决公式(24)，测绘获得目标的 D0As〇
[0139] 实施例1
[0140]步骤一、建立未知互耦误差条件下单基地ΜΜ0雷达的接收信号模型，根据互耦矩 阵结构特点，利用线性变换消除未知互耦影响。
[0141] 考虑一个窄带单基地ΜΜ0雷达系统，发射阵列和接收阵列分别由Μ和N个天线组 成，阵列均为半波长空间距的均匀线性阵列。在单基地ΜΜ0雷达中发射和接收阵列紧密共 置，因此对于一个远场目标可以看作它们具有相同的角度（即波达方向（D0A))。发射阵列利 用Μ个天线发射Μ个具有相同带宽和中心频率的不同窄带正交波，假设目标数目为Ρ，θ ρ表示 第Ρ个目标的D0A。在有K = k+1个非零互耦系数和min{M，N}>2k的情况下，发射阵列和接收 阵列都考虑互耦的影响。于是，在某个快拍下接收端匹配滤波器的输出表示为
[0142] λλ(7) = As{t) + n(t) (30)
[0143] 其中，冲)=[咕),為(0,_..々(0]、€^1是非高斯信号，_)￡1：胍1是 具有零均值的高斯白噪声或色噪声向量。1_~] Τ和 ?及,）=[l,，in(V，…，(Λ?-分别是发射和接收导向矢量，p = 1，2，. . . Ρ。
_f其中Ct和Cr分别是发射和接收阵列 的互耦矩阵，互耦矩阵为带状对称Toeplitz矩阵并且表示为
[0144]
[0145] 其中cij(i = t，r; j = 0，l, · · ·，k)是互親系数，0〈|cik|〈. · .〈|cii|〈|ci〇| =1，与天线 间距离有关。基于j结构和Kronecker积特性，互親条件下ΜΙΜΟ雷达接收数据为
[0146]
[0147] 其中 C = (C;#Cr：)eCMVxMV，/? = [α网的)]eC麻ρ 是导向矩阵，并且 α(€) = ?,(&)?α,.(4)，Ρ = 152，···/>。关于信号和噪声作如下假设:信号是零均值平稳非高斯 信号，彼此相互独立;噪声是零均值高斯噪声，可能为白噪声或空间相关色噪声；噪声和信 号互相独立。
[0148] 稀疏表示方法由于未知互耦矩阵的存在而失效。为构造有效的四阶累积量稀疏表 示框架，首先需要消除发射和接收阵列的互耦影响。利用公式（31)中互耦矩阵带状对称 Toeplitz结构特点，对于Ctat(0p)，p = l，2, . . .，Ρ，定义狀χΜ维选择矩阵 Λ = ..? U以选择Ct的中间M = ikf - 2Α：行。令z.= ，.有
[0149]
[0150] 其中 9(4> = /[1 + 1>,(乂+7，]是标量，且运(0 = (1 z-r*-1)1。定义另一个选 /=1 择矩阵>/，=[0& 4 0&JeC#XjV，# 类似公式(33)，可以得到
[0151] - m、-;" -,.'-m …，
[0152] 其弓，因此进一步推导出 1=1
[0153]
= (\σρ)α{σρ)，μ = ΚΔ"·.，Γ
[0154] 其中c(0p)=Ct(0p)Cr(0p)，如0,) = 0，⑷"）?? (~Κ令，二义?々 eC丽遍，将J乘 到X(t)的左边，得到
[0155]
[0156]其中y(t)是新接收数据矢量，对角矩阵D = diag[ I c%) I，· · ·，I c(0p) I ]包含未知 互耦信息。i = R⑷),...，々'说.,)]e C 是新导向矩阵和噪声。互耦矩阵转变成了对角 矩阵，因此互耦不再对导向矩阵存在影响，公式(36)中互耦误差得到补偿。
[0157] 步骤二、构造降维转换矩阵，对消除互耦后的数据进行降维处理。
[0158] 在消除未知互耦影响后，为了降低下面四阶累积量计算的复杂度，基于单基地 Μ頂0雷达的特殊结构构造降维转换。】的第ρ列以&)衡,.(&) e 为
[0159]
[0160] 可以发现和4)中包含很多重复项，并且满足
[0161] 祕),,）= 61^,,)(38)
[0162] 其中G和b(0p)的具体表达式如下所示
[0163]
[0164] 其中%>C(i7+A^，, = e ,： -:1 :。重复项在下面的四阶累积量计算中将导致大量的冗杂信息，进而在稀疏信 号恢复中导致很大的计算量。为了解决这个问题，构造降维转换如下
[0165]
[0166] 因此，将降维转换矩阵J3与数据矢量y(t)相乘，得到新数据Κ?)如下
[0167]
[0168] U: 111 = [6(私…,A((V)] e C(A7+n)xi3，_:)是新噪声向量。同时，F计算为
[0169]
[0170] 其中J 。降维转换不仅将接收数据矩阵的维数从Wx J显著减少 到(M + # -1) X /，而且避免了附加空间色噪声，J为快拍数。
[0171] 步骤三、基于新的数据矩阵构造具有特殊形式的四阶累积量矩阵。
[0172] 在补偿了未知互耦误差影响和降低了接收数据维数之后，通过收集J，⑴的J个快 拍，数据，的四阶累积量定义为
[0173]
[0174] 其中氣,是，的第h个元尜，丨+ ，/ = 1,2,3,4收集J个快拍条件下， 卜(&，心n)的样本估计中的每一项按如下公式获得
[0175] ^ ?=1
[0176] 根据公式(43)，设计四阶累积量矩阵C,如下所示
[0177]
[0178] 其中表示.中第k2行且kM的元素。基于信号和噪声的假设，推导出
[0179]
[0180] 其中和分别是b(M中第k2和第lu个元素，急=cum(\,<,w;)是信 号8[)的四阶累积量。构造四阶累积量矩阵A后，(泣+#-1) 2个值被合理地排列在；中。由 于CV中包含1^和1?的所有可能的排列置换，推导出q如下所示
[0181]
[0182] 其中Cy e C(A'+Ar-1)x(iU+A-" .《Dmc = diag( | c(9〇 |4, · · ·，| c(9p) | 4)和Cs = diag(Pi，· · ·， fr)为对角矩阵。同时，SFB中的第i行且第j列元素为
[0183] SFB(i,j)=|[Fb(0j)]i|2[Fb(0j)]i(48)
[0184] 其中丨3/5^4_? +々-1)，1彡_]_彡卩，[卩13(00]1是卩13(0〇中的第1个值。
[0185] 步骤四、对四阶累积量观测矩阵进行降维处理。
[0186] 为了进一步降低下面D0A估计稀疏表示方法的计算复杂度，四阶累积量观测矩阵 G能够从(丑+及-1) X (及+没-1)维转换成为(M + W -1) X P维，如下所示
[0187]
[0188] 其中q =(：,_%包含绝大部分信号能量，可以用来代替&对D0A进行精确估计。^由 对应着大奇异值的右奇异向量组成，此外，(^ = tMFH，U和V分别由左和右奇异向量组成， /1 = diag(斗6?,…,)是对角矩阵，岣之…> 6^+以是奇异值。令Γ =凡￡ ， Dm4PCs均为对角矩阵。当不考虑互耦时原始恢复矩阵模型中的某些行是零向量时，T中的对 应行也必然为零，即使Dmc和的组成完全未知。这表明和Vs的引入不会改变稀疏 解，故通过公式(36)消除互耦后，在公式(49)现有形式中，稀疏表示方法能够成功地对D0A 进行估计。
[0189] 步骤五、获得稀疏表示框架下的相应模型，利用导向矢量和相应噪声子空间的正 交性，设计权值矩阵以加强稀疏解。
[0190] 为了应用稀疏表示方法估计D0A，令彳堯为所有感兴趣DOAs的离散样本网格。可能 的DOAs数目远远大于目标数目，即L>>P，进而为可能DOAs构造完备字典 馬=，啊4)，％4)，"彳(4)]?^_#-1>￡。令々满足之=馬7^当空域某点有真实目标存在 时，％中对应行为非零向量而其他为零向量。因此与原始矩阵T具有相同行支持，我们能 够将DOA估计问题转换成为找到完备字典中对应％非零行的位置。为估计非零行的最小数 目，直接的稀疏度量是1〇_范数处罚，但1〇_范数最小化是非凸优化且ΝΡ困难而无法解决。h-范数处罚能够用来合理地解决这个问题。为了加强h-范数处罚以更好地接近Ιο-范数处罚， 我们提出如下加权h-范数最小化稀疏表示方法。将完备字典g按列分为两部分，一部分包 含可能对应目标真实DOAs的P个导向矢量，另一部分由剩下的对应着恢复矩阵％中零行的 导向矢量组成。因此，先能够表示为毛=[與 1jf1]，于是选择U中从(p+n到(反+及-1)的 共(麗+ #-1-户)列作为噪声子空间1??表示为
[0191]
[0192]其中当时，已经被证明Wi中的每一项趋于零，因此基于此特性，设个权 值矩阵如下所示
[0193]
[0194] 其中％丨(/) / max(灰p)比 31 (./) / max(Wf =1)小，且 0 < )⑴ / max(d < 1， 同时当①时，(/) / max(?^'>) - 0，1 (Z_)和dy)分别是1和iff1的第i项和第 j项。為潜在地包含两部分，故(句％)(/; 1潜在地包含R(/: 1和2S实现了大权值处罚那些 稀疏向量Tfj中可能为零的元素，小权值保留值更大的项。换言之，加权矩阵的引入能够加 强稀疏解，提高DOA估计精确度。
[0195] 步骤六、设计加权1:范数约束最小化稀疏表示框架，获得恢复矩阵，寻找恢复矩阵 中的非零行，获得目标精确的DOA。
[0196] 构造加权稀疏表示框架如下所示
[0197]
[0198] 其中[0ll W，:)ll2，且#i ，η是正则化参数，通过收集J个快拍秦是真实 值^的样本估计，从公式(43)(44)(45)和(49)获得。τι是设置误差量的正则化参数，τι数值的 选择在最后的DOA估计中起到很重要的作用。
[0199] 令为四阶累积量估计误差，为了获得να(Δ^)的分布特性，我们首 先推导vec(ACf：)的分布，并且Μ；-.=之-C\.為是真实值；的样本估计，因此Δ& 。 基于 的元素结构，可以认为vec(G)是一系列J的四阶累积量估计值。根据文献(Signals， SystemsandComputers，1993，2: 1186-1190)，对于统计过程,， Λ，〇,)] = 6 (?為為為），且四阶累积量估计vec(g)中的所有项& Λ Λ)具 有力[弋(u,，κ)-_(a2，w\)]满足渐近正态分布的特性，因此多快拍情况下， VeC(Aq)是具有零均值的渐近正态分布，即
[0200] νβο(Λ^) ~ AsN(0,F) (53)
[0201] 其中AsN(y，V)表示均值为μ，协方差矩阵为Fe€@+#-1)2χ(?+^) 2的渐近正态分布。 并且根据文献可知四阶累积量估计误差的协方差如下所示
[0202]
[0203] 其中1 =让1，1?，1?，1^4，0 = 1^1'，1^'2，1?'，1^'4，缩写表示〇44 = 〇44(丁，0)，1112(入）=1112(1^1， 入八且肥以'卜肥仏'，）。根据
可以获得估 计值，并且乂 ⑴=(，卜· Λ (，)，A，S 4，⑴= ' ⑴X.少 k3，kl3 = k4 ; k21 = k3，k22 = k2，k23 = k4 ; k31 = k4，k32 = k3，k33 = k2，同理k ' 具有相同定义。同时， 利用公式(44)获得k阶样本矩％。
[0204] 当让2 =让3 =让4且让'2 =让3'=让'4，￥由1^〇*.，力.=切外'(々|為為，々:)，《(々|為 成，％来自于。推导出i和lu，1?的关系如下所示
[0205]
[0206] 其中[· ]rem表不值i./(Af+_iy-_1)的余数，[· ]ceil表不不小于值Ι'/:0^ + ?ν"-1)的最 小的整数。」和1沽'2的关系与公式(55)相同。由于八<^.=八&匕，利用矢量化操作的特性， 可以推导出
[0207]
[0208] 其中五v e 。基于公式(56)和vec(ACf.)中渐近正态分布的线性转换 不变特性，推导出vecWG)中关于均值#和协方差矩阵P的结论如下所示
[0209]
[0210] 其中f e 。即vec(A^V)满足均值为零且协方差矩阵为EvV(Ev)H的渐 近正态分布。因此，根据文献（IEEETransactionsonSignalProcessing，2005，53(8): 3010-3022)，本发明中正则化参数τι选择为在四阶累积量估计误差的方差均一化基础上，具有 (Μ +及-1)P自由度高概率l-ε置信区间的卡方分布的上限值，ε = Ο . 〇〇 1足够。利用Mat lab 软件，公式M/2/V;v(l - &(W - 1)P)可以用来计算η。另一方面，公式(52)可以通过SOC(二 阶锥)程序软件包计算，例如SeDuMi和CVX。最后，通过解决公式(52)，测绘If1获得目标的 D0As〇
[0211] 本发明的效果可通过以下内容说明：
[0212] ( - )运算复杂度分析：
[0213] 本发明的主要计算量是四阶累积量的计算和通过S0C程序包解决公式(52)获得恢 复矩阵的稀疏解。前者由于所设计的四阶累积量矩阵的特殊形式只需要 后者需要〇(/-?，其中U是离散样本网格总数，因此本发明的主要计算复杂度是 Ο[(及+ ?? - 1):: J + /」? />]"然而，在FOC-MUSIC方法中，由于所有的（lu，k2，k3，k4)在矩阵C x中都 需要被任意地排列置换，因而需要的计算量为^^^^ +斤及乂成^-巧:^降维转换 的应用和具有优势的四阶累积量矩阵形式使得计算量大大减少，故本发明的计算复杂度比 F0C-MUSIC方法更合理。
[0214](二)仿真条件与内容：
[0215] 我们呈现一些仿真结果来证实本发明的有效性和优点。用相似ESPRIT，FOC-MUSIC 和h-SVD方法与本发明进行对比。角度估计的均方根误差(RMSE)用来评估角度估计性能， 其定义为
[0216]
[0217] 其中是第i次MonteCarlo试验D0A真实值θρ的估计，Q是MonteCarlo试验总次数， 下面的仿真中选择Q = 500。考虑具有Μ个发射天线和N个接收天线的窄带单基地ΜΜ0雷达系 统，发射和接收阵列均为半波长空间距的均匀线性阵列。此外，仿真中的接收数据是混有高 斯噪声和阵元间互耦影响的二进制相移键控(BPSK)信号，大多数情况下，假设不相关目标 的数目Ρ是已知的，SNR定义为SNR=101 Ogl?(||CΛSM&/|l，l?),S'和N'是J快拍下信号和噪 声矩阵。此外，对于本发明以及F0C-MUSIC方法和h-SVD方法，离散样本网格均从-90°到90° 变换，均匀间隔为〇.〇5°。
[0218](三)仿真结果：
[0219] 1、本发明互耦K = 2和K = 3时，SNR = -10dB和SNR = 0dB不同情况下三个目标角度估 计的空间谱
[0220] 图二展示了本发明对于非零互耦系数K = 2和K = 3时，SNR=-10dB和SNR = 0dB不同 情况下的空间谱，其中M=N=7，J = 2000,三个不相关目标的DOAs分别为0: = -20°，θ2 = 〇°和 03 = 2〇°。当1( = 3时，发射阵列的非零互耦系数为[(^，(^1，(^2] = [1，0.6+」0.2,0.02+」0.1]， 接收阵列为[~0，(^1，(^2] = [1，0.5+]_0.3,0.01 + ]_0.2]。当1( = 2时，它们为[(^0，(^1] = [1， 0.0174+j0.0377]和[crQ，crl] = [l，0.0521-j0.1029]，并且将这种情况应用到下面的仿真 中。如图二所示，本发明的空间峰值很尖锐，旁瓣压制很低，即使在k增加，SNR降低的情况 下。这表明本发明能够提供极好的估计性能。
[0221] 2、不同方法互耦白噪声时三个目标角度估计均方根误差和信噪比关系
[0222]图三展示了高斯白噪声条件下，不同方法D0A估计的RMSE和SNR的关系，其中三个 不相关目标为θι = _20°，θ2 = 〇°和θ3 = 20°，同时M = N = 7, J = 4000。从图三可以看出，在低 SNR区域h-SVD方法比其他方法具有更好的性能。然而，本发明由于利用了所设计的四阶累 积量矩阵的特殊形式和加权稀疏表示框架，当SNR超过大约OdB时，本发明与相似ESPRIT， F0C-MUSIC和h-SVD方法相比，提供最好的角度估计性能。
[0223] 3、不同方法互耦色噪声时三个目标角度估计均方根误差和信噪比关系
[0224]图四展示了高斯色噪声情况下，D0A估计的RMSE和SNR的关系，其他条件与图三保 持一致。与li-SVD，相似ESPRIT和F0C-MUSIC方法相比，图四很明显说明了本发明的估计性 能在所有SNR区域都是最好的。因为利用了四阶累积量的特殊性质，本发明成功地抑制了色 噪声，然而在li-SVD和相似ESPRIT方法中，色噪声导致了性能的下降，尤其在低SNR区域。
[0225] 4、不同方法互耦白或色噪声时两个目标角度估计均方根误差和角度间隔关系
[0226] 图五和图六展示了RMSE和角度间隔的关系，其中M=N=7，J = 4000,图五中考虑高 斯白噪声以及SNR = 5dB，而图六考虑的是高斯色噪声且SNR = 0dB。每种方法中有两个不相 关目标，D0A分别为θ1 = 〇°和02 = 〇° + Δ Θ，并且Δ Θ从2°到14°变化。由图五可知，当Δ Θ在某个 值之内时，本发明的估计性能和相似ESPRIT方法相近，并且都比F0C-MUSIC和h-SVD方法性 能好。当角度间隔超过大约6°时，本发明在所有方法中性能最好。如图六所示，色噪声下本 发明对于空间相近目标提供最好的角度估计性能，即本发明具有最高的空间角度分辨率。
[0227] 5、不同方法互耦白或色噪声时三个目标角度估计均方根误差和快拍数关系
[0228] 图七和图八展示了不同方法D0A估计的RMSE和快拍数的关系，M = N = 7。图七中噪 声是高斯白噪声且SNR = 5dB，而图八为高斯色噪声且SNR=0dB。考虑三个具有不同D0A的不 相关目标为θι = _11.5°，θ2 = 〇°和θ3 = 11.5°。从图七和图八可以看出，对于高斯白噪声和色 噪声，本发明在所有快拍区域都提供很好的估计性能，此外，本发明相比于相似ESPRIT， FOC-MUSIC和1 i-SVD方法的估计优势随着J的增加而越来越突出。
[0229] 6、不同方法互耦白或色噪声时三个目标角度估计分辨概率和信噪比关系
[0230]图九和图十分别展示了高斯白噪声和高斯色噪声条件下，不同方法的目标分辨概 率和SNR的关系，其中M=N=7，J = 2000。三个不相关目标的DOAs分别为0: = -11°，θ2 = 〇°和 θ3 = 11°。此外，当三个目标所有D0A估计误差的绝对值都在0.1°之内时，目标被认为成功分 辨。如图所示，当SNR足够高时，所有方法提供100 %的目标分辨概率，然而，每种方法的分辨 概率在某一点开始下降，定义为SNR临界值。从图九和图十可以清楚看到，本发明与相似 ESPRIT，F0C-MUSIC和h-SVD方法相比具有最低的SNR临界值，同时它们表明无论噪声是高 斯白噪声还是色噪声，本发明都提供更高的目标分辨概率。
【主权项】
1. 互禪条件下基于四阶累积量稀疏表示的ΜΙΜΟ雷达波达方向估计方法，其特征在于： 包括如下步骤： (1) 发射阵列发射相互正交的相位编码信号，接收端进行匹配滤波处理后获得接收数 据，并利用发射和接收阵列都具有的互禪矩阵带状对称Toeplitz结构特点，通过线性变换 消除未知互禪的影响； (2) 构造降维转换矩阵，对消除互禪后的数据进行降维处理，进而基于新的数据矩阵构 造具有特殊形式的四阶累积量矩阵； (3) 对四阶累积量观测矩阵进行降维处理，获得稀疏表示框架下的相应模型，并利用导 向矢量和相应噪声子空间的正交性，设计权值矩阵W加强稀疏解； (4) 设计加权h范数约束最小化的稀疏表示框架，利用编程软件包S0C计算方法，获得恢 复矩阵，寻找恢复矩阵中的非零行，实现对高斯白噪声和色噪声条件下ΜΙΜΟ雷达系统中目 标D0A的精确估计。2. 根据权利要求1所述的互禪条件下基于四阶累积量稀疏表示的ΜΙΜΟ雷达波达方向估 计方法，其特征在于：步骤（1)中利用单基地ΜΙΜΟ雷达互禪矩阵带状对称Toepl itz结构特 点，通过线性变换消除未知互禪的影响按如下步骤： (1.1) 在快拍下，接收数据为 x(t)二CAs(t)+n(t) 其中^ = 和Cr分别是发射和接收阵列的互禪矩阵，?是Kronecker积操作dA二 [曰（目1)，. . .，a(0p)]且a巧，）= ",巧，），巧，），尸=1，2，...尸，P是远场目标总数。分别是发射和接收导 向矢量，θρ是第P个目标的D0A，M和N分别是发射和接收天线数;s(t) = ki(t)，...，sp(t)]T 是非高斯信号，n(t)是零均值高斯白噪声或色噪声； (1.2) 消除未知互禪影响，构造选择矩阵J，新的数据矢量为其中 J = J|0，/2，二[〇ι,χΑ 与壬 〇iw.]，J:'二[0'、χ/, /.\ 〇、J，Iq表示QXQ维单位矩 阵，府=Μ - 2/t，凉=iV - ，k+1为非零互禪系数个数。W)是新噪声，^ =脚6 )，...诚马,)]， a稱）=坏0,巧,）=[l，e曲虹0"，...，ew南-I仙叫，6.巧,）=p典 包含未知互禪信息的对角矩阵。3. 根据权利要求1所述的互禪条件下基于四阶累积量稀疏表示的ΜΙΜΟ雷达波达方向估 计方法，其特征在于:所述的步骤(2)中构造降维转换矩阵，对消除互禪后的数据进行降维 处理，进而基于新的数据矩阵构造具有特殊形式的四阶累积量矩阵按如下步骤： (2.1) 对消除互禪后的数据矩阵进行降维处理，构造降维转换矩阵J3如下其中巧=[与，片，…，每-ι]τ，且尤讯二恥X,"，/巧，〇趴(办。,_,)]，w = Q，l，."，原一1; (2.2) 将降维转换矩阵J3与数据矢量y(t)相乘，新的数据矢量为其中B=[b(0i)，...，b(0p)] i，而句是去除互禪和冗杂行后的新噪声向量； (2.3)构造具有特殊形式的四阶累积量矩阵：其中表示中第k2行且第ki列的元素，通过收集j：，的的J个快拍，获得来的四阶 累积量矩阵Cf的样本估计括f中每一项其中私是采的第ki个元素，1如(^M +南-1，/ = 1，2,3,4，在收集J个快拍条件下，4. 根据权利要求1所述的互禪条件下基于四阶累积量稀疏表示的ΜΙΜΟ雷达波达方向估 计方法，其特征在于:所述的步骤(3)中对四阶累积量观测矩阵进行降维处理，获得稀疏表 示框架下的相应模型，并利用导向矢量和相应噪声子空间的正交性，设计权值矩阵W加强 稀疏解按如下步骤： (3.1) 对四阶累积量观测矩阵进行降维处理：其中Vs由对应Ρ个大奇异值的右奇异向量组成，U和V分别由奇异值分解的左和右奇异向 量组成，; /1二化Ig(巧，…，似。、|)，Dmc和Cs分别为包含奇异值，未知互禪和信号 四阶累积量的对角矩阵，巧^扣;Sfb中第i行第j列元素为SFB(i，j)= I [化 (ej)]i|2[讯（目J)]1，[即（目J)]1是即（目J)第i个值。 (3.2) 利用导向矢量和噪声子空间的正交性，设计权值矩阵加强稀疏解：其中巧=(鸣。)Η[/。，馬=巧6(如...，&(每X]是稀疏表示框架下的完备字典，曲;;，为所 有感兴趣D0A的离散样本网格，L> >Ρ，钱i包含完备字典马中可能对应目标真实DOAs的Ρ个 导向矢量。眠2=(钱))Η。;,，帮)由馬中剩下导向矢量组成;Un为噪声子空间，由U中第(P+1) 到(麻+办-_1)列组成；^似三(钱巧,)似潜在地包含哺时的第i个值是W中第i行2_ 范数，[W 似],=11 ^(/，:化。5. 根据权利要求1所述的互禪条件下基于四阶累积量稀疏表示的ΜΙΜΟ雷达波达方向估 计方法，其特征在于:所述步骤(4)中设计加权h范数约束最小化的稀疏表示框架，利用编 程软件包S0C计算方法，获得恢复矩阵，寻找恢复矩阵中的非零行，实现对高斯白噪声和色 噪声条件下ΜΙΜΟ雷达系统中目标D0A的精确估计按如下步骤： 设计加权h范数约束最小化的稀疏表示框架，其中ξ是真实值的样本估计;臀堪LXl维向量，峨勺,二?Ι写化:化，稀疏矩阵攝满足 =玄决；η正则化参数选择为在四阶累积量估计误差的方差均一化基础上，具有 (府+历-1)户自由度高概率1 -ε置信区间的卡方分布的上限值，ε = 0.001足够；利用公式 C月这im(l-6，（Ι/+本U Ρ计算η;最后，通过SOC程序计算，测绘巧1寻找Ρ个峰值获得目标 的DOAs。
【文档编号】G01S7/02GK105974366SQ201610280200
【公开日】2016年9月28日
【申请日】2016年4月29日
【发明人】周卫东, 刘婧, 朱鹏翔, 宫文贺, 刘可
【申请人】哈尔滨工程大学