pxxgra.m

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

function pxxgra

global Apsd;
global npsdgraf;
global nnpsdgraf;
global psdaxes;
global psddeltaf;
global archivo;
global letras;
global numeropsd;
global psdFmax;
global psdFmin;
global psdAmax;
global psdAmin;
global deltat;
global psdtipo;

color1=['y' 'g' 'c' 'r' 'b' 'm' 'w' 'y' 'g' 'c' 'r' 'b' 'm' 'w' 'y' 'g' 'c' 'r' 'b' 'm' 'w'];
color2=['y--' 'g--' 'c--' 'r--' 'b--' 'm--' 'w--' 'y--' 'g--' 'c--' 'r--' 'b--' 'm--' 'w--' 'y--' 'g--' 'c--' 'r--' 'b--' 'm--' 'w--'];

axes(psdaxes);
for k=nnpsdgraf;
    j=0;
    for i=npsdgraf;
        j=j+1;
        if i==k;
           fult=psddeltaf(j)*(numeropsd(j+1)-numeropsd(j)-1);
           X = [0:psddeltaf(j):fult];
           if psdtipo~=1;
              loglog(X,Apsd(numeropsd(j):(numeropsd(j+1)-1),[1]).*deltat(i),color1(i));
           else;
              plot(X,Apsd(numeropsd(j):(numeropsd(j+1)-1),[1]).*deltat(i),color1(i));
           end;
           hold on;
%           loglog(X,(Apsd(numeropsd(j):(numeropsd(j+1)-1),[1])+...
%                     Apsd(numeropsd(j):(numeropsd(j+1)-1),[2])).*deltat(i),color2(3*i-2:3*i));
%           loglog(X,(Apsd(numeropsd(j):(numeropsd(j+1)-1),[1])-...
%                     Apsd(numeropsd(j):(numeropsd(j+1)-1),[2])).*deltat(i),color2(3*i-2:3*i));
        end;
    end;
end; 
 
if psdFmin~=0 | psdFmax~=0 | psdAmin~=0 | psdAmax~=0;
   axis([psdFmin psdFmax psdAmin psdAmax]);
end;

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

xlabel('frecuencia (Hz)');
ylabel('Amplitud (cm2/seg3)');
hold off;