matlab及多元函数-李志强.pptVIP

1 多元函数的偏导数与全微分 1、符号求导指令diff 多元函数的偏导数与全微分可以通过Matlab中的diff( )指令直接求解，现以二元函数 为例，现要求出 ，则可以用下面的函数求出 f=diff(diff(f, x, m), y, n) 或 f=diff(diff(f, y, n), x, m) 实际上，当m, n取不同值时，该指令可以完成以下功能：  2 空间图形的绘制 画动态图像 函数comet3专门用于x=f(t), y=g(t), z=h(t)所示的曲线的动画。 例:圆锥的顶点为原点，轴线为Z轴，设圆锥的顶角为 ，旋转角度为 ，直线速度为 。动点的起始位置在原点，可得圆锥螺线的方程为 M文件： t=0:0.01:60 =pi/60; =2; =3; x= *sin( )*cos( .*t).*t; y= *sin( )*sin( .*t).*t; z= *t.*cos( ); comet3(x,y,z); 练习：绘制圆柱螺线 空间曲线图的绘制 基本命令： 1 plot3——根据向量x,y,z绘制空间曲线图 例1：绘制参数方程组 的空间曲线图。 M文件： t=1:0.01:10*pi; x=sin(t); y=cos(t); z=t; plot3(x,y,z); grid on;%加网格线 2 box on ——显示边框 t=1:0.01:10*pi; x=sin(t); y=cos(t); z=t; plot3(x,y,z,r-o); box on; %加边框 空间网线图的绘制 基本命令：mesh（x,y,z）%绘制由二元函数z=f(x,y) 所确定的网格线图。 例2：绘制函数z=x^2+y^2在定义域 的 空间网线图。 x=linspace(0,4,30); （x=0:0.2:4;） y=linspace(0,3,20); ( y=0:0.15:3;) [X,Y]=meshgrid(x,y); %生成数据点阵 Z=X.^2+Y.^2; mesh(X,Y,Z); title(‘网线图’)； box on;%(加边框) 空间曲面图的绘制 基本命令：surf（x,y,z）%绘制由二元函数 z=f(x,y)所确定的空间曲面图。 例3：x=linspace(0,4,30); （x=0:0.2:4;） y=linspace(0,3,20); ( y=0:0.15:3;) [X,Y]=meshgrid(x,y); %生成数据点阵 Z=X.^2+Y.^2; surf(X,Y,Z); title(‘曲面图’)； box on； 3 重积分的计算 1、矩形区域上的二重积分的数值计算函数dblquad 设二重积分 ， 其中积分区域 是矩形域. 则二重积分函数dblquad的调用格式如下： 1）q = dblquad(fun,a,b,c,d)：输入量fun是被积函数 ，默认精度 ； 2）q = dblquad(fun,a,b,c,d,tol)：同上，但限定精度tol； 3）q = dblquad(fun,a,b,c,d,tol, @quadl)：指定用quadl代替quad计算二重积分. 4）q = dblquad(fun,a,b,c,d,tol, @myquadf)：指定使用用户自己编写得求二重积分的Matlab程序myquadf.m代替quad计算二重积分. 2、长方体上三重积分数值计算函数triplequad 设三重积分 ，其中积分区域 是长方体域，则三重 积分函数triplequad的调用格式如下： 1）q = tripleq