clear?all;
clc;
C=3e8;???????????????%光速
ima=sqrt（-1）;????????%虛數(shù)單位
%%場景中心參數(shù)
R0=2e4;??????????????%場景中心斜距
%%雷達參數(shù)
fc=10e9;?????????????%工作頻率10GHz
lamda=C/fc;??????????%波長
SNR=-10;
%%天線參數(shù)
D=4;?????????????????%天線孔徑長度
Lsar=lamda*R0/D;?????%雷達合成孔徑長度
%%距離向參數(shù)
Tr=10e-6;????????????%發(fā)射脈沖時寬
B=180e6;?????????????%帶寬
Kr=B/Tr;?????????????%調(diào)頻率
Fsr=200e6;???????????%距離采樣率
nrn_Tr=Fsr*Tr;???????%脈沖區(qū)間采樣數(shù)
nrn=2048;????????????%距離向采樣數(shù)
tnrn=2*R0/C+[-nrn/2:nrn/2-1]‘/Fsr;??%距離向快時間
%%方位向參數(shù)
Fa=100;??????????????%方位采樣率
PRF=500;?????????????%脈沖重復頻率
V=100;???????????????%飛機速度
ar=5;????????????????%飛機加速度
nan=512;?????????????%方位向采樣數(shù)
tnan=[-nan/2:nan/2-1]/PRF;??%方位向慢時間
%%分辨率參數(shù)
DY=C/2/B;???????????%距離向分辨率
DX=DY;??????????????%為保持距離分辨率和方位分辨率相同，圖像比例一致
wt=lamda/（2*DX*（nan/PRF））;???%方位向分辨率計算出應該轉(zhuǎn)多大的角度wt，wt等于角速度*積累時間
%%目標參數(shù)，飛機模樣
Ntarget=11;???????	%點目標個數(shù)
deltaX=15;
deltaY=15;
points=zeros（Ntarget3）;??????????????
points（:1）=[2*deltaX;deltaX;deltaX;0;0;0;0;0;-deltaX;-deltaX;-2*deltaX];
points（:2）=[0;0;-2*deltaY;2*deltaY;deltaY;0;-deltaY;-2*deltaY;0;-2*deltaY;0];
points（:3）=ones（Ntarget1）;
figure;
h=stem（points（:1）points（:2）‘fill‘）;
grid?on;
set（h‘Linestyle‘‘none‘）;
axis（[-50?50?-50?50]）
title（‘點目標真實圖‘）;xlabel（‘x‘）;ylabel（‘y‘）;

%雷達回波
s=zeros（nrnnan）;
for?n=1:nan
Rs=R0+V*tnan（n）+1/2*ar*tnan（n）^2;??
for?m=1:Ntarget
???rp=sqrt（points（m1）^2+points（m2）^2）;??
???phi=atan2（points（m2）points（m1））+wt*tnan（n）;
???R=sqrt（Rs^2+rp^2+2*Rs*rp*sin（phi））;
???td=2*R/C;
???win_r=（abs（tnrn-td）<=Tr/2）;?????%距離包絡（矩形窗函數(shù)）
???s_mn=points（m3）.*exp（ima*pi*Kr*（tnrn-td）.^2）.*win_r.*exp（-ima*4*pi*R/lamda）;
???s（:n）=s（:n）+s_mn;
end
s（: