基于频移迭代的多普勒谱中心频率估计方法
【技术领域】
[0001] 本发明属于雷达信号处理领域，涉及一种基于频移迭代的多普勒谱中心频率估计 方法。本发明适用于各种相干机制的岸基微波多普勒雷达系统和高频地波雷达。
【背景技术】
[0002] 海洋观测是人类认识海洋自然属性和环境特征的基本手段，是海洋事业发展的基 石。近年来，美、英、法、德、俄、挪威、新西兰等国家都在积极开展岸基雷达海洋遥感技术的 研究，从而获取无向浪高谱、海浪统计（如有效浪高、浪周期、浪向）及海洋表面流等海洋动 力学要素。
[0003] 随着无线电海洋探测技术的发展，岸基微波多普勒雷达已被广泛应用于海洋表面 动力学参数遥测。该雷达通过获取照度区内的回波多普勒谱，估计其中心频率得到海面水 质点的径向速度，然后利用径向速度与海面浪高的转换关系，直接测量海浪参数如有效浪 高、平均浪周期等。该方法无需校准即可得到准确的海浪谱等结果，是一种"直接"测量海 浪的方法，而其中准确估计海洋回波多普勒谱中心频率是准确计算海浪谱的前提。因此，如 何有效估计多普勒谱的中心频率是微波多普勒雷达探测海浪的重要问题，其估计精度直接 决定了雷达海浪探测的准确性。
[0004] 从原理上说，频率是指在单位时间内完整波动的次数。在很多实际应用中，都可以 把运动中的传输波简化为正弦信号，在某一固定点的波质点在单位时间内走过的波动的个 数，就是波动频率。传统意义上的频率是指在一段时间上信号的傅里叶频率，由于时间参数 在傅里叶变换中已被去除，所以傅里叶频率与时间无关。与之相反，瞬时频率却指示的是信 号在某一时间点或一小段时间上的频率，其值是时间的函数。如果信号的相关函数存在，那 么一定时间内的傅里叶频率和瞬时频率的均值就完全一致。雷达回波的多普勒谱是某一时 间段内信号的频谱，其中心频率就是多普勒获取时间内的均值。
[0005] 目前，多普勒谱中心频率估计方法可分为时域和频域两大类。其中时域上具有代 表性的是协方差矩估计法。该方法利用相关函数来估计谱中心频率，其不足之处在于中心 频率的估计偏差会随谱歪斜程度增大而增大，且要求采样周期尽量短。在频域上具有代表 性的是谱矩法，该方法是将多普勒谱的能量中心所在位置作为多普勒谱中心频率。由于谱 矩法计算量小且使用方便，它已成为当前最广泛的中心频率估计方法。
[0006] 然而，实际接收的雷达回波无可避免受到噪声"污染"，较小的信噪比降低了中心 频率的估计精度。为了克服这一难题，国内外学者做了大量的研究工作，目前，对于多普勒 谱中心频率的估计已经取得了一定的进展，但是估计效果依然没有达到最好，主要是受限 于如何确定噪声水平，并将信号和噪声分离，再提取谱参数。
【发明内容】
[0007] 本发明利用了频移迭代的方法，通过多次频移估计直至收敛从而得到中心频率， 可以有效地提高多普勒谱中心频率估计的精度。
[0008] 本发明的目的在于：基于实际的岸基多普勒雷达系统，提供一种具有更高估计精 度的基于频移迭代的多普勒谱中心频率估计方法，从而更加准确地获得海洋动力学参数， 如有效波高和平均浪周期等。
[0009] 为实现上述目的，本发明提供的中心频率估计方法如下：
[0010] -种基于频移迭代的多普勒谱中心频率估计方法，包括以下步骤：
[0011] 步骤1，通过实测多普勒谱初步确定噪声基底，并将该噪声基底作为门限，利用截 断法实现信号和噪声的初步分离；转步骤2 ;
[0012] 步骤2,对截断后剩余的多普勒谱采用谱矩法第一次估计多普勒频移；转步骤3 ;
[0013] 步骤3,以该多普勒频移为对称中心，两边同时增加相同的带宽构成新的积分区 间，对新积分区间内的多普勒谱再次利用谱矩法估计多普勒频移；转步骤4 ;
[0014] 步骤4,重复步骤3,经过数次迭代后，进行收敛性判断，得到多普勒频移的收敛结 果，并将该结果作为多普勒谱中心频率的最佳估计值。
[0015] 步骤1所述的截断法实现信号和噪声的初步分离的具体方法包括以下子步骤：
[0016] 步骤1. 1，在实测多普勒谱的两边各选取5%的点作为左右噪声区间，分别计算左 右噪声区间的噪声平均值NoiselefJP Noise "ght，并将两者之中的较大值作为噪声区间的门 限 Noisethrehold，艮P Noisethrehold= max (Noise left，Noiseright);
[0017] 步骤I. 2,从多普勒谱最大幅值开始，向左依次查找，直到多普勒谱幅值等于门限 值，将该幅值所对应的频率作为信号区间的左边界同理，从多普勒谱最大幅值向右查 找获得信号区间的右边界 fright ' 且有 fleft< f right;
[0018] 步骤I. 3,记积分区间B。= [f left，f"ght]，将该区间作为回波多普勒谱的信号区间。
[0019] 步骤2与步骤3中所述的谱矩法估计多普勒频移是依照如下公式实现的：
[0021] 其中，匕表示多普勒谱的频点，S(f J表示多普勒谱每个频点所对应的幅值，仁表 示估计的多普勒谱频移，η表示迭代次数。
[0022] 所述步骤3中，进行第一次迭代时是以第一次利用谱矩法得到的初始频移f。为对 称中心，两边同时增加相同的带宽ΛΒ。，且ΔΒ。满足AB Q=min(|fQ_flrft|，|f；lght-f Q|)，此 时的积分区间变为B1= [f ^ΔΒ。，&+ΔΒ。];所得到的新的积分区间队内再次使用谱矩法估 计多普勒谱频移可得到经过第一次迭代后的频移A，接着以A为对称中心，两边同时增加 相同的带宽AB 1，其中从而得到新的积分区间B2;在积分 区间B2内再次使用谱矩法估计多普勒谱频移得到第二次迭代的频移f 2，以此类推，一直重 复η次得到第η次迭代后的频移fn。
[0023] 所述步骤4中，经过了 η次迭代，得到了第η次谱矩法估计频移后，进行收敛性判 断，若满足fn= fn :，则将fn作为该多普勒谱中心频率的最佳估计值；若不满足，重复步骤3 直至两者相等。
[0024] 因此，本发明具有如下优点：
[0025] 1.该估计方法由于进行多次迭代直至频移收敛，因此在第一次确定积分边界的时 候无需过于严格，从而降低了初始边界对于多普勒谱中心频率估计的影响。
[0026] 2.该估计方法可提高多普勒谱中心频率估计的精度，为后续处理获取无向浪高 谱、海浪统计参数（如有效波高、浪周期、浪向）及海洋表面流等海洋动力学要素提供了更 加准确的信息。
[0027] 3.该估计方法可应用于多种不同谱形的中心频率估计，应用范围广，实用性强。
【附图说明】
[0028] 图1是本发明的算法流程图。
[0029] 图2是基于实测数据的微波雷达回波多普勒谱。
[0030] 图3是本发明基于实测数据所估计的频移的收敛过程。
[0031] 图4是在雷达信号处理过程中其他部分处理方法完全相同的情况下，运用谱矩法 与本发明所得到的中心频率估计值对比图。
【具体实施方式】
[0032] 特定距离元的雷达回波是该距离元内多种随机波动海浪反射电磁波的结果。雷达 接收到来自某一特定距离元的