%%========================================================
clear;clc;close?all;
%%========================================================
%%參數(shù)——常量
C=3e8;???????????????????????????%電磁波速度

%%參數(shù)——雷達(dá)規(guī)格
Fc=3e9;??????????????????????????%載波頻率
lambda=C/Fc;?????????????????????%波長?

%%參數(shù)——目標(biāo)區(qū)域
X0=500;??????????????????????????%方位向范圍：[-X0X0]
Y=10000;?????????????????????????%測繪中心Y坐標(biāo)
Y0=500;??????????????????????????%成像寬度：2*Y0

%%參數(shù)
theta=0;?????????????????????????%平臺波束斜視角
theta=theta*pi/180;??????????????%rad
V=100;???????????????????????????%平臺速度：100?m/s
H=4000;??????????????????????????%平臺高度：4?km
Rb=sqrt（H^2+Y^2）;
Rc=Rb/cos（theta）;????????????????%目標(biāo)區(qū)域距離向范圍：[Rc-Y0Rc+Y0]

%%參數(shù)——雷達(dá)天線
D=4;?????????????????????????????%天線方位向尺寸
Lsar=lambda*Rc/D;????????????????%SAR合成孔徑長度
Tsar=Lsar/V;?????????????????????%SAR合成孔徑時間

%%參數(shù)——慢時間域
Ka=-2*（V*cos（theta））^2/lambda/Rc;????????????%doppler調(diào)頻斜率
Ba=abs（Ka*Tsar）;????????????????%doppler調(diào)頻頻寬
PRF=Ba;?????????????????????????%脈沖重復(fù)頻率
PRT=1/PRF;??????????????????????%脈沖重復(fù)時間
ds=PRT;?????????????????????????%慢時域采樣間距
????????????????????????????????%為了進(jìn)行FFT，對PRT、PRF、ds進(jìn)行運(yùn)算并更新
Nslow=ceil（（2*X0+Lsar）/V/ds）;???%慢時域采樣點(diǎn)數(shù)
Nslow=2^nextpow2（Nslow）;????????%使采樣點(diǎn)數(shù)為2的冪；N點(diǎn)

%Nslow=1024;

sn=linspace（-（X0+Lsar/2）/V（X0+Lsar/2）/VNslow）;%慢時域離散時間序列，1*N矩陣
PRT=（2*X0+Lsar）/V/Nslow;????????%更新PRT、PRF、ds
PRF=1/PRT;
ds=PRT;
%%參數(shù)——快時域
Tr=5e-6;????????????????????????%Chirp脈沖持續(xù)時間：5?us
Br=50e6;????????????????????????%chirp信號調(diào)頻帶寬：50?MHz
Kr=Br/Tr;???????????????????????%chirp信號調(diào)頻斜率
Fsr=2*Br;???????????????????????%快時域采樣頻率
dt=1/Fsr;???????????????????????%快時域采樣間距
Rmin=sqrt（（Y-Y0）^2+H^2）;????????%最小斜距
Rmax=sqrt（（Y+Y0）^2+H^2+（Lsar/2）^2）;?????%最大斜距
Nfast=ceil（2*（Rmax-Rmin）/C/dt+Tr/dt）;???%快時域采樣點(diǎn)數(shù)
Nfast=2^nextpow2（Nfast）;????????????????%使采樣點(diǎn)數(shù)為2的冪；M點(diǎn)
tm=linspace（2*Rmin/C2*Rmax/C+TrNfast）;%快時域離散時間序列，1*M矩陣
dt=（2*Rmax/C+Tr-2*Rmin/C）/Nfast;????????%更新Fsr、dt
Fsr=1/dt;

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

%%參數(shù)——點(diǎn)目標(biāo)
Ntarget=3;??????????????????????????%點(diǎn)目標(biāo)個數(shù)

%格式[方位向?距離向水平投影?目標(biāo)散射系數(shù)]
%標(biāo)準(zhǔn)點(diǎn)位置
Ptarget=[0Y1
?????????62*DXY1
?????????75Y+25*DY1?];
%分辨率點(diǎn)位置
%Ptarget=[0Y1
%?????????3*DXY1
%?????????0Y+2*DY1?];

disp（‘參數(shù):‘）
disp（‘快時域采樣數(shù)：‘）;disp（Nfast）
disp（‘慢時域采樣數(shù)：‘）;disp（Nslow）
disp（‘距離向分辨率：‘）;disp（DY）
disp（‘方位向分辨率：‘）;disp（DX）?????
disp（‘合成孔徑長度：‘）;disp（Lsar）?????
%%========================================================
%%產(chǎn)生單、多點(diǎn)回?fù)芑夭ㄐ盘?，一個N*M矩陣
K=Ntarget;???????????????????????
N=Nslow;?????????????????????????
M=Nfast;????????????????????????
T=Ptarget;??????????????????????
Srnm=zeros（NM）;
for?k=1:1:K??????????????????????
????sigma=T（k3）;?????????????????%獲取點(diǎn)目標(biāo)散射系數(shù)
????Dslow=sn*V-T（k1）;????????????%慢時域中雷達(dá)與目標(biāo)方位向距離，1*N矩陣
????R=sqrt（Dslow.^2+T（k2）^2+H^2）;%斜距，1*N矩陣
????tau=2*R/C;????????????????????%時間延遲，1*N矩陣
????Dfast=ones（N1）*tm-tau‘*ones（1M）;??????????????????%[t-2*R/C]項(xiàng)，N*M

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

?????文件???????5366??2009-05-15?18:57??pangbanbi1.m

?????文件???????6315??2009-05-21?13:50??chengxiang.m

?????文件???????5324??2009-05-21?14:37??dingweixianshi.m

?????文件?????????68??2009-04-10?13:32??ftx.m

?????文件?????????72??2009-04-10?13:32??fty.m

?????文件?????????71??2009-04-10?13:32??iftx.m

?????文件?????????75??2009-04-10?13:32??ifty.m

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

????????????????17291????????????????????7