filgra.m

德国人开发的地震处理分析软件

function filgra

global Afil;
global nfilgraf;
global nnfilgraf;
global filaxes;
global fildeltaf;
global archivo;
global letras;
global numerofil;
global filFmax;
global filFmin;
global filAmax;
global filAmin;
global tipofilG;
       
color=['y' 'g' 'c' 'r' 'b' 'm' 'w' 'y' 'g' 'c' 'r' 'b' 'm' 'w' 'y' 'g' 'c' 'r' 'b' 'm' 'w'];

axes(filaxes);
for k=nnfilgraf;
    j=0;
    for i=nfilgraf;
        j=j+1;
        if i==k;
           Nf=numerofil(j+1)-numerofil(j);
           filff=fildeltaf(j)*Nf;
           X = [fildeltaf(j):fildeltaf(j):filff];
           if tipofilG==1;
              plot(X,Afil(numerofil(j):numerofil(j+1)-1),color(i));
           else;
              loglog(X,Afil(numerofil(j):numerofil(j+1)-1),color(i));
           end;
           hold on;
        end;
    end;
end;  

if filFmin~=0 | filFmax~=0 | filAmin~=0 | filAmax~=0;
   axis([filFmin filFmax filAmin filAmax]);
end;

j=0;
for i=nnfilgraf;
    xtext = 0.05+j*0.25;
    text('Position',[xtext 1.1],'Units','normalized','Color',color(i),'String',archivo(i,1:letras(i)));
    j=j+1;
end;

xlabel('frecuencia (Hz)');
ylabel('Amplitud');
hold off;