详细说明梯度散度旋度Word文件下载.docx

详细说明梯度散度旋度Word文件下载.docx

《详细说明梯度散度旋度Word文件下载.docx》由会员分享，可在线阅读，更多相关《详细说明梯度散度旋度Word文件下载.docx（11页珍藏版）》请在冰豆网上搜索。

%计算函数值,zz也是25x25的矩阵

mesh（xx,yy,zz）:

%画岀立体网状图

■2.2

x=linspace（-2,2,25）;

%在x轴上取25点y=linspace（-2,2,25）;

%在y轴上取25点[xx,yy]=meshgrid（x,y）;

%xx和yy都是25x25的矩阵

zz=xx.*exp（-XX."

%计算函数值,zz也是25x25的矩阵

surf（xx,yy,zz）;

%画出立体曲面图

x=linspace（-2,2,25）;

%在x轴上取25点

y=linspace（-2,2,25）;

%在y轴上取25点

[xx,yy]=meshgrid（x,y）:

zz二sqrt（xx."

2+yy."

%函数变化，zz也是25x25的矩阵

mesh（xx,yy,zz）;

%画出立体网状图

%在x轴上取25点

y=linspace（-2,2,25）;

%在y轴上

取25点

[xx,yy_=meshgrid（x,y）;

%xx和

yy都是25x25的矩阵

zz=sqrt（xx."

2）;

%计算函数

值，zz也是25x25的矩阵

%画出立体曲面图梯度

%在x轴上取25点y=linspace（-2,2,25）:

%在y轴上取25点[xx,yy]二meshgrid（x,y）;

zz二xx・*exp（一xx・"

2一yy."

h二contour（zz,12）;

%画等高线clabel（h）:

%标注等高线值

[dx,dy]=gradient（zz,.2,.2）;

%求梯度holdon;

quiver（dx,dy）:

%画矢量箭头

y轴上取25点［xx,yy］=meshgrid（x,y）;

%xx和yy都是25x25的矩阵zz二xx・*exp（一xx・"

contour（zz）;

%特定范围内的等高线［dx,dy］二gradient（zz,.2,.2）;

holdon;

%在x轴上取25点y=linspace（-2,2,25）;

%在y轴上取25点［xx,yy］二meshgrid（x,y）;

zz=sqrt（xx・"

［dx,dy］二gradient（zz,・2八2）;

quiver（dx,dy）;

%在x轴上取25点y二1inspace（-2,2,25）:

y轴上取25点［xx,yy］=meshgrid（x,y）;

%特定范围内的等高线［dx,dy］二gradient（zz八2八2）;

holdon;

quiver（dx,dy）

[x,y,z]=meshgrid（~3:

1:

3）;

%构建7*7矩阵u二3.*x."

2+2.*y."

2+z."

2;

%构建三维矢量场（U,V,W）v=x."

2+3.*y."

2+2.*z."

w=2.*x."

2+乙"

div=divergence（x,y,z,u,v,w）;

%求（U,V,W）的散度figure%显示散度图像

fori=l：

7%根据x,y,z值系统自动分成7个7*7矩阵的散度值

%figure%7个图层单独显示

mesh（div（:

:

i））;

%以散度值为点做三维图

%不被下一个图覆盖end

旋度

10

S

to

io

■

v=x."

2+3・*y."

w=2.*x."

[curlx,curly,curlz]=curl（x,y,z,u,v,w）;

figure%旋度有3个矢量分别有7个7*7矩阵

fori二1:

7

mesh（cur1x（:

:

i））;

end

figure

fori=l:

mesh（curly（:

mesh（curlz（: