




已阅读5页,还剩8页未读, 继续免费阅读
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
一、函数的极限 二、函数的导数 三、函数的极值 四、函数的积分 3/13 MatlabMatlab实验实验(五)(五) 一. 函数极限的实现 格式:limit(F,x,a) 计算当xa时,F(x)的极限值, limit(F,x,a,right) 计算当xa+时, F的右极限, limit(F,x,a,left) 计算当xa-时, F的左极限, 特别地,当a=0时有: 解: syms x %定义变量 limit(1-cos(x)*x(-2) 注意:求极限时,先要定义自变量,然后直接将函 数放入limit的括号内,不用引号. ans =1/2 省略了自变量 的变化过程 4/13 MatlabMatlab实验实验(五)(五) 1.一元函数的导数:计算y = f(x) 导数的命令为 :diff(y) 例2.计算下列函数的导数 y=sym(1+x)*log(1+x+sqrt(2*x+x2)-sqrt(2*x+x2); dy=diff(y); b=simplify(dy); 解:syms x 结果为: 二. 函数导数的实现 5/13 MatlabMatlab实验实验(五)(五) 例2.计算下列函数的导数 y=sym(asin(x)/sqrt(1-x2)+0.5*log(1-x)/(1+x); dy=diff(y);b=simplify(dy); 解:syms x 高阶导数可直接计算:diff(S,v,n) 求S对v的n阶导数 6/13 MatlabMatlab实验实验(五)(五) 2. 偏导数的计算 计算 z=f(x,y) 的偏导数的方法为: 首先定义自变量: syms x y; 然后建立函数:z=sym(f(x,y) 用diff求导:dzdx=diff(z,x) ,dzdy=diff(z,y) 例3. 求 的一阶偏导数 解:syms x y; z=sym(exp(x/y); dzdx=diff(z,x) ,dzdy=diff(z,y) 7/13 MatlabMatlab实验实验(五)(五) 三. 求函数的极大值与极小值 在Matlab中有求函数极小值的命令: 计算F在a, b之间取极小值时的x与y(即fval). 命令:x,fval = fminbnd(F,a,b) 解:f=inline(2*x.3-6*x.2-18*x+7) 例4. 求 在区间(-2,4)内极小值 x,fval = fminbnd(f, -2, 4) 故 函数在x=3时,有极小值-47 输出结果为:x = 3.0000 fval = -47.0000 8/13 MatlabMatlab实验实验(五)(五) 注意:如果计算极大值,可将f(x)前面添负号, 则-f(x)的极小值点,即f(x)的极大值点.极大值为-fval 例5. 求 在区间(-2,4)内极大值 解:f=-2*x.3+6*x.2+18*x-7 ; x,fval = fminbnd(f, -2, 4) x = -1.0000 fval = -17.0000 故f(x)在x= -1时有极大值17 注意:计算函数极值时,不能用sym(f(x)表示法 但是可以用:y=f(x) 注意符号! 9/13 MatlabMatlab实验实验(五)(五) 四、不定积分、定积分与广义积分的计算 1.符号函数的积分 格式 : int(s,v,a,b) 其中,s积分表达式; v积分变量; a积分下限,b积分上限 如果求不定积分,无穷积分请大家猜想格式如何? 例6. 计算 解:s=x*exp(-x) g=int(s,x) ans =-x*exp(-x)-exp(-x) 注意:计算结果只给出一个原函数,没有任意常数C 10/13 MatlabMatlab实验实验(五)(五) 2.梯形法数值积分 格式 : I=trapz(x,y) 其中,x是积分区间a,b的取值(向量),y是相 应的函数值 3.辛普森法 格式 : I=quadl(fun,a,b) 注意:quadl最后是字母l, 不是数字1 例7. 计算 方法1:输入 y=(1+sin(x)*exp(x)/(1+cos(x); I1=int(y,x,0,pi/2) 符号运算,不 要点乘除 11/13 MatlabMatlab实验实验(五)(五) 例7. 计算 方法2:输入 x=0:0.01:pi/2; y=(1+sin(x).*exp(x)./(1+cos(x); I2=trapz(x,y) 方法3:输入 x=0:0.01:pi/2; I3=quadl(1+sin(x).*exp(x)./(1+cos(x),0,pi/2) 结果为:I1=exp(1/2*pi) I2=4.8030 I3 =4.8105 12/13 MatlabMatlab实验实验(五)(五) 五 . 函数的级数展开式 格式:taylor(F,a) 功能:F在x=a处的泰勒级数前5项 格式:taylor(F,v) 功能:F对变量v的泰勒展式前5项 格式:taylor(F,v,n) 功能:求F的n 阶泰勒展式,且 (n缺省时默认为 5) 例8. 求 的麦克劳林级数 解:syms x, y=x*exp(-x), taylor(y,9) ans=x-x2+1/2*x3+1/6*x4+1/24*x5+1/120*x6 +1/720*x7+1/5040*x8 13/13 MatlabMatlab实验实验(五)(五) simplify(y),simple(y)化简函数y=f(x) diff(y), diff(z,x)计算y = f(x) 导数, 偏导数 diff(y,n)计算y = f(x) n阶导数 x,fval = fminbnd(F,a,b)曲线F在a,b内极小值点 (x,y) limit(F,x,a) 计算x
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- GB/T 45906.1-2025变电站二次系统第1部分:通用要求
- GB/T 45905.8-2025电力现货市场运营第8部分:多周期电力电量平衡计算方法
- 考评员考试题库及答案
- 森林运动会教学课件
- 森林抚育规程培训课件
- 森林写意课件
- 关于风湿病试题及答案
- (2025)70周岁以上老年人换长久驾照三力测试题库(附答案)
- 2025年外贸营销人员技能水平测试试题集
- 2025年市政工程项目成本控制与管理考试预测试题
- 浙江仙琚制药股份有限公司年产2.5亿粒性激素软胶囊生产线技术改造项目环评报告
- DB37/T 3658-2019地质灾害治理工程施工技术规范
- 个体户餐饮股东协议书
- 挖机台班合同协议书
- 2025年中国智慧养殖行业市场占有率及投资前景预测分析报告
- 电影院安全生产与安全管理规定制度
- 废气处理合同协议
- 2025-2030中国余热回收行业市场现状供需分析及投资评估规划分析研究报告
- 无人机物流配送服务手册
- 二年级上册语文课内阅读理解每日一练(含答案)
- 基层管理培训课程
评论
0/150
提交评论