paingra.m

地震资料处理的matlab源码软件包

function paingra

global Aine;
global inegraf;
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 ineUmon;
global ineBeta;
global Umax;
global UUmax;
global R;
global RR;
global HE;
global HHE;
global Xmax;
global XXmax;

j=0;
for i=ninegraf;
    j=j+1;
    if i==inegraf;
       ineini=10*deltat(j);
       NNine=(numeroine(j+1)-numeroine(j))/6;
       inedt=inedeltat(j);
       inetf=ineini+(NNine-1)*inedt;
    end;
end;
  
k=0;
for i=ineini:inedt:inetf;
    k=k+1;
    DPA  =[DPA  ,(UUmax(:,k)/ineUmon+ineBeta*HHE(:,k)./RR(:,k)/ineUmon./XXmax(:,k).*UUmax(:,k))];
end;

DPA(1,:)=DPA(1,:)*0;

[DDPA,I]=sort(DPA);

k=0;
for i=ineini:inedt:inetf;
    k=k+1;
    l=0;
    for j=[I(:,k)'];
        l=l+1;
        RRR(l,k)=RR(j,k);
    end;
end;

XPA=[0 0.10 0.20 0.30 0.40];

k=0;
for i=ineini:inedt:inetf;
    k=k+1;
    IDPA  =[IDPA  ;(interp1(DDPA(:,k),RRR(:,k),XPA))'];
end;

RIDPA=[R(1,:)' R(1,:)' R(1,:)' R(1,:)' R(1,:)']./IDPA;

j=0;
for i=ninegraf;
    j=j+1;
    if i==inegraf;
       Aine=[Aine(1:numeroine(j)-1+4*NNine,:);[IDPA;RIDPA];...
             Aine(numeroine(j+1):length(Aine),:)];
    end;
end;

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)+4*NNine:(numeroine(j)+5*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;