clear?;close?all;
clc;
%???參數(shù)說明

Fs=548.571428571429;?????????????????%MHz	采樣率
Fr=-2e14;???????????????%MHz/us	距離調(diào)頻率
start=6.5956e-003;??????????????%ms	數(shù)據(jù)窗開始時間
Tr=4.175e-05?;??????????????????%us	脈寬
R0=23306.25??;??
f0=5.3e+09??;???????????????????%GHz	雷達(dá)頻率
lamda=0.03125;??????????????????%m	雷達(dá)波長
Fa=533.330793?;????????????????%Hz	脈沖重復(fù)頻率
Vr=154.195864;????????????????????????%m/s	有效雷達(dá)速率
Kr=-2e14;
Ka=0;????????????????????????%Hz/s	方位調(diào)頻率
Fc=0;???????????????????????%Hz	多普勒中心頻率
c=299790000;????????????????????%光速
range_size=16384;
%???讀取數(shù)據(jù)
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）;
%-----定義外部參數(shù)

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