matlab中提取wav文件的多普勒频移，提供matlab代码

可以使用matlab的fftn函数对wav文件进行FFT变换，得到频谱图后通过计算频率偏移即可获得多普勒频移。

以下是一个简单的代码示例：

% 读入wav文件
[y,fs] = audioread('example.wav');

% 取其中一段信号进行处理，这里假设信号为50s到60s
windowSize = 10 * fs;
startIndex = 50 * fs + 1;
endIndex = startIndex + windowSize - 1;
y = y(startIndex:endIndex);

% 进行FFT
N = length(y);
Y = fft(y, N);
f = (0:N-1) * fs/N;
Y = Y .* conj(Y) / N;
Y = 20*log10(abs(Y));

% 找到频谱峰值，假设多普勒频移小于1kHz
[~, maxIndex] = max(Y(1:N/2));
if f(maxIndex) < 1000
    dopplerShift = f(maxIndex);
else
    dopplerShift = 0;
end

disp(['多普勒频移为：', num2str(dopplerShift), ' Hz']);

需要注意的是，这个代码示例只是提供了一个简单的实现方法，实际上需要根据具体的需求进行调整和优化。