wigb函数matlab中画地震剖面用

function wigb （ascalxzamx）
%WIGB: Plot seismic data using wiggles
%
%  WIGB（ascalxzamx） 
%
%  IN    a: seismic data
%        scale: multiple data by scale
%        x: x-axis;
%        z: vertical axis （time or depth）
%	 x and z are vectors with offset and time.
%
%	 If only ‘a‘ is enter ‘scalxzamnamx‘ are decided automatically; 
%	 otherwise ‘scal‘ is a scalar; ‘x z‘ are vectors for annotation in 
%	 offset and time amx are the amplitude range.


if nargin == 0 nx=10;nz=10; a = rand（nznx）-0.5; end;

[nznx]=size（a）;

trmx= max（abs（a））;
if （nargin <= 4）; amx=mean（trmx）;  end;
if （nargin <= 2）; x=[1:nx]; z=[1:nz]; end;
if （nargin <= 1）; scal =1; end;

if nx <= 1; disp（‘ ERR:PlotWig: nx has to be more than 1‘）;return;end;

 % take the average as dx