wiggle(1).m

这是matlab在地球物理数据处理方面的源码

function handle = wiggle(z,taxis,zaxis,sc,ornt)
%WIGGLET1	[ A.K.Booer 24-Mar-1993, Modified JHaverstock March 20, 1996 ]
%
%   wiggle(x,taxis,zaxis) - plots a shaded wiggle plot
%
%   x - matrix of waveform columns
%   taxis - time axis vector
%   zaxis - space axis vector
% 
%   wiggle(x,t,z,mag)   - scales waveforms by given magnification
%   wiggle(x,t,z,m,'d') - changes time axis orientation to 'downwards'
%
%   WIGGLE returns handle to created graphics objects, so that
%   set(wiggle(x,t,z), 'face',c1, 'edge',c2') will specify colours.
% 

% Modifications:
% 24-Mar-1993   akb   Based on "wp" of 2-Oct-1992
%
[n,m] = size(z);
n1 = n:-1:1;
if nargin < 2, taxis = 1:n; end
if nargin < 3, zaxis = 1:m; end
if nargin < 4, sc = 1; end
if nargin < 5, ornt = 'a'; end	% default orientation "across"
colour = 'w'; 				% default colour
zaxis = zaxis(:)';
% 
z = z * (2*mean(diff(zaxis))) * (sc / max(max(z) - min(z))); 	% scale
%
if ornt(1) == 'd'
  h = fill(zaxis(ones(2*n,1),:)+[z;min(z(n1,:),0)],taxis([1:n n1]),colour);...
      set(h,'edge',colour), set(gca,'Ydir','r')
  else
  h = fill(taxis([1:n n1]),zaxis(ones(2*n,1),:)-[z;min(z(n1,:),0)],colour);...
      set(h,'edge',colour), set(gca,'Ydir','r')
end
%
if nargout > 0, handle = h; end
%