matlab计算旋度散度程序

以下是Matlab计算旋度和散度的程序示例：

% 定义三个矢量场
x = linspace(-1, 1, 10);
y = linspace(-1, 1, 10);
[X,Y] = meshgrid(x,y);
u = Y.*sin(X);
v = -X.*cos(Y);
w = X.*Y;

% 计算旋度和散度
[dx, dy] = gradient(u);
[dxx, ~] = gradient(dx);
[~, dyy] = gradient(dy);
curl = dyy - dxx;

[dx, dy] = gradient(v);
[dxx, ~] = gradient(dx);
[~, dyy] = gradient(dy);
curl = curl + dyy - dxx;

[dx, dy] = gradient(w);
[~, dxx] = gradient(dx);
[dyx, ~] = gradient(dy);
div = dx + dyx;

% 绘制矢量场、旋度和散度
figure;
quiver(X, Y, u, v); % 矢量场
title('Vector Field');

figure;
surf(X, Y, curl); % 旋度
title('Curl');

figure;
surf(X, Y, div); % 散度
title('Divergence');