ampgra.m

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

function ampgra

global Aamp;
global nampgraf;
global nnampgraf;
global ampaxes;
global ampdeltaf;
global archivo;
global letras;
global numeroamp;
global ampFmax;
global ampFmin;
global ampAmax;
global ampAmin;
global amptipoE;
global amptipoG;

color=['y' 'g' 'c' 'r' 'b' 'm' 'w' 'y' 'g' 'c' 'r' 'b' 'm' 'w' 'y' 'g' 'c' 'r' 'b' 'm' 'w'];

axes(ampaxes);hold off;
for k=nnampgraf;
    j=0;
    for i=nampgraf;
        j=j+1;
        if i==k;
           fult=ampdeltaf(j)*(numeroamp(j+1)-numeroamp(j)-1);
           X = [0:ampdeltaf(j):fult];
           if amptipoG==1;
              loglog(X,Aamp(numeroamp(j):(numeroamp(j+1)-1)),color(i));
           else;
              plot(X,Aamp(numeroamp(j):(numeroamp(j+1)-1)),color(i));
           end;
           hold on;
        end;
    end;
end;

if ampFmin~=0 | ampFmax~=0 | ampAmin~=0 | ampAmax~=0 ;
   axis([ampFmin ampFmax ampAmin ampAmax]);
end;

j=0;
for i=nnampgraf;
    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;

if amptipoE==1;
   xlabel('frecuencia (Hz)');
else;
   xlabel('periodo (seg)');
end;
ylabel('Amplitud');

hold off;