%%========================================================
%%RTPC方法：SAR原理的直觀解釋
%%========================================================
clear;clc;close?all;
%%========================================================
%%參量--常數(shù)
C=3e8;???????????????????????????%傳播速度（光速）
%%參量--雷達特性
Fc=1e9;??????????????????????????%載頻
lambda=C/Fc;?????????????????????%波長
%%參數(shù)--目標區(qū)域
Xmin=0;??????????????????????????%方位向坐標范圍[XminXmax]
Xmax=50;??????????????????
Yc=10000;????????????????????????%測繪帶中心
Y0=500;??????????????????????????%距離向坐標范圍[Yc-Y0Yc+Y0]
?????????????????????????????????%測繪帶寬2*Y0
%%參數(shù)--平臺信息
V=100;???????????????????????????%SAR?飛行速度?100?m/s
H=5000;??????????????????????????%高度?5000?m
R0=sqrt（Yc^2+H^2）;???????????????%SAR距測繪帶中心距離
%%參數(shù)--天線
D=4;?????????????????????????????%天線孔徑最大尺寸
Lsar=lambda*R0/D;????????????????%合成孔徑長度：Lsar決定于天線運動過程中所能接收到的來自同一目標單元的回波信號的最大作用范圍，
?????????????????????????????????%?????????????它等于真實天線所能覆蓋的最大范圍?????????????????????????????????????????????
Tsar=Lsar/V;?????????????????????%合成孔徑時間
%%參數(shù)--慢時間域
Ka=-2*V^2/lambda/R0;????????????%多普勒調(diào)頻斜率
Ba=abs（Ka*Tsar）;????????????????%多普勒調(diào)頻帶寬
PRF=Ba;????????????????????????%PRF：（PRF的取值范圍：1.1BdPRT=1/PRF;?????????????????????%PRT
ds=PRT;?????????????????????????%慢時間域抽樣間隔
Nslow=ceil（（Xmax-Xmin+Lsar）/V/ds）;?%慢時間域抽樣數(shù)
Nslow=2^nextpow2（Nslow）;??????????%取2的階乘數(shù)便于FFT計算
sn=linspace（（Xmin-Lsar/2）/V（Xmax+Lsar/2）/VNslow）;%慢時間域離散時間陣（SAR對所有點有作用的時間）
PRT=（Xmax-Xmin+Lsar）/V/Nslow;????%更新
PRF=1/PRT;
ds=PRT;
%%參數(shù)--快時間域
Tr=5e-6;?????????????????????????%脈沖持續(xù)?10us
Br=30e6;????????????????????????%LFM信號帶寬?30MHz
Kr=Br/Tr;????????????????????????%LFM斜率
Fsr=3*Br;????????????????????????%快時間抽樣頻率：理論上1.5倍以上的抽樣頻率即可，但是在matlab中，3倍Br時波形仍成鋸齒狀
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）;%快時間域抽樣個數(shù)
Nfast=2^nextpow2（Nfast）;???????????????????%取2的階乘便于FFT計算
tm=linspace（2*Rmin/C2*Rmax/C+TrNfast）;?%快時間域離散時間陣
dt=（2*Rmax/C+Tr-2*Rmin/C）/Nfast;????%更新
Fsr=1/dt;
%%參數(shù)--分辨率
DY=C/2/Br;??????????????????????????%距離分辨率
DX=D/2;????????????????????????????%方向分辨率
%%參數(shù)--點目標個數(shù)
Ntarget=1;???????????????????????????%點目標個數(shù)
%目標參數(shù)格式?[x?y?散射系數(shù)]
Ptarget=[XminYc1????????????????????%目標的坐標
?????????????XminYc+DY1
??????????????XminYc+2*DY1
???????????????XminYc+3*DY1
??????????????XminYc+4*DY1
???????????????XminYc+5*DY1
????????????????XminYc+6*DY1
?????????????????XminYc+7*DY1
??????????????????XminYc+8*DY1
???????????????????XminYc+9*DY1
???????????????????XminYc+10*DY1
????????????????????XminYc+11*DY1
?????????????????????XminYc+12*DY1
??????????????????????XminYc+13*DY1
???????????????????????XminYc+14*DY1
????????????????????????XminYc+15*DY1
???????????Xmin+DXYc1????????????????????%目標的坐標
?????????????Xmin+DXYc+DY1
??????????????Xmin+DXYc+2*DY1
???????????????Xmin+DXYc+3*DY1
??????????????Xmin+DXYc+4*DY1
???

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????9738??2010-01-14?14:08??sar_p.m

-----------?---------??----------?-----??----

?????????????????9738????????????????????1