- 1、本文档共37页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
3.4.3函数的数值积分与微分 1.函数的数值积分 (1) 低阶数值积分函数quad 格式一:q=quad(‘fun’,a,b) 功能:采用自适应的Simpson积分方法,返回函数‘fun’ 在上限 a 和下限 b之间的数值积分。当给定一个输入值向量,‘fun’ 必须返回一个输出向量。函数‘fun’可以是函数名、函数句柄或字符串。 格式二:q=quad(‘fud’,a,b,tol) 功能:按指定绝对误差tol返回数值积分值, tol 缺省值为 1e-6。 (2) 高阶数值积分函数quadl ,具有与quad同样的调用方法,在处理软奇异函数时比quad更好。 例3-50:求y=e-x+x2在[0,1]区间上的定积分。 (3) 梯形面积法的积分函数trapz( ) 格式一:T=trapz(Y) 功能:以单位间隔,采用计算若干梯形面积的和来计算某函数的近似积分。如果Y为向量,计算Y的积分;如果Y是矩阵,得一个每列积分的行向量;如果Y为多维数组,则沿第一个非单元素维计算。 格式二:T=trapz(X,Y) 功能:用梯形积分法,依据X计算Y的积分。如果X为矢量,则Y必须是同大小的矢量;如果X是一列向量,并且数组Y第一非单元素维长度为length(X),则在该维中计算。 例3-51: x1=-1:0.17:2;y1=humps(x1); Areal=trapz(x1,y1) x2=-1:0.07:2;y2=humps(x2);Aarea=trapz(x2,y2) (4) 双重积分函数dblquad MATLAB提供了一个求双重积分的函数dblquad,其基本调用格式为: 格式:Q=dblquad(fun,xmin,xmax,ymin,ymax,tol) 功能:按指定精度tol,对指定函数 f(x, y)在[xmin, xmax]范围和[ymin, ymax]范围进行双重积分。精度tol缺省时默认精度为1e-6。 例3-52:在0xπ, πy2π区域内计算ysinx+xcosy的二重积分 (5) 不定积分的计算 对于函数不定积分的计算,可以采用定积分函数来求不定积分的数值解。方法是: 固定积分下限,用for循环,把积分上限逐步增加即可实现。 例3-53:试求humps函数以x=0为下限的不定积分。 2.函数的数值微分 (1)通过计算数组中元素间的差分函数diff来粗略计算微分函数。 例3-54:试利用diff函数计算y=f(x)的微分,已知该函数的数据如下: X: 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 Y: -0.447 1.978 3.28 6.16 7.08 7.34 7.66 9.56 9.48 9.30 11.2 (2)对数据拟合后再利用polyder函数微分 为了实现对数值较准确的微分,最好先用最小二乘曲线拟合这种数据,然后对所 得到的多项式进行微分。 例3-54:试利用polyder函数计算y=f(x)的微分,已知该函数的数据如下: X: 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 Y: -0.447 1.978 3.28 6.16 7.08 7.34 7.66 9.56 9.48 9.30 11.2 3.4.4 常微分方程的数值求解 1.初始值的常微分方程求解 格式:[T,Y] = ode45(‘odefun’,tspan,y0) 功能:返回由文件‘odefun’所定义的具有初始条件为y0、时间t变化范围为[t0, tfinal] 的微分方程y’= f(t,y)的解,其中tspan = [t0,tfinal]。向量T中的每一列对应着矩阵Y的每一列。 例3-55:用数值积分法求解微分方程:y’’=-y+1-t2/2π。设初始时间t0=0,终止时间tf=3π,初始条件y(0)=0,y’(0)=0。 任何高阶常微分方程都可以变换成一阶微分方程,即表示成右函数形式,这是利用龙格-库塔法求解微分方程的前提。 2.边界条件的常微分方程求解 对于如下的微分方程: 用于边界条件的常微分方程求解问题: 函数bvp4c 格式:sol = bvp4c(‘odefun’,bcfun,solinit) 其中‘odefun’为常微分方程函数,bcfun为边界条件函数,solinit为求解的初始值。输出sol是一个结构,它
文档评论(0)