请给出多普勒fft的matlab代码

在MATLAB中，要计算多普勒频谱通常涉及到对连续信号进行快速傅里叶变换（FFT），然后处理频率轴上因多普勒效应引起的偏移。下面是一个简单的示例，假设我们有一个包含加速度信号的数组`signal`，并且已知目标的速度`v`（单位可能是米/秒），用于计算多普勒频移：

% 假设信号已经存在，例如从文件读取或者生成
signal = ...; % 你的信号数据

% 获取采样率和信号长度
fs = ...; % 采样率
T = length(signal)/fs; % 时间间隔

% 计算多普勒频移 (Hz)
f_doppler = v / T;

% 对信号进行FFT
fft_result = fft(signal);

% 计算频域范围
N = length(fft_result);
frequencies = (-NyquistFrequency:NyquistFrequency-N+1) * fs/N; % Nyquist Frequency = fs/2

% 应用多普勒频移到频率轴
shifted_frequencies = frequencies + f_doppler;

% 获取频谱并只保留正半部分（因为负频部分是对称的）
shifted_fft = fft_result(1:N/2+1) ./ N;
shifted_freqs_positive = shifted_frequencies(1:N/2+1);

% 绘制多普勒频谱
plot(shifted_freqs_positive, abs(shifted_fft).^2, 'b');
xlabel('Frequency [Hz]');
ylabel('Magnitude Spectrum');
title(['Doppler FFT with velocity ' num2str(v) ' m/s']);

% 如果有实际的多普勒峰值，你可以找到最大值的位置作为估计值
[~, idx] = max(abs(shifted_fft));
estimated_velocity = shifted_freqs_positive(idx) * T;

%