给出脉冲雷达多普勒测速的matlab程序，要求显示目标速度并绘制发射波和回波的频域波形

以下是一个简单的脉冲雷达多普勒测速的MATLAB程序示例：

% 生成发送信号
fs = 1000; % 采样频率
T = 1/fs; % 采样时间
fc = 100; % 载波频率
t = 0:T:1; % 时间序列
x = cos(2*pi*fc*t) + 0.5*cos(2*pi*2*fc*t); % 发送信号

% 生成回波信号
v = 10; % 目标速度
f0 = 2*fc*v/3e8; % 多普勒频移
tau = 2*v/3e8; % 时延
y = [zeros(1,round(tau/T)),cos(2*pi*(fc+f0)*t) + 0.5*cos(2*pi*(2*fc+2*f0)*t)]; % 回波信号

% 频谱分析
N = length(x); % 采样点数
X = fftshift(fft(x)); % 发送信号频谱
Y = fftshift(fft(y)); % 回波信号频谱
f = (-N/2:N/2-1)*fs/N; % 频率序列

% 显示结果
subplot(3,1,1)
plot(t,x)
title('发送信号')
subplot(3,1,2)
plot(t,y)
title('回波信号')
subplot(3,1,3)
plot(f,abs(Y))
title('回波信号频谱')
xlabel('频率（Hz）')

运行程序后，会在图形界面中显示发送信号、回波信号以及回波信号的频谱。其中，目标速度为$v=10$，对应的多普勒频移为$f_0=2fcv/c=2\times100\times10/3\times10^8=0.0067\text{ Hz}$，时延为$\tau=2v/c=2\times10/3\times10^8=6.67\times10^{-8}\text{ s}$。