txygra.m

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

function txygra

global Atf;
global ntfgraf;
global nntfgraf;
global psdaxes;
global tfdeltaf;
global archivo;
global letras;
global numerotf;
global psdFmax;
global psdFmin;
global psdAmax;
global psdAmin;
global psdtipo;
       
color=['y' 'g' 'c' 'r' 'b' 'm' 'w' 'y' 'g' 'c' 'r' 'b' 'm' 'w' 'y' 'g' 'c' 'r' 'b' 'm' 'w'];

axes(psdaxes);
l=0
for k=nntfgraf(1,:);
    l=l+1;
    j=0;
    for i=ntfgraf(1,:);
        j=j+1;
        if k==i & nntfgraf(2,l)==ntfgraf(2,j); 
           fult=tfdeltaf(j)*(numerotf(j+1)-numerotf(j)-1);
           X = [0:tfdeltaf(j):fult];
           if psdtipo~=1;
              loglog(X,abs(Atf(numerotf(j):(numerotf(j+1)-1),[2])),color(i));
           else;
              plot(X,abs(Atf(numerotf(j):(numerotf(j+1)-1),[2])),color(i));
           end;
           hold on;
        end;
    end;
end;  

if psdFmin~=0 | psdFmax~=0 | psdAmin~=0 | psdAmax~=0;
   axis([psdFmin psdFmax psdAmin psdAmax]);
end;

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

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