内容简介作者简介 - 科学与工程计算系

nx =
6.2968e+000
nxx =
6.2918e+000
nxxx =
6.3356e+000
(3)
e=norm(y-A*x),ee=norm(y-A*xx),eee=norm(y-A*xxx)
e =
6.9863e-001
ee =
4.7424e-002
eee =
4.7424e-002
4.5 函数的数值导数和切平面
4.5.1 法线
【例 4.5.1-1】曲面法线演示。
y=-1:0.1:1;x=2*cos(asin(y));
[X,Y,Z]=cylinder(x,20);
surfnorm(X(:,11:21),Y(:,11:21),Z(:,11:21));
view([120,18])
1.2
0.8
0.6
0.4
0.2
-0.2
-3
4.5.2 偏导数和梯度
4.5.2.1 理论定义
4.5.2.2 数值计算指令
1
0
-2
-1
0
图 4.5.1-1
【例 4.5.2.2-1】用一个简单矩阵表现 diff 和 gradient 指令计算方式。
F=[1,2,3;4,5,6;7,8,9]
Dx=diff(F)
Dx_2=diff(F,1,2)
[FX,FY]=gradient(F)
[FX_2,FY_2]=gradient(F,0.5)
F =
1 2 3
4 5 6
7 8 9
7
1
5
0
-5