function?[y]?=?stretch（nscattaupf0bscat_rangerrecscat_rcswinid）
eps?=?1e-16;
htau?=?taup/2;
c?=?3e8;
trec?=?2.*rrec/c;
n?=?fix（2.*trec*b）;%fix向0取整
m?=?power_integer_2（n）;
nfft?=?2.^m;
%初始化矩陣
x（nscat1:n）?=?0;
y（1:n）?=?0;
%確定適當的窗口
if?（winid?==?0）
????win（1:n）?=?1;
????win?=?win‘;????%不加窗
else
????if（winid?==?1）
????win?=?hamming（n）;???%加漢明
????else
????????if（winid?==?2）
????????????win?=?kaiser（npi）;??%參數為pi的kaiser窗
????????else
????????????if（winid?==?3）
????????????????win?=?chebwin（n60）;???%旁瓣為-60db的切比雪夫窗
????????????end
????????end
????end
end

deltar?=?c/2/b;
max_rrec?=?deltar*nfft/2;
maxr?=?max（scat_range）;

if（rrec>max_rrec|maxr>=rrec）
????‘Error.Receive?window?is?too?large;or?scatterers?fall?outside?window‘;
????return
end

t?=?linspace（0taupn）;
for?j?=?1:1:nscat
????range?=?scat_range（j）;
????psi1?=?4.*pi*range*f0/c-4.*pi*b*range*range/c/c/taup;
????psi2?=?（2*4.*pi*b*range/c/taup）.*t;
????x（j:）?=?scat_rcs（j）.*exp（i*psi1+i.*psi2）;
????y?=?y?+?x（j:）;
end

figure
plot（treal（y）‘k‘）;
xlabel（‘時間延遲‘）;
ylabel（‘未壓縮的回波信號‘）;

ywin?=?y.*win‘;
yfft?=?fft（yn）./n;
out?=?fftshift（abs（yfft））;
figure
delinc?=?rrec/n;
dist?=?linspace（（-rrec/2）rrec/2n）;
plot（distout‘k‘）;
xlabel（‘相對距離/m‘）;
ylabel（‘壓縮回波信號‘）;

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

?????文件???????1355??2018-03-21?21:20??stretch.m

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

?????????????????1355????????????????????1