www.pudn.com > 87361026FEA.rar > contura.m, change:2002-10-01,size:2344b

function [c,h]=contoura(x,y,z,a,b),
% CONTOURA CONTOURA plot.
%       CONTOURA(Z) is a contour plot of matrix Z treating the values in Z
%       as heights above a plane.
%       CONTOURA(X,Y,Z), where X and Y are vectors and Z is a matrix, specifies
%       the X- and Y-axes used on the plot.
%       CONTOURA(X,Y,Z), where X, Y, and Z are vectors, interpolates
%       the data onto a rectangular grid.
%       CONTOURA(Z,N) and CONTOURA(X,Y,Z,N) draw N contour lines,
%       overriding the default automatic value.
%       CONTOURA(Z,V) and CONTOURA(X,Y,Z,V) draw LENGTH(V) contour lines
%       at the values specified in vector V.
%       CONTOURA(...,'linetype') draws with the color and linetype specified,
%       as in the PLOT command.
%
%       C = CONTOURA(...) returns contour matrix C as described in CONTOURC
%       and used by CLABEL.
%       [C,H] = CONTOURA(...) returns a column vector H of handles to LINE
%       objects, one handle per line.
%
%       NOTE: Do not rename this CONTOUR, it calls CONTOUR.

%       D. Thomas 11/17/94
%       Copyright (c) 1984-94 by The MathWorks, Inc.
A = 'xyzab';
C = ',,,, ';

args = [A(1:nargin);[C(1:nargin-1) ' ']];
args = [args(:)' ' '];

if (nargout == 1),
argo = 'c=';
elseif (nargout == 2),
argo = '[c,h]=';
else
argo = '';
end

if (nargin > 2)
if (min(size(x)) == 1 & min(size(y)) == 1 & min(size(z)) == 1),

if  ~(size(x(:),1) == size(y(:),1) & size(y(:),1) == size(z(:),1)),
error('Sizes of X, Y, and Z must be equal');
end
xmin=min(x(:));
ymin=min(y(:));
xmax=max(x(:));
ymax=max(y(:));
N = ceil(length(x(:))^(1/2));
Xl = xmin:((xmax-xmin)/(N-1)):xmax;
Yl = ymin:((ymax-ymin)/(N-1)):ymax;
[Xi,Yi]=meshgrid(Xl,Yl);
Z = griddata(x,y,z,Xi,Yi);
if (length(args) > 7),
args = [',',args];
end
eval([argo,'contour(Xl,Yl,Z',args(7:length(args)),');']);
else,
eval([argo,'contour(' args ');']);
end
else,
eval([argo,'contour(' args ');']);
end