【水声通信】基于Matlab的水声信号Doppler谱分析与目标速度估计【附源码】

1. 从“火车汽笛”到“水下听诊”：理解水声多普勒效应

想象一下，你站在铁路边，一列火车鸣着汽笛呼啸而过。当火车朝你开来时，汽笛声听起来又尖又高；而当它远离你时，声音立刻变得低沉。这个我们生活中都体验过的现象，就是著名的多普勒效应。声音的频率（也就是我们感觉到的音调高低）会随着声源和听者之间的相对运动而改变。朝你运动，频率变高，声音变尖；离你而去，频率变低，声音变沉。

现在，让我们把这个场景搬到水下。水下世界并非一片寂静，充满了各种声音：潜艇的螺旋桨声、鱼群的游动声、海底地质活动声，甚至是海洋哺乳动物的“歌声”。水声通信和探测，本质上就是利用声波在水下“看”和“听”的技术。当我们需要知道一个水下目标（比如一艘潜航器或一群鱼）的速度时，上面提到的“火车汽笛”原理就派上了大用场，这就是水声信号的多普勒谱分析。

简单来说，我们向水中发射一个已知频率的声波信号（就像一声固定的哨音），这个信号遇到运动的目标后会被反射回来。由于目标在运动，反射回来的声波频率就发生了变化，产生了多普勒频移。通过精确测量这个频移量，我们就能像侦探一样，反推出目标相对于我们的运动速度。这就像是给水下世界装上了一台“速度雷达”，只不过用的是声波而不是电磁波。

我刚开始接触这个领域时，觉得理论很抽象。后来在海上实验，亲眼看到屏幕上一条清晰的频移谱线随着目标潜艇的加速而移动，瞬间就明白了多普勒分析的威力。它不仅是海洋科学研究中追踪洋流、监测生物迁徙的利器，在工程应用上，比如水下导航、避碰、资源勘探等领域更是不可或缺的核心技术。今天，我就带你用Matlab这把“瑞士军刀”，亲手实现一次完整的水声信号多普勒分析，从原理到代码，一步步算出水下目标的速度。

2. 动手前的准备：理解核心公式与Matlab工具箱

在打开Matlab写第一行代码之前，我们得先把核心的“换算公式”搞清楚。多普勒速度估计的基石是下面这个公式：

fd = (2 * v * f0 * cosθ) / c

看起来有点复杂？别怕，我们拆开看：

• fd：这就是我们最终要从信号里找出来的多普勒频移，单位是赫兹（Hz）。简单理解，它就是回声频率比发射频率高了多少或者低了多少。
• v：这是我们要求解的目标径向速度，也就是目标沿着声波方向运动的速度分量。单位通常是米/秒（m/s）。如果目标正对着我们开来，v是正的；远离我们，v就是负的。
• f0：这是我们发射的声波信号的原始频率，单位也是赫兹（Hz）。这是已知的，比如一个10kHz的声呐脉冲。
• θ：目标运动方向与声波传播方向之间的夹角。这是个关键参数。如果目标正好沿着声波方向运动（迎面而来或背向而去），cosθ等于1或-1，此时频移最大。如果目标是横着跑（垂直方向），cosθ=0，那就没有多普勒频移，我们也就测不到它的径向速度了。
• c：水中的声速，单位米/秒（m/s）。这不是一个固定值！它会随着海水的温度、盐度和深度（压力）变化，通常在1450米/秒到1550米/秒之间浮动。在实际应用中，用一个粗略的常数值（比如1500 m/s）会带来误差，所以高精度的测量需要根据实测水文数据对声速进行修正。

我们的任务流程很清晰：采集到回声信号 -> 分析信号频谱，揪出那个微小的频率变化fd -> 把fd、已知的f0、c和估算的θ代入上面的公式 -> 解算出目标速度v。

接下来是工具准备。Matlab之所以是信号处理的神器，是因为它有强大的内置函数和工具箱。我们这次用到的核心工具是信号处理工具箱。你可以用 ver 命令查看是否已安装。我们主要会用到以下几个函数：

• fft：快速傅里叶变换的“主将”，负责把时域信号（幅度随时间变化）转换成频域信号（不同频率成分的强度），这是我们观察频谱和寻找频移的必经之路。
• pwelch：用于估计信号的功率谱密度。相比于对整段信号做一次FFT，它采用分段加窗平均的方法，能得到更平滑、更稳定的频谱图，特别适合分析像水声信号这种带有噪声的数据。
• findpeaks：频谱图上找到峰值点的“探测器”。在多普勒谱分析中，回声信号的能量会集中在“发射频率+多普勒频移”处形成一个尖峰，这个函数能帮我们自动定位这个峰值的频率，非常方便。
• 基本的绘图函数如plot、xlabel、title等，用于可视化我们的每一步结果。

我建议你在动手编码前，先在Matlab命令窗口里简单试试这几个函数，比如对一个正弦波做FFT看看它的频谱线在哪里，熟悉一下它们的输入输出格式。磨刀不误砍柴工，理解这些工具，后面的路会顺畅很多。

3. 实战第一步：模拟生成带多普勒效应的水声信号

在真实的海上试验之前，我们通常先用计算机模拟一个信号来练手和验证算法。这能让我们在完全可控的环境下，知道“标准答案”应该是什么。我们来模拟一个最常见的情况：发射一个单频脉冲信号，然后一个目标以恒定速度运动，产生回声。

首先，我们设定一些基本参数，这些参数模拟了一个典型的中频水声通信场景：

% 基本参数设置
fs = 100000;          % 采样频率 100 kHz，决定了我们能分析的最高频率
T = 0.1;              % 信号持续时间 0.1秒
t = 0:1/fs:T-1/fs;    % 时间向量
f0 = 10000;           % 发射信号频率 10 kHz
c = 1500;             % 假设水中声速 1500 m/s
v_target = 5;         % 目标径向速度 5 m/s (正值为靠近)
theta = 0;            % 目标运动方向与声波方向夹角0度（迎面而来）

接下来，根据多普勒公式，我们先计算出理论上的多普勒频移 fd：

fd = (2 * v_target * f0 * cosd(theta)) / c; % cosd用于角度输入
f_echo = f0