Documentation

divergence

Compute divergence of vector field

Syntax

div = divergence(X,Y,Z,U,V,W)
div = divergence(U,V,W)
div = divergence(X,Y,U,V)
div = divergence(U,V)

Description

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

The arrays X, Y, and Z, which define the coordinates for the vector components 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.

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

collapse all

Display the divergence of vector volume data as slice planes. Use color to indicate divergence.

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