function  seiplot(trace,sample,disoffset,dt,amp);

%%%%%%%%%%%定义画地震图的函数seiplot%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%trace 是地震道数，sample是一道采样点数，
%disoffset是检波距，dt采样时间,amp 是地震振幅扩大的倍数
%trace=601;
%sample=1001;
%disoffset=10;
%dt=0.004;
%%%%%%%%%%读取十进制数据%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

fid=fopen('output.dat','r');%%%%%%%output.dat是十进制无道头数据%
temp_data=fscanf(fid,'%f');
for n=1:trace
    for m=1:sample
        md2(m,n)=temp_data((n-1)*sample+m);
    end
end

%%%%%%%%%%%%%%%%%%%%%%%画图部分%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
xx=(0:(m-1))*dt*0.001;
for n=1:trace
    md(:,n)=md2(:,n)*amp+(n-1)*disoffset;
    fill(md(:,n),xx,'k');
    hold on;
end
axis ij;
set(gca,'FontName','times','FontSize',14);
ylabel('Time (s)','FontName','times','FontSize',14);
xlabel('Position (m)','FontName','times','FontSize',14);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%注，在使用此函数是只需在MATLAB命令窗口中输入%%%%%%%%%%%%%

%%%%seiplot(601,1001,10,0.004,30)即可%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%运行结果%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

http://s5/bmiddle/674956e7h73f91a0606e4&690