%%?&&&&&&&&&&
%?date:2018.3.6
%?點(diǎn)目標(biāo)RDA算法
%%?&&&&&&&&&&
clear;clc;close?all;
C=3e8;????????????%光速
Fc=1e9;???????????%載波頻率1GHz
lambda=C/Fc;??????%波長(zhǎng)

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

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

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

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

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

%%?合成孔徑回波仿真
Srnm=zeros（NM）;????????????????????????%生成零矩陣存儲(chǔ)回波信號(hào)
for?k=1:1:K?????????????????????????????%總共K個(gè)目標(biāo)
????sigma=T（k3）;???????????????????????%得到目標(biāo)的反射系數(shù)
????Dslow=sn*V-T（k1）;??????????????????%方位向距離，投影到方位向的距離
????R=sqrt（Dslow.^2+T（k2）^2+H^2）;??????%實(shí)際距離矩陣
????tau=2*R/C;??????????????????????????%信號(hào)延時(shí)
????Dfast=ones（N1）*tm-tau‘*ones（1M）;??%時(shí)間矩陣
????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））））;???%?轉(zhuǎn)換到頻域，并做復(fù)共軛
%?S_rf?=?fftshift（fft（fftshift（Srnm2）[]2）2）;?????%?頻域
%?S_rcf?=?S_rf.*（ones（