用自定义函数colorgrad给出向量梯度图

你可以使用 MATLAB 自带的 `gradient` 函数来计算向量场的梯度，然后使用 `quiver` 函数绘制向量场。以下是一个示例代码，使用自定义函数 `colorgrad` 来绘制向量梯度图：

% 定义向量场
[x, y] = meshgrid(-2:.2:2, -2:.2:2);
u = y;
v = -x;

% 计算向量场的梯度
[dx, dy] = gradient(u, v);

% 计算梯度大小和方向
mag = sqrt(dx.^2 + dy.^2);
angle = atan2(dy, dx);

% 绘制向量梯度图
figure;
color_gradient = colorgrad(angle);
quiver(x, y, u, v);
colormap(color_gradient);
colorbar;

在上面的代码中，我们首先定义了一个向量场 `u` 和 `v`，然后使用 `gradient` 函数计算了向量场的梯度 `dx` 和 `dy`。接下来，我们使用 `sqrt` 和 `atan2` 函数计算了梯度大小和方向，并使用自定义函数 `colorgrad` 生成了一个颜色映射。最后，我们使用 `quiver` 函数绘制了向量场，并使用 `colormap` 和 `colorbar` 函数将颜色映射添加到图像中。当你运行代码后，你将会看到一个向量梯度图的图像。