clear?all;
clc;
close?all;
C=3e8;??????????????%光速
ima=sqrt（-1）;?????????%虛數單位
%%場景中心參數
R0=2e4;?????????????%場景中心斜距
Wr=0.12e3;??????????%場景寬度
%%雷達參數
fc=10e9;?????????????%工作頻率10GHz
lamda=C/fc;??????????%波長
H=5000;?????????????%高度?5000m
SNR=-10;
%%天線參數
D=4;????????????????%天線孔徑長度
Lsar=lamda*R0/D;?????%雷達合成孔徑長度
%%距離向參數
Tr=10e-6;????????????%發射脈沖時寬
B=180e6;????????????%帶寬
Kr=B/Tr;????????????%調頻率
Fsr=200e6;??????????%距離采樣率
dt=1/Fsr;????????????%距離采樣時間
nrn=round（（Tr+2*Wr/C）*Fsr/2）*2;????????%距離向的取樣數
nrn_Tr=Fsr*Tr;????????????????????????%脈沖區間的取樣數
Nfast_signal=exp（ima*pi*Kr*（[-nrn_Tr/2:nrn_Tr/2-1].‘/Fsr）.^2）;???%參考信號
Nfast=2.^ceil（log2（nrn+nrn_Tr））;??%為了達到線性卷積的要求而取兩者之和
r_signal=[zeros（round（Nfast/2-nrn_Tr/2）1）;Nfast_signal;zeros（round（Nfast/2-nrn_Tr/2）1）];
%%方位向參數
Fa=100;??????????????%方位采樣率
PRF=500;????????????%脈沖重復頻率
PRT=1/PRF;??????????%脈沖重復時間
V=100;
ar=5;
wt=0.05;
nan=round（PRF*lamda*R0/D/V/2）*2;?????%方位向取樣數
%%分辨率參數
DY=C/2/B;???????????%距離向分辨率
DX=D/2;?????????????%方位向分辨率
%%目標參數
Ntarget=11;???????	???%點目標個數
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‘）;
%雷達回波
nrn=2048;
nan=512;
tnan=[-nan/2:nan/2-1]*PRT;
tnrn=2*R0/C+[-nrn/2:nrn/2-1]‘/Fsr;
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