psdascii.m

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

function psdascii

global Apsd;
global Atf;
global Acor;
global Acorxy;
global psdgrafX;
global psdgrafY
global psddeltaf;
global tfdeltaf;
global deltat;
global npsdgraf;
global ntfgraf;
global ncorgraf;
global ncorxygraf;
global numeropsd;
global numerotf;
global numerocor;
global numerocorxy;
global exten;
global psdX;
global corX;
global psdXY;
global corXY
global tfXY;

exten='.psd'; 

if psdgrafY==[];
   j=0;
   for i=npsdgraf;
       j=j+1;
       if i==psdgrafX;
          psdX=[psddeltaf(j) psddeltaf(j) psddeltaf(j);...
                Apsd(numeropsd(j):numeropsd(j+1)-1,[1]).*deltat(i),...
                (Apsd(numeropsd(j):numeropsd(j+1)-1,[1])+Apsd(numeropsd(j):numeropsd(j+1)-1,[2])).*deltat(i),...
                (Apsd(numeropsd(j):numeropsd(j+1)-1,[1])-Apsd(numeropsd(j):numeropsd(j+1)-1,[2])).*deltat(i)]';
       end;
   end;

   j=0;
   for i=ncorgraf;
       j=j+1;
       if i==psdgrafX;
          corX=[deltat(i);Acor(numerocor(j):numerocor(j+1)-1)]';
       end;
   end;

else;

   j=0;
   for i=npsdgraf;
       j=j+1;
       if i==psdgrafX;
          psdXY=[psddeltaf(j) psddeltaf(j) psddeltaf(j);...
                Apsd(numeropsd(j):numeropsd(j+1)-1,[1]),...
                Apsd(numeropsd(j):numeropsd(j+1)-1,[1])+Apsd(numeropsd(j):numeropsd(j+1)-1,[2]),...
                Apsd(numeropsd(j):numeropsd(j+1)-1,[1])-Apsd(numeropsd(j):numeropsd(j+1)-1,[2])];
       end;
   end;

   j=0;
   for i=npsdgraf;
       j=j+1;
       if i==psdgrafY;
          psdXY=[[psdXY] [psddeltaf(j) psddeltaf(j) psddeltaf(j);...
                Apsd(numeropsd(j):numeropsd(j+1)-1,[1]),...
                Apsd(numeropsd(j):numeropsd(j+1)-1,[1])+Apsd(numeropsd(j):numeropsd(j+1)-1,[2]),...
                Apsd(numeropsd(j):numeropsd(j+1)-1,[1])-Apsd(numeropsd(j):numeropsd(j+1)-1,[2])]]';
       end;
   end;

   j=0;
   for i=ntfgraf(1,:);
       j=j+1;
       k=ntfgraf(2,j);
       if i==psdgrafX & k==psdgrafY;
          tfXY=[tfdeltaf(j) tfdeltaf(j) tfdeltaf(j);...
                abs(Atf(numerotf(j):numerotf(j+1)-1,[2])),...
                angle(Atf(numerotf(j):numerotf(j+1)-1,[2])).*180/pi,...
                Atf(numerotf(j):numerotf(j+1)-1,[3])]';
       end;
   end;

   j=0;
   for i=ncorgraf;
       j=j+1;
       if i==psdgrafX;
          corXY=[deltat(i);Acor(numerocor(j):numerocor(j+1)-1)];
       end;
   end;

   j=0;
   for i=ncorgraf;
       j=j+1;
       if i==psdgrafY;
          corXY=[[corXY] [deltat(i);Acor(numerocor(j):numerocor(j+1)-1)]];
       end;
   end;

   j=0;
   for i=ncorxygraf(1,:);
       j=j+1;
       k=ncorxygraf(2,j);
       if i==psdgrafX & k==psdgrafY;
          corXY=[[corXY] [deltat(i);Acorxy(numerocorxy(j):numerocorxy(j+1)-1)]]';
       end;
   end;
end;

ascii;

clc;