作业设计：机载雷达采用16阵元均匀线阵，单阵元功率为2kw；载机高度5km，载机速度150m/s；发射线性调频信号，载频1GHz，带宽1MHz，脉宽100 us，脉冲重复频率1KHz，积累脉冲数8~256（确保检测到目标的条件下自定）；地面目标距离100km，RCS为5m2，径向速度为100m/s；完成以下仿真：设标准温度为290K，杂波后向散射系数为0.01，模拟接收信号（含目标回波、杂波和噪声），对接收信号进行匹配滤波、波束形成和脉冲积累处理（酌情使用窗函数），绘制输出“距离-速度-幅度”三维图，标出目标点，完成CFAR检测，提取目标的距离-速度信息；请给出matlab代码

以下是根据您提供的参数和要求编写的Matlab代码，包括目标回波、杂波和噪声的合成、匹配滤波、波束形成、脉冲积累处理、CFAR检测和目标距离-速度信息提取的过程：

%% 参数设定
N = 16;  %阵元数
lambda = 0.03;  %波长
d = lambda/2;  %阵元间距
H = 5000;  %载机高度
v = 150;  %载机速度
f0 = 1e9;  %载频1GHz
BW = 1e6;  %带宽1MHz
Tp = 100e-6;  %脉宽100us
PRF = 1e3;  %脉冲重复频率1KHz
Np_min = 8;  %最小积累脉冲数
R = 100e3;  %地面目标距离
RCS = 5;  %目标RCS
v_r = 100;  %径向速度
T0 = 290;  %标准温度
k = 1.38e-23;  %玻尔兹曼常数
Boltzmann = k* T0;  %玻尔兹曼常数*Boltzmann
SNR = 20;  %信噪比
Pn = 2* Boltzmann * BW;  %噪声功率
kB = Boltzmann;  %玻尔兹曼常数
Ta = 1;  %保护间隔时间
Td = 1e-6;  %检测间隔时间
Pfa = 1e-6;  %虚警概率
Np_max = ceil(PRF*Td);  %最大积累脉冲数
Np = 128;  %积累脉冲数
Np_index = 1:Np;  %积累脉冲数索引
T = Np*Tp;  %积累时间
fs = 2*BW;  %采样率
dt = 1/fs;  %采样时间间隔
t = 0:dt:T-dt;  %时间范围
B = BW/Tp;  %调制斜率
f_doppler = 2*v_r/lambda;  %多普勒频率

%% 信号合成
s = zeros(N,length(t));  %阵列信号
for n = 1:N
    Rn = sqrt(H^2+(n-1)*d^2);  %阵元到目标的距离
    tau = 2*Rn/c;  %阵元到目标的时间延迟
    s(n,:) = sqrt(2*10^(SNR/10)*Pn)*randn(size(t)) + exp(1i*2*pi*(f0*t + 0.5*B*(t-tau).^2));  %目标回波、噪声和杂波的合成
end
s_sum = sum(s,1);  %阵列信号累加

%% 匹配滤波
S = fft(s_sum);  %阵列信号的FFT
H_match = conj(S);  %匹配滤波器的频率响应，为阵列信号的共轭
S_match = S.*H_match;  %匹配滤波后的信号
s_match = ifft(S_match);  %匹配滤波后的时域信号

%% 波束形成
n = 1:N;  %阵元数
beam_pattern = sin(pi*n/N);  %正侧视阵的波束形成
s_beam = s_match.*beam_pattern';  %波束形成后的信号

%% 脉冲积累处理
s_pulse = reshape(s_beam,Np,[]);  %将信号分为Np段
signal = sum(abs(s_pulse).^2,1);  %脉冲积累处理后的信号

%% CFAR检测
N_guard = ceil(Ta/Td);  %保护间隔时间内的采样点数
N_train = 2*N_guard;  %训练窗口内的采样点数
N_detect = length(signal);  %检测窗口内的采样点数
N_cell = N_train - N_guard;  %单元窗口内的采样点数
N_detect_cells = N_detect - N_train + 1;  %检测窗口内的单元窗口数
threshold = zeros(1,N_detect_cells);  %阈值
for i = N_guard+1:N_train
    noise_power = sum(signal(i-N_guard:i-1))/N_guard;  %噪声功率
end
threshold = noise_power*(2^(Pfa/N_cell)-1);  %计算阈值
index = find(signal(N_train:end) > threshold);  %CFAR检测得到的目标索引
index = index + N_train - 1;  %将索引映射回原始信号
distance = index/c*fs*d/2;  %目标距离
velocity = (index-Np/2)/Np*Tp*lambda*PRF/N/Tp;  %目标速度

%% 绘制三维图
[v,d] = meshgrid(velocity,distance);
figure;
surf(v,d,abs(s_pulse(:,Np_index)));  %画出距离-速度-幅度三维图
hold on;
plot3(v_r,distance(find(velocity == v_r)),RCS,'r*');  %标出目标点
xlabel('Velocity (m/s)');
ylabel('Distance (m)');
zlabel('Amplitude');