putcmp.m

这是matlab在地球物理数据处理方面的源码

function []=PutCMP(str,CMPNumber,traces,offset)
%
% str=string of 7 character ????000
% example nCMP000
% CMPNumber: integer CMP number
% traces should contain a CMP gather
% offset should contain offsets of the CMP gather
% To run  []=PutCMP(str,CMPNumber,traces,offset)
% Example [] = PutCMP('stuf000',29,traces,offset) saves the traces and 
% their offset into a file called stuf029.dat.
% These data can then be read back into matlab in via the program getcmp.m 
%
ntraces=length(offset);
filename=str;
filename(8:11)='.dat';
s=num2str(CMPNumber);
if(CMPNumber<10);
   str(7:7)=s
elseif (CMPNumber<100);
   str(6:7)=s
elseif(CMPNumber<999)
   str(5:7)=s;
else
  'error'
  break
end
fid=fopen(filename,'w','l');
for i=1:ntraces
  fwrite(fid,offset(i),'float');
  fwrite(fid,traces(:,i),'float');
end  
fclose(fid);
return