地震数据的MATLAB中读取程序

function s=segread（fn）
%s=segread（fn）
echo all off
fr=fopen（fn）;
if fr==0
    disp（‘文件打开失败！‘）;
end
fseek（fr4‘bof‘）;
fbs=fread（fr1‘int16‘）;
s=sprintf（‘文件描述块大小：%d （Byte）‘fbs）;
disp（s）;
tnum=fread（fr1‘int16‘）;
s=sprintf（‘共有：%d 道地震记录‘tnum）;
disp（s）;
fseek（fr32‘bof‘）;
tp=fread（frtnum‘int32‘）;
disp（‘