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

函数的数值导数和切平面 函数的数值导数和切平面 数学实验第五次课 微积分实验 4.1 函数的数值导数和切平面 4.1.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]) 图 4.5.1-1 4.1.2 偏导数和梯度 4.1.2.1 理论定义 4.1.2.2 数值计算指令 【例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）。设在(a , b ) 处有电荷+q ，在(-a , -b ) 处有电荷-q 。那么在电荷所在平面上任何一点的电势和场强分别为 ，E =-?V 。其中