函数的数值导数和切平面.doc

数学实验第五次课 微积分实验 函数的数值导数和切平面 法线 【例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]) 图 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) 【例 4.5.2.2-2】研究偶极子(Dipole)的电势（Electric potential）和电场强度（Electric field density）。设在处有电荷，在处有电荷。那么在电荷所在平面上任何一点的电势和场强分别为，。其中。。又设电荷，，。 clear;clf;q=2e-6;k=9e9;a=1.5;b=-1.5;x=-6:0.6:6;y=x; [X,Y]=meshgrid(x,y); rp=sqrt((X-a).^2+(Y-b).^2);rm=sqrt((X+a).^2+(Y+b).^2); V=q*k*(1./rp-1./rm); [Ex,Ey]=gradient(-V); AE=sqrt(Ex.^2+Ey.^2);Ex=Ex./AE;Ey=Ey./AE; cv=linspace(min(min(V)),max(max(V)),49); contourf(X,Y,V,cv,k-) %axis(square) title(\fontname{隶书}\fontsize{22}偶极子的场),hold on quiver(X,Y,Ex,Ey,0.7) plot(a,b,wo,a,b,w+) plot(-a,-b,wo,-a,-b,w-) xlabel(x);ylabel(y),hold off 图 4.5.2.2-1 函数的零点 多项式的根 一元函数的零点 利用MATLAB作图指令获取初步近似解 任意一元函数零点的精确解 【例 4.6.2.2-1】通过求的零点，综合叙述相关指令的用法。 （1） y=inline(sin(t)^2*exp(-a*t)-b*abs(t),t,a,b); %1 （2） a=0.1;b=0.5;t=-10:0.01:10; y_char=vectorize(y); % 3 Y=feval(y_char,t,a,b); clf,plot(t,Y,r);hold on,plot(t,zeros(size(t)),k); xlabel(t);ylabel(y(t)),hold off 图4.6-1 （3） 由于Notebook中无法实现zoom、ginput指令涉及的图形和鼠标交互操作，因此下面指令必须在MATLAB指令窗中运行，并得到如图4.6-2所示的局部放大图及鼠标操作线。 zoom on [tt,yy]=ginput(5);zoom off 图 4.6-2 tt （4） [t4,y4,exitflag]=fzero(y,tt(4),[],a,b) %11 （5） [t3,y3,exitflag]=fzero(y,tt(3),[],a,b) （6） op=optimset(fzero) op=optimset(tolx,0.01); op.TolX ans = 0.0100 （7） [t4n,y4n,exitflag]=fzero(y,tt(4),op,a,b) 多元函数的零点 【例 4.6.3-1】求解二元函数方程组的零点。 （1） x=-2:0.5:2;y=x;[X,Y]=meshgrid(x,y); F1=sin(X-Y);F2=cos(X+Y); v=[-0.2, 0, 0.2]; contour(X,Y,F1,v) hold on,contour(X,Y,F2,v),hold off 图4.6-3 （2） [x0,y0]=ginput(2); disp([x0,y0]) -0.7926 -0.7843 0.7926 0.7843 （3） fun=[sin(x(1)-x(2)),cos(x(1)+x(2))]; %12 [xy,f,exit]=fsolve(fun,[x0(2),y0(2)]) %13 Optimization terminated successfully