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