x=linspace(-2, 2, 10); %在-2到2取25个点
y=linspace(-2, 2, 10); %在-2到2取25个点
z=linspace(-2, 2, 10);
[xx, yy, zz]=meshgrid(x, y, z); %生成网格采样点
r = sqrt( xx.^2 + yy.^2 +zz.^2); %生成矩阵Z
% h=contour(xx,yy,zz, 12); %以12个等高线层级绘制矩阵Z的等高线图
% clabel(h); %写等高线的值
[dx, dy, dz]=gradient(r,0.2,0.2,0.2); %求梯度
grid on
hold on %做下一幅图时保持原来图像
quiver3(xx,yy,zz,dx,dy,dz); %画矢量图箭头
axis equal; %等比例显示
散度(二维)的仿真
syms x y z real %定义符号变量
F = [ cos(x+2*y), sin(x-2*y) ]; %定义函数F
g = divergence(F,[x y]); %求函数F的散度,符号形式
divF=matlabFunction(g); %将散度转换为函数形式
x = linspace(-2.5,2.5,20);
y = linspace(-2.5,2.5,20);
[X