svgra.m

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

function svgra

global Ares;
global nresgraf;
global nnresgraf;
global resaxes;
global resdeltat;
global archivo;
global letras;
global numerores;
global resTmax;
global resTmin;
global resAmax;
global resAmin;
global restipo;
       
color=['y' 'g' 'c' 'r' 'b' 'm' 'w' 'y' 'g' 'c' 'r' 'b' 'm' 'w' 'y' 'g' 'c' 'r' 'b' 'm' 'w'];

axes(resaxes);
for k=nnresgraf;
    j=0;
    for i=nresgraf;
        j=j+1;
        if i==k;
           restf=resdeltat(j)*(numerores(j+1)-numerores(j));
           X = [resdeltat(j):resdeltat(j):restf];
           if restipo==1;
              loglog(X,Ares(numerores(j):(numerores(j+1)-1),[2]),color(i));
           else;
              plot(X,Ares(numerores(j):(numerores(j+1)-1),[2]),color(i));
           end;
           hold on;
        end;
    end;
end;  

if resTmin~=0 | resTmax~=0 | resAmin~=0 | resAmax~= 0;
   axis([resTmin resTmax resAmin resAmax]);
end;

j=0;
for i=nnresgraf;
    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)');
hold off;