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

好的，以下是一个简单的MATLAB R2017代码，用于模拟脉冲多普勒雷达的辐射和目标回波信号：

% 雷达参数设置
fc = 10e9;        % 中心频率
bw = 1e9;         % 带宽
prf = 5e3;        % 脉冲重复频率
pulse_width = 10e-6; % 脉冲宽度
tx_power = 1;     % 发射功率
antenna_gain = 30; % 天线增益
losses = 2;       % 系统损耗

% 目标参数设置
target_rcs = 1;   % 回波截面
target_range = 1000; % 距离
target_velocity = 100; % 速度

% 确定采样时间和采样点数
sampling_time = 1/prf;
sampling_points = round(pulse_width*fc);

% 产生脉冲信号
pulse = cos(pi*(0:sampling_points-1)/(sampling_points-1));
pulse = sqrt(tx_power)*pulse/max(abs(pulse)); 

% 产生频率步进信号
freq_step = bw/pulse_width;
freqs = fc + (-bw/2:freq_step:bw/2);

% 产生回波信号
target_delay = 2*target_range/(3e8);
target_doppler = 2*target_velocity/fc;
target_phase_shift = exp(1i*2*pi*target_doppler*freqs*target_delay);
echo = target_rcs*target_phase_shift;

% 产生接收信号
rx_signal = zeros(1,length(pulse)+length(echo)-1);
for i = 1:length(freqs)
    rx_signal = rx_signal + ...
        fftfilt(pulse.*exp(1i*2*pi*freqs(i)*sampling_time*(0:length(pulse)-1)),echo(i)*pulse);
end

% 添加噪声
snr_db = 20;
snr_linear = 10^(snr_db/10);
noise_power = norm(rx_signal)^2/(length(rx_signal)*snr_linear);
noise = sqrt(noise_power/2)*(randn(1,length(rx_signal))+1i*randn(1,length(rx_signal)));
rx_signal = rx_signal + noise;

% 计算功率
tx_power_dbm = 10*log10(tx_power)+30;
rx_power_dbm = 10*log10(mean(abs(rx_signal).^2))+30;
rx_power_dbm = rx_power_dbm + antenna_gain - losses;
fprintf('发射功率（dBm）：%f, 接收功率（dBm）：%f\n',tx_power_dbm,rx_power_dbm);

% 可视化结果
figure;
subplot(211);plot(abs(pulse));title('脉冲信号');
subplot(212);plot(abs(rx_signal));title('接收信号');

这个代码主要分为以下几部分：

1. 雷达参数设置：设置了雷达的中心频率、带宽、脉冲重复频率、脉冲宽度、发射功率、天线增益和系统损耗等参数；
2. 目标参数设置：设置了目标的回波截面、距离和速度等参数；
3. 产生脉冲信号：根据采样时间和采样点数，产生了一个宽度为脉冲宽度的脉冲信号；
4. 产生频率步进信号：根据带宽和脉冲宽度，产生了一个频率步进信号；
5. 产生回波信号：根据目标的距离、速度和回波截面，产生了一个回波信号；
6. 产生接收信号：通过对脉冲信号和回波信号进行滤波和相乘，得到了接收信号；
7. 添加噪声：根据信噪比，添加了一定程度的高斯白噪声；
8. 计算功率：计算了发射功率和接收功率，并输出结果；
9. 可视化结果：将脉冲信号和接收信号可视化。

这个代码只是一个简单的示例，实际应用中可能需要更加复杂的处理，比如考虑多径效应、杂波、多目标等问题。