SAR成像的CS算法

clear ;close all;
clc;
%   参数说明

Fs=548.571428571429;                 %MHz	采样率
Fr=-2e14;               %MHz/us	距离调频率
start=6.5956e-003;              %ms	数据窗开始时间
Tr=4.175e-05 ;                  %us	脉宽
R0=23306.25  ;  
f0=5.3e+09  ;                   %GHz	雷达频率
lamda=0.03125;                  %m	雷达波长
Fa=533.330793 ;                %Hz	脉冲重复频率
Vr=154.195864;                        %m/s	有效雷达速率
Kr=-2e14;
Ka=0;                        %Hz/s	方位调频率
Fc=0;                       %Hz	多普勒中心频率
c=299790000;                    %光速
range_size=16384;
%   读取数据
chemin=‘.\image\‘;
[nomchemin]=uigetfile（[chemin‘*.*‘][‘Input the picture‘]100100）;
file=fopen（[cheminnom]‘r‘）;
data1=fread（file[range_size*2 2048]‘float‘）;
data=zeros（range_size2048）;
for i=1:range_size
    data（i:）=data1（2*i-1:）+1i*data1（2*i:）;
end
data=data‘;
% load data.mat
% data=double（data）;
[length_alength_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=（2*T_start+length_r/Fs）/2/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