ruingra.m

来自「地震资料处理的matlab源码软件包」· M 代码 · 共 57 行

function ruingra

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;
       
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];
           V2 = polyfit(X',Aine(numeroine(j)+3*NNine:(numeroine(j)+4*NNine-1),2),3);
           V3 = polyfit(X',Aine(numeroine(j)+3*NNine:(numeroine(j)+4*NNine-1),3),3);
           V4 = polyfit(X',Aine(numeroine(j)+3*NNine:(numeroine(j)+4*NNine-1),4),3);
           V5 = polyfit(X',Aine(numeroine(j)+3*NNine:(numeroine(j)+4*NNine-1),5),3);
           plot(X,Aine(numeroine(j)+3*NNine:(numeroine(j)+4*NNine-1),:),color(i));
           hold on;
           plot(X,polyval(V2,X),'w');
           plot(X,polyval(V3,X),'w');
           plot(X,polyval(V4,X),'w');
           plot(X,polyval(V5,X),'w');
        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');
hold off;