设计警戒雷达系统的参数如下： 工作频率：L波段，频率范围1-2GHz，频率中心1.5GHz，对应波长0.2m 天线增益：54.2 dB 方向性系数：77.43 方位波束宽度：0.31° 俯仰波束宽度：30° 信号带宽B=2.38 MHz PRI：2.8 ms 脉冲宽度：0.42 us 脉冲压缩增益：2.4e6 扫过目标时间0.0086s 脉冲积累数1.853 多普勒滤波积累增益0.0035dB 雷达天线发射功率6千兆瓦 中心频率1.5Ghz，中频信号采样率4Ghz，基带数据采样率100Mhz 最大探测距离：216.8 km，给出其仿真代码

由于警戒雷达系统的参数较多，这里只能给出一部分仿真代码，供参考：

% 雷达参数
f0 = 1.5e9; % 中心频率
lambda = 0.2; % 波长
G = 10^(54.2/10); % 天线增益
K = 77.43; % 方向性系数
BW_h = 0.31; % 方位波束宽度
BW_v = 30; % 俯仰波束宽度
B = 2.38e6; % 信号带宽
PRI = 2.8e-3; % 脉冲重复间隔
T = 0.42e-6; % 脉冲宽度
Np = 1.853; % 脉冲积累数
Nd = 0.0035; % 多普勒滤波积累增益
Ptx = 6e9; % 雷达天线发射功率
T_scan = 0.0086; % 扫过目标时间

% 目标参数
R_max = 216.8e3; % 最大探测距离

% 仿真时间
T_sim = 10*PRI;

% 仿真步长
dt = 1/(4*f0);

% 生成基带信号
fs = 100e6; % 基带数据采样率
t = 0:dt:T_sim-dt;
f_IF = 4e9; % 中频信号采样率
t_IF = 0:1/f_IF:T_sim-1/f_IF;
x = sin(2*pi*10e6*t) + sin(2*pi*20e6*t); % 基带信号

% 生成发射信号
A = sqrt(2*Ptx*G*lambda^2/(4*pi)^3); % 发射信号振幅
fc = f0 + f_IF; % 发射信号中心频率
st = zeros(size(t));
for ii = 1:length(t)
    st(ii) = A*cos(2*pi*fc*t(ii) + 2*pi*B/T*(t(ii)-PRI/2)*rectpuls(t(ii)-PRI/2,T)); % 调制
end

% 天线扫描
theta_scan = (-BW_h/2:0.001:BW_h/2)*pi/180; % 方位角扫描范围
phi_scan = (-BW_v/2:0.1:BW_v/2)*pi/180; % 俯仰角扫描范围
P = zeros(length(phi_scan),length(theta_scan));
for ii = 1:length(phi_scan)
    for jj = 1:length(theta_scan)
        P(ii,jj) = beam_pattern(theta_scan(jj),phi_scan(ii),f0,lambda,G,K,BW_h,BW_v); % 计算天线指向系数
    end
end

% 雷达回波
R = zeros(length(phi_scan),length(theta_scan),length(t)); % 回波矩阵
for ii = 1:length(phi_scan)
    for jj = 1:length(theta_scan)
        r = sqrt(R_max^2 + (R_max*tan(phi_scan(ii)))^2 + (R_max*tan(theta_scan(jj)))^2); % 目标距离
        delay = 2*r/c; % 时延
        signal = st.*exp(-1i*2*pi*fc*(t-delay)); % 目标回波信号
        signal_rx = signal.*P(ii,jj)*cosd(90-phi_scan(ii))*cosd(theta_scan(jj)); % 经过天线接收的信号
        R(ii,jj,:) = signal_rx.*conj(signal_rx); % 回波功率
    end
end

% 处理回波信号
Nfft = 2^nextpow2(length(t)); % FFT点数
f = (-Nfft/2:Nfft/2-1)*fs/Nfft; % 频率轴
R_fft = fftshift(fft(R,Nfft,3),3); % FFT
R_range = abs(R_fft); % 距离FFT
R_doppler = abs(squeeze(sum(R_fft,2))); % 多普勒FFT

% 绘图
figure
imagesc(theta_scan*180/pi,phi_scan*180/pi,sum(R_range,3)) % 距离像
xlabel('方位角/度')
ylabel('俯仰角/度')
title('距离像')

figure
imagesc(theta_scan*180/pi,phi_scan*180/pi,R_doppler) % 多普勒像
xlabel('方位角/度')
ylabel('俯仰角/度')
title('多普勒像')