Global Index (short | long) | Local contents | Local Index (short | long)
[scale,maxlen] = quiver ( arg1, arg2, arg3, arg4, arg5, arg6 )
QUIVER Quiver plot.
function [scale,maxlen] = quiver ( arg1, arg2, arg3, arg4, arg5, arg6 )
QUIVER(X,Y,DX,DY) draws little arrows at every (X,Y) pair in
matrices X and Y. The (DX,DY) pairs in matrices DX and DY
determine the direction and magnitude of the arrows.
QUIVER(x,y,DX,DY), with two vector arguments replacing the first
two matrix arguments, must have length(x) = n and
length(y) = m where [m,n] = size(DX) = size(DY). In this case, the
arrows are the quads (x(j), y(i), DX(i,j), DY(i,j)).
Note that x corresponds to the columns of DX and DY and y corresponds
to the rows.
QUIVER(DX,DY) uses x = 1:n and y = 1:m. In this case DX and DY
are defined over a geometrically rectangular grid.
QUIVER(X,Y,DX,DY,S) and QUIVER(DX,DY,S) apply scalar S as a scale
factor to the lengths of the arrow. For example, S = 2 doubles
their relative length and S = 0.5 halves them.
A final trailing string argument specifies linetype and color using
any legal line specification as described under the plot command.
For example, try
xord = -2:.2:2;
yord = -2:.2:2;
[x,y] = meshgrid(xord,yord);
z = x .* exp(-x.^2 - y.^2);
[px,py] = gradient(z,.2,.2);
contour(x,y,z),hold on, quiver(x,y,px,py), hold off
See also GRADIENT, COMPASS, FEATHER, ROSE.
------------------------------------------------------------------------------
Enhancement by Alexis Lau 7/12/93:
*) Arrows centered on the grid points
*) Arrow head closer to the shaft
*) Extend meaning of S:
positive or unspecified - same as in factory setting
negative - absolute scaling
e.g. (consider you are plotting a wind field, a wind
vector with 1 unit amplitude will be mapped to abs(S)
units of the axes in the plot).
See also CY_QUIVER, QUIVER2
------------------------------------------------------------------------------
8/31/93
*) allow to use global variable QUIVERWIDTH to control shaft width
similar to COASTWIDTH in 'drawcoasts'
------------------------------------------------------------------------------
function [scale,maxlen] = quiver ( arg1, arg2, arg3, arg4, arg5, arg6 )
global QUIVERWIDTH
if isempty(QUIVERWIDTH);
QUIVERWIDTH=get(gca,'defaultlinelinewidth');
end
if exist('skip')~=1; skip = 1.e-18; end;
% Charles R. Denham, MathWorks 3-20-89
% Modified 12-19-91, LS.
% Copyright (c) 1984-92 by The MathWorks, Inc.
xx = [0 1 .70 1 0.70].' - 0.5;
yy = [0 0 .15 0 -.15].';
arrow = xx + yy.*sqrt(-1);
eval(['last = arg' int2str(nargin) ';']);
if isstr(last)
eval(['style = arg' int2str(nargin), ';']);
narg = nargin - 1;
else
style = 'r-';
narg = nargin;
eval(['last = arg' int2str(narg-1) ';']);
if isstr(last)
error('Only the final argument can be a string.');
end
end
if narg == 0
error('First 2 or 4 arguments must be numeric.')
end
eval(['lastdim = min(size(arg' int2str(narg) '));']);
if lastdim == 1
if isstr(eval(['arg' int2str(narg)]))
error('Scalar scale argument expected.')
end
eval(['scale = arg' int2str(narg) ';']);
narg = narg - 1;
else
scale = 1;
end
if narg == 0
error('First 2 or 4 arguments must be numeric.')
end
if min(size(arg1)) > 1
[m,n] = size(arg1);
else
m = max(size(arg1));
n = max(size(arg2));
end
if narg == 2
if isstr(arg1) | isstr(arg2)
error('Input must be numeric.')
end
[xx,yy] = meshgrid(1:n, 1:m);
px = arg1;
py = arg2;
else
if isstr(arg1) | isstr(arg2) | isstr(arg3) | isstr(arg4)
error('Input must be numeric.')
end
if min(size(arg1)) == 1 & min(size(arg2)) == 1
[xx,yy] = meshgrid(arg1,arg2);
else
xx = arg1;
yy = arg2;
end
px = arg3;
py = arg4;
end
px = px(:); py = py(:);
z = (px + py.*sqrt(-1)).';
len = abs(z); phi = angle(z);
z = len.*cos(phi) + sqrt(-1)*len.*sin(phi);
grid = xx + yy.*sqrt(-1); grid = grid(:);
equal_len = 0;
if scale > 0
% figure out delx and dely so spacing is accounted for in z
delx = xx(1,2)-xx(1,1);
dely = yy(2,1)-yy(1,1);
if dely == 0; dely = delx; end;
maxlen = max(sqrt((px/delx).^2+(py/dely).^2));
scale = scale*0.90 ./ maxlen;
elseif scale < 0;
% absolute scaling
scale = abs(scale);
maxlen = 0.9 / scale;
else
% same length for all vectors
maxlen = 3; equal_len = 1;
end
%%%%%
pa=sqrt(px.^2 + py.^2); pa=(find(pa>=skip));
z=z(pa); grid=grid(pa);
if equal_len;
z = z ./ abs(z) * maxlen;
end
%%%%%
if scale~=0;
a = scale * arrow * z + ones(5,1) * grid.';
else
a = arrow * z + ones(5,1) * grid.';
end
% append nan's so we get one handle
a = [a; nan*ones(1,size(a,2))];
a = a(:);
cax = newplot;
curlinewidth = get(gca,'defaultlinelinewidth');
set(gca,'defaultlinelinewidth',QUIVERWIDTH);
plot(real(a), imag(a), style);
set(gca,'defaultlinelinewidth',curlinewidth);
next = lower(get(cax,'NextPlot'));
if ~ishold
minx = min(min(xx));
miny = min(min(yy));
maxx = max(max(xx));
maxy = max(max(yy));
axis([minx maxx miny maxy]);
view(0,90);
end