clear?all;

Tr=4e-6;?????????????????????????????????????
Br=50e6;?????????????????????????????????????
Kr=Br/Tr;???????????????????????????????????????
Nr=512;
fs=100e6;
dt=1/fs;??
c=3e8;????????????????????????????????????
fc=75e9;??????????????????????????????????
lamda=c/fc;?
Rc=865000;
?

v=7098.2;????????????????????????????????????????????
Lsar=lamda*Rc/6;??????????????????????????????????????
Ta=Lsar/v;
Na=512;
%?theta=28/180*pi;
fa=2*v/6*2;
du=1/fa;????????????????????????????????????????
?
%%%%%%%%%%%%%%%%%%%%%%%%%載入回波數據1%%%%%%%%%%%%%%%%%%%%%%%%%%%
load?planechoh1.mat;
%?load?zhui1.mat;
%?load?ping1.mat;
?

tr=-Tr/2:dt:Tr/2;
pr_t=exp（i*2*pi*1/2*Kr*tr.^2）;
pr_t=[zeros（156）pr_t];??????????????????
pr_f=fft（pr_tNr）;
wr=fftshift（hamming（Nr））.‘;
for?flaga=1:Na
????echo（flaga:）=fft（fftshift（echo（flaga:）））;
????echo（flaga:）=echo（flaga:）.*conj（pr_f）.*wr;
????echo（flaga:）=ifft（echo（flaga:））;
end
?

Ka=2*v^2/lamda/Rc;
ta=-Ta/2:du:Ta/2;
ta=ta+0.8/6*du;
pa_t=exp（i*2*pi*1/2*Ka*ta.^2+i*0*2*pi*ta）;
pa_t=[zeros（175）pa_t];????????????????
pa_f=fft（pa_tNa）;
wa=fftshift（hamming（Na））;
for?flagr=1:Nr
????echo（:flagr）=fft（fftshift（echo（:flagr）））;
????echo（:flagr）=echo（:flagr）.*（pa_f）.‘.*wa;
????echo（:flagr）=ifft（echo（:flagr））;
end
%?fig1=echo;
?
%?a1=echo（55:105180:280）;
%?figure
%?imagesc（abs（a1））;

I=mat2gray（abs（echo））;
S=medfilt2（I）;
figureimshow（S）

save?slc1.mat?echo?


%%%%%%%%%%%%%%%%%%%%%%%%%?第二組成像??%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%載入回波數據2%%%%%%%%%%%%%%%%%%%%%%%%%%%
load?planechoh2.mat;
%?load?zhui2.mat;
%?load?ping2.mat;
?

tr=-Tr/2:dt:Tr/2;
pr_t=exp（i*2*pi*1/2*Kr*tr.^2）;
pr_t=[zeros（156）pr_t];??????????????????
pr_f=fft（pr_tNr）;
wr=fftshift（hamming（Nr））.‘;
for?flaga=1:Na
????echo（flaga:）=fft（fftshift（echo（flaga:）））;
????echo（flaga:）=echo（flaga:）.*conj（pr_f）.*wr;
????echo（flaga:）=ifft（echo（flaga:））;
end
?

Ka=2*v^2/lamda/Rc;
ta=-Ta/2:du:Ta/2;
ta=ta+0.8/6*du;
pa_t=exp（i*2*pi*1/2*Ka*ta.^2+i*0*2*pi*ta）;
pa_t=[zeros（175）pa_t];???????????????
pa_f=fft（pa_tNa）;
wa=fftshift（hanning（Na））;
for?flagr=1:Nr
????echo（:flagr）=fft（fftshift（echo（:flagr）））;
????echo（:flagr）=echo（:flagr）.*（pa_f）.‘.*wa;
????echo（:flagr）=ifft（echo（:flagr））;
end


I=mat2gray（abs（echo））;
S=medfilt2（I）;
figureimshow（S）
save?slc2.mat?echo

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

?????文件????2106836??2012-06-19?00:33??成像\ping1.mat

?????文件????1908754??2012-06-19?00:34??成像\ping2.mat

?????文件????3286893??2012-06-15?20:34??成像\planechoh1.mat

?????文件????2989756??2012-06-15?20:38??成像\planechoh2.mat

?????文件???????2430??2012-06-30?16:17??成像\slc.m

?????文件????2106923??2012-06-19?01:13??成像\zhui1.mat

?????文件????1923488??2012-06-19?01:14??成像\zhui2.mat

?????目錄??????????0??2013-03-11?11:03??成像

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

?????????????14325080????????????????????8