获取多普勒频谱的matlab代码

下面是一个简单的 MATLAB 代码，用于获取多普勒频谱：

% 设置参数
fs = 1000; % 采样率
f0 = 100; % 信号频率
v = 10; % 目标速度
c = 343; % 声速

% 生成信号
t = 0:1/fs:1; % 时间向量
x = sin(2*pi*f0*t); % 信号

% 生成回波信号
lambda = c/f0; % 波长
tau = 2*v/lambda; % 延迟时间
delay = round(tau*fs); % 延迟采样点数
y = [zeros(1,delay) x]; % 回波信号

% 计算多普勒频谱
nfft = 1024; % FFT点数
freq = (-fs/2:fs/nfft:fs/2-fs/nfft); % 频率向量
z = fftshift(fft(y,nfft)); % FFT
Doppler = abs(z).^2; % 多普勒频谱

% 绘制结果
figure;
plot(freq,Doppler);
xlabel('多普勒频率(Hz)');
ylabel('功率谱密度');

这个代码可以生成一个多普勒频谱图，其中 x 轴表示多普勒频率，y 轴表示功率谱密度。你可以根据需要修改代码中的参数，例如信号频率、目标速度、采样率等。