一种声学多普勒长期测速精度估计方法技术

一种声学多普勒长期测速精度估计方法，属于水声导航及组合导航技术领域。解决了声学多普勒测速长期精度估计准确性差的问题。本发明专利技术先计算换能器波束照射的海底散射区域内所有散射体的指向性、散射强度、传播损失和多普勒相位；将换能器接收到散射体反射信号时延相同的N个散射体作为一个微元，利用所述N个散射体的散射强度、传播损失和多普勒相位，计算微元的脉冲响应函数的二阶统计量，再计算当散射体个数趋近无穷时，二阶统计量的极限值，进而计算散射区域的多普勒谱的一阶谱距，通过一阶谱距建立声学多普勒长期测速精度的解析公式，对声学多普勒长期测速精度值估计。本发明专利技术适用于声学多普勒长期测速精度估计。声学多普勒长期测速精度估计。声学多普勒长期测速精度估计。

【技术实现步骤摘要】
一种声学多普勒长期测速精度估计方法


[0001]本专利技术属于水声导航及组合导航


技术介绍

[0002]声学多普勒测速通过倾斜波束照向海底，并依靠提取散射区域的回波信号多普勒信息进行速度测量的，其特点是测速精度高，是现代舰船导航系统的重要组成。声学多普勒测速长期精度是指测量速度的平均值偏离真值的大小，它是排除测速方差影响后的残余精度。换句话说，当测速数据被累加或平均时，由于测速长期精度的影响，结果将产生误差。在一个简单的航位推算应用中，声学多普勒测速数据可以积分一次以获得位置，其中测速长期精度引起的位置误差将随着时间线性增加。测速长期精度是声学多普勒测速声纳设计者和使用者常用的评价指标。然而，目前评价测速长期精度的方法较单一，准确性差，通常采用模型仿真和曲线拟合方式获得量化测速长期精度的公式。类似这样的测速长期精度估计方法物理机制不够明确，不利于在工程上广泛应用。

技术实现思路

[0003]本专利技术是为了解决声学多普勒测速长期精度估计准确性差的问题，现提供一种声学多普勒长期测速精度估计方法。
[0004]本专利技术所述一种声学多普勒长期测速精度估计方法，包括：
[0005]步骤一、以换能器到散射体的矢量r与换能器的运动速度矢量v方向的夹角为α角，以换能器到散射体的矢量r在垂面投影与垂直方向的夹角为η角，建立(α,η)角度坐标系；
[0006]步骤二、在(α,η)角度坐标系内，计算换能器波束照射的海底散射区域内所有散射体的指向性、散射强度、传播损失和多普勒相位；/>[0007]步骤三、将换能器接收到散射体反射信号时延相同的N个散射体作为一个微元，利用所述N个散射体的散射强度、传播损失和多普勒相位，计算微元的脉冲响应函数其中，τ为时延变量，t为时间变量，为微元的时延；
[0008]步骤四、计算微元脉冲响应函数的二阶统计量其中，Δt为相关时间，τ'＝τ+Δτ，Δτ时间延迟差；
[0009]步骤五、计算当微元中散射体个数N
→
∞时，二阶统计量的极限值；
[0010]步骤六、利用所述二阶统计量的极限值，计算散射区域的系统相关函数A
h
(τ,Δt)；
[0011]步骤七、利用所述散射区域的系统相关函数A
h
(τ,Δt)，计算散射区域的多普勒谱S(ρ)，再计算多普勒谱S(ρ)的一阶谱距M1；
[0012]步骤八、计算多普勒谱S(ρ)的一阶谱距M1相对波束方向上的理论多普勒频率ω
d
的误差，利用所述误差建立声学多普勒长期测速精度的解析公式，获取声学多普勒长期测速精度值。
[0013]进一步地，本专利技术中，步骤二中，计算换能器能采集到的所有散射体的指向性的方法均相同，具体为：
[0014]利用公式：
[0015][0016]计算获得指向性G(α,η)，其中，Δ1表示指向性G(α,η)在α角度方向上的等效波束开角，Δ2表示指向性G(α,η)在η角度方向上的等效波束开角，α0表示换能器散射波束中心角，
[0017][0018]其中，G(θ)表示声线与波束中心轴线夹角θ的指向性函数，α(x,z)表示用x和z变量表示的α变量，η(x,z)表示用x和z变量表示的η变量；
[0019]η(x,z)＝arctan(z/H)
[0020][0021]式中，H表示换能器到达海底自身投影点的垂直距离。
[0022]进一步地，本专利技术中，步骤二中，计算换能器能采集到的所有散射体的散射强度的方法相同，具体为：
[0023][0024]其中，表示散射体的强度，表示波束中心散射体的散射强度，表示散射体散射波的入射角，指波束中心处散射体的入射角，ο表示α相对波束中心角α0的差角，P
MSP
‑
in
、P
MSP
‑
sact
分别表示入射和散射均平方声压，ε是散射强度变化率。
[0025]进一步地，本专利技术中，步骤二中，计算换能器能采集到的所有散射体的传播损失的方法相同，具体为：
[0026][0027]其中，A为散射体的传播损失，r0指参考距离，r1表示散射体距换能器的距离，β是声吸收系数，P
MSP
(r1)、P
MSP
(r0)分别指与散射体相距r1、r0位置处的均平方声压。
[0028]进一步地，本专利技术中，步骤二中，计算换能器能采集到的所有散射体的多普勒相位的方法相同，具体为：
[0029]ΔΦ≈ω
d
Δt
‑
poΔt,
[0030]ΔΦ为散射体的多普勒相位，Δt为相关时延，ω
d
为理论多普勒频率，c为换能器
表层声速，v为换能器的运动速度；v表示速度矢量，|v|＝v。
[0031]进一步地，本专利技术中，步骤三中，计算微元的脉冲响应函数计算微元的脉冲响应函数的方法为：
[0032][0033]式中，φ
n
分别表示微元内第n散射体的幅度，多普勒频率和随机相位。
[0034]进一步地，本专利技术中，步骤四中，计算微元脉冲响应函数的二阶统计量的方法为：
[0035][0036]其中，E{
·
}表示统计期望运算，δ(
·
)表示狄利克雷函数，表示在时延为τ
′
和相关时间为Δt的微元脉冲响应函数。
[0037]进一步地，本专利技术中，步骤五中，二阶统计量的极限值为：
[0038][0039]进一步地，本专利技术中，步骤六中，计算散射区域的系统相关函数A
h
(τ,Δt)的方法为：
[0040][0041]其中，E{|a(α,η)|2}表示(α,η)位置的散射体的功率增益的期望，E{e
j(ρ(α,η))Δt
}表示(α,η)位置的散射体的多普勒相位函数的期望。
[0042]进一步地，本专利技术中，步骤七中，计算散射区域的多普勒谱S(ρ)的具体方法为：
[0043][0044]其中，dΩ表示任意散射体的面积微元；
[0045][0046]其中，r表示散射体距离换能器的斜距，r
‑
1≈r；
[0047]S(ρ)是散射区域所有散射体E{|a(α,η)|2}和E{e
j(ρ(α,η)Δt)
}乘积后累加的傅里叶变换获得；
[0048][0049]其中，参考距离r0已用1m替代，下角标T和R分别代表发射和接收过程，G
T
(α,η)表示发射过程的指向性，G
R
(α,η)表示接收过程的指向性，r
T1
和r
R1
分别表示发射时和接收时的散射体到换能器的斜距，β是声吸收系数；
[0050][0051]其中，r＝r
T1
≈r
R1
，r
‑
1≈r；
[0052]本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种声学多普勒长期测速精度估计方法，其特征在于，包括：步骤一、以换能器到散射体的矢量r与换能器的运动速度矢量v方向的夹角为α角，以换能器到散射体的矢量r在垂面投影与垂直方向的夹角为η角，建立(α,η)角度坐标系；步骤二、在(α,η)角度坐标系内，计算换能器波束照射的海底散射区域内所有散射体的指向性、散射强度、传播损失和多普勒相位；步骤三、将换能器接收到散射体反射信号时延相同的N个散射体作为一个微元，利用所述N个散射体的散射强度、传播损失和多普勒相位，计算微元的脉冲响应函数其中，τ为时延变量，t为时间变量，为微元的时延；步骤四、计算微元脉冲响应函数的二阶统计量其中，Δt为相关时间，τ'＝τ+Δτ，Δτ时间延迟差；步骤五、计算当微元中散射体个数N
→
∞时，二阶统计量的极限值；步骤六、利用所述二阶统计量的极限值，计算散射区域的系统相关函数A
h
(τ,Δt)；步骤七、利用所述散射区域的系统相关函数A
h
(τ,Δt)，计算散射区域的多普勒谱S(ρ)，再计算多普勒谱S(ρ)的一阶谱距M1；步骤八、计算多普勒谱S(ρ)的一阶谱距M1相对波束方向上的理论多普勒频率ω
d
的误差，利用所述误差建立声学多普勒长期测速精度的解析公式，获取声学多普勒长期测速精度值。2.根据权利要求1所述的一种声学多普勒长期测速精度估计方法，其特征在于，步骤二中，计算换能器能采集到的所有散射体的指向性的方法均相同，具体为：利用公式：计算获得指向性G(α,η)，其中，Δ1表示指向性G(α,η)在α角度方向上的等效波束开角，Δ2表示指向性G(α,η)在η角度方向上的等效波束开角，α0表示换能器散射波束中心角，其中，G(θ)表示声线与波束中心轴线夹角θ的指向性函数，α(x,z)表示用x和z变量表示的α变量，η(x,z)表示用x和z变量表示的η变量；所述x和z表示直角坐标系的两个坐标轴；η(x,z)＝arctan(z/H)式中，H表示换能器到达海底自身投影点的垂直距离。3.根据权利要求1所述的一种声学多普勒长期测速精度估计方法，其特征在于，步骤二中，计算换能器能采集到的所有散射体的散射强度的方法相同，具体为：
其中，表示散射体的强度，表示波束中心散射体的散射强度，表示散射体散射波的入射角，指波束中心处散射体的入射角，ο表示α相对波束中心角α0的差角，P
MSP
‑
in
、P
MSP
‑
sa...

【专利技术属性】
技术研发人员：曹忠义，孙大军，李雪松，张殿伦，勇俊，刘鑫，
申请(专利权)人：哈尔滨工程大学，
类型：发明
国别省市：