极限导数微分与符号计算.ppt

数学实验补充补充 极限极限/导数导数/积分积分/FUNTOOL/MATSERV 内容：本讲主要针对高等数学补充极限、导数、内容：本讲主要针对高等数学补充极限、导数、 积分相关运算积分相关运算;补充重要的符号计算器；补充重要的符号计算器； 补充基于补充基于MATLAB计算服务器搭建知识计算服务器搭建知识目的：学习极限目的：学习极限/导数导数/积分相关函数的指令实积分相关函数的指令实现现 及工具实现；了解及工具实现；了解MATSERV的搭建过的搭建过程程 要求：能够解决高等数学中的极限要求：能够解决高等数学中的极限/导数导数/积分积分问问 题；能够动手搭建题；能够动手搭建MATSERV计算服务器计算服务器掌握极限掌握极限(左、右极限左、右极限)计算函数计算函数 limit掌握导数掌握导数(1阶导、高阶导、偏导阶导、高阶导、偏导)计算函数计算函数 diff掌握积分掌握积分(不定积分、定积分、数值积分不定积分、定积分、数值积分)计算计算函数函数 int trapz quad quad8了解了解MATLAB计算服务器搭建知识计算服务器搭建知识补充补充 引入与导言引入与导言01.数学实验关于专题引入和概述~ 极限极限,导数导数,积分是我们在高等数学学习中积分是我们在高等数学学习中接触过的最基本也是最重要的概念接触过的最基本也是最重要的概念.一方面它一方面它们是很多数学工具的基础们是很多数学工具的基础(比如微分方程比如微分方程);另另一方面它们又是工程计算和科学研究直接面一方面它们又是工程计算和科学研究直接面对的问题对的问题. 微分微分(导数导数)运算比较简单运算比较简单,任何一个由基本任何一个由基本初等函数经过四则及复合运算构成的函数初等函数经过四则及复合运算构成的函数,都都可以用导数公式和求导法则算出它们的导数可以用导数公式和求导法则算出它们的导数. 积分运算则相对复杂得多积分运算则相对复杂得多,仍有许多函数仍有许多函数“积不出来积不出来”,由于它们的原函数无法由基本初由于它们的原函数无法由基本初等函数经过四则及复合运算构成等函数经过四则及复合运算构成,计算这类定计算这类定积分问题我们也只能采用数值方法积分问题我们也只能采用数值方法. 借助借助 MATLAB 我们得以快速解决这些问题我们得以快速解决这些问题!补充补充引入与导言引入与导言02.数学实验基本调用格式：基本调用格式：limit(f)功能功能:计算计算limit(f,x,a)功能功能:计算计算limit(f,x,inf)功能功能:计算计算limit(f,x,a,'right')功能功能:计算计算limit(f,x,a,'left')功能功能:计算计算求极限运算1补充补充求极限运算求极限运算03注意：注意：默认默认x趋于趋于0；；在左，右极限在左，右极限不相等，或有不相等，或有一个不存在时，一个不存在时，默认为求右极默认为求右极限；限；.数学实验求极限运算2应用示例〔熟悉应用类型）：应用示例〔熟悉应用类型）：1、求极限、求极限syms x;y=((1+tan(x))/(1+sin(x)))^(1/x^3);limit(y)2、求极限、求极限syms n;y=(1+1/n)^n;limit(y,n,inf) 3、求极限、求极限补充补充求极限运算求极限运算04syms x;y=5*x+log(sin(x)+exp(sin(x)));limit(y,x,3,'left').数学实验求导数运算1补充补充求导数运算求导数运算[1]一元函数求导一元函数求导基本调用格式：基本调用格式：diff(f) 功能功能-求函数求函数f的一阶导数的一阶导数diff(f,n) 功能功能-求函数求函数f的的n阶导数阶导数应用示例：应用示例：4，求，求的一阶、二阶导数的一阶、二阶导数syms a b x;y=(a*x+tan(3*x))^(1/2)+sin(x)*cos(b*x);d1y=diff(y), disp('***'),pretty(d1y), disp('***')d2y=diff(y,2), disp('***') ,pretty(d2y), 05.数学实验求导数运算2[2]多项式拟合求导多项式拟合求导方法说明：方法说明：先利用先利用polyfit将函数拟合成多项式函数，然后利将函数拟合成多项式函数，然后利用多项式函数求导命令用多项式函数求导命令polyder求导求导应用示例：应用示例：4，用，用4阶多项式拟合函数阶多项式拟合函数并求并求x=2处的一阶与二阶导数处的一阶与二阶导数x=0:.1:8;y=cos(x).*log(3+x.^2+exp(x.^2));p=polyfit(x,y,4);plot(x,y,'b',x,polyval(p,x),'r');poly2str(p,'x'),p1=polyder(p);p2= polyder(p1);y=polyval(p,2),y1=polyval(p1,2), y2=polyval(p2,2)补充补充求导数运算求导数运算06.数学实验求导数运算3[3]参数方程求导参数方程求导方法说明：方法说明：对参数方程对参数方程x=x(t);y=y(t);先求出先求出dy/dt和和dx/dt然后代入公式然后代入公式dy/dx= dy/dt / dx/dt 即可即可应用示例：应用示例：5，求参数方程，求参数方程syms t;x=t*(1-sin(t));y=t*cos(t);dx=diff(x,t);dy=diff(y,t);pretty(dy/dx)补充补充求导数运算求导数运算07.数学实验求导数运算4补充补充求导数运算求导数运算08[4]多元函数求导多元函数求导方法说明：方法说明：对指定变量求导，此处涉及偏导对指定变量求导，此处涉及偏导…应用示例：应用示例：6，求，求 对对z 的偏导数的偏导数syms a b x y z;u=a*exp(b*x+y+z^2);pretty(diff(u,z))7，对，对syms x y;z=x^3*y^2+sin(x*y);diff(z,x,3).数学实验求导数运算5补充补充求导数运算求导数运算098，以，以 为例验证罗必塔法则：为例验证罗必塔法则：syms a b xf=a^x-b^x;g=x;l1=limit(f/g,x,0)df=diff(f,x);dg=diff(g,x);l2=limit(df/dg,x,0)l1==l2.数学实验求积分运算1补充补充求积分运算求积分运算[1]不定积分不定积分方法说明：方法说明：int(f)对默认变量积分；对默认变量积分；int(f,v)对指定变量积分对指定变量积分应用示例：应用示例：8，计算，计算syms x;y=1/(sin(x)^2*cos(x)^2);pretty(int(y))9，计算，计算syms a x;y=1/(a^2-x^2);pretty(int(y,x))10.数学实验求积分运算2补充补充求积分运算求积分运算11[2] 定积分定积分-符号解法符号解法方法说明：方法说明：int(f,x,a,b) 依据微积分基本公式计算依据微积分基本公式计算应用示例：应用示例：10，计算，计算syms a x;f=sqrt(x^2+a);pretty(int(f,x,-2,2))11，对变上限函数，对变上限函数 求导求导 syms t x;f= sqrt(1-t^2);pretty(diff(int(f,t,0,x^2))).数学实验求积分运算3补充补充求积分运算求积分运算12[3] 定积分定积分-数值解法数值解法方法说明：方法说明：当定积分当定积分-符号解法失效时，必须用定积分符号解法失效时，必须用定积分-数数值解法来近似计算定积分的值。

矩形公式值解法来近似计算定积分的值矩形公式sum，，复合梯形公式复合梯形公式trapz，复合辛普森公式，复合辛普森公式quad/quad8的区别在于替代等距曲边梯形的方的区别在于替代等距曲边梯形的方式不同：式不同：.数学实验求积分运算4补充补充求积分运算求积分运算13应用示例：应用示例：sum使用一次用于求向量或矩阵每一列的和，使用一次用于求向量或矩阵每一列的和，若使用两次则先按列求和再按行求和若使用两次则先按列求和再按行求和(行列总和行列总和)12，矩形法计算，矩形法计算 在在x=0与与x=10之之间所围面积间所围面积dx=0.1;x=0:dx:10;y=-x.^2+115;sum(y(1:length(x)-1))*dx( 的近似值的近似值).数学实验求积分运算5补充补充求积分运算求积分运算14trapz(x,y)用复合梯形公式计算定积分，用复合梯形公式计算定积分，x为积分变量分为积分变量分点向量，点向量，y为被积函数分点函数值向量为被积函数分点函数值向量quad('fun',a,b,tol,trace)用复合辛普森公式计算定积分，用复合辛普森公式计算定积分，fun为被积函数为被积函数表达式字符串或表达式字符串或m函数文件名，函数文件名，a，，b是积分下是积分下上限，上限，tol表示精度表示精度(缺省缺省0.001)，，trace=1图示积图示积分过程分过程(默认默认=0不显示不显示) %quad8与与quad类似，精度更高类似，精度更高.数学实验求积分运算6补充补充求积分运算求积分运算1513，用两种方法求，用两种方法求x=2:.1:5;y=log(x)./(x.^2);tt=trapz(x,y) %复合梯形公式复合梯形公式 fun=inline(' log(x)./(x.^2) ','x');ss=quad(fun,2,5) %复合辛普森公式复合辛普森公式.数学实验符号计算器的使用1补充补充求积分运算求积分运算16.数学实验符号计算器的使用2补充补充求积分运算求积分运算17图形化符号函数计算器的使用：f=为图形窗口1的控制函数，其缺省值为x；g=为图形窗口2的控制函数，其缺省值为1；x=为两窗口函数的自变量取值范围，缺省a=为常数，缺省值为1/2。

df/dx 计算函数f对x的导法式，并赋给fint f 计算函数f的积分函数，并赋给fsimple f 计算函数f的最简表达式，并赋给fsyms x) simplify(cos(x)^2+sin(x)^2); simplify((x^2+5*x+6)/(x+2)); expand(cos(x+y)); expand((x-2)*(x-4)); syms x y; factor(x^3-y^3); factor(x^3+3*x^2+3*x+1); num f 取表达式f的分子，并赋给fden f 取表达式f的分母，并赋给f1/f 求f的倒数函数，并赋给ffinv 求f的反函数，并赋给f数学实验符号计算器的使用3补充补充求积分运算求积分运算18f±a 计算f(x) ±a，并赋给ff*a 计算f(x)*a，并赋给ff/a 计算f(x)/a，并赋给ff^a 计算 ，并赋给ff(x+a) 计算f(x+a)，并赋给ff(x*a) 计算f(ax) ，并赋给ff+a 计算f(x)+a，并赋给ff±g 计算两函数之和/差，并赋给f。

f*g 计算两函数之积，并赋给ff/g 计算两函数之比，并赋给ff(g) 计算复合函数f(g(x)) g=f 将f的函数值赋给gswap 交换f与g的函数表达式数学实验MATLAB WEB SERVER19应用前景：应用前景：通过架构通过架构 Matlab Web Server，为公司或个人提，为公司或个人提供供 Matlab 计算服务〔商业计算服务〔商业/免费），这是免费），这是Matlab 很有前景的一个应用方向很有前景的一个应用方向实现这项技术，我们需要一个实现这项技术，我们需要一个WWW服务程序，服务程序，比如比如IIS或者或者APACHE，推荐，推荐APACHE，因为它，因为它是使用最广泛的服务器平台，且是开源的程序，是使用最广泛的服务器平台，且是开源的程序，下面我们以下面我们以APACHE为例详细讲解为例详细讲解Matlab 计算服务的实现，另外我们还需要计算服务的实现，另外我们还需要Matlab 的内的内联服务器联服务器 WEBSERVER 沟通沟通APACHE 和和 Matlab 计算内核计算内核补充补充搭建计算服务器搭建计算服务器.数学实验MATLAB WEB SERVER201.安装并配置安装并配置APACHE服务器程序服务器程序:1.1 下载最新的下载最新的Apache win32安装程序安装程序1.2 配置域名配置域名/主机信息〔随意输入）主机信息〔随意输入）1.3 更改安装目录为更改安装目录为C:\，测试页面，测试页面index. html.en1.4 设置设置CGI目录执行权限：目录执行权限：在在C:\Apache\conf\httpd.conf 中〔中〔618行）行）修改修改 Options None 为为 Options ExecCGI1.5 设置演示图片保存目录设置演示图片保存目录icons：： 在在C:\Apache\conf\httpd.conf 中〔中〔580行）行）用用#号注释掉号注释掉580行到行到587行，在行，在588行新增行新增Alias /icons/ "C:/Apache/htdocs/"1.6 修改完毕重启修改完毕重启Apache〔通过开始菜单〔通过开始菜单restart））补充补充搭建计算服务器搭建计算服务器.数学实验MATLAB WEB SERVER212. 配置配置MATLAB WEB SERVER服务器服务器(需启动服务需启动服务!):2.1 在在C:\MATLAB6p5\webserver\matlabserver.conf 中中修改修改 -m 1 为为 -m 5 （提高同时运行的（提高同时运行的Matlab进程数）进程数）2.2 在在C:\MATLAB6p5\toolbox\webserver\wsdemos\下下复制复制 matweb.conf 和和 matweb.exe 到到 C:\Apache\cgi-bin\复制复制 wsdemos\ 目录下其他文件到目录下其他文件到 C:\Apache\htdocs\2.3 在在C:\Apache\cgi-bin\matweb.conf中中修改所有修改所有mlserver= 为为mlserver=localhost修改所有修改所有mldir=/toolbox/webserver/wsdemos 为为mldir= C:/Apache/htdocs3.修改完毕重启修改完毕重启Apache，验证，验证localhost/index.html补充补充搭建计算服务器搭建计算服务器.数学实验That’s all~3Q!下次课内容下次课内容微分方程数值解微分方程数值解22.。