帮我写一个用MATLAB R2017来实现脉冲多普勒雷达产生雷达辐射和目标回波仿真信号的代码

好的，下面是MATLAB R2017版本的脉冲多普勒雷达产生雷达辐射和目标回波仿真信号的代码：

% 设置参数
fc = 10e9; % 雷达中心频率
c = 3e8; % 光速
lambda = c/fc; % 波长
bw = 2e6; % 带宽
prf = 1e3; % 脉冲重复频率
t = -1e-5:1/prf:1e-5; % 时间范围
n = length(t); % 采样点数
f = linspace(-bw/2,bw/2,n); % 频率范围

% 产生脉冲信号
pulse = rectpuls(t,2e-6);
s_pulse = fft(pulse);

% 产生高斯白噪声
noise_power = 1e-10;
noise = sqrt(noise_power)*randn(size(pulse));

% 产生回波信号
range = 5000; % 目标距离
target_rcs = 1; % 目标雷达反射截面积
target_amp = sqrt(target_rcs)*exp(-1j*4*pi*fc*range/c);
echo = conv(pulse,target_amp,'same') + noise;
s_echo = fft(echo);

% 绘制时域和频域脉冲信号和回波信号
figure(1);
subplot(2,2,1);
plot(t,pulse);
title('Pulse Signal');
xlabel('Time (s)');
ylabel('Amplitude');

subplot(2,2,2);
plot(f,abs(fftshift(s_pulse)));
title('Pulse Signal Spectrum');
xlabel('Frequency (Hz)');
ylabel('Amplitude');

subplot(2,2,3);
plot(t,real(echo));
title('Echo Signal');
xlabel('Time (s)');
ylabel('Amplitude');

subplot(2,2,4);
plot(f,abs(fftshift(s_echo)));
title('Echo Signal Spectrum');
xlabel('Frequency (Hz)');
ylabel('Amplitude');

这段代码会产生一个脉冲多普勒雷达的信号，并产生一个距离为5000米、反射截面积为1平方米的目标回波信号，最后绘制出时域和频域的脉冲信号和回波信号。你可以根据需要修改参数来模拟不同的场景。