airborne_radar_clutter_simulation.m

%机载雷达地杂波模拟
%V1.0 By 牛志永2018-5-12
%本程序完全按照现实场景模拟机载雷达录取数据的过程，目前考虑一个接收阵元
%在该程序框架下，可以非常方便地扩展到支持线阵，面阵场景
clc;
close all;
clear all;
%雷达参数
BW = 10e6;%信号带宽
Fs = 1*BW;
C  = 3e8;
fc = 1e9;
lambda = C/fc;
PRF = 100;
PRT = 1/PRF;
Pulse_num = 100;
%快时间维距离门数在第3部分获得
dr = C/2*1/Fs;
% 平台及环境参数
Va = 200;%平台速度，沿X轴正向飞行
H  = 5000;%平台高度，
X0 = 0;  %0时刻，平台坐标（00H）
azimuth_angle = 80; %波束方位角，波束照射方向在X-Y平面的投影与X轴正向的夹角
depression_angle = 10; %波束俯视角，波束照射方向与水平面的的夹角，也被称为擦地角
azimuth_angle_rad = azimuth_angle/180*pi;
depression_angle_rad = depression_angle/180*pi;
beam_width = 2.5; %波束宽度，圆锥波束。
beam_width_rad = beam_width/180*pi; %波束宽度，单位为弧度圆锥波束。
R_center = H/sin（depression_angle/180*pi）;%波束照射中心点到0时刻雷达距离
X_center = R_center * cos（depression_angle/180*pi） *cos（azimuth_angle/180*pi）;%波束照射中心对应的X轴坐标
Y_center = R_center * cos（depression_angle/180*pi） *sin（azimuth_angle/180*pi）;%波束照射中心对应的Y轴坐标
Z_center = 0;

%% 1、生成待探测目标区域网格认为波束照射中心即为网格中心。
net_Nr = 1000;%Y
net_Na = 1000;%X
dest_net = zeros（net_Nrnet_Na）;
net_L = C/（2*BW）;%网格长度
radar_coordinate = [X0 0 H]; %雷达坐标
% %确定被照射到的网格
irradiated_table=zeros（net_Nrnet_Na）;
[irradiated_table r_range]= irradiate_region（radar_coordinateazimuth_angle_raddepression_angle_rad...
                   beam_width_rad[net_Nr net_Na net_L X_center Y_center]1）;
% figure;imagesc（X_center+（（0:net_Na-1）-round（net_Na/2））*net_LY_center+（（0:net_Nr-1）-round（net_Nr/2））*net_Labs（irradiated_table））;
% xlabel（‘X‘）;ylabel（‘Y‘）;
        RI = imref2d（size（irradiated_table））;
        RI.XWorldLimits = X_center + （[0 net_Na-1]-round（net_Na/2）） *net_L;
        RI.YWorldLimits = Y_center + （[0 net_Nr-1]-round（net_Nr/2）） *net_L;;
figure;imshow（abs（irradiated_table）RI）;
% figure;mesh（abs（irradiated_table））;
%% 2、确定网格散射