%%?&&&&&&&&&&
%?date:2018.3.6
%?點目標RDA算法
%%?&&&&&&&&&&
clear;clc;close?all;
C=3e8;????????????%光速
Fc=1e9;???????????%載波頻率1GHz
lambda=C/Fc;??????%波長

%%目標區域參數
Xmin=0;?????%目標區域方位向范圍
Xmax=50;????%[XminXmax]
Yc=10000;???%成像區域中線
Y0=500;?????%目標區域距離向范圍[Yc-Y0Yc+Y0]成像寬度為2*Y0
%%軌道參數
V=100;??????????????%雷達平臺速度100m/s
H=0;????????????????%雷達平臺高度
R0=sqrt（Yc^2+H^2）;??%最短距離
%%天線參數
D=4;????????????????%方位向天線孔徑長度
Lsar=lambda*R0/D;???%SAR合成孔徑長度
Tsar=Lsar/V;????????%合成孔徑時間

%%?慢時間域參數（方位向）
Ka=-2*V^2/lambda/R0;%線性調頻率
Ba=abs（Ka*Tsar）;%調制帶寬
PRF=2*Ba;???????%脈沖重復頻率
PRT=1/PRF;??????%脈沖重復時間
ds=PRT;?????????%脈沖重復周期?
Nslow=ceil（（Xmax-Xmin+Lsar）/V/ds）;%脈沖數，ceil為取整函數
Nslow=2^nextpow2（Nslow）;?%量化為2的指數
sn=linspace（（Xmin-Lsar/2）/V（Xmax+Lsar/2）/VNslow）;???%慢時域的時間矩陣
PRT=（Xmax-Xmin+Lsar）/V/Nslow;????%更新脈沖重復時間???
PRF=1/PRT;?%?更新脈沖重復頻率
fa=linspace（-0.5*PRF0.5*PRFNslow）;%慢時間域頻率（方位頻率）

%%?快時間域參數（距離向）
Tr=5e-6;????%脈沖寬度（脈沖持續時間5us）
Br=30e6;????%chirp頻率調制帶寬為30MHz
Kr=Br/Tr;???%chirp調頻率
Fsr=2*Br;?????????????????????????%快時間域采樣頻率，為3倍帶寬
dt=1/Fsr;?????????????????????????%快時間域采樣間隔
Rmin=sqrt（（Yc-Y0）^2+H^2）;???????????????%近距離點
Rmax=sqrt（（Yc+Y0）^2+H^2+（Lsar/2）^2）;????%遠距離點
Nfast=ceil（2*（Rmax-Rmin）/C/dt+Tr/dt）;???%快時域采樣數量
Nfast=2^nextpow2（Nfast）;????????????????%更新為2的冪次
tm=linspace（2*Rmin/C2*Rmax/C+TrNfast）;%快時域的離散時間矩陣
dt=（2*Rmax/C+Tr-2*Rmin/C）/Nfast;???%更新間隔
Fsr=1/dt;
fr=linspace（-0.5*Fsr0.5*FsrNfast）;%快時間域頻率（距離頻率）

%%?分辨率參數設置
DY=C/2/Br;????%距離向分辨率
DX=D/2;???????%方位向分辨率

%%?目標設定
Ntarget=3;???%目標數目
Ptarget=[XminYc1???????????????%目標位置
?????????XminYc+10*DY1
?????????Xmin+20*DXYc+50*DY1];?
K=Ntarget;????%目標數目
N=Nslow;????%慢時間采樣數
M=Nfast;?????%快時間采樣數
T=Ptarget;???%目標位置

%%?合成孔徑回波仿真
Srnm=zeros（NM）;????????????????????????%生成零矩陣存儲回波信號
for?k=1:1:K?????????????????????????????%總共K個目標
????sigma=T（k3）;???????????????????????%得到目標的反射系數
????Dslow=sn*V-T（k1）;??????????????????%方位向距離，投影到方位向的距離
????R=sqrt（Dslow.^2+T（k2）^2+H^2）;??????%實際距離矩陣
????tau=2*R/C;??????????????????????????%信號延時
????Dfast=ones（N1）*tm-tau‘*ones（1M）;??%時間矩陣
????phase=pi*Kr*Dfast.^2-（4*pi/lambda）*（R‘*ones（1M））;??%相位
????Srnm=Srnm+sigma*exp（1i*phase）.*（0end
%%?圖形一
figure（1）?
subplot（211）?
imagesc（abs（Srnm））;title（‘SAR?data‘）?
subplot（212）?
imagesc（angle（Srnm））?

row=tm*C/2;??%%?行
col=sn*V;????%%?列
%%?距離壓縮
%?tr=tm-2*R/C;

%?ta=sn-（Xmin+Xmax）/2/V;

%?ht_rc?=?exp（1i*pi*Kr*tr.^2）.*（abs（tr）
%ht_rc?=?kaiser（N?2.5）‘.*exp（1i*pi*Kr*tr.^2）.*（abs（tr）%ht_rc?=?exp（1i*pi*Kr*（tm-2*Rmin/C）.^2）.*（0%（abs（tau-2*Rc/c）%?%?beta?=?2.5
%?
%?Hf_rc?=?conj（fftshift（fft（fftshift（ht_rc））））;???%?轉換到頻域，并做復共軛
%?S_rf?=?fftshift（fft（fftshift（Srnm2）[]2）2）;?????%?頻域
%?S_rcf?=?S_rf.*（ones（