matlab在微积分中的应用.ppt

2.求函数的极值点 (1)求一元函数的极值问题 命令形式1：fmin(fun,x1,x2) 功能:在区间[x1，x2]内求函数fun的极小值点 命令形式2：fminbnd(fun,x1,x2) 功能:在区间[x1，x2]内求函数fun的极小值点 两个函数功能相同,命令1是matlab早期版本中使用的. 注意：fun可以是字符串，内联函数或M函数文件名。 (2)求多元函数极值问题 命令形式1：fminsearch(fun,x0) 功能:用单纯形法求多元函数fun在x0附近的极值点。 命令形式2：fminunc(fun,x0) 功能:用拟牛顿法求多元函数fun在x0附近的极值点。 例27.求函数 的极小值点。 例28.求函数 在点(0,5,4)附近的极小值。 function f=myfun(p) x=p(1); y=p(2); f=100*(y-x^2)^2+(1-x)^2; x0=[-1.2,1] x=fminunc(myfun,x0) function f=myfun(p) x=p(1); y=p(2); z=p(3); f=x^4+sin(y)-cos(z); x0=[0,5,4] [xmin,fval]=fminsearch(myfun,x0) 2.（1） 1.已知 ，求 。 （2） * 湖北师范学院数学与统计学院数学建模实验电子教案 * * 微积分的基础知识及其在Matlab中的实现 明巍 数学与统计学院 数学建模种常用的微积分知识在Matlab中的实现 1. 极限运算 2. 求导运算 3. 积分运算 4. 函数的Taylor展开 5. 数值积分 6. 线性方程和非线性方程的求解 7. 求和及求极值方法 计算单侧极限 limit(f,x,a,’left’) 计算单侧极限 limit(f,x,a,’right’) 计算 limit(f,x,inf) 计算 limit(f,x,a) 功能 命令 注意：在左右极限不相等或左右极限有一个不存在时，Matlab的默认状态是求右极限。 例1.求极限 与极限 解： syms x; y1=(1+4*x)^(1/x); y2=(exp(x)-1)/x; a=limit(y1,x,0) b=limit(y2,x,0) a=exp(4) b=1 例2.求极限 解：syms x; y=sqrt(x)-2^(-1/x); limit(y,x,0,right) ans=0； 1．一元函数求导 命令形式1：diff（f） 功能：求函数f的一阶倒数，其中f为符号函数。 命令形式2：diff（f,n） 功能：求函数f的n阶倒数，其中f为符号函数。 例3.求函数 的二阶倒数 例4.设 ，求 解：syms x; f=3*x^3+5*x+1; diff(f,2) ans=18*x 解：syms x; y=3*x^2-2*x+1 B=diff(y),x=1; eval(B) B=6*x-2 ans=4 2.多元函数的偏导数 命令形式1：diff（f,xi） 功能：求多元函数f对变量xi的一阶偏导。 命令形式2：diff（f,xi,n） 功能：求多元函数f对变量xi的n阶偏导。 例5.求 关于x的偏导。 解：syms x y; z=x^2*sin(2*y); B=diff(z,x) B=2*x*sin(2*y) 3.全微分、参数方程求导及隐函数求导 (1)若函数 在某点 的两个偏导数存在且连续，则函数在该点的全微分为 （2）对参数方程 所确定的函数 ， Matlab中求函数全微分的命令为： 根据公式 ， 连续两次利用 命令就可以求出结果。 (3)隐函数求导 方程 所确定的隐函数 ，其导数为 方程 确定的隐函数 ，其导数为 在Matlab中按照上述公式，分别求出函数的偏导数再相除就可以得到隐函数的导数。 1．一元函数的不定