基于稀疏分解理论的声矢量阵宽带测向方法
【专利摘要】本发明公开了基于稀疏分解理论的声矢量阵宽带测向方法，它涉及水声矢量传感器阵探测领域中针对目标宽带连续谱噪声的一种测向技术及实现方法。本发明将稀疏化思想引入声矢量阵测向系统，根据宽带连续谱噪声信号带宽内各频点分量建立一种联合稀疏约束，最终得到统一的空域稀疏分解形式，从而实现对宽带连续谱噪声的测向。本方法能够形成较为尖锐的谱峰和幅度较低的噪底；在阵列快拍数较低时性能稳定，且能分辨相干信号源；在目标来波方向接近阵列轴向时，也不存在常规的波束形成(CBF)法和最小方差无失真响应(MVDR)法在低信噪比下出现的噪底起伏现象。
【专利说明】
基于稀疏分解理论的声矢量阵宽带测向方法
技术领域
[0001] 本发明属于水声信号处理领域，特别涉及一种基于稀疏分解理论的声矢量阵宽带 连续谱噪声测向方法，可用于被动水声矢量传感器阵列探测目标辐射的宽带连续谱噪声信 号。
【背景技术】
[0002] 对水面水下各类目标自身辐射的噪声包括窄带线谱噪声和宽带连续谱噪声。在被 动声纳窄带线谱探测方面，传统方法主要通过空域的常规波束合成(CBF)处理获得空间增 益，从而完成测向、检测等任务。而对于宽带连续谱噪声信号，CBF方法不能直接使用，通常 的做法是将宽带信号拆分成窄带信号，再利用CBF进行测向。这种处理思路没有充分利用宽 带信号的信息，它的方位分辨能力受瑞利限限制，不能有效分辨位于主瓣内的两个目标。最 小方差无失真响应(MVDR)法在阵列快拍数较低时的目标方位估计偏差较大。上述两种方法 在低信噪比下会出现噪底起伏现象。
[0003] 基于稀疏分解理论的空间谱估计是空间谱估计的另一个方向。信号的稀疏分解理 论于上世纪90年代提出，广泛应用于信号压缩、识别、提取等领域，其特有的信号稀疏性理 论适用于信号在空间内的分布情况，因此，信号的稀疏分解理论及其算法为声矢量阵列的 空间谱估计提供了一个新的解决思路和方向。
【发明内容】
[0004] 本发明所要解决的技术问题在于提供一种基于稀疏分解理论的声矢量阵宽带测 向方法。本发明将稀疏化思想引入声矢量阵测向系统中，并基于多重分片（Multiple Measurement Slice)的正交匹配追踪(MMS-0MP)算法，将宽带信号带宽内各频点分量建立 一种联合稀疏约束，最终得到统一的空域稀疏分解形式，从而实现对宽带连续谱噪声信号 的测向。
[0005] 本发明所要解决的技术问题是由以下技术方案实现的：
[0006] 基于稀疏分解理论的声矢量阵宽带测向方法，采用分片稀疏表示方法，实现对声 矢量阵接收的数据分片稀疏表示，通过联合稀疏约束求解该稀疏表示模型，对求解得到的 稀疏矩阵进行空间谱计算，从而实现目标信号的测向。具体包括以下步骤：
[0007] (1)将声矢量阵在预设时间段内的接收信号进行频域建模，得到接收信号在多个 预定频率处的频域快拍；
[0008] (2)计算每个频域快拍对应的协方差矩阵；
[0009] (3)根据空域扫描范围和宽带连续谱噪声的带宽建立过完备原子库，并对协方差 矩阵进行稀疏表示；
[0010] (4)将稀疏表示后的所有协方差矩阵的同一列按频率前后顺序堆叠形成三维张量 协方差矩阵，同时，过完备原子库也按照频率前后顺序堆叠形成三维张量过完备原子库后， 将三维张量协方差矩阵用三维张量过完备原子库和三维系数矩阵来表示;所述的三维系数 矩阵为未知量；
[0011] (5)分别初始化信号支撑集合、已选原子集合和迭代次数；
[0012] (6)计算过完备原子库中每个原子与当前残差信号的内积，利用准则公式对内积 进行处理得到过完备原子库中每个原子与当前残差信号的近似程度值;所述的残差信号的 初始值为三维张量协方差矩阵；
[0013] ⑴选择最大的近似程度值，根据最大的近似程度值找出对应原子的索引值，并更 新信号支撑集合和已选原子集合；
[0014] (8)利用最小二乘法计算更新后的已选原子集合在当前残差信号方向上的正交投 影后，根据正交投影和当前残差信号计算下一次残差信号；
[0015] (9)判定更新后的信号支撑集合个数是否大于迭代次数或者当前残差信号的2-范 数是否小于预设阈值，如果更新后的信号支撑集合个数大于迭代次数或者当前残差信号的 二范数小于预设阈值，则执行步骤（10);否则，将下一次残差信号更新为当前残差信号，跳 转到步骤(6);
[0016] (10)根据最终的正交投影和最终的已选原子集合计算得到三维系数矩阵;根据最 终的信号支撑集合和宽带连续谱噪声的来波方向的对应关系，计算三维系数矩阵的F-范 数，并根据F-范数得到宽带连续谱噪声的空间谱估计；
[0017] (11)对空间谱估计进行谱峰搜索得到前K个空间角度值;所述的K为迭代次数。
[0018] 其中，步骤(3)所述的根据空域扫描范围和宽带连续谱噪声的带宽建立过完备原 子库，具体包括步骤：
[0019] (301)将协方差矩阵中的每一列表示成
[0021] 其中，Rx(fm)为协方差矩阵;fm为频率，且匕已[心，5]，1为宽带连续谱噪声带宽内 的频点数;A(fm)为导向矢量矩阵，其每一列对应一个来波方向；s(fm)为宽带连续谱噪声;A h (fm)为导向矢量矩阵的共辄转置;rn(fm)为协方差矩阵中的某一列，n=l，2,…，N，N为声矢 量阵的阵元数；
[0022] (302)将-180度~180度的空间进行均匀划分，构造过完备化的阵列流型矩阵，满 足
[0023] rn(fm) =G(fm)fn(fm)；
[0024] 其中，fn(fm)为系数矩阵，其非零行对应真实的信号来波方向，其它非真实来波方 向的数据均为〇;G(f m)为过完备化的阵列流型矩阵即过完备原子库。
[0025] 其中，所述的步骤(4)具体包括以下步骤：
[0026] (401)将稀疏表示后的不同频率的协方差矩阵，按频率前后顺序堆叠成一个三维 张量协方差矩阵；
[0027] (402)将协方差矩阵中的每一列放入稀疏表不后的协方差矩阵对应的位置，同时， 将过完备原子库和系数矩阵也按频率前后顺序堆叠存放形成三维张量过完备原子库与三 维系数矩阵，并且满足
[0028] R=G?F
[0029] 其中，R为三维张量协方差矩阵；G为三维张量过完备原子库；F为三维系数矩阵， 它的每一水平切片都具有相同的稀疏结构。
[0030]其中，步骤(6)所述的利用准则公式对内积进行处理得到过完备原子库中每个原 子与当前残差信号的近似程度值;具体为：
[0031 ]第k次迭代时，原子gl与当前残差信号股的近似程度值为：
[0033] 其中，Ai,k为近似程度值;N为声矢量阵的阵元数;Μ为宽带连续谱噪声带宽内的频 点数；Μ · I U为2-范数;gi(fm)为三维张量过完备原子库G中对应角度的垂直分片中对应 频率f m的列向量;为当前残差信号股^1>中第n个垂直切分片中对应频率fm的列向 量;i为三维张量过完备原子库中原子的索引值，iel，I为过完备原子库的原子的索引值个 数。
[0034] 其中，所述的步骤(7)具体为:挑选符合条件的原子g补充已选原子集合，选择条件 为g = afgg. max(4,4);其中，Q(k)为已选原子集合；
[0035]更新后的信号支撑集合Ω (k)= Ω U {ik}，更新后的已选原子集合 QW ;其中，ik为第k次迭代时，三维张量过完备原子库中原子的索引值;8;(:为 第k次迭代时，三维张量过完备原子库的第ik列数据。
[0036]其中，所述的步骤(8)具体包括以下步骤：
[0037] (801)利用最小二乘法计算更新后的已选原子集合在当前残差信号方向上的正交 投影矩阵，计算公式为：
[0038] pq(/)=Q/Q；
[0039] 其中，PQ(f)为正交投影矩阵;Qf为更新后的已选原子集合，Q f=[gu gl2 ···];<?}为 更新后的已选原子集合的伪逆；
[0040] (802)根据正交投影矩阵和当前残差信号，计算下一次残差信号，计算方法为：
[0041] RLk,(fJ-Rr,(fm)(l-PQ(f)) ；
[0042] 其中，Rti/,,)为当前残差信号股中第X个水平切分片中对应频率fm的列向 量，m=l,2,…,Μ; I为单位矩阵。
[0043] 其中，步骤（10)所述的根据最终的信号支撑集合和宽带连续谱噪声来波方向的对 应关系，计算三维系数矩阵的F-范数，并根据F-范数得到宽带连续谱噪声的空间谱估计;具 体包括以下步骤：
[0044] (111)将三维系数矩阵按三维张量过完备原子库中原子的索引值划分成列片结 构；
[0045] (112)根据列片结构划分，计算三维系数矩阵的功率谱估计;计算公式为：
[0046]
[0047]其中，为三维系数矩阵的功率谱估计；i为三维张量过完备原子库中原 子的索引值；0,为索引值i对应的来波方向，当I个原子组成的过完备原子库对应的来波方 向在整个空域均勾分布时：；I I · I |f表不F-范数;Fe(i)为每一列片结构的 数据矩阵。
[0048] 本发明与现有技术相比具有如下优点：
[0049] 1、本发明与传统方法相比，具有尖锐的谱峰，具有更低的旁瓣水平，更有利于信号 的分辨和检测。
[0050] 2、本发明单次快拍即可完成对目标辐射的宽带连续谱信号进行快速测向，测向精 度高。
[0051] 3、本发明能分辨相干信号源，在目标来波方向接近阵列轴向时，也不存在传统常 规方法的噪底起伏现象。
【附图说明】
[0052]图1是本发明的声矢量阵列接收信号原理图；
[0053] 图2是本发明的声矢量阵宽带测向方法流程图；
[0054] 图3是本发明的三维张量协方差矩阵R、三维系数矩阵E和三维G的构造方式示意 图；
[0055] 图4是本发明三维张量过完备原子库(G的列片划分示意图；
[0056] 图5是本发明三维张量协方差矩阵徽的稀疏分解示意图；
[0057] 图6是本发明三维张量协方差矩阵:K的列片分解示意图；
[0058] 图7是本发明三维系数矩阵F的稀疏分片结构示意图；
[0059]图8是使用本发明方法声矢量阵宽带处理结果；
[0060]图9是本发明方法声矢量阵宽带处理与其它方法的处理结果对比；
[00611图10是本发明方法估计的角度方向的RMSE (均方根误差)与其它方法的对比。
【具体实施方式】
[0062]下面结合附图，对本发明的技术方案和效果作进一步详细说明。
[0063] 由于舰艇辐射噪声的宽带特性，在将稀疏分解理论应用到矢量水听器阵列测向 中，必然会产生宽带信号的稀疏分解问题。宽带信号频率能量分布的不均匀性，将导致频带 内各频点的空域稀疏特性不相同，若对各频点进行独立的窄带空域稀疏分解，再将分解结 果简单相加，不仅损失了信噪比，也不能处理相干信号。为了得到带宽内统一的空域稀疏分 解，需要一种方式能够将各频点的稀疏结构进行综合。多重分片处理是能够对多个平行子 分解进行联合统一约束的一种处理方式。基于MMS-0MP算法将窄带信号的稀疏分解扩展为 多个窄带信号的联合统一稀疏分解。通过MMS-0MP算法，多个窄带信号可以得到一个统一的 稀疏表示形式，从而得到宽带信号的稀疏分解。在基于MMS-0MP的矢量阵宽带测向算法中， 根据宽带信号带宽内各频点分量建立带宽内的联合约束，最终得到统一的空域稀疏分解形 式，从而确定宽带信号的来波方向。
[0064] 图1是本发明的声矢量阵列接收信号原理图。
[0065]图中声矢量阵列共有N个阵元，按照等间距排列，阵元间距为d，N个阵元对应有N个 接收通道。每个接收通道接收到水声宽带连续谱噪声信号以后经过滤波放大调理和AD变换 之后，输出一路数字信号以供后续处理。
[0066]本发明的基于稀疏分解理论的声矢量阵宽带测向方法，具体流程如图2所示，实现 步骤如下：
[0067]步骤1:参考图1，对声矢量阵在预设时间段内的接收信号进行频域模型建模。
[0068] (101)考虑由K个目标发出或反射的宽带连续谱噪声信号以波速c经水声信道传 播，在目标距离声矢量阵足够远的条件下，宽带连续谱噪声信号近似为远场平面波，并以来 波方向9i(i = l，2，. . .，K)入射至声矢量阵。若上述宽带连续谱噪声信号分布在频率fi-fM 之间，则声矢量阵在时间段T内的接收信号在频域的数学模型可表示为：
[0069] X=[x(fi) x(f2) ... x(fM)]
[0070] 其中，x(fm)(m=l，2, . . .，M)为声矢量阵在频率匕处的频域快拍。
[0071] (102)在基于扩展法的声矢量阵处理中，将x(fm)表示成如下形式：
[0073]其中，s(fm)为声矢量阵接收的目标声压信号在频率fm处的分量;a(0，f m)为目标在 频率fm处对应的导向矢量，对于第i个目标，有a(H) = [l ... ，其中 = = sin0 1]τ为声矢量传感器方向向量；符号發为 Kronecker积;N(fm) = [ni(fm) n2(fm) ... n3N(fm)]T为空间加性噪声矢量。
[0074] (103)将上式改写为矩阵乘积的形式为：
[0075] x(fm) =A(fm)S(fm)+N(fm)
[0076] 其中，A(./:J=[a(A./:J0u(q) a(A，./;")?u(6y …a的丄 量矩阵;S(fm) = [S1(fm) s2(fm) ... sK(fm)]T为信号在频率fm处的分量。
[0077] 步骤2:计算频率^对应的频域快拍的协方差矩阵。
[0078] x(fm)的协方差矩阵Rx(fm)为：
[0079] Rx(fm) =E[x(fm)xH(fm)]
[0080] =A(fm)E[S(fm)SH(fm)]AH(fm)+Rn(fm)
[0081] =A(fm)Rs(fm)AH(fm)+Rn(fm)
[0082] 由前述数学模型可知，矩阵A(fm)中含有信号在频率匕处分量的来波方向信息，对 宽带信号的来波方向估计即是对宽带信号所在的各频点处分量来波方向的综合估计。
[0083] 步骤3:根据空域扫描范围和宽带连续谱噪声的带宽建立过完备原子库G(fm)，并 对协方差矩阵进行稀疏表示。
[0084]建立过完备原子库G(fm)的过程：
[0085] (301)首先，在忽略声矢量阵列接收噪声的基础上，将Rx(fm)中的每一列表示成形 式：
[0086] rn(fm) =A(fm)Rx(fm)AH(fm)
[0087] =A(fm)s(fm)
[0088] 其中，A(fm)为导向矢量矩阵，它中的每一列对应一个来波方向。
[0089] (302)将{θ〇。~θ18〇。}角度空间完备化，构造一个完备化的阵列流型矩阵G(f m)(即 过完备原子库），满足
[0090] rn(fm)=G(fm)fn(fm)n=l,2,··· ,Ν
[0091] 其中，fn(fm)中非零行对应真实的信号来波方向，其它非真实来波方向的数据均为 0〇
[0092]在尽可能充分利用信号带宽的前提下，对宽带信号的稀疏表示需在频带[fifM]内 进行。
[0093]步骤4:对稀疏表示后的协方差矩阵Rx(fm)和过完备原子库分别按频率前后顺序堆 叠形成三维张量协方差矩阵和三维张量过完备原子库后，将三维张量协方差矩阵用三维张 量过完备原子库和三维系数矩阵来表示。
[0094] (401)将稀疏表示后的不同频率的协方差矩阵Rx(fm)，f me [hfM]按频率顺序堆叠 成一个三维张量协方差矩阵E ;
[0095] (402)将化(匕)放入Rx(fm)对应的位置时，过完备原子库G(f m)和系数矩阵f(fm)也 按频率顺序堆叠存放形成三维张量过完备原子库G与三维系数矩阵F，并且满足
[0096]
[0097] 其中，三维系数矩阵F为未知量，它的每一水平切片F(fm)都具有相同的稀疏结构。 [0098]参考图3、4完成数据构造，参考图5、6对三维张量协方差矩阵进行稀疏表示，建立 分片稀疏表不模型。
[0099]步骤5:初始化参数设置，包括信号支撑集合、已选原子集合和迭代次数。
[0100]初始化信号支撑集合Ω = 0，已选原子集合Q = 0，迭代次数K等参数。
[0101] 步骤6:寻找与当前残差信号最为匹配的原子:计算过完备原子库G(fm)中每个原 子与当前残差信号趿 (^的内积（股@=?)，并利用准则公式得到过完备原子库G(fm)中每 个原子与当前残差信号的近似程度值A 1>k。
[0102] 设三维张量过完备原子库为第k次迭代时，定义原子gl与当前残差 信号K<M)的近似程度值为：
[0104] 其中，| | · | |2为2-范数;N为阵元数;Μ为带宽内频点数;Rt'/,,)为当前残差信号 K(m中第n个垂直切分片中对应频率fm的列向量; gi (fm)为三维张量过完备原子库&中对应 角度的垂直分片中对应频率f m的列向量。
[0105] 步骤7:选择最大值对应原子§4的索引值，根据索引值更新信号支撑集合和已 选原子集合。
[0106] 挑选符合条件的原子g补充原子集，选择条件为：
[0108] 更新信号支撑集合Ω ?= Ω m U {ik}，更新已选原子集合其 中，ik为第k次迭代时，三维张量过完备原子库中原子的索引值泳为第k次迭代时，三维张 量过完备原子库的第ik列数据。
[0109] 步骤8:通过最小二乘法计算更新残差信号:利用最小二乘法计算更新后的已选原 子集合在当前残差信号方向上的正交投影PQ(f)，进而重新组合得到残差矩阵股(i)，完成更新 残差矩阵解算。
[0110] (801)利用最小二乘法计算更新后的已选原子集合在当前残差信号方向上的正交 投影矩阵PQW，计算公式为 [mu] Pq(/)=Q/Q}
[0112]其中，Qf为已选原子集合Qf=[gii gi2 ...]。
[0113] (802)根据正交投影矩阵和当前残差信号，计算下一次残差信号Rf(/")，计算方 法为：
[0115] 其中，(/",)为当前残差信号政夂υ中第X个水平切分片中对应频率frn的列向量， m=l，2,…，M;I为单位矩阵。
[0116] 步骤9:判定是否终止迭代:根据更新后的信号支撑集合个数是否满足设定的迭代 次数K或者通过判断当前残差信号的2-范数是否小于预设阈值，来判定是否终止步骤(6)~ (8)的迭代。
[0117]步骤10:根据最终的正交投影和最终的已选原子集合计算得到三维系数矩阵;根 据最终的信号支撑集合和宽带连续谱噪声的来波方向的对应关系，计算三维系数矩阵的 Frobenius范数(简称F-范数），并根据F-范数得到宽带连续谱噪声的空间谱估计；
[0118] (111)由于三维系数矩阵1F的每一水平分片F(fm)都具有相同的行稀疏结构，因此 将三维系数矩阵F按三维张量过完备原子库G中原子索引值i划分成列片结构，每一列片 数据矩阵表示为内（1);
[0119] (112)根据列片划分，计算三维系数矩阵F的功率谱估计，计算公式为：
[0120] p_-。-涓)-|聊|;
[0121] 其中，I I · I |F表示F-范数;01为索引值i对应的来波方向。当I个原子组成的过完备 原子库对应的来波方向在整个空域均匀分布时，
[0122] 参考图7,计算阵列接收数据的空间谱。经过K次分解后得到三维系数矩阵F，利用 原子集索引值与噪声信号来波方向具有对应关系，通过计算相应的原子集索引对应系数矩 阵的F-范数得到空间谱估计P? S-?P (Θi)。
[0123] 步骤11:对空间谱估计进行谱峰搜索得到前K个空间角度值;所述的K为迭代次数。
[0124] 本发明的效果可以通过以下仿真说明：
[0125] 1.仿真条件与方法
[0126] 设有一 16元声矢量阵接收带宽分别为150Hz~170Hz和180Hz~200Hz的两个目标 信号，以阵列轴向为0°，其来波方向分别为78°和-146°，接收机带宽为140Hz~210Hz，带内 信噪比为OdB。
[0127] 2.仿真内容与结果
[0128] 图8是将迭代次数限制为2次的稀疏分解结果，迭代次数等价于已知宽带信号包含 的目标数，当目标数未知时，可采用残差能量阈值来迭代次数，但是有可能会增加不必要的 计算量。
[0129] 在不改变前述仿真参数的条件下，将本文提出的基于MMS-0MP的矢量阵宽带测向 算法与常用的宽带非相干MVDR、非相干CBF和STMV算法进行比较，可得如图9所示的结果。本 发明的方法能够形成较为尖锐的谱峰和幅度较低的噪底，这也是基于稀疏分解的阵列测向 算法的一个特点，在目标来波方向接近阵列轴向时，也不存在上述三种算法的噪底起伏现 象。
[0130] 不改变其余仿真参数，对本文提出的算法在不同带内信噪比下的测向误差进行仿 真，可得如图10所示的结果。在不同带内信噪比下，本发明的矢量阵宽带测向方法的测向均 方根误差与STMV算法相近，其相对较窄的目标谱峰是与常规空间谱估计方法相比的优势所 在。综上所述，本发明的声矢量阵宽带信号测向方法能够形成较为尖锐的谱峰和幅度较低 的噪底。相比传统方法增强了目标信号的检测和分辨能力，在阵列快拍数较低时性能较为 稳定。在目标来波方向接近阵列轴向时，有效解决了上述传统方法的噪底起伏现象。
【主权项】
1. 基于稀疏分解理论的声矢量阵宽带测向方法，采用空域稀疏分解理论，实现对水下 目标福射的宽带连续谱噪声的测向，其特征在于:包括W下步骤： (1) 将声矢量阵在预设时间段内的接收信号进行频域建模，得到接收信号在多个预定 频率处的频域快拍； (2) 计算每个频域快拍对应的协方差矩阵； (3) 根据空域扫描范围和宽带连续谱噪声的带宽建立过完备原子库，并对协方差矩阵 进行稀疏表示； (4) 将稀疏表示后的所有协方差矩阵的同一列按频率前后顺序堆叠形成Ξ维张量协方 差矩阵，同时，过完备原子库也按照频率前后顺序堆叠形成Ξ维张量过完备原子库后，将Ξ 维张量协方差矩阵用Ξ维张量过完备原子库和Ξ维系数矩阵来表示;所述的Ξ维系数矩阵 为未知量； (5) 分别初始化信号支撑集合、已选原子集合和迭代次数； (6) 计算过完备原子库中每个原子与当前残差信号的内积，利用准则公式对内积进行 处理得到过完备原子库中每个原子与当前残差信号的近似程度值;所述的残差信号的初始 值为Ξ维张量协方差矩阵； (7) 选择最大的近似程度值，根据最大的近似程度值找出对应原子的索引值，并更新信 号支撑集合和已选原子集合； (8) 利用最小二乘法计算更新后的已选原子集合在当前残差信号方向上的正交投影 后，根据正交投影和当前残差信号计算下一次残差信号； (9) 判定更新后的信号支撑集合个数是否大于迭代次数或者当前残差信号的2-范数是 否小于预设阔值，如果更新后的信号支撑集合个数大于迭代次数或者当前残差信号的二范 数小于预设阔值，则执行步骤（10);否则，将下一次残差信号更新为当前残差信号，跳转到 步骤(6); (10) 根据最终的正交投影和最终的已选原子集合计算得到Ξ维系数矩阵;根据最终的 信号支撑集合和宽带连续谱噪声的来波方向的对应关系，计算Ξ维系数矩阵的F-范数，并 根据F-范数得到宽带连续谱噪声的空间谱估计； (11) 对空间谱估计进行谱峰捜索得到前K个空间角度值;所述的K为迭代次数。2. 根据权利要求1所述的基于稀疏分解理论的声矢量阵宽带测向方法，其特征在于:步 骤(3)所述的根据空域扫描范围和宽带连续谱噪声的带宽建立过完备原子库，具体包括步 骤： (301) 将协方差矩阵中的每一列表示成其中，Rx(fm)为协方差矩阵；fm为频率，且fme[fl，fM]，M为宽带连续谱噪声带宽内的频 点数;A(fm)为导向矢量矩阵，其每一列对应一个来波方向；S(fm)为宽带连续谱噪声;AH(fm) 为导向矢量矩阵的共辆转置;rn(fm)为协方差矩阵中的某一列，n = l，2,…，N，N为声矢量阵 的阵元数； (302) 将-180度~180度的空间进行均匀划分，构造过完备化的阵列流型矩阵，满足 rn(fm) =G(fm)fn(fm)； 其中，fn(f")为系数矩阵，其非零行对应真实的信号来波方向，其它非真实来波方向的 数据均为〇;G(fm)为过完备化的阵列流型矩阵即过完备原子库。3. 根据权利要求2所述的基于稀疏分解理论的声矢量阵宽带测向方法，其特征在于:所 述的步骤(4)具体包括W下步骤： (401) 将稀疏表示后的不同频率的协方差矩阵，按频率前后顺序堆叠成一个Ξ维张量 协方差矩阵； (402) 将协方差矩阵中的每一列放入稀疏表示后的协方差矩阵对应的位置，同时，将过 完备原子库和系数矩阵也按频率前后顺序堆叠存放形成Ξ维张量过完备原子库与Ξ维系 数矩阵，并且满足 股=枝.F 其中，股为Ξ维张量协方差矩阵；G为Ξ维张量过完备原子库;F为Ξ维系数矩阵，它的 每一水平切片都具有相同的稀疏结构。4. 根据权利要求3所述的基于稀疏分解理论的声矢量阵宽带测向方法，其特征在于:步 骤(6)所述的利用准则公式对内积进行处理得到过完备原子库中每个原子与当前残差信号 的近似程度值;具体为： 第k次迭代时，原子gi与当前残差信号盼的近似程度值为：其中，A i,k为近似程度值;N为声矢量阵的阵元数;Μ为宽带连续谱噪声带宽内的频点 数；II · I U为2-范数;gi(fm)为Ξ维张量过完备原子库拉中对应角度目i的垂直分片中对应 频率fm的列向量；Rl4 -1乂_4)为当前残差信号炒A-n中第η个垂直切分片中对应频率fm的列向 量;i为Ξ维张量过完备原子库中原子的索引值，iei，I为过完备原子库的原子的索引值个 数。5. 根据权利要求4所述的基于稀疏分解理论的声矢量阵宽带测向方法，其特征在于:所 述的步骤（7)具体为：挑选符合条件的原子g补充已选原子集合，选择条件为;其中，qw为已选原子集合. 更新后的信号支撑集合Ω (kl二Q (1-U U { ik}，更新后的已选原子集合;其中，ik为第k次迭代时，Ξ维张量过完备原子库中原子的索引值;馬为 第k次迭代时，Ξ维张量过完备原子库的第ik列数据。6. 根据权利要求5所述的基于稀疏分解理论的声矢量阵宽带测向方法，其特征在于:所 述的步骤(8)具体包括W下步骤： (801)利用最小二乘法计算更新后的已选原子集合在当前残差信号方向上的正交投影 矩阵，计算公式为：其中，PQ(f)为正交投影矩阵;Qf为更新后的已选原子集合，Qf=[gil gi2 为更新 后的已选原子集合的伪逆； (802)根据正交投影矩阵和当前残差信号，计算下一次残差信号，计算方法为：其中，R产1乂。为当前残差信号胁*-11中第X个水平切分片中对应频率fm的列向量，m二 1，2,…，M;I为单位矩阵。7.根据权利要求6所述的基于稀疏分解理论的声矢量阵宽带测向方法，其特征在于:步 骤(10)所述的根据最终的信号支撑集合和宽带连续谱噪声来波方向的对应关系，计算Ξ维 系数矩阵的F-范数，并根据F-范数得到宽带连续谱噪声的空间谱估计;具体包括W下步骤： (111) 将Ξ维系数矩阵按Ξ维张量过完备原子库中原子的索引值划分成列片结构； (112) 根据列片结构划分，计算Ξ维系数矩阵的功率谱估计;计算公式为：其中，P"s-D"p(0i)为Ξ维系数矩阵的功率谱估计；i为Ξ维张量过完备原子库中原子的 索引值；θι为索引值i对应的来波方向，当I个原子组成的过完备原子库对应的来波方向在 整个空域均匀分布时，有表示F-范数;Fe(i)为每一列片结构的数 据矩阵。
【文档编号】G01S7/539GK106093921SQ201610585976
【公开日】2016年11月9日
【申请日】2016年7月25日 公开号201610585976.9, CN 106093921 A, CN 106093921A, CN 201610585976, CN-A-106093921, CN106093921 A, CN106093921A, CN201610585976, CN201610585976.9
【发明人】王晓庆, 李晋, 王大宇
【申请人】中国电子科技集团公司第五十四研究所