当前位置: 首页 > 专利查询>浙江大学专利>正文

基于近场测量数据的阵列诊断方法技术

技术编号:27058002 阅读:45 留言:0更新日期:2021-01-15 14:36
本发明专利技术公开了一种基于近场测量数据的阵列诊断方法,包括以下步骤:(1)平面阵列阵元接收声波入射信号,基于接收的声波入射信号构建参考阵列和待测阵列的近场波束图;(2)计算参考阵列的近场波束图和待测阵列的近场波束图之差得到差分波束图,并将差分波束图看作参考阵列与待测阵列的阵元系数差分向量与观测矩阵相乘的结果;(3)利用贝叶斯压缩感知算法求解参考阵列与待测阵列的阵元系数差分向量,从而定位到故障阵元的位置和系数,实现阵列诊断的诊断。该阵列诊断方法基于贝叶斯压缩感知的方法采用近场测量数据,将故障阵元诊断问题转化为贝叶斯概率框架下的稀疏向量求解,有效地提高了诊断的准确率和计算效率。

【技术实现步骤摘要】
基于近场测量数据的阵列诊断方法
本专利技术涉及声纳阵列信号处理领域,具体涉及一种使用近场测量数据的阵列诊断方法。
技术介绍
近年来,相控阵三维成像声纳技术由于其在水下物理、生物、地质等方面的应用而得到了迅速发展。为了获得高分辨率和高性能的波束图,声呐阵列需要大量的阵元组成,因此其中不可避免地存在故障阵元。为了保持系统成像的性能,必须准确检测故障阵元的位置和系数。因此,研究故障阵元诊断技术具有重要意义。现有的故障诊断技术有论文ACompressedSensingApproachforArrayDiagnosisFromaSmallSetofNear-FieldMeasurements提出的基于l1-范数法的故障阵元诊断技术。论文Accuratediagnosisofconformalarraysfromnear-fielddatausingthematrixmethod提出的基于矩阵方法的故障阵元诊断技术。
技术实现思路
本专利技术的目的是提供一种基于近场测量数据的阵列诊断方法,该阵列诊断方法基于贝叶斯压缩感知的方法使用近场测量数据,将故障阵元诊断问题转化为贝叶斯概率框架下的稀疏向量求解,有效地提高了诊断的准确率和计算效率。为实现上述专利技术目的,本专利技术提供以下技术方案:一种基于近场测量数据的阵列诊断方法,包括以下步骤:(1)平面阵列阵元接收声波入射信号,基于接收的声波入射信号构建参考阵列和待测阵列的近场波束图;(2)计算参考阵列的近场波束图和待测阵列的近场波束图之差得到差分波束图,并将差分波束图看作参考阵列与待测阵列的阵元系数差分向量与观测矩阵相乘的结果;(3)利用贝叶斯压缩感知算法求解参考阵列与待测阵列的阵元系数差分向量,从而定位到故障阵元的位置和系数,实现阵列诊断。步骤(1)中,对于一个N个阵元均匀分布且不存在故障阵元的平面参考阵列,其近场波束图F如下所示:其中,wn表示阵元系数,w表示阵元系数矩阵,λ表示信号波长,表示第p个近场测量点位置,p=1,2,…,P,P为自然数,表示第n个阵元的位置,Φ表示观测矩阵。在近场波束图F时,通过构建声呐阵列的近场球面波传播模型,计算出近场测量点的观测矩阵Φ。由公式(1)可知近场波束图综合的关键在于计算出不同阵元与测量点的路程差,即的值,的计算具体如下所示:其中,D是阵列的孔径大小,αp,βp表示第p个测量点的仰角和方位角,rp表示测量点到阵列中心的距离,xn,yn表示第n个阵元的坐标。由于阵元故障的位置的随机性,即在实际系统中每一个阵元出现故障的概率是相等的,因此为保证算法的普适性,近场测量点的位置应该满足均匀分布,本专利技术中,当测量点数P为平方数时,测量点的均匀分布是易于实现的,此时其中,αmin和βmin表示仰角和方位角的最小值,αmax和βmax表示仰角和方位角的最大值。然而,当P不是平方数时,上述策略便无法实现。为此,采用泊松盘采样方法对近场测量点进行采样,实现了任意数量的均匀测量点分布。首先确定测量点之间的最小间距r为如下所示:然后随机生成一个测量点,以该测量点为圆心在距离为r到2r的圆环内随机生成k个测量点,判断其与已确定的测量点的距离,若小于r则清除,大于r则保留下来。对这k个点操作完毕之后将圆心的点标记为非活跃,再对之前保留的点重复上述操作直到生成P个测量点。存在故障阵元的待测阵列的近场波束图如下所示:其中,表示待测阵列的阵元系数矩阵,Φ表示观测矩阵,为测量过程中产生的均值为零的复数高斯噪声向量。步骤(2)中,计算待测阵列与参考阵列的近场波束图之差得到差分波束图:其中,ΔF表示待测阵列与参考阵列的近场波束图之差,Δw表示待测阵列与参考阵列的阵元系数之差,步骤(3)中,阵元故障诊断可视为在待测阵列与参考阵列的近场波束图之差已知的情况下寻找其阵元系数之差,问题可描述为:minΔw||Δw||0,subjectto∫∫||ΔF-ΦΔw||2dαdβ=e(6)基于该问题描述,采用贝叶斯压缩感知算法对问题描述求解,得到阵元系数差分向量Δw,从而定位到故障阵元的位置和系数,实现阵列诊断的诊断。具体地,采用贝叶斯压缩感知算法对问题描述求解时,由于Δw是稀疏的,则公式(6)的问题描述转化为贝叶斯后验概率问题,如下所示:其中,G=R,I,P(ΔwG|ΔFG)表示后验概率,R,I分别表示实部和虚部,令ΦR和ΦI分别表示Φ的实部和虚部,则ΔF的计算如下所示:同时,P(ΔwG|ΔFG)可以进一步表示为:P(ΔwG|ΔFG)=∫P(ΔwG|ΔFG,γG)P(γG|ΔFG)dγGG=R,I(11)其中,γG为控制后验概率分布的超参数向量,其值可以通过求解最大似然函数得到,如下所示:其中,a,b为用户自定义比例控制参数,diag(γG)表示以γG为对角元素的对角矩阵,上标T表示矩阵的转置,G=R,I,最终,阵元系数差分向量ΔwG求解如下所示:ΔwG=(diag(γG)+ΩGTΩG)-1ΩGTΔFG(13)与现有技术相比,本专利技术具有以下有效效果:本专利技术提供的基于近场测量数据的阵列诊断方法中,通过接收的声波入射信号构建参考阵列和待测阵列的近场波束图,将参考阵列与待测阵列的差分波束图看作参考阵列与待测阵列的阵元系数差分向量与观测矩阵相乘的结果,利用贝叶斯压缩感知算法求解参考阵列与待测阵列的阵元系数差分向量,从而定位到故障阵元的位置和系数以实现阵列诊断。相对于其他阵列诊断方法,该方法具有更低的诊断误差和更高的计算效率。附图说明为了更清楚地说明本专利技术实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图做简单地介绍,显而易见地,下面描述中的附图仅仅是本专利技术的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动前提下,还可以根据这些附图获得其他附图。图1为本专利技术的基于近场测量数据的阵列诊断方法中近场测量点示意图;图2为本专利技术的基于近场测量数据的阵列诊断方法示例中的待测阵列阵元系数;图3为本专利技术的基于近场测量数据的阵列诊断方法的诊断误差;图4为作为对比的l1范数法诊断误差。具体实施方式为使本专利技术的目的、技术方案及优点更加清楚明白,以下结合附图及实施例对本专利技术进行进一步的详细说明。应当理解,此处所描述的具体实施方式仅仅用以解释本专利技术,并不限定本专利技术的保护范围。本实施例中,设计所考虑的是一个33×33的二维换能器阵列。换能器按半波长间距均匀分布在一个矩形平面内,换能器的水平和垂直间距相等,系数服从切比雪夫分布,载波频率为f=300kHz,声速为1500m/s,仰角和方位角的取值范围分别是(-90°,90°)和(-90°,90°)。近场聚焦距离为D2/4λ,测量波束数为15×15,假设阵列的本文档来自技高网...

【技术保护点】
1.一种基于近场测量数据的阵列诊断方法,其特征在于,包括以下步骤:/n(1)平面阵列阵元接收声波入射信号,基于接收的声波入射信号构建参考阵列和待测阵列的近场波束图;/n(2)计算参考阵列的近场波束图和待测阵列的近场波束图之差得到差分波束图,并将差分波束图看作参考阵列与待测阵列的阵元系数差分向量与观测矩阵相乘的结果;/n(3)利用贝叶斯压缩感知算法求解参考阵列与待测阵列的阵元系数差分向量,从而定位到故障阵元的位置和系数,实现阵列诊断。/n

【技术特征摘要】
1.一种基于近场测量数据的阵列诊断方法,其特征在于,包括以下步骤:
(1)平面阵列阵元接收声波入射信号,基于接收的声波入射信号构建参考阵列和待测阵列的近场波束图;
(2)计算参考阵列的近场波束图和待测阵列的近场波束图之差得到差分波束图,并将差分波束图看作参考阵列与待测阵列的阵元系数差分向量与观测矩阵相乘的结果;
(3)利用贝叶斯压缩感知算法求解参考阵列与待测阵列的阵元系数差分向量,从而定位到故障阵元的位置和系数,实现阵列诊断。


2.如权利要求1所述的基于近场测量数据的阵列诊断方法,其特征在于,步骤(1)中,参考阵列的近场波束图F为:



其中,wn表示第n个阵元系数,w表示参考阵列的阵元系数矩阵,λ表示信号波长,表示第p个近场测量点位置,p=1,2,…,P,P为自然数,表示第n个阵元的位置,n=1,2,…,N,N为自然数,Φ表示观测矩阵;
其中,阵元与近场测量点的路程差具体计算为:



其中,D是阵列的孔径大小,αp,βp表示第p个测量点的仰角和方位角,rp表示测量点到阵列中心的距离,xn,yn表示第n个阵元的坐标;
待测阵列的近场波束图为:



其中,表示待测阵列的阵元系数矩阵,Φ表示观测矩阵,为测量过程中产生的均值为零的复数高斯噪声向量。


3.如权利要求2所述的基于近场测量数据的阵列诊断方法,其特征在于,步骤(2)中,差分波束图ΔF为:



其中,Δw表示待测阵列与参考阵列的阵元系数之差,


4.如权利要求3所述的基于近场测量数据的阵列诊断方法,其特征在于,将阵元故障诊断视为在待测阵列与参考阵列的差分波束图ΔF已知的情况下寻找其阵元系数差分向量Δw,问题描述为:
minΔw‖Δw‖0,subjectto∫∫‖ΔF-ΦΔw‖2dαdβ=e(5...

【专利技术属性】
技术研发人员:陈耀武林振伟蒋荣欣刘雪松
申请(专利权)人:浙江大学
类型:发明
国别省市:浙江;33

网友询问留言 已有0条评论
  • 还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。

1