本发明专利技术涉及一种基于柯西-高斯模型的新型高分辨方位估计方法，以满足水下宽带高频三维成像掩埋雷探测声纳对目标方位分辨力的要求。利用柯西-高斯信号模型，将方位谱表示成一个正则化的空域傅里叶变换约束优化问题，提出一种对模型参数宽容的谱估计器，通过寻求源信号在空间中的稀疏分布特征，在单频域快拍内实现高分辨方位谱估计，且无需解相干即可分辨相干源。本发明专利技术有益的效果是：通过代价函数的约束优化实现了单快拍内的方位高分辨。且本方法对信号概率分布模型参数的宽容性高，因而具有工程实用性。数值仿真和海试数据分析结果说明本发明专利技术方法适用于欠快拍条件下的小孔径基阵的目标方位高分辨。
【技术实现步骤摘要】
本专利技术涉及传感器阵列信号处理领域，主要是。
技术介绍
信号方位估计问题一直以来在声纳、雷达和通信等领域引起了广泛的研究兴趣。 常规波束形成通过空间匹配滤波来检测和估计目标方位，但目标方位分辨力取决于基阵波长单位的孔径大小。自适应波束形成[1]和子空间高分辨波达方向(direction of arrival, D0A)估计算法，如MUSIC，能够获取高分辨性能，但需要足够的信噪比条件，并且需要较长的时间平均以提高协方差矩阵估计的精度。信号子空间分解算法需要预估计信号源数，当存在相干源时，需要通过子阵空间平滑来解相干。因而，在基阵孔径受限制的场合和高动态条件下，常规的高分辨处理方法难以有效应用。空域约束优化处 (spatial processing optimized and constrained, SP0C) 是近几年出现的一种新型单快拍高分辨算法，SPOC算法基于MAP准则，通过信号模型的约束，最优估计出在每一时刻，信号源在空间中的稀疏功率分布函数，因而突破了快拍数和相干源的限制。韩国Kam等人将SPOC算法成功应用于海底合成孔径声纳布放目标的方位高分辨处理。本专利技术利用波束形成与空域傅里叶变换在数学上的等效性，研究一种迭代非参数化方位谱估计技术。利用柯西(Cauchy)分布来近似空域DFT信号样本，并对信号最小模约束代价函数进行正则化，对正则化代价函数进行优化得到空域高分辨DFT。提出一种次优的、对信号模型参数宽容的方位谱估计器，以提高算法的稳健性。
技术实现思路
本专利技术的目的是克服现有技术的缺点和不足，提供，以满足水下宽带高频三维成像掩埋雷探测声纳对目标方位分辨力的要求。本专利技术解决其技术问题采用的技术方案这种基于柯西-高斯模型的新型高分辨方位估计方法，包括如下步骤(1)将基阵接收时域信号变换到频域X (k) = [X0(k), ... , Xsh (k) ]τ (1)将第k个频点第η个波束输出信号表示为空域DFT形式tool ο] Y(n,k)=Yjm^XSkyj2xfkmd^ejctoon] O揪—'Ν,Ν>Μ(2)N为空域DFT点数，第η个离散空间方位满足θ n = asin {cn/ (Nfkd)} (3)(2)利用DFT的性质，将空间DFT谱估计表示成线性反演求解形式X (k) = FY (k) (6)其中，F为Rmrier变换矩阵，Y(k)第k个频点的波束输出信号向量；(3)构造正则代价函数J = Σ 二 log (I Y{n, k) I2 )+\\X-FY Il^⑷(4) Y = ( λ Q^+FhF) ^1FhX(H)其中，义/σ#，Q为对角矩阵Q = diag{IN+YYii/a2Y}(12)对式(11)中的调节参数λ取固定的经验值，再按以下方式迭代求解高分辨空间谱1)以Y = DFT {X}作为初始值，并代入式(17)得到Q(°)；0 = diag{ITH}(17)2)第i次迭代，有Y⑴=Q(B)Fh ( λ IJFQ(H)Fh) ^1X (5)利用上式更新Q矩阵Q⑴=diag{Y⑴Y⑴H}；3)判断收敛条件I J⑴—J—1) IzlJ(H)I < ε (6)ε为一个小量，若不满足条件，令i = i+Ι，返回第2)步；若满足条件，则输出当前迭代周期的空间谱Y(i)；(5)输出满足收敛条件的高分辨空间谱到后续的图像显示处理系统。本专利技术有益的效果是本专利技术所研究的一种基于柯西-高斯信号模型的高分辨方位估计技术，通过代价函数的约束优化实现了单快拍内的方位高分辨。且本方法对信号概率分布模型参数的宽容性高，因而具有工程实用性。数值仿真和海试数据分析结果说明本专利技术方法适用于欠快拍条件下的小孔径基阵的目标方位高分辨。类似地，本专利技术方法可用于舰艇辐射噪声的高分辨线谱估计。附图说明图1为CBF频率-波数谱(无法分辨归一化频率为0. 2处的两个信号)；图2为RCG-SPEC算法频率-波数谱(提高了算法的稳健性，同时保证了高分辨性能)；具体实施例方式下面结合附图和实施例对本专利技术作进一步说明本专利技术所述的这种基于柯西-高斯模型的新型高分辨方位估计方法，包括以下步骤步骤一将空间谱估计表示成空域DFT的线性反演求解形式X (k) = FY (k) (6)其中，X(k)为第k个频点的基阵频域信号向量，F为Fourier变换矩阵，Y(k)第k 个频点的波束输出信号向量。步骤二 基于柯西-高斯信号模型，求空域DFT线性反演的约束优化正则解首先构造代价函数本文档来自技高网...
【技术保护点】
１．一种基于柯西－高斯模型的新型高分辨方位估计方法，其特征是：包括如下步骤：（１）将基阵接收时域信号变换到频域：Ｘ（ｋ）＝［Ｘ０（ｋ），．．．，ＸＭ－１（ｋ）］Ｔ （１）将第ｋ个频点第ｎ个波束输出信号表示为空域ＤＦＴ形式：（ｍａｔｈ）??（ｍｒｏｗ）?（ｍｉ）Ｙ（／ｍｉ）?（ｍｒｏｗ）?（ｍｏ）（（／ｍｏ）?（ｍｉ）ｎ（／ｍｉ）?（ｍｏ），（／ｍｏ）?（ｍｉ）ｋ（／ｍｉ）?（ｍｏ））（／ｍｏ）?（／ｍｒｏｗ）?（ｍｏ）＝（／ｍｏ）?（ｍｓｕｂｓｕｐ）?（ｍｉ）＆Ｓｉｇｍａ；（／ｍｉ）?（ｍｒｏｗ）?（ｍｉ）ｍ（／ｍｉ）?（ｍｏ）＝（／ｍｏ）?（ｍｎ）０（／ｍｎ）?（／ｍｒｏｗ）?（ｍｒｏｗ）?（ｍｉ）Ｍ（／ｍｉ）?（ｍｏ）－（／ｍｏ）?（ｍｎ）１（／ｍｎ）?（／ｍｒｏｗ）?（／ｍｓｕｂｓｕｐ）?（ｍｓｕｂ）?（ｍｉ）Ｘ（／ｍｉ）?（ｍｉ）ｍ（／ｍｉ）?（／ｍｓｕｂ）?（ｍｒｏｗ）?（ｍｏ）（（／ｍｏ）?（ｍｉ）ｋ（／ｍｉ）?（ｍｏ））（／ｍｏ）?（／ｍｒｏｗ）?（ｍｓｕｐ）?（ｍｉ）ｅ（／ｍｉ）?（ｍｒｏｗ）?（ｍｏ）－（／ｍｏ）?（ｍｉ）ｊ（／ｍｉ）?（ｍｎ）２（／ｍｎ）?（ｍｉ）＆ｐｉ；（／ｍｉ）?（ｍｓｕｂ）?（ｍｉ）ｆ（／ｍｉ）?（ｍｉ）ｋ（／ｍｉ）?（／ｍｓｕｂ）?（ｍｉ）ｍｄ（／ｍｉ）?（ｍｉ）ｓｉｎ（／ｍｉ）?（ｍｓｕｂ）?（ｍｉ）＆ｔｈｅｔａ；（／ｍｉ）?（ｍｉ）ｎ（／ｍｉ）?（／ｍｓｕｂ）?（ｍｏ）／（／ｍｏ）?（ｍｉ）ｃ（／ｍｉ）?（／ｍｒｏｗ）?（／ｍｓｕｐ）?（／ｍｒｏｗ）?（／ｍａｔｈ）（ｍａｔｈ）??（ｍｒｏｗ）?（ｍｏ）＝（／ｍｏ）?（ｍｓｕｂｓｕｐ）?（ｍｉ）＆Ｓｉｇｍａ；（／ｍｉ）?（ｍｒｏｗ）?（ｍｉ）ｍ（／ｍｉ）?（ｍｏ）＝（／ｍｏ）?（ｍｎ）０（／ｍｎ）?（／ｍｒｏｗ）?（ｍｒｏｗ）?（ｍｉ）Ｎ（／ｍｉ）?（ｍｏ）－（／ｍｏ）?（ｍｎ）１（／ｍｎ）?（／ｍｒｏｗ）?（／ｍｓｕｂｓｕｐ）?（ｍｓｕｂ）?（ｍｉ）Ｘ（／ｍｉ）?（ｍｉ）ｍ（／ｍｉ）?（／ｍｓｕｂ）?（ｍｒｏｗ）?（ｍｏ）（（／ｍｏ）?（ｍｉ）ｋ（／ｍｉ）?（ｍｏ））（／ｍｏ）?（／ｍｒｏｗ）?（ｍｓｕｐ）?（ｍｉ）ｅ（／ｍｉ）?（ｍｒｏｗ）?（ｍｏ）－（／ｍｏ）?（ｍｉ）ｊ（／ｍｉ）?（ｍｎ）２（／ｍｎ）?（ｍｉ）＆ｐｉ；ｍｎ（／ｍｉ）?（ｍｏ）／（／ｍｏ）?（ｍｉ）Ｎ（／ｍｉ）?（／ｍｒｏｗ）?（／ｍｓｕｐ）?（ｍｏ），（／ｍｏ）?（ｍｉ）Ｎ（／ｍｉ）?（ｍｏ））（／ｍｏ）?（ｍｉ）Ｍ（／ｍｉ）?（ｍｏ）－（／ｍｏ）?（ｍｏ）－（／ｍｏ）?（ｍｏ）－（／ｍｏ）?（ｍｒｏｗ）?（ｍｏ）（（／ｍｏ）?（ｍｎ）２（／ｍｎ）?（ｍｏ））（／ｍｏ）?（／ｍｒｏｗ）?（／ｍｒｏｗ）?（／ｍａｔｈ）Ｎ为空域ＤＦＴ点数，第ｎ个离散空间方位满足θｎ＝ａｓｉｎ｛ｃｎ／（Ｎｆｋｄ）｝ （３）（２）利用ＤＦＴ的性质，将空间ＤＦＴ谱估计表示成线性反演求解形式：Ｘ（ｋ）＝ＦＹ（ｋ） （６）其中，Ｆ为Ｆｏｕｒｉｅｒ变换矩阵，Ｙ（ｋ）第ｋ个频点的波束输出信号向量；（３）构造正则代价函数：（ｍａｔｈ）??（ｍｒｏｗ）?（ｍｏｖｅｒ）?（ｍｉ）Ｊ（／ｍｉ）?（ｍｏ）～（／ｍｏ）?（／ｍｏｖｅｒ）?（ｍｏ）＝（／ｍｏ）?（ｍｓｕｂｓｕｐ）?（ｍｉ）＆Ｓｉｇｍａ；（／ｍｉ）?（ｍｒｏｗ）?（ｍｉ）ｎ（／ｍｉ）?（ｍｏ）＝（／ｍｏ）?（ｍｎ）０（／ｍｎ）?（／ｍｒｏｗ）?（ｍｒｏｗ）?（ｍｉ）Ｎ（／ｍｉ）?（ｍｏ）－（／ｍｏ）?（ｍｎ）１（／ｍｎ）?（／ｍｒｏｗ）?（／ｍｓｕｂｓｕｐ）?（ｍｉ）ｌｏｇ（／ｍｉ）?（ｍｒｏｗ）?（ｍｏ）（（／ｍｏ）?（ｍｓｕｐ）?（ｍｒｏｗ）?（ｍｏ）｜（／ｍｏ）?（ｍｉ）Ｙ（／ｍｉ）?（ｍｒｏｗ）?（ｍｏ）（（／ｍｏ）?（ｍｉ）ｎ（／ｍｉ）?（ｍｏ），（／ｍｏ）?（ｍｉ）ｋ（／ｍｉ）?（ｍｏ））（／ｍｏ）?（／ｍｒｏｗ）?（ｍｏ）｜（／ｍｏ）?（／ｍｒｏｗ）?（ｍｎ）２（／ｍｎ）?（／ｍｓｕｐ）?（ｍｏ））（／ｍｏ）?（／ｍｒｏｗ）?（ｍｏ）＋（／ｍｏ）?（ｍｓｕｂｓｕｐ）?（ｍｒｏｗ）?（ｍｏ）｜（／ｍｏ）?（ｍｏ）｜（／ｍｏ）?（ｍｉ）Ｘ（／ｍｉ）?（ｍｏ）－（／ｍｏ）?（ｍｉ）ＦＹ（／ｍｉ）?（ｍｏ）｜（／ｍｏ）?（ｍｏ）｜（／ｍｏ）?（／ｍｒｏｗ）?（ｍｎ）２（／ｍｎ）?（ｍｎ）２（／ｍｎ）?（／ｍｓｕｂｓｕｐ）?（ｍｏ）－（／ｍｏ）?（ｍｏ）－（／ｍｏ）?（ｍｏ）－（／ｍｏ）?（ｍｒｏｗ）?（ｍｏ）（（／ｍｏ）?（ｍｎ）４（／ｍｎ）?（ｍｏ））（／ｍｏ）?（／ｍｒｏｗ）?（／ｍｒｏｗ）?（／ｍａｔｈ）（４）Ｙ＝（λＱ－１＋ＦＨＦ）－１ＦＨＸ ...
【技术特征摘要】
1. 一种基于柯西-高斯模型的新型高分辨方位估计方法，其特征是包括如下步骤(1)将基阵接收时域信号变换到频域 X(k) = D(0(k)，···，ΧηΟΟΓ (ι)将第k个频点第η个波束输出信号表示为空域DFT形式= Ztlx^j2xfimdsiaejc= Σ :1οΧΛ^2πη,/Ν,Ν>Μ⑵N为空域DFT点数，第η个离散空间方位满足 θ n = asin{cn/(Nfkd)}(3)(2)利用DFT的性质，将空间DFT谱估计表示成线性反演求解形式 X(k) = FY (k)(6)其中，F为Rmrier变换矩阵，Y(k)第k个频点的波束输出信号向量；(3)构造正则代价函数J 二 Σ二 l0§ (I nn, A) I2 )+11 X-...
【专利技术属性】
技术研发人员：蒋飚，从卫华，朱必波，郑乐一，
申请(专利权)人：中国船舶重工集团公司第七一五研究所，
类型：发明
国别省市：86
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。