SAR成像与CS算法实现MATLAB代码解析

SAR-CSA算法成像（matlab）
(2010-12-28 11:05:56) 
 
clear all;   close all;
t1=cputime;                                                  
Fs=32.317e+006;                            
Fr=7.2135e+11   ;                     
start=6.5956e-003;                     
Tr=4.175e-05 ;                               
R0=9.88646462e+05   ;                        
f0=5.3e+09  ;                          
lamda=0.05667;                    
Fa=1.25698e+03 ;                    
Vr=7062;                             
Kr=0.72135e+012;
Ka=1733;                        
Fc=-6900;                             
c=299790000;                    
load CDdata1.mat
data=double(data);                        
[length_a,length_r]=size(data);                
T_start=6.5956e-003;                   %数据窗开始时间
tau=T_start:double(1/Fs):T_start+double(length_r/Fs)-double(1/Fs);  
R_ref=(T_start+length_r/Fs)/2*c;            
f_a=(-Fa/2+Fc):(Fa/length_a):(Fa/2+Fc-Fa/length_a);      
f_r=0:Fs/length_r:Fs-Fs/length_r;            
D = (1 - (f_a*lamda/2/Vr).^2).^0.5;         %距离徙动因子
alpha = 1./D - 1;                    
R = R_ref./D;                           %距离多普勒域中更精确的双曲线距离等式
Z=(R0*c*f_a.^2)./(2*Vr^2*f0^3.*D.^3);