linbase.m

来自「德国人开发的地震处理分析软件」· M 代码 · 共 61 行

function linbase

global Limiteizq;
global Limiteder;
global A;
global deltat;
global archivo;
global deltat;
global numero;
global nlbgraf;
global MaxA;
global ngraf;
global grafArange;
global fc;

%if Limiteizq==Limiteder;

   izq=numero(nlbgraf);
   der=numero(nlbgraf+1)-1;
   Nn=der-izq+1;
%elseif Limiteizq < Limiteder;
%   posizq = fix(Limiteizq/deltat(nlbgraf));
%   posder = fix(Limiteder/deltat(nlbgraf));
%
%   izq=numero(nlbgraf)+posizq;
%   der=numero(nlbgraf)+posder;
%
%   if der > numero(nlbgraf+1)-1; der=numero(nlbgraf+1)-1; end
%else
%   posizq = fix(Limiteder/deltat(nlbgraf));
%   posder = fix(Limiteizq/deltat(nlbgraf));
%
%   izq=numero(nlbgraf)+posizq;
%   der=numero(nlbgraf)+posder;
%
%   if der > numero(nlbgraf+1)-1; der=numero(nlbgraf+1)-1; end
%
%end

Nfft=2^(fix(log(der-izq+1)/log(2))+1);
df=1/deltat(nlbgraf)/Nfft;
Nceros=fix(fc/df)+1;
filtro=[zeros(Nceros+1,1);ones(Nfft-2*Nceros-1,1);zeros(Nceros,1)];

x = A(izq:der);
y = detrend(x);
yfft = fft(y,Nfft);

yfft=yfft.*filtro;
z=ifft(yfft);
A(izq:der)=real(z(1:Nn));

Amax = max(max(A(numero(nlbgraf):numero(nlbgraf+1)-1)),abs(min(A(numero(nlbgraf):numero(nlbgraf+1)-1))));
MaxA(nlbgraf) = Amax;
Amax=max(MaxA(ngraf));
grafArange=[-Amax Amax];
ngraf=nlbgraf;

grafica;