heingra.m

来自「德国人开发的地震处理分析软件」· M 代码 · 共 50 行

function heingra

global Aine;
global ninegraf;
global nninegraf;
global ineaxes;
global inedeltat;
global deltat;
global archivo;
global letras;
global numeroine;
global ineTmax;
global ineTmin;
global ineAmax;
global ineAmin;
global NNine;
       
color=['y' 'g' 'c' 'r' 'b' 'm' 'w' 'y' 'g' 'c' 'r' 'b' 'm' 'w' 'y' 'g' 'c' 'r' 'b' 'm' 'w'];

axes(ineaxes);
for k=nninegraf;
    j=0;
    for i=ninegraf;
        j=j+1;
        if i==k;
           ineini=10*deltat(j);
           NNine=(numeroine(j+1)-numeroine(j))/6;
           inetf=ineini+(NNine-1)*inedeltat(j);
           X = [ineini:inedeltat(j):inetf];
           plot(X,Aine(numeroine(j)+2*NNine:(numeroine(j)+3*NNine-1),:),color(i));
           hold on;
        end;
    end;
end;  

if ineTmin~=0 | ineTmax~=0 | ineAmin~=0 | ineAmax~=0;
   axis([ineTmin ineTmax ineAmin ineAmax]);
end;

j=0;
for i=nninegraf;
    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('periodo (seg)');
ylabel('Amplitud ((cm/seg)^2)');
hold off;