合成孔径雷达成像csChirp Scaling算法Matlab程序含星载观测数据，整副实际成像

合成孔径雷达（SAR）成像算法——CS（Chirp Scaling）算法，真实场景成像，包括读取数据，变量初始化，算法流程，结果输出。并包含仿真报告文件（内含基本理论，算法流程，数据结果和具体代码）。是初学者上手，学习的好资料！！！

C=3.0e8;                    %光速
PRF=1395.0;                 %脉冲重复频率
Fs=1.12492800e+7;           %采样频率
Tp=4.00025602e-5;           %发射脉冲宽度
fC=9.600000000000000e+009;  %发射信号中心频率
gama=-2.3724394e+11;        % Chirp rate
B=-gama*Tp;                  %发射带宽
lambda=C/fC;                %中心频率波长
f0=4.7346377e+6;            %采样中频 

u=3.986008*1e5*1e9;         %宇宙常量
Re=6378.135e3;              %地球高度
h=215e3;                    %卫星高度
Vs=sqrt（u/（h+Re））;          %天速
Vg=Re/（h+Re）*Vs;         %载机速度
v=sqrt（Vs*Vg）;

fdC=1397.0;               %提多普勒中心
kr=-9357.0;               %中心调频率

Rs=-2*v.^2/（kr*lambda）;
%fai=aCos（（（Re+h）.^2-Rs.^2-Re.^2）/（2*Re*Rs））;        %入射角
%beta=aCos（（（Re+h）.^2+Rs.^2-Re.^2）/（2*（Re+h）*Rs））;   %视角

Da=12;                      %天线口径
azibeamwidth=lambda/Da;     %方位波束宽度 
Ta=Rs*azibeamwidth/v;       %合成孔径时间
DeltaF=2*v*sin（azibeamwidth）/lambda;
Band_Ratio=DeltaF/PRF;      %多普勒带宽


tp_num=fix（Tp*Fs）;          %距离上的脉宽
tp_nan=fix（Ta*PRF）;         %方位上的脉宽


fd1=fopen（‘raws.dat‘‘rb‘）;
s=fread（fd1‘uchar‘）;
%fClose（fd1）
s=s（1:2:2097152）+1j*s（2:2:2097152）;
s=reshape（s10241024）;
s=s-mean（mean（s））;
[ran_numazi_num]=size（s）;

%%==CS算法
faM=2*v/lambda;
fa=[-azi_num/2:azi_num/2-1]/azi_num*PRF;
fr=[-ran_num/2:ran_num/2-1]‘/ran_num*Fs;
t=[-ran_num/2:ran_num/2-1]‘/Fs;
sin_t=fa./faM;
Cos_t=（1-（fa./faM）.^2）.^0.5;
a_fa=1./Cos_t-1;
gama_e=1./（1/gama-Rs*2*lambda*sin_t.^2./C^2./Cos_t.^3）;
R_fa=Rs*（1+a_fa）;





for k=1:ran_num
    s1（k:）=fftshift（fft（s（k:）））;
end

for k=1:azi_num
    H1=exp（-1j*pi*（gama_e（k）.*a_fa（k）*（t-2*R_fa（k）/C）.^2））;
    s11（:k）=s1（:k）.*H1;
end

for n=1:azi_num
    s2（:n）=fftshift（fft（s11（:n）））;
end

for n=1:ran_num
    H2=exp（-1j*pi*（1./（gama_e.*（1+a_fa）））.*fr（n）.^2+1j*4*pi*Rs*a_fa.*fr（n）/C）;
    s22（n:）=s2（n:）.*H2;
end

s3=ifft（s22[]1）;

deltaR=（-ran_num/2:ran_num/2-1）‘/Fs*C;
Rb=deltaR+Rs;

for n=1:azi_num
    H3=exp（1j*2*pi/v*Rb*（faM^2-fa（n）.^2）.^0.5）.*exp（1j*4*pi/C^2*gama_e（n）*a_fa（n）*（1+a_fa（n））*（Rb-Rs）.^2）;
    s33（:n）=（s3（:n）.*H3）;
end

s4=ifft（s33[]2）;

imwrite（abs（s4）/200‘piCture1.jpg‘‘jpg‘）;
A=imread（‘piCture1.jpg‘‘jpg‘）;
figure（1）;
imshow（A）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        2314  2018-07-13 17:17  cs0.m
     文件     2097152  1999-05-06 22:03  raws.dat
     文件      675840  2018-09-04 17:03  雷达探测与成像大作业（cs算法）.doc