divergence - Compute divergence of vector field

Syntax

Description

div = divergence(X,Y,Z,U,V,W) computes the divergence of a 3-D vector field U, V, W.

The arrays X, Y, and Z, which define the coordinates for U, V, and W, must be monotonic, but do not need to be uniformly spaced. X, Y, and Z must have the same number of elements, as if produced by meshgrid.

div = divergence(U,V,W) assumes X, Y, and Z are determined by the expression

[X Y Z] = meshgrid(1:n,1:m,1:p)

where [m,n,p] = size(U).

div = divergence(X,Y,U,V) computes the divergence of a 2-D vector field U, V.

The arrays X and Y, which define the coordinates for U and V, must be monotonic, but do not need to be uniformly spaced. X and Y must have the same number of elements, as if produced by meshgrid.

div = divergence(U,V) assumes X and Y are determined by the expression

[X Y] = meshgrid(1:n,1:m)

where [m,n] = size(U).

Examples

This example displays the divergence of vector volume data as slice planes, using color to indicate divergence.

load wind
div = divergence(x,y,z,u,v,w);
h = slice(x,y,z,div,[90 134],59,0);
shading interp
daspect([1 1 1])
axis tight
camlight
set([h(1),h(2)],'ambientstrength',.6)

See Also

curl | isosurface | streamtube

How To

• Example – Displaying Divergence with Stream Tubes

© 1984-2012- The MathWorks, Inc.    -   Site Help   -   Patents   -   Trademarks   -   Privacy Policy   -   Preventing Piracy   -   RSS