fxygra.m

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

function fxygra

global Atf;
global ntfgraf;
global nntfgraf;
global psdaxes;
global tfdeltaf;
global archivo;
global letras;
global numerotf;
global psdFmax;
global psdFmin;
       
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];
           semilogx(X,angle(Atf(numerotf(j):(numerotf(j+1)-1),[2])).*180/pi,color(i));
           hold on;
        end;
    end;
end;  

if psdFmin~=0 | psdFmax~=0 ;
   axis([[psdFmin psdFmax] -200 200]);
end;

j=0;
for i=ntfgraf(1,:);
    j=j+1;
    k=ntfgraf(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('Fase (