matlab中计算一段wav信号的多普勒频移，提供matlab代码

以下是matlab代码示例：

% 读取wav文件
[x, Fs] = audioread('audio.wav');
% 配置多普勒频移分析参数
DopplerFFTSize = 2048; % FFT分析点数
DopplerWindow = hamming(DopplerFFTSize); % 窗函数
DopplerOverlap = DopplerFFTSize / 2; % 窗口重叠
DopplerMaxShift = Fs / 2; % 最大多普勒频移
% 计算多普勒频移
[s, f, t] = spectrogram(x, DopplerWindow, DopplerOverlap, DopplerFFTSize, Fs);
[v, i] = max(abs(s)); % 找到最大值点
dopplerFrequency = (i - 1) * Fs / DopplerFFTSize; % 计算多普勒频移
% 输出结果
fprintf('多普勒频移为 %fHz\n', dopplerFrequency);