注册 |登录
已有 11183 次阅读 2017-9-30 11:20 |个人分类:教学辅导|系统分类:科研笔记
先定义函数
f = @(x,y) x.^2 + y.^2;
求偏导数
[xx,yy] = meshgrid(-5:0.25:5); [fx,fy] = gradient(f(xx,yy),0.25);
切点(x0,y0) = (1,2). 函数值 f(1,2) = 5.
(x0,y0) = (1,2)
f(1,2) = 5
x0 = 1; y0 = 2; t = (xx == x0) & (yy == y0); indt = find(t); fx0 = fx(indt); fy0 = fy(indt);
创建切平面方程
z = @(x,y) f(x0,y0) + fx0*(x-x0) + fy0*(y-y0);
画出曲面和切平面
surf(xx,yy,f(xx,yy),'EdgeAlpha',0.7,'FaceAlpha',0.9) hold on surf(xx,yy,z(xx,yy)) plot3(1,2,f(1,2),'r*')
换个角度看一下
view(-135,9)
推荐到博客首页
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2024-9-16 08:25
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社