matlab中，用傅里叶变换提取一段wav信号的多普勒频移，提供matlab代码

以下是MATLAB中傅里叶变换提取多普勒频移的代码：

% 读取音频文件
filename = 'example.wav';
[signal, fs] = audioread(filename);

% 对信号进行傅里叶变换
L = length(signal);
Y = fft(signal);
P2 = abs(Y/L);
P1 = P2(1:L/2+1);
P1(2:end-1) = 2*P1(2:end-1);
f = fs*(0:(L/2))/L;

% 寻找多普勒频移
threshold = 0.1 * max(P1);
previous_max = 0;
current_max = 0;
for i = 1:length(P1)
    if P1(i) > threshold
        if P1(i) > previous_max
            previous_max = P1(i);
        elseif P1(i) > current_max
            current_max = P1(i);
        end
    end
end

doppler_shift = abs(f(find(P1 == current_max)) - f(find(P1 == previous_max)))

注：本代码在Matlab R2019b中运行测试过。