function?speed=SNR（speed）
%信噪比法
%?clear;clc
%?[datpathname]=?uigetfile?（‘*.out‘?‘Load?Ramac?format?radar?file‘）;
%?[HeaderFields]=gprmax（[pathnamedat]）;

%?dt=1.0e9*Header.dt;
%?nlen=Header.iterations;
%?nt=Header.nrx;
%?dx=1;
%?
%?DAT13=Fields.ez‘;
%?DAT=DAT13‘+mean（mean（DAT13‘））;
%?
%?am1=0.25*max（max（abs（DAT13）））;
%?for?jj=1:nt
%?DAT13（:jj）=（1+jj*0.3）*DAT13（:jj）;
%?end
A=speed.A+mean（mean（speed.A））;


m=length（A（1:））;
n=length（A（:1））;
R（1:m）=0;
for?i=1:m
????d=0;p=0;
????r=0;
????for?j=11:n-11
????????for?t=j-5:j
????????????temp1（t）=（A（ti）.^2）*（exp（-0.5*（t-j））.^2）;
????????????d=d+temp1（t）;
????????end
????????for?t=j:j+5
????????????temp2（t）=（A（ti）.^2）*（exp（-0.5*（j-t））.^2）;
????????????p=p+temp2（t）;
????????end
????????r（j）=（p-d）./d;
????end
???
????y=max（r）;
????k=find（r==y）;
????if（size（k）==1）
????????R（i）=k;
????else
????????R（i）=k（1）;
????end
end
figure

imagesc（A）;
%?for?ii=1:nt
%?????sig（1:nlen）=A13（ii1:nlen）;
%?????am（ii）=2*max（sig）;
%?????h=plot（（ii）*dx+sig（1:nlen）/（am1*3）[1:nlen]*dt‘k-‘）;hold?on;
%?????set（h‘Linewidth‘1.5）
%?end
%?hold?off
axis?ij
axis?tight
hold?on;
%?plot（R*dt）;
plot（R）;
%?colorbar;
title（‘信噪比最大法‘）
xlabel（‘Trace?number‘）;
ylabel（‘time?[ns]‘）;
speed.Cmin=0;
speed.Cmax=0;
speed.amp=0;
p=length（R（:1））;
q=length（R（1:））;
if?p>q
%?????speed.time=C.*speed.dt+163.4;
????speed.time=R.*speed.dt;
else
%?????speed.time=C‘.*speed.dt+163.4;
????speed.time=R‘.*speed.dt;
end

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

?????文件???????1518??2010-04-08?09:47??SNR.m

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

?????????????????1518????????????????????1