clear?all
%%?頻域脈壓&方位匹配
C=3.0e8;????????????????????%光速
c=C;
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）+j*s（2:2:2097152）;
s=reshape（s10241024）;
s=s-mean（mean（s））;
[ran_numazi_num]=size（s）;
%%?成像
x=s;
[nrnnan]=size（x）;
tnrn=[-nrn/2:nrn/2-1].‘/Fs;
fr=[-nrn/2:nrn/2-1].‘/nrn*Fs;%距離頻率
H_match=exp（-1j*pi/gama*fr.^2）;
%?H=exp（-1j*pi/gama*fr.^2）;
H1=exp（1j*4*pi*Rs/c*fr）;
x1=fftshift（fft（fftshift（x1）[]1）1）;
x2=zeros（nrnnan）;
%距離脈壓
for?n=1:nan
???x2（:n）=x1（:n）.*H_match.*H1;
end
x3=ifftshift（ifft（ifftshift（x21）[]1）1）;
figure;imagesc（abs（x3））
%方位脈壓
deltaR=c/2/Fs;
R=Rs+[-nrn/2:nrn/2-1]*deltaR;
faM=2*v/lambda;
fa=[-nan/2:nan/2-1]/nan*PRF;
x4=fftshift（fft（fftshift（x32）[]2）2）;
for?nn=1:nrn
%????Ha=exp（1j*2*pi*Rs/v*sqrt（faM.^2-fa.^2））?;
???Ha=exp（1j*2*pi*R（nn）/v*sqrt（faM.^2-fa.^2））?;%場景很小，直接用場景中心匹配都不會散焦
???x5（nn:）=x4（nn:）.*Ha;
end
x6=ifftshift（ifft（ifftshift（x52）[]2）2）;
figure;imagesc（abs（x6）[0?600]）

imwrite（abs（x6）/150‘YD_result1.jpg‘‘jpg‘）
imshow（‘YD_result.jpg‘）

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

?????文件????2097152??1999-05-06?22:03??raws.dat

?????文件???????1975??2017-06-17?22:01??SAR_image_simple.m

-----------?---------??----------?-----??----

??????????????2099127????????????????????2