matlab求散度的代码

在MATLAB中，可以使用gradient函数来计算向量场的散度。下面是一个简单的示例代码：

% 定义向量场
x = linspace(-2, 2, 20);
y = linspace(-2, 2, 20);
[X, Y] = meshgrid(x, y);
U = -Y;
V = X;

% 计算散度
[DX, DY] = gradient(U, x, y);
divergence = DX + DY;

% 绘制向量场和散度图
figure;
quiver(X, Y, U, V);
hold on;
contour(X, Y, divergence);
xlabel('x');
ylabel('y');
title('Vector Field and Divergence');
legend('Vector Field', 'Divergence');

这段代码首先定义了一个二维向量场，然后使用gradient函数计算了向量场的散度。最后，通过quiver函数绘制了向量场，通过contour函数绘制了散度图。