?%%?正側視點目標成像--CSA
?clc;
?clear;
?close?all;
?%%?參數設置
?c=3e8;?????????????????????????
?fc=3e9;????????????????????????
?lambda=c/fc;??????????????????%載波頻率（Hz）、波長0.1（m）
?Ralpha=1.2;??Aalpha=1.3;??????%距離向、方位向過采樣率
?%----------觀察場景相關參數----------%
?H=8000;???????????????????????%平臺運行高度
?Yc=6000;??????????????????????%場景中心線長度
?R0=sqrt（H^2+Yc^2）;
?Xmin=-200;????????????????????%[Xmin?Xmax]
?Xmax=200;?????????????????????%在合成孔徑長度Ls基礎上，額外方位向范圍
?Yw=300;??????????????????????
??%----------天線孔徑設置----------%
?La=1.5;
?theta=0.886*lambda/La;
?Ls=R0*theta;
?%----------慢時間參數設置----------%
?v=150;????????????????????????%平臺運行速度
?Ts=Ls/v;
?Xwid=Ls+Xmax-Xmin;????????????%慢時間域時間窗長度
?Twid=Xwid/v;
?Ka=2*v^2/lambda/R0;???????????%方位向調頻率
?Ba=abs（Ka*Ts）;????????????????%方位向帶寬
?PRF=Aalpha*Ba;
?PRT=1/PRF;
?dx=PRT;???????????????????????%方位向采樣間隔
?N=ceil（Twid/dx）;
?N=2^nextpow2（N）;??????????????%提高fft效率
?x=linspace（（Xmin-Ls/2）/v（Xmax+Ls/2）/vN）;
???????????????????????????????%慢時間域時間序列
?X=v*x;????????????????????????%慢時間域時間序列對應方位向距離
?PRT=Twid/N;???????????????????%更新
?PRF=1/PRT;
?dx=PRT;
?%----------快時間參數設置----------%
?Tr=5e-6;??????????????????????%脈沖寬度
?Br=100e6;?????????????????????%帶寬100MHz
?Kr=Br/Tr;?????????????????????%調頻率
?Fs=Ralpha*Br;?????????????????%距離向采樣率
?dt=1/Fs;
?Rmin=sqrt（（Yc-Yw）^2+H^2）;
?Rmax=sqrt（（Yc+Yw）^2+H^2）;
?Rm=Rmax-Rmin+c*Tr/2;??????????%斜距測繪帶寬
?M=ceil（2*Rm/c/dt）;????????????%采樣點數?
?M=2^nextpow2（M）;???????
?t=linspace（2*Rmin/c-Tr/22*Rmax/c+Tr/2M）;
???????????????????????????????%快時間域時間序列
?r=c*t/2;??????????????????????%快時間域時間序列對應斜距
?dt=（2*Rmax/c+Tr-2*Rmin/c）/M;??%更新
?Fs=1/dt;
?%----------目標參數設置----------%
?Ptarget=[0R01;50R0+1000.8;700R0+1000.8;];
?Ntarget=length（Ptarget）;
?%%?生成SAR回波
?s0=zeros（NM）;
?h=waitbar（0‘SAR回波生成‘）;
?for?k=1:Ntarget
????R=sqrt（Ptarget（k2）^2+（X-Ptarget（k1））.^2）;?
????delay=2*R/c;???????????????%距離徙動帶來的距離時間延時
????Delay=ones（N1）*t-delay‘*ones（1M）;???%生成N*M的矩陣
????Phase=1j*pi*Kr*Delay.^2-1j*4*pi*fc*（R‘*ones（1M））/c;
????s0=s0+Ptarget（k3）*rectpuls（Delay/Tr）.*...
??????rectpuls（（X-Ptarget（k1））‘*ones（1M）/Ls）.*exp（Phase）;
???%s0=s0+Ptarget（k3）*exp（Phase）.*（abs（Delay）<=Tr/2）.*（abs（（X-Ptarget（k1））‘*ones（1M））<=Lsar/2）;
????waitbar（k/Ntarget）;
?end
?close（h）;
?%----------繪制原始信號----------%
?h=waitbar（0‘繪制原始信號‘）;
?figure（1）;
?imagesc（r*1e-3Xabs（s0））;
?axis?tight;
?xlabel（‘距離向（km）\rightarrow‘）;
?ylabel（‘\leftarrow方位向（m）‘）;
?title（‘原始信號‘）;
?close（h）;
?clear?R?delay?Delay?Phase;
?%%?方位向fft、補余RCMC（在距離多普勒域中進行補余RCMC）
?h=waitbar（0‘方位向fft、補余RCMC‘）;
?fr=（-M/2:M/2-1）*Fs/M;?????????%距離向頻率
?fa=（-N/2:N/2-1）*PRF/N;????????%方位向頻率
?Srd=ftx（s0）;??????????????????%距離多普勒域
?Km=Kr;????????????????????????%這里我們近似等于
?D_fa_v=sqrt（1-（c^2*fa.^2/4/v^2/fc^2））.‘*ones（1M）;
?Ssc=exp（1j*pi*Km*（1./D_fa_v-1）.*（（ones（N1）*t-2*R0/c./D_fa_v）.^2））;
?S1=Srd.*Ssc;
?figure（2）;
?imagesc（r*1e-3faabs（iftx（S1）））;
?axis?tight;
?xlabel（‘距離向（km）\rightarrow‘）;
?ylabel

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

?????文件???????5801??2017-12-28?20:45??CS\CSA.m

?????文件???????5871??2014-03-28?16:30??CS\CSA_basic.m

?????文件????????966??2014-03-25?14:22??CS\CSA_tu7_3.m

?????文件????????956??2014-04-06?15:16??CS\CSA_tu7_5.m

?????文件???????2812??2014-11-18?20:19??CS\CS_scsj.m

?????文件???????2808??2014-11-18?20:24??CS\cs_scsj2.m

?????文件???????4643??2014-11-14?09:25??CS\CS_theta.m

?????文件???????4718??2014-11-16?10:44??CS\CS_theta1.m

?????文件????????115??2014-03-23?14:42??CS\ftx.m

?????文件????????116??2014-03-23?14:42??CS\fty.m

?????文件????????119??2014-03-23?14:42??CS\iftx.m

?????文件????????120??2014-03-23?14:43??CS\ifty.m

?????文件????????481??2014-03-23?20:14??CS\islrfunc.m

?????文件????????270??2006-06-30?12:06??CS\plot_img.m

?????文件????????277??2014-03-23?20:14??CS\pslrfunc.m

?????文件????????115??2014-03-23?14:42??RD\ftx.m

?????文件????????116??2014-03-23?14:42??RD\fty.m

?????文件????????205??2014-11-18?14:57??RD\Hist.m

?????文件????????713??2014-11-18?14:58??RD\Histeq_f.m

?????文件????????119??2014-03-23?14:42??RD\iftx.m

?????文件????????120??2014-03-23?14:43??RD\ifty.m

?????文件????????481??2014-03-23?20:14??RD\islrfunc.m

?????文件????????666??2014-11-18?14:59??RD\medfilt_f.m

?????文件???????1045??2014-11-18?14:59??RD\Pinghua_f.m

?????文件????????335??2014-11-12?20:12??RD\plot_img.m

?????文件???????1651??2014-11-20?10:52??RD\poumian.m

?????文件????????277??2014-03-23?20:14??RD\pslrfunc.m

?????文件???????8287??2017-05-11?10:06??RD\RDA.m

?????文件???????7515??2014-10-28?10:35??RD\RDA_basic.m

?????文件???????2584??2014-11-20?09:14??RD\RDA_scsj.m

............此處省略31個文件信息