Вычислите расхождение векторного поля
вычисляет числовое расхождение 3-D векторного поля с векторными компонентами div
= divergence(X
,Y
,Z
,Fx
,Fy
,Fz
)Fx
, Fy
, и Fz
.
Массивы X
Y
, и Z
, которые задают координаты для векторных компонентов Fx
, Fy
, и Fz
, должно быть монотонным, но не должны быть расположенными равными интервалами. X
Y
, и Z
должны быть трехмерные массивы, одного размера, который может быть произведен meshgrid
.
принимает сетку по умолчанию точек выборки. Узлы решетки по умолчанию div
= divergence(Fx
,Fy
,Fz
)X
Y
, и Z
определяются выражением [X,Y,Z] = meshgrid(1:n,1:m,1:p)
, где [m,n,p] = size(Fx)
. Используйте этот синтаксис, когда это необходимо, чтобы сохранить память, и не касаются абсолютных расстояний между точками.
вычисляет числовое расхождение 2D векторного поля с векторными компонентами div
= divergence(X
,Y
,Fx
,Fy
)Fx
и Fy
.
Матрицы X
и Y
, которые задают координаты для Fx
и Fy
, должно быть монотонным, но не должны быть расположенными равными интервалами. X
и Y
должны быть 2D матрицы, одного размера, который может быть произведен meshgrid
.
divergence
вычисляет частные производные в его определении при помощи конечных разностей. Для внутренних точек данных частные производные вычислены с помощью центральной разности. Для точек данных вдоль ребер частные производные вычисляются с помощью одностороннего (прямого) различия.
Например, рассмотрите 2D векторное поле F, который представлен матрицами Fx
и Fy
в местоположениях X
и Y
с размером m
- n
. Местоположения являются 2D сетками, созданными [X,Y] = meshgrid(x,y)
, где x
вектор из длины n
и y
вектор из длины m
расхождение
затем вычисляет частные производные ∂Fx / ∂x и ∂Fy / ∂y как
dFx(:,i) = (Fx(:,i+1) - Fx(:,i-1))/(x(i+1) - x(i-1))
и
dFy(j,:) = (Fy(j+1,:) - Fy(j-1,:))/(y(j+1) - y(j-1))
для внутренних точек данных.
dFx(:,1) = (Fx(:,2) - Fx(:,1))/(x(2) - x(1))
и
dFx(:,n) = (Fx(:,n) - Fx(:,n-1))/(x(n) - x(n-1))
для точек данных на левых краях и правых краях.
dFy(1,:) = (Fy(2,:) - Fy(1,:))/(y(2) - y(1))
и
dFy(m,:) = (Fy(m,:) - Fy(m-1,:))/(y(m) - y(m-1))
для точек данных в верхних краях и базовых краях.
Числовое расхождение векторного поля равно div = dFx + dFy
.
streamtube
| gradient
| curl
| isosurface