MATLAB第4讲-MATLAB 数值计算.ppt

  1. 1、本文档共58页,可阅读全部内容。
  2. 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
MATLAB第4讲-MATLAB 数值计算

数值线性代数问题 4.3 函数和数值积分 一、函数的绘图及分析 1、绘制函数曲线 格式: fplot(‘函数名’,lims,’s’) 功能: 绘制指定函数的曲线, lims为x,y 轴的最小最大值, s可指定线形 ■函数和数值积分库(funfun) ■特殊函数库(specfun) 例题 4.3 函数和数值积分 2、函数的简易绘图 格式:ezplot(‘函数名’,lims) 例如:ezplot(‘sin’) ezplot(x^2 - x - 1) ; ezplot(sin(t),cos(t)); (画一椭圆) 除 ezplot外,还有: ezpolar, ezplot3, ezmesh , ezsurf等。 例题 4.3 函数和数值积分 3、求函数极小值 格式:fmin(‘函数名’,x1,x2) 求函数在x=[x1 x2]之间的极小值 4、求函数零点 格式:fzero(‘函数名’,x0)。 求函数在x0附近的过零点 。 例如:fzero(‘humps’,1)=1.2995 4.3 函数和数值积分 二、特殊函数 特殊的数学函数, 如:贝塞尔函数、误差函数等。 4.3 函数和数值积分 三、函数的数值积分 1、定积分(一维数值积分) 格式:quad(‘函数名’,x1,x2) (或 quad8 ,高阶方法) 对函数在区间[x1,x2]内的定积分 例如: quad(‘humps’,0,1)=29.8583 利用定积分可以求不定积分的数值解 数值积分 [y, n]=quad( F, a, b, tol); [y, n]=quad8(F, a, b, tol); [y, n]=quadl(F, a, b, tol); 其中F为被积函数, a, b为积分上下限, tol 为误差限 y 为积分结果,n 为被积分函数的调用次数。 tol 的缺省值对 quad 函数为 1e-3 对 quad8 和 quadl 函数为 1e-6 quad8 和 quadl 函数使用的算法比 quad 高得多 例子:求无穷定积分 f=inline(1/sqrt(2*pi)*exp(-x.^2/2),x); [y,kk]=quad(f,-8,8) y = 1.00000197533430 kk = 81 [y1,kk1]=quad8(f,-8,8) y1 = 1.00000000000023 kk1 = 161 [y,kk]=quad(f,-15,15) y = 0.99999920879563 kk = 89 [y1,kk1]=quad8(f,-15,15) y1 = 0.99999999999999 kk1 = 769 该无穷定积分的理论值为 1 4.3 函数和数值积分 三、函数的数值积分 2、二重数值积分 格式: dblquad(‘函数名’,x1,x2,y1,y2) 例子:双重积分 用函数 dblquad(F, xmin, xmax, ymin, ymax, tol) 求双重积分 f=inline(exp(-x.^2/2).*sin(x.^2+y),x,y); y=dblquad(f,-2,2,-1,1) y = 1.57450259485169 4.3 函数和数值积分 四、常微分方程的求解 (ODE ,Ordinary Differential Equation ) 格式:[x,y]=ode23(‘函数名’,x0,xn,y0) (2,3阶) 或 ode45 (4,5阶) 函数名为微分方程名(含x,y) x0,xn为x区间,y0为y初值。 4.3 函数和数值积分 四、常微分方程的求解 例如:求微分方程: dy/dx=-3y+2x, y(1)=2 区间为[1,3]的解 首先建立函数:f=myde(x,y) (myde.m) f= -3*y+2*x [x,y]=

文档评论(0)

yaocen + 关注
实名认证
内容提供者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档