clear?all;clc;close?all;
fc=3e9;?????????????????%載波頻率
PRF=2000;???????
Br=5e6;?????????????????%帶寬
fs=10*Br;???????????????%采樣頻率
Tp=5e-6;????????????????%脈寬
Kr=Br/Tp;???????????????%頻率變化率
c=3e8;??????????????????%光速
lamda=c/fc;?????????????%波長
Tr=1/PRF;???????????????%脈沖重復(fù)周期
N_mc=1.5/60*PRF;????????%脈沖個數(shù)
t=0:1/fs:15*Tp+Tp;??????%采樣時間
N_r=length（t）;??????????%采樣點(diǎn)數(shù)
N_target=5;?????????????%目標(biāo)個數(shù)
Rmax=c/2*15*Tp;?????????????????????????????%目標(biāo)最大距離
R_t=Rmax*abs（rand（1N_target））;?????????????%目標(biāo)的距離
RCS_t=10*（exp（i*2*pi*rand（1N_target）））;????%目標(biāo)RCS，幅度為10，相位在（02pi）之間隨機(jī)分布
Vmax=lamda*PRF/2;???????????????????????????%目標(biāo)最大速度
v=Vmax*（（1+rand（1N_target））/2）;????????????%目標(biāo)速度
%%?生成目標(biāo)矩陣
sr=zeros（N_mcN_r）;
for?i=1:N_mc
????ta=（i-1）*Tr;
????sri=0;
????for?k=1:N_target
????????tao=2*（R_t（k）-v（k）.*（ta+t））/c;
????????srj=RCS_t（k）.*rectpuls（t-tao-Tp/2Tp）.*exp（-1j*2*pi*fc*tao+1j*pi*Kr.*（t-tao-Tp/2）.^2）;
????????sri=sri+srj;
????end
????sr（i:）=sri;
end
%%?距離壓縮前的回波
tm=（1:N_mc）/PRF;
R=c*t/2;
figure（1）;
image（Rtm255*abs（sr）/max（max（abs（sr））））??
figure（2）;
plot（t*c/2abs（sr（1:）））?????????
%%?距離壓縮
st=rectpuls（t-Tp/2Tp）.*exp（1i*pi*Kr*（t-Tp/2）.^2）;
stf=conj（fft（st））;
for?i=1:N_mc
????sr（i:）=ifft（fft（sr（i:））.*stf）;????????????
end
figure（3）;
image（Rtm255*abs（sr）/max（max（abs（sr））））????????????????
figure;
plot（t*c/2abs（sr（1:）））???????????????????????
?????????????
sr=fft（sr[]1）;
V=linspace（0PRF50）*lamda/2;
figure;image（RV255*abs（sr）/max（max（abs（sr））））??

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

?????文件???????1624??2017-05-22?11:46??echo.m

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

?????????????????1624????????????????????1