




版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、从积分和式到求积公式 插值型求积公式 求积公式的代数精度 复合梯形公式 MATLAB求积分命令,第五章 数值积分与数值微分,椭圆周长计算:,x =a cos t y =b sin t,0 t 2,椭圆积分,x = a sin cos y = b sin sin z = c cos ,D= (,) |0 2, 0 ,思考题: 椭球面的面积计算,椭球面积的积分表达式? 对二重积分的计算问题? 三维体积的离散数据计算?,积分和式的计算: f (x)Ca, b单增, 令 h = (b a)/n, xk = a + kh, (k = 0,1,2,n),近似计算,4.8612e+004 4.8660e+0
2、04 4.8683e+004 4.8695e+004 4.8701e+004 4.8704e+004 4.8706e+004 4.8707e+004 4.8707e+004,4.8803e+004 4.8755e+004 4.8731e+004 4.8719e+004 4.8713e+004 4.8710e+004 4.8709e+004 4.8708e+004 4.8708e+004 4.8708e+004,a=7782.5 c=972.5,P.170人造卫星的轨道长度计算,数值求积公式的一般形式(机械求积公式),Rf 为数值求积公式余项, x0, x1, , xn为求积 结点; A0, A
3、1, , An为求积系数.,矩形公式: 取A0 =(b a ),特别地, 时分别称为左矩公式,中矩公式,右矩公式。,梯形公式: 取A0 =A1 =(b a )/2,Simpson 公式 取A0 = A3 =(b a )/6, A1 =2(b a )/3,,插值型求积公式: 在 a,b上取 a x0 x1 x2 xnb 作Lagrange插值,令,插值求积法,插值型求积公式的余项,例2 梯形公式的误差余项,即,例3取 x0 =a, x1 =0.5(a+b), x2 = b ,则 h=0.5(b a ),A0= (b-a)/6 A1=2(b-a)/3 A2= (b-a)/6,著名的 Simpson
4、 公式,定义: 若一个求积公式对f(x)= xi(i=0,1,.,m)能精确成立,但对f(x)= xm+1不精确成立,则称该公式具有m次代数精度。 令机械求积公式对f(x)= xi(i=0,1,.,n)精确成立,那么得线性方程组 当节点xk (k=0,1,.,n)给定且互异时,系数Ak可由上式确定。,定理:(n+1)个节点的求积公式为插值型的充要条件是该公式至少有n次代数精度.,类似有: Simpson公式具有3阶代数精度,例. 矩形公式 代数精度为0,解: 取f(x)= 1, x, x2 令求积公式准确成立,容易验证, 对f (x) = x3 求积公式式不能准确成立. 因此这一公式只具有2次
5、代数精度,取等距结点xj = a + jh时,插值型求积公式称为Newton-Cotes公式,定理: 当n为偶数时, n阶Newton-Cotes公式至少有(n+1)阶代数精确度。,Newton-Cotes公式代数精度至少为n,复合梯形求积公式,将积分区间a,b n 等分.令h=(b-a)/n . xj=a+jh,取,递推,得,给定允许误差界0,当,时,结束计算并以T2n作为定积分的近似值.,f=inline(sqrt(7782.5*sin(x).2+59621550*cos(x).2); t=0.25*pi*(f(0)+f(pi/2); n=1;h=pi/2;e=1;k=0; while e0.01 s=0.5*(t+h*sum(f(.5*h:h:pi/2); e=abs(s-t);t=s; n=2*n;h=h/2;k=k+1 end 4*t,ans=4.8707e+004, (循环次数k= 2),复合梯形公式计算,将程序第二行改为t=0.5*pi*f(pi/2)用积分部分和式计算,循环次数k=13,f=inline(sqrt(7782.5*sin(x).2+59621550*cos(x).2); T=4*quad(f, 0, pi/2) Ans=4.8707e+004,MATLAB求定
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 人教版七年级道法下册 第三单元第六课 传承核心思想理念(上课、学习课件)
- 招商与销售宠物食品的技巧试题及答案
- 2025全年红家具买卖合同样本
- 2025中外合资经营企业合同汽车零部件生产
- 2025路灯采购合同范本
- 2025设备采购合同
- 2025长途汽车停车场租赁合同
- 美容师行业创新与实践案例分析试题及答案
- 政治经济学选择题
- 车辆保养记录在评估中的运用试题及答案
- 钢筋机械连接技术规程 JGJ 107-2016
- 德育主题班会 《遵义会议》教学课件
- 马克思主义基本原理期末复习考试要点
- MTI学位论文写作要求课件(PPT 70页)
- 行业标准:TSG T7007-2016 电梯型式试验规则
- 生产计划与物料管理PMC
- 2022版义务教育(语文)课程标准(含2022年修订和新增部分)
- ADAScog(老年痴呆量表—认知)
- 药品经营中存在的问题及监管对策
- ZN73A(X)永磁式户内中压真空断路器
- 水利工程供水定价成本监审办法
评论
0/150
提交评论