一种基于最优空时聚焦技术的声学覆盖层回声降低测量方法技术

一种基于最优空时聚焦技术的声学覆盖层回声降低测量方法，包括以下步骤：1)生成最优空时聚焦发射信号；2)有试样情况下直达信号及回波信号的采集；3)回声降低测量值的计算。本发明专利技术利用各发射阵元及接收阵元之间的信道相干信息，通过目标函数最优化方法，获取各发射阵元最优聚焦发射信号，经同步发射，实现预设水听器处的最优空时聚焦，达到回声降低测量中低频高分辨率及混响抑制的目的，尤其适用于低频段及多层壳体模型回声降低的测量；同时，压力罐测量实验验证了本发明专利技术在声学覆盖层回声降低测量中的有效性。

【技术实现步骤摘要】
一种基于最优空时聚焦技术的声学覆盖层回声降低测量方法
水下声学覆盖层是水声工程中使用广泛且至关重要的水下部件，为评估声学覆盖层试样的回声降低性能，实验室条件下大件样品的回声降低测量是一种重要的手段，本专利技术涉及一种水下声学覆盖层大样回声降低的低频测量方法。
技术介绍
声学覆盖层指敷设在水下部件上的专用功能声学材料与结构，主要由消声瓦、隔声瓦、抑振瓦、去耦瓦、阵列消音器等具有不同声学功能的系列产品综合组成的声学防护体系。声学覆盖层一方面可吸收主动探测声波，降低水下结构的声目标强度，另一方面可作为抑制结构自身辐射噪声的材料。回声信号的强弱与目标的反射特性密切相关，它是衡量水下声学覆盖层吸声性能的重要指标，工程上用“目标强度”来描述目标声反射本领的大小。有限空间近场条件下，无法直接测量试样的目标强度，但可在相同的入射声波条件下，通过测量敷设声学覆盖层前/后模型回声的变化量(即回声降低)来衡量声学覆盖层对模型目标强度的抑制效果。现有的声学覆盖层回声降低的测量方法包括空间傅里叶变换方法，宽带脉冲压缩法，时间反转聚焦方法、多通道空时逆滤波法等。有限空间环境下，混响及多径效应严重，无法实现平面波的精确分解，从而限制了空间傅里叶变换方法的应用；宽带脉冲压缩法可实现时域脉冲信号波形聚焦，但低频条件下，其单通道发射的特点导致指向性变差，衍射及混响干扰较严重；时间反转聚焦方法可实现空时聚焦，但由于其没有抵消信道幅度的影响，聚焦的时域信号脉宽受限，无法实现最优聚焦效果，不利于低频情况下回声降低的测量；多通道空时逆滤波技术没有用到接收阵及发射阵各阵元信道传递函数之间的相干信息，对于具有多个接收阵元的测量系统，没有实现最佳空时聚焦。此外，对于多层壳体模型，通常需分别评估每一层试样的回声性能以实现最优的综合吸声性能评估，因此不但须对直达波和回波进行分离，还须对各层壳体回波进行分离，尤其随着频率的降低，测量对回声信号的时域波形聚焦具有更高的要求。
技术实现思路
为克服现有技术的上述缺点，实现低频、有限空间条件下，声学覆盖层大样回声降低测量中各吸声层反射回波与直达波的分离并进行多径干扰抑制，进一步提高现有的声学覆盖层回声降低测量技术在低频段及多层壳体模型测试中测量结果的精度，本专利技术提出一种基于最优空时聚焦技术的声学覆盖层大样回声降低测量方法，该方法应用多通道目标函数最优化方法，通过目标函数最优解控制多通道发射信号，实现在目标位置处的最优空时聚焦，获取空时域高分辨聚焦信号，从而达到最优聚焦入射、混响抑制及多层壳体回波分离的目的。该方法适用于低频段，在有效减小测量误差、提高测量精度的同时，降低了测量实验对试样尺寸和实验空间的要求。压力罐测量实验验证了本专利技术在声学覆盖层回声降低测量中的有效性。本专利技术为解决上述技术问题所采用的技术方案如下。一种基于最优空时聚焦技术的声学覆盖层回声降低测量方法包括如下步骤：1)生成最优空时聚焦发射信号：无试样情况下，各发射换能器依次发射初始导引信号，各水听器依次接收信号，利用最小二乘方法获取各个通道的信道响应函数，通过目标函数最优化方法获得各发射通道的二次发射信号；2)有试样情况下直达信号及回波信号的采集：将敷设声学覆盖层的试样放入测试环境，各换能器同步发射1)中计算得到的发射信号，则根据目标函数最优化方法的原理，发射信号将在预设水听器位置处产生空时域均为狄拉克函数的波形，水听器记录试样的直达信号pi及经附近试样反射的回波信号pr；3)回声降低测量值的计算：利用回声降低的计算公式可得其中Er表示回声降低值。本专利技术的技术构思为：通过初始导引信号发射及水听器逐一接收，得到包含电路信道和水声信道信息的接收信号，进而通过目标函数最优化方法，获得换能器阵每一阵元的二次发射信号，然后各换能器同步发射，获得预设水听器处的最优聚焦信号，从而最优化地实现测量信号的空间聚焦及时域脉冲压缩，达到抑制低频混响的目的。与现有的回声降低测量方法相比，本专利技术的技术优势主要表现在：利用各发射阵元及接收阵元之间的信道相干信息，通过目标函数最优化方法，获取各发射阵元最优聚焦发射信号，经同步发射，实现预设水听器处的最优空时聚焦，达到回声降低测量中低频高分辨率及混响抑制的目的，尤其适用于低频段及多层壳体模型回声降低的测量。附图说明图1是本专利技术方法的测量系统整体示意图。图2是压力罐环境下5mm厚钢板试样回声降低理论及测量试验结果对比图。具体实施方式下面结合附图对本专利技术作进一步描述。参照图1～图2，一种基于最优空时聚焦技术的声学覆盖层回声降低测量方法，用于有限空间声学覆盖层大样回声降低的测量。整套测量方法的技术方案如下：1)生成最优空时聚焦发射信号；首先，利用最小二乘算法估计信道响应；无试样情况下，各发射换能器分别发射导引信号，频域表示为en(ω)，n＝1,2,…,N，经电路信道和水声信道传递后被第m个水听器接收，输出信号表示为y(ω,rn；rm)，m＝1,2,…,M；其中N、M分别为发射换能器和水听器的个数，ω表示角频率，rn、rm分别表示第n个发射换能器和第m个水听器的位置；则频域上该过程可表示为y(ω,rn；rm)＝en(ω)·G(ω,rn；rm)(2)其中G(ω,rn；rm)为第n个发射通道与第m个接收通道之间的信道响应函数；根据最小二乘原理，定义目标函数P(G)＝(Y-EG)H(Y-EG)(3)其中上标H表示共轭转置，G表示不同收发通道响应函数组成的矩阵Y是输出信号组成的矩阵E是由输入信号组成的对角矩阵，即E＝diag[e1(ω),e2(ω),…,eN(ω)](6)其中diag表示以向量[e1(ω),e2(ω),…,eN(ω)]为对角元的矩阵；为求P(G)的最小值，令其对G的偏导数等于零，即计算得到信道响应G＝(EHE)-1EHY(8)得到信道响应后，为求解最优化的二次发射信号，建立如下目标函数J(s1(ω),s2(ω),…,sN(ω))＝∫W|Φ(ω,r)-exp(-iωT)δ(r-rc)|2dr(9)其中，W为水听器的位置集合。另设sn(ω)对应的时域信号为sn(t)，n＝1,2,…,N，t表示时间；当sn(t)在t＝0时刻发射，经信道传递后设其在预设聚焦位置rc处和T时刻聚焦，要实现空时最佳聚焦，预设位置处的聚焦信号在频域上的理论值可表示为exp(-iωT)δ(r-rc)，即空时均为狄拉克函数；设Φ(ω,rm)，m＝1,2,…,M，为第m个水听器振元处的接收信号，可表示为其中，xn(ω,rm)，m＝1,2,…,M，为第m个水听器位置处的接收信号，表示为xn(ω,rm)＝sn(ω)·G(ω,rn；rm)(11)将Φ(ω,rm)带入目标函数并进行最小化求解，可得最优二次发射信号为Γ(ω)s(ω)＝exp(-iωT)g*(ω)(12)其中g*(ω)表示函数g(ω)的共轭函数，Γ矩阵的组成元素为Γnm(ω)＝∫WG(ω,rn；r)G*(ω,rm；r)dr(13)向量g(ω)的组成元素为gn(ω)＝G(ω,rn；rc),n＝1,2,…,N(14)即Γ为各信道响应函数的互相关表示，g(ω)为各发射阵元到预设聚焦位置之间的响应函数组成的向量；则最优二次发射信号可表示为其中表示矩阵Γ的伪逆矩阵，对该信号做逆傅里叶变换得到时域最优空时聚焦的发射信号本文档来自技高网...

【技术保护点】
1.一种基于最优空时聚焦技术的声学覆盖层回声降低测量方法，包括以下步骤：1)生成最优空时聚焦发射信号：无试样情况下，各发射换能器依次发射初始导引信号，各水听器依次接收信号，利用最小二乘方法获取各个通道的信道响应函数，通过目标函数最优化方法获得各发射通道的二次发射信号；2)有试样情况下直达信号及回波信号的采集：将敷设声学覆盖层的试样放入测试环境，各换能器同步发射1)中计算得到的发射信号，则根据目标函数最优化方法的原理，发射信号将在预设水听器位置处产生空时域均为狄拉克函数的波形，水听器记录试样的直达信号pi及经附近试样反射的回波信号pr；3)回声降低测量值的计算：利用回声降低的计算公式可得

【技术特征摘要】
1.一种基于最优空时聚焦技术的声学覆盖层回声降低测量方法，包括以下步骤：1)生成最优空时聚焦发射信号：无试样情况下，各发射换能器依次发射初始导引信号，各水听器依次接收信号，利用最小二乘方法获取各个通道的信道响应函数，通过目标函数最优化方法获得各发射通道的二次发射信号；2)有试样情况下直达信号及回波信号的采集：将敷设声学覆盖层的试样放入测试环境，各换能器同步发射1)中计算得到的发射信号，则根据目标函数最优化方法的原理，发射信号将在预设水听器位置处产生空时域均为狄拉克函数的波形，水听器记录试样的直达信号pi及经附近试样反射的回波信号pr；3)回声降低测量值的计算：利用回声降低的计算公式可得其中Er表示回声降低值。2.如权利要求1所述的一种基于最优空时聚焦技术的声学覆盖层回声降低测量方法，其特征在于：所述步骤1)中，生成最优空时聚焦发射信号的过程如下：首先，利用最小二乘算法估计信道响应；无试样情况下，各发射换能器分别发射导引信号，频域表示为en(ω)，n＝1,2,…,N，经电路信道和水声信道传递后被第m个水听器接收，输出信号表示为y(ω,rn；rm)，m＝1,2,…,M，其中N、M分别为发射换能器和水听器的个数，ω表示角频率，rn、rm分别表示第n个发射换能器和第m个水听器的位置；则频域上该过程可表示为y(ω,rn；rm)＝en(ω)·G(ω,rn；rm)(2)其中G(ω,rn；rm)为第n个发射通道与第m个接收通道之间的信道响应函数；根据最小二乘原理，定义目标函数P(G)＝(Y-EG)H(Y-EG)(3)其中上标H表示共轭转置，G表示不同收发通道响应函数组成的矩阵Y是输出信号组成的矩阵E是由输入信号组成的对角矩阵，即E＝diag[e1(ω),e2...

【专利技术属性】
技术研发人员：李建龙，马骁晨，李素旋，
申请(专利权)人：浙江大学，
类型：发明
国别省市：浙江,33