%機(jī)載雷達(dá)地雜波模擬
%V1.0?By?牛志永2018-5-12
%本程序完全按照現(xiàn)實場景模擬機(jī)載雷達(dá)錄取數(shù)據(jù)的過程，目前考慮一個接收陣元
%在該程序框架下，可以非常方便地擴(kuò)展到支持線陣，面陣場景
clc;
close?all;
clear?all;
%雷達(dá)參數(shù)
BW?=?10e6;%信號帶寬
Fs?=?1*BW;
C??=?3e8;
fc?=?1e9;
lambda?=?C/fc;
PRF?=?100;
PRT?=?1/PRF;
Pulse_num?=?100;
%快時間維距離門數(shù)在第3部分獲得
dr?=?C/2*1/Fs;
%?平臺及環(huán)境參數(shù)
Va?=?200;%平臺速度，沿X軸正向飛行
H??=?5000;%平臺高度，
X0?=?0;??%0時刻，平臺坐標(biāo)（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時刻雷達(dá)距離
X_center?=?R_center?*?cos（depression_angle/180*pi）?*cos（azimuth_angle/180*pi）;%波束照射中心對應(yīng)的X軸坐標(biāo)
Y_center?=?R_center?*?cos（depression_angle/180*pi）?*sin（azimuth_angle/180*pi）;%波束照射中心對應(yīng)的Y軸坐標(biāo)
Z_center?=?0;

%%?1、生成待探測目標(biāo)區(qū)域網(wǎng)格認(rèn)為波束照射中心即為網(wǎng)格中心。
net_Nr?=?1000;%Y
net_Na?=?1000;%X
dest_net?=?zeros（net_Nrnet_Na）;
net_L?=?C/（2*BW）;%網(wǎng)格長度
radar_coordinate?=?[X0?0?H];?%雷達(dá)坐標(biāo)
%?%確定被照射到的網(wǎng)格
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、確定網(wǎng)格散射