ariacal.m

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

function ariacal

global numero;
global A;
global Aria;
global ariagraf;
global nariagraf;
global nnariagraf;
global deltat;
global numeroaria;
global directorio;
global directorio1;

nnariagraf=ariagraf;

No=numero(ariagraf);
Nf=numero(ariagraf+1)-1;
Nn=Nf-No+1;

%Sum=0;
%Yi=A(No);
%for i=No:Nf;
%    Yj=A(i)*A(i);
%    Sum=Sum+(Yi+Yj);
%    Iaria=[Iaria;Sum];
%    Yi=Yj;
%end;

%cd d:\matlab\toolbox\prima
cd (directorio1);

fid=fopen('Acel.dat','wt');
fprintf(fid,'%12.4e\n',A(No:Nf));
fclose(fid);

fid=fopen('Sdat.dat','wt');
fprintf(fid,'%8i',Nn);
fclose(fid);

!ariacal

load Iaria.sal;
delete Sdat.dat;
delete Iaria.sal;
delete Acel.dat;

Iaria=Iaria/2*deltat(ariagraf);
       
ind=0;
j=0;
for i=nariagraf;
    j=j+1;
    if i==ariagraf;
       Aria=[Aria(1:numeroaria(j)-1);Iaria;...
             Aria(numeroaria(j+1):length(Aria))];
       dif=Nn-numeroaria(j+1)+numeroaria(j);
       if j==length(nariagraf);
          numeroaria=[numeroaria(1:j);numeroaria(j)+Nn];
       else; 
          numeroaria=[numeroaria(1:j);numeroaria(j)+Nn;...
                    (numeroaria(j+2:length(numeroaria))+dif)];
       end;
       ind=1;
       break; 
    end;
end;

if ind==0;
   Aria=[Aria;Iaria];
   nariagraf=[nariagraf ariagraf];
   numeroaria=[numeroaria;numeroaria(length(numeroaria))+Nn];
end;

cd (directorio);

ariagra;