Matlab画地球剖面图,分享用matlab显示地震记录的波形变面积图

wigb函数是网上找的，自己添了个主程序方便使用。

主程序

function PlotRecord_Wigb(filename,nx,nt,dx,dt,zy)

%波形变面积显示地震记录

% clear; clf;clc;

% nx=200; %道数

% nt=2001; %采样点数

% dx=10.;     %道间距

% dt=0.001;   %时间采样间隔

% zy=1.0;     %增益

% filename='r200-2001-1ms.dat';

clf;

fid=fopen(filename,'r');

[A,COUNT]=fread(fid,[nt,nx],'float32');

fclose(fid);

x=0:dx: (nx-1)*dx;

z=0:dt: (nt-1)*dt;

wigb(A,zy,x,z);

end

子函数

function wigb (a,scal,x,z,amx)

%WIGB: Plot seismic data using wiggles.

%

%  WIGB(a,scal,x,z,amx)

%

%  IN    a:     seismic data

%        scale: multiple data by scale

%        x:     x-axis;

%        z:     vertical axis (time or depth)

%

%