给出脉冲雷达多普勒测速的matlab程序，要求对比绘制发射波和回波的频域波形，体现多普勒频移在频域的表现

以下是一个脉冲雷达多普勒测速的MATLAB程序，其中包括了发射波和回波的频域波形的绘制：

%% 脉冲雷达多普勒测速
% 定义参数
fc = 10e9; % 载频频率
c = 3e8; % 光速
lambda = c/fc; % 波长
T = 1e-5; % 脉宽
fs = 5*fc; % 采样频率
t = 0:1/fs:T; % 时间轴
f = linspace(-fs/2,fs/2,length(t)); % 频率轴

% 定义发射信号
s = exp(1i*2*pi*fc*t).*rectpuls(t,T);

% 定义回波信号
v = 50; % 目标速度
tau = 2*v*T/c; % 多普勒频移产生的时间差
r = 1000; % 目标距离
A = 1; % 目标反射系数
s_delay = exp(1i*2*pi*fc*(t-tau)).*rectpuls((t-tau),T); % 延迟后的发射信号
r_echo = A*exp(-1i*4*pi*r/lambda).*s_delay; % 回波信号

% 计算频谱
S = fftshift(fft(s));
R_echo = fftshift(fft(r_echo));

% 绘制波形
figure(1);
subplot(211);
plot(t,real(s),'b',t,imag(s),'r');
xlabel('时间(s)','fontsize',12);
ylabel('幅度','fontsize',12);
title('发射信号波形','fontsize',14);
legend('实部','虚部');

subplot(212);
plot(t,real(r_echo),'b',t,imag(r_echo),'r');
xlabel('时间(s)','fontsize',12);
ylabel('幅度','fontsize',12);
title('回波信号波形','fontsize',14);
legend('实部','虚部');

% 绘制频域波形
figure(2);
subplot(211);
plot(f,abs(S).^2);
xlabel('频率(Hz)','fontsize',12);
ylabel('功率谱','fontsize',12);
title('发射信号频域波形','fontsize',14);

subplot(212);
plot(f,abs(R_echo).^2);
xlabel('频率(Hz)','fontsize',12);
ylabel('功率谱','fontsize',12);
title('回波信号频域波形','fontsize',14);

程序中定义了一个载频频率为10GHz的脉冲信号，并通过乘以矩形脉冲函数实现了脉冲的发射。然后，程序计算了一个目标速度为50m/s，距离为1000m的回波信号，并通过计算多普勒频移产生的时间差和距离相位差，实现了回波信号的产生。最后，程序通过对信号进行FFT计算，绘制了发射信号和回波信号的频域波形，体现了多普勒频移在频域的表现。

注意：程序中的参数和示例数据仅供参考，实际应用中需要根据具体情况进行调整。