一种多普勒频偏测量方法.pdf

(19)中华人民共和国国家知识产权局(12)发明专利申请(10)申请公布号CN109405955A(43)申请公布日2019.03.01(21)申请号201811217212.X(22)申请日2018.10.18(71)申请人哈尔滨工程大学地址150001黑龙江省哈尔滨市南岗区南通大街145号哈尔滨工程大学科技处知识产权办公室(72)发明人殷敬伟葛威杨光韩笑(51)Int.Cl.G01H3/04(2006.01)权利要求书1页说明书4页附图3页(54)发明名称一种多普勒频偏测量方法(57)摘要本发明属于水声信号处理领域，具体涉及一种多普勒频偏测量方法。本发明属于水声信号处理领域。涉及一种水声信道中的多普勒频偏测量方法。本发明中的多普勒频偏测量方法是在传统的FFT运算基础上进行的。先对接收的信号进行FFT运算，得到的结果仅为频率的粗略估计值。为进一步提高估计精度，在粗略估计值周围进行插值迭代运算，迭代两次即可得到频率的精细估计值。这样，将频率估计结果与原信号频率相减便可得到高精度的多普勒频偏测量值。本发明的优点在于估计性能受多途干扰影响小，能在低信噪比下进行多普勒频偏估计；频率估计精度高，理论上估计方差接近于克拉美罗界。CN109405955ACN109405955A权利要求书1/1页1.一种多普勒频偏测量方法，其特征在于，该方法包括以下步骤：(1)发射端发射单频信号；(2)接收端将采集到的信号进行粗略频率估计；(3)利用粗略频率估计结果进行精细频率估计；(4)将接收信号频率估计值与发射信号频率估计值相减得到多普勒频偏测量值。2.根据权利要求1所述的一种多普勒频偏测量方法，其特征在于：首先发射端发射单频信号，接收端接收信号并采集完毕，在所述步骤(2)中的粗略频率估计时设接收信号为s(k)，采样点数为N，对输入信号作N点FFT计算，并寻找计算结果的模值平方的最大值即：其中第一参数k表示第k个采样点。3.根据权利要求2所述的一种多普勒频偏测量方法，其特征在于：所述步骤(3)中的频率精细估计时设定量化频率偏差初值设定第二参数i从1到2进行迭代计算：其中，k表示第k个采样点，s(k)表示第k个采样点的接收信号，第i次迭代后量化频率偏差为粗略频率估计模值平方的最大值，p为第三参数，β为第四参数，β按下式计算：得到接收信号s(k)的频率估计值：其中，fs为接收信号的采样频率，N为采样点数，为粗估计频率模值平方的最大值，为第二次迭代后量化频率偏差值。4.根据权利要求3所述的一种多普勒频偏测量方法，其特征在于：所述步骤(4)得到多普勒频偏测量值Δf为：其中，f0为发射信号的频率，为接收信号s(k)的频率估计值。5.根据权利要求1所述的一种多普勒频偏测量方法，其特征在于：所述的步骤(2)粗略频率估计通过对接收信号进行FFT计算实现。6.根据权利要求1所述的多普勒频偏测量方法，其特征在于：所述步骤(3)精细频率估计通过在粗略频率估计值周围进行迭代插值实现。2CN109405955A说明书1/4页一种多普勒频偏测量方法技术领域[0001]本发明属于水声信号处理领域，具体涉及一种多普勒频偏测量方法。背景技术[0002]近些年来，水下移动节点，如水下无人潜器，间的声学通信正在受到更多的关注。由于水中声信号的传播速度小，约1500m/s，与其它通信系统相比，水声通信更易受多普勒效应的影响。即使接收机/发射机之间的慢速运动或者固有的电流波运动都会对发射信号产生显著的多普勒效应。目前水声通信中最常用的抑制多普勒效应的方法分两步进行：第一步先对接收信号中的多普勒频偏进行估计，计算出多普勒因子；第二步利用估计出的多普勒因子对接收信号进行重采样，以此达到多普勒补偿的目的。此方法中，只有频偏估计的足够准确，多普勒补偿才能达到理想的效果。因此，多普勒频偏估计在抑制多普勒效应的过程中是十分关键的。[0003]频率估计法是多普勒频偏估计的重要方法。这种方法通过发射单频信号，在接收端对接收到的信号进行频率估计，接收信号频率与发射信号频率之差即为多普勒频偏。然而，在复杂的水声信道环境中，频率估计的性能易受噪声和多途干扰的影响，从而影响了多普勒频偏估计的精度。发明内容[0004]本发明的目的在于提供一种多普勒频偏测量方法。[0005]本发明的目的是这样实现的：[0006]一种多普勒频偏测量方法，包括以下步骤：[0007](1)发射端发射单频信号；[0008](2)接收端将采集到的信号进行粗略频率估计；[0009](3)利用粗略频率估计结果进行精细频率估计；[0010](4)将接收信号频率估计值与发射信号频率估计值相减得到多普勒频偏测量值。[0011]首先发射端发射单频信号，接收端接收信号并采集完毕，在所述步骤(2)中的粗略频率估计时设接收信号为s(k)，采样点数为N，对输入信号作