第7章 函数的导数和积分PPT讲稿.ppt

《第7章 函数的导数和积分PPT讲稿.ppt》由会员分享，可在线阅读，更多相关《第7章 函数的导数和积分PPT讲稿.ppt（42页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、第7章 函数的导数和积分第1页，共42页，编辑于2022年，星期一7.1 7.1 函数的导数函数的导数n7.1.1 7.1.1 函数导数的解析解函数导数的解析解n使用符号求导函数的调用格式是：n diff(fun,x,n)n其中，fun是函数符号表达式；x是符号自变量；n是求导的阶数（n为1时可以省略）。n例例7-1 7-1 计算函数的一阶和二阶导数。n%求函数导数解析解nsyms x%定义表达式中的符号变量nf=sqrt(cos(x)-x*sin(x);%定义函数表达式ndisp(函数的1阶导数：),f1=diff(f)ndisp(函数的2阶导数：),f2=diff(f,2)第2页，共42页

2、，编辑于2022年，星期一nM文件运行结果：n函数的1阶导数：nf1=n-1/2/cos(x)(1/2)*sin(x)-sin(x)-x*cos(x)n函数的2阶导数：nf2=n-1/4/cos(x)(3/2)*sin(x)2-1/2*cos(x)(1/2)-2*cos(x)+x*sin(x)第3页，共42页，编辑于2022年，星期一n7.1.2 7.1.2 二维函数和参数方程的偏导数二维函数和参数方程的偏导数n1、二维函数的偏导数n使用符号求二维函数fun(x,y)偏导数的函数调用格式是：n diff(diff(fun,x),y)n diff(diff(fun,y),x)n其中，fun是函数

3、符号表达式；x和y是符号自变量n例例7-27-2 计算二维函数的2阶偏导数：第4页，共42页，编辑于2022年，星期一n%求二维函数2阶偏导数解析解nsyms x y%定义表达式中的符号变量nf=x2*y/(x+y)3;%定义二维函数表达式ndisp(对变量x的2阶偏导数：)nd2x=diff(f,x,2)ndisp(对变量y的2阶偏导数：)nd2y=diff(f,y,2)ndisp(函数的2阶偏导数：)ndxy=diff(diff(f,x),y)ndisp(函数的2阶偏导数的简化符号表达式：)nsdxy=simplify(dxy)第5页，共42页，编辑于2022年，星期一nM文件运行结果：n

4、对变量x的2阶偏导数：nd2x=n2*y/(x+y)3-12*x*y/(x+y)4+12*x2*y/(x+y)5n对变量y的2阶偏导数：nd2y=n-6*x2/(x+y)4+12*x2*y/(x+y)5n函数的2阶偏导数：ndxy=n2*x/(x+y)3-6*x*y/(x+y)4-3*x2/(x+y)4+12*x2*y/(x+y)5n函数的2阶偏导数的简化符号表达式：nsdxy=n-x*(x2-7*x*y+4*y2)/(x+y)5第6页，共42页，编辑于2022年，星期一第7页，共42页，编辑于2022年，星期一n2、参数方程的偏导数n设参数方程为 和 ，计算参数方程k阶导数 的函数调用格式是

5、：n diff(f,t,k)/diff(g,t,k)n例例7-37-3 计算参数方程导数 解析解，并计算当 时的数值解。n%求参数方程导数的解析解nsyms t%定义参数方程中的符号变量nx=log(cos(t);%定义参数方程1ny=cos(t)-t*sin(t);%定义参数方程2第8页，共42页，编辑于2022年，星期一ndydx1=diff(y,t)/diff(x,t);ndisp(参数方程的1阶导数的简化符号表达式：)nsdydx1=simplify(dydx1)ndydx2=diff(y,t,2)/diff(x,t,2);ndisp(参数方程的2阶导数的简化符号表达式：)nsdydx

6、2=simplify(dydx2)n%计算t=pi/3时的二阶导数值(有效数字6位)ns2=vpa(subs(dydx2,t,pi/3),6);ndisp(t=pi/3时的二阶导数值：),s2nM文件运行后得到的计算结果：n参数方程的1阶导数的简化符号表达式：nsdydx1=n(2*sin(t)+t*cos(t)*cos(t)/sin(t)第9页，共42页，编辑于2022年，星期一n参数方程的2阶导数的简化符号表达式：nsdydx2=n-(-3*cos(t)+t*sin(t)*cos(t)2nt=pi/3时的二阶导数值：ns2=n.148275n因此，参数方程一阶和二阶导数的解析解是：第10页

7、，共42页，编辑于2022年，星期一n7.1.3 n7.1.3 n维函数的偏导数维函数的偏导数n1、n维隐函数的偏导数n设n维隐函数为 ，由于变量之间偏导数的关系是n则计算变量 对变量 的偏导数的函数调用格式是：n -diff(f,xj)/diff(f,xi)第11页，共42页，编辑于2022年，星期一n例例7-47-4 计算二维隐函数n偏导数的解析解。n%求二维隐函数偏导数的解析解nsyms x y%定义函数表达式中的符号变量nf=(x2-2*x)*exp(-x2-y2-x*y);%定义二维隐函数表达式ndydx=-diff(f,x)/diff(f,y);%计算二维隐函数的偏导数dy/dxn

8、disp(二维隐函数偏导数的简化符号表达式)nsdydx=simplify(dydx)第12页，共42页，编辑于2022年，星期一nM文件运行后得到的计算结果：n二维隐函数偏导数的简化符号表达式：nsdydx=n-(-2*x+2+2*x3+x2*y-4*x2-2*x*y)/x/(x-2)/(2*y+x)n因此，二维隐函数偏导数的解析解是：第13页，共42页，编辑于2022年，星期一n2、jacobion矩阵n设有n个自变量 的m个函数n相应的偏导数构成jacobion矩阵为第14页，共42页，编辑于2022年，星期一n运用MATLAB符号工具箱函数可以直接求出jacobion矩阵n jacob

9、ion(f,x)n其中，f是n维函数，x是n维向量。n该函数用于计算n维函数f对n维向量x的jacobion矩阵，其第i行、第j列的值为。当f为数量时，所得的值为f的梯度。当v为数量时，该函数相当于diff(f,x)。n例例7-57-5 已知三维函数 试计算它的偏导数。第15页，共42页，编辑于2022年，星期一n%计算多维函数的导数nsyms x y znf1=x*y*z;f2=x2-y;f3=x+z;ndisp(多维函数表达式向量：)nf=f1;f2;f3ndisp(多维函数的偏导数矩阵：)nJ=jacobian(f,x y z)nM文件运行结果：n多维函数表达式向量：nf=x*y*zn

10、x2-yn x+z第16页，共42页，编辑于2022年，星期一n多维函数的偏导数矩阵：nJ=n y*z,x*z,x*yn 2*x,-1,0n 1,0,1n因此，三维函数的偏导数矩阵是：第17页，共42页，编辑于2022年，星期一n7.1.4 7.1.4 数值微分数值微分n根据函数在一些离散点的函数值，推算它在某点的导数或某高阶导数的近似值。通常用差商代替微商，或用一能近似代替该函数的较简单的函数（如多项式、样条函数）的相应导数作为所求导数的近似值。针对实验数据X（向量或多维数组）数值微分，由于函数表达式是未知的，需要采用数值解法（中心差分方法的微分算法）。n例例7-67-6 已知一组实验数据如

11、表7-1所示，试计算它的4阶微分。表7-1 实验数据x0.00000.20000.40000.60000.80001.0000y0.39270.56720.69820.79410.86140.9053第18页，共42页，编辑于2022年，星期一n%用中心差分方法的数值微分算法nfunction dy,dx=diff_ctr(y,dt,n)nyx1=y 0 0 0 0 0;yx2=0 y 0 0 0 0;yx3=0 0 y 0 0 0;nyx4=0 0 0 y 0 0;yx5=0 0 0 0 y 0;yx6=0 0 0 0 0 y;nswitch nn case 1n dy=(-diff(yx1

12、)+7*diff(yx2)n+7*diff(yx3)-diff(yx4)/(12*dt);l0=3;n case 2n dy=(-diff(yx1)+15*diff(yx2)-15*diff(yx3)+diff(yx4)/(12*dt2);l0=3;第19页，共42页，编辑于2022年，星期一n case 3n dy=(-diff(yx1)+7*diff(yx2)-6*diff(yx3)-6*diff(yx4)+7*diff(yx5)-diff(yx6)/(8*dt3);l0=5;n case 4n dy=(-diff(yx1)+11*diff(yx2)-28*diff(yx3)+28*dif

13、f(yx4)-11*diff(yx5)+diff(yx6)/(6*dt3);l0=5;nEndndy=dy(l0+1:end-l0);ndx=(1:length(dy)+l0-2-(n2)*dt;第20页，共42页，编辑于2022年，星期一n%用中心差分方法的数值微分算法调用M文件nsx=0.0000 0.2000 0.4000 0.6000 0.8000 1.0000;nsy=0.3927 0.5672 0.6982 0.7941 0.8614 0.9053;ny=sy,sx;ndt=0.2;nn=4;ndy,dx=diff_ctr(y,dt,n)nM文件运行结果：ndy=19.2187-2

14、17.2292 581.9417-611.7500 250.0458 -23.0271ndx=0.6000 0.8000 1.0000 1.2000 1.4000 1.6000第21页，共42页，编辑于2022年，星期一n7.1.5 7.1.5 函数的梯度和梯度的模函数的梯度和梯度的模n在向量微积分中，标量场的梯度是一个向量场，标量场中某一点上的梯度指向标量场增长最快的方向。梯度是一个向量，当某一函数在某点处沿着该方向的方向导数取得该点出的最大值，即函数在该点处沿方向变化最快，变化率最大，称为为该梯度的模。nn维函数 ，在点 的梯度，是函数各维一阶偏导数组成的向量第22页，共42页，编辑于20

15、22年，星期一n梯度的模是函数各维一阶偏导数平方和的开方n1、计算梯度的符号函数 jacobian(fun,X)n其中，fun是符号函数表达式，X是定义求导变量组成的向量（可以省略）。n2、计算矩阵或向量范数的符号函数n norm(fun,X)n其中，fun是符号型函数表达式，X是定义求导变量组成的符号型向量。第23页，共42页，编辑于2022年，星期一n例例7-77-7 计算函数在 梯度和梯度模:n%计算多维函数的梯度和梯度的模nsyms x1 x2 x3%定义符号变量nf=2*x12+5*x22+x32+2*x2*x3+2*x1*x3-6*x2+3ndisp(三维函数的梯度：)ngf=ja

16、cobian(f)nxk=1,-1,2;ndisp(函数在xk点的梯度值：)ngfk=subs(subs(subs(gf,xk(1),xk(2),xk(3)ndisp(函数在xk点的梯度的模：)ngmk=norm(gfk)第24页，共42页，编辑于2022年，星期一nM文件运行结果：n三维函数的梯度：ngf=n 4*x1+2*x3,10*x2+2*x3-6,2*x3+2*x2+2*x1n函数在xk点的梯度值：ngfk=8 -12 4n函数在xk点的梯度的模：ngmk=14.9666n因此，三维函数的梯度，以及在的梯度和梯度的模是：第25页，共42页，编辑于2022年，星期一第26页，共42页，

17、编辑于2022年，星期一7.2 7.2 函数的积分函数的积分n设 是函数 的一个原函数，则将函数 的所有原函数 （C为任意常数）叫做函数的不定积分，记作 。n定积分 就是求函数 在区间a,b 中函数图线下包围的面积。第27页，共42页，编辑于2022年，星期一n7.2.1 7.2.1 不定积分的解析解不定积分的解析解n计算不定积分解析解的函数调用格式是：n F=int(fun,x)n其中，fun是被积函数；x是自变量，如果只有一个变量，x可以省略。n应当指出，计算结果F(x)只是积分原函数，实际的不定积分应该是F(x)+C构成的曲线族，C是任意常数。n例例7-87-8 计算不定积分n%计算不定

18、积分nsyms x a b;%定义表达式中的符号变量ny=sin(a*x)*cos(b*x);%定义被积函数第28页，共42页，编辑于2022年，星期一ndisp(被积函数的不定积分：)nF=int(y,x)nM文件运行结果：n被积函数的不定积分：nF=n-1/2/(a+b)*cos(a+b)*x)-1/2/(a-b)*cos(a-b)*x)n因此，不定积分 第29页，共42页，编辑于2022年，星期一n7.2.2 7.2.2 定积分定积分n计算定积分的函数调用格式是：n I=int(fun,x,a,b)n其中，fun是被积函数；x是自变量；(a,b)是定积分的积分区间，计算无穷积分时，需要将

19、积分区间设置为(-inf,inf)。n例例7-97-9 计算定积分 和 的解析解和数值解。n%计算定积分nsyms x;%定义表达式中的符号变量ny1=x*exp(x)/(1+x)2;%定义被积函数1第30页，共42页，编辑于2022年，星期一na1=0;b1=1;%定义积分区间1ndisp(被积函数1定积分的解析解：)nF1=int(y1,x,a1,b1)%计算y1的定积分ndisp(被积函数1定积分的数值解：)nF1v=double(F1)ny2=cos(x)/sqrt(x);%定义被积函数2na2=0;b2=inf;%定义积分区间2ndisp(被积函数2定积分的解析解：)nF2=int(

20、y2,x,a2,b2)ndisp(被积函数2定积分的数值解：)nF2v=double(F2)nM文件运行结果：第31页，共42页，编辑于2022年，星期一n 被积函数1定积分的解析解：nF1=1/2*exp(1)-1n 被积函数1定积分的数值解：nF1v=.3591n 被积函数2定积分的解析解：nF2=1/2*2(1/2)*pi(1/2)n 被积函数2定积分的数值解：nF2v=1.2533n因此，两个定积分的解析解和数值解是：第32页，共42页，编辑于2022年，星期一n7.2.3 7.2.3 数值积分数值积分n求某函数的定积分时，在多数情况下，被积函数的原函数很难用初等函数表达出来；另外，许

21、多实际问题中的被积函数往往是列表函数或其他形式的非连续函数。对这类函数只能使用数值积分计算函数的近似值。n1、梯形法求解数值积分n梯形法求解数值积分的：将积分区间分割成许多足够小的分区间，取各个矩形的面积分别代替积分函数的曲线在各小段区间上围出来的曲边梯形的面积，再将所有这样的矩形面积加起来的总和，近似地等于函数在这个区间上的定积分。n梯形法求解数值积分的函数调用格式是：n s=trapz(x,y)n其中，x是向量；y的行数应该等于x向量的元素数 第33页，共42页，编辑于2022年，星期一n例例7-10 7-10 用梯形法计算表7-2数据的定积分。n%梯形法计算定积分nx=0:0.1:1.2

22、;ny=0 2.2077 3.2058 3.4435 3.241 2.8164 2.311 1.8101 1.3602 0.9817 0.6791 0.4473 0.2768;ns=trapz(x,y)nM文件运行结果：ns=n 2.2642第34页，共42页，编辑于2022年，星期一n2、变步长辛普生法求解数值积分n变步长辛普生法求解数值积分的基本思想是：在梯形法的基础上，将积分区间逐次分半，计算出每个子区间的定积分近似值，并且求和。n变步长辛普生法求解数值积分的函数调用格式是n s,n=quadl(fun,a,b,t0l)n其中，s是返回的定积分值；n n是返回被积函数调用次数；n fun

23、是被积函数；n a和b分别是定积分的下限和上限;n tol是控制积分精度（默认 ）。n该函数适用于求解给定函数的数值积分。第35页，共42页，编辑于2022年，星期一n例例7-11 7-11 用变步长辛普生法计算定积分：n%变步长辛普生法计算定积分nf=inline(exp(-0.5.*x).*sin(x+pi/6),x);%定义被积函数和自变量na=0;b=3*pi;%定义积分区间nformat long;%数值按长格式显示ns,n=quadl(f,a,b,1e-15)%计算y的定积分nM文件运行结果：ns=0.90084078781889nn=828第36页，共42页，编辑于2022年，星

24、期一n%绘制积分图形nfplot(exp(-0.5.*x).*sin(x+pi/6),a,b);%符号函数绘图nhold on;nx=a:pi/100:b;ny=exp(-0.5.*x).*sin(x+pi/6);nfill(a,x,b,0,y,0,y);%填充积分区域nxlabel(bfit x);ylabel(bfit y);ntitle(bf 变步长辛普生法计算定积分);ngtext(y=e-0.5xsin(x+pi/6);n输出图形（如图7-2所示）第37页，共42页，编辑于2022年，星期一第38页，共42页，编辑于2022年，星期一n7.2.4 7.2.4 函数的重积分函数的重积分

25、n1 1、二维函数的数值积分、二维函数的数值积分（双重积分）n对于二维函数的双重定积分问题n可以直接使用函数进行求解矩形区域的数值积分，其调用格式是：n y=dblquad(fun,xm,xM,ym,yM,tol)n其中，输入参数fun是二维被积函数，可以用函数文件或函数inline()定义；n (xm,xM,ym,yM)是指定的矩形区域；n tol是指定的计算精度（默认精度是 ）。第39页，共42页，编辑于2022年，星期一n例例7-12 7-12 试计算二维被积函数在矩形区域内的双重定积分：n%二维函数矩形区域的双重定积分nfun=inline(exp(-x.2/2).*sin(x.2+y

26、),x,y);nepsilon=1e-10;%计算精度nI=dblquad(fun,-2,2,-1,1,epsilon);%计算双重积分nformat long ndisp(二维函数矩形区域的双重定积分：),InM文件运行结果：n二维函数矩形区域的双重定积分：nI=1.57449815921952第40页，共42页，编辑于2022年，星期一n2 2、三维函数的数值积分、三维函数的数值积分（三重积分）n对于三维函数的三重定积分问题n可以直接使用函数进行求解长方体区域的数值积分，其调用格式是：n y=triplequad(fun,xm,xM,ym,yM,zm,zM,tol)n其中，输入参数fun是

27、三维被积函数，可以用函数文件或函数inline()定义；n(xm,xM,ym,yM,zm,zM)是指定的长方体区域；ntol是指定的计算精度（默认精度是 ）。第41页，共42页，编辑于2022年，星期一n例例7-13 7-13 试计算三维被积函数在长方体区域内的双重定积分：nfun=inline(y.*sin(x)+z.*cos(x),x,y,z);%定义被积函数和变量nepsilon=1e-10;%计算精度nI=triplequad(fun,0,pi,0,1,-1,1,epsilon);%计算双重积分，积分区间 nformat long ndisp(三维函数矩形区域的三重定积分：),InM文件运行结果：n三维函数矩形区域的三重定积分：nI=1.99999999999975第42页，共42页，编辑于2022年，星期一