gettrace.m

来自「这是matlab在地球物理数据处理方面的源码」· M 代码 · 共 20 行

function [trace]=gettrace(tracenum,shotnum)
if(shotnum>0&shotnum<0)
  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');
  if(tracenum>1)
    status=fseek(fid1,512*4*(tracenum-1),'cof');
  end;
  trace=fread(fid1,512,'float');
  status=fclose(fid1);
else
  'Wrong trace number'
end