fftgra.m

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

function fftgra

global Afft;
global nfftgraf;
global nnfftgraf;
global fftaxes;
global fftdeltaf;
global archivo;
global letras;
global numerofft;
global fftFmax;
global fftFmin;
global fftAmax;
global fftAmin;
global ffttipo;

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

axes(fftaxes);hold off;
for k=nnfftgraf;
    j=0;
    for i=nfftgraf;
        j=j+1;
        if i==k;
           fult=fftdeltaf(j)*(numerofft(j+1)-numerofft(j)-1);
           X = [0:fftdeltaf(j):fult];
           if ffttipo~=1;
              loglog(X,Afft(numerofft(j):(numerofft(j+1)-1)),color(i));
           else;
              plot(X,Afft(numerofft(j):(numerofft(j+1)-1)),color(i));
           end;
           hold on;
        end;
    end;
end;

if fftFmin~=0 | fftFmax~=0 | fftAmin~=0 | fftAmax~=0 ;
   axis([fftFmin fftFmax fftAmin fftAmax]);
end;

j=0;
for i=nnfftgraf;
    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 (cm/seg)');
hold off;