getshot.m

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

function [shot]=getshot(shotnum)
if(shotnum>0&shotnum<86)
  s=num2str(shotnum);
  s1=s;
  if(shotnum<10)
    s1='0';
    s1(2:2)=s;
  end
  filen='shot00.dat';
  filen(5:6) = s1;
  fid1=fopen(filen,'r','b');
  shot1=fread(fid1,512*48,'float');
  shot=reshape(shot1,512,48);
  clear shot1;
  status=fclose(fid1);
else
  'Wrong shot number'
   break
end;