计算机数值方法第四章-数值积分与微分PPT.ppt

延安大学 数学与计算机科学学院;本章要点:;第4章 数值积分 ;一、 引言;以上这些现象,牛顿-莱布尼兹很难发挥作用;如果我们用两端点“高度”;如果改用区间中点;对定积分;式中;定义1. 若求积公式 ;梯形公式和矩形公式均具有一次代数精度。;例：;三、插值型的求积公式;作为积分;由插值型的求积公式的余项可推得;四、求积公式的收敛性与稳定性; 牛顿－柯特斯公式是指等距节点下使用 拉格朗日插值多项式建立的数值求积公式。;而;有;即有;等距节点时;; n 1 1/2 1/2 2 1/6 4/6 1/6 3 1/8 3/8 3/8 1/8 4 7/90 16/45 2/15 16/45 7/90 5 …　　…　　　…　　　…　　…　　…; 二、 偶数阶求积公式的代数精度; 三、低阶牛顿－柯特斯公式及其余项;上式称为梯形求积公式,也称两点公式，记为;梯形公式的余项为;2.辛普森公式及其余项;称为辛普森求积公式，也称三点公式或抛物线公式;3.柯特斯公式及其余项;求积公式为;上式称为柯特斯求积公式，也称五点公式;例: 用n=2和n=3的牛顿-柯特斯公式;§2 复合求积法;Newton-Cotes积分公式;低阶Newton-Cotes积分公式;低阶Newton-Cotes积分公式;低阶Newton-Cotes积分公式;Newton-Cotes积分法的稳定性;b;1. 复合求积公式;1.1 复合求积公式;x0;1.2 复合梯形求积公式;1.3 复合辛普森求积公式;复合辛普森公式;xi+1/2;1.4 复合柯特斯求积公式;2. 复合求积法的余项及收敛阶;1, 2, 4 阶牛顿-柯特斯积分公式的余项分别为：;;R1(f ) 是h的2阶无穷小量，所以;因此复合辛普森公式是4阶收敛的。;因此复合柯特斯公式是6阶收敛的。;2.5 三种复合求积公式比较;例1.;0 1 0.125 00.25 00.375 00.5 00.625 00.75 00.875 0 1 0分别由复合梯形、辛普森、柯特斯公式得：;比较三个复合公式的计算结果：;3. 步长的自动选择; 1.变步长公式 在实际计算中，借助于计算机来完成积分步长 h 的自动选择，即采用变步长求积公式。具体地讲，就是将步长逐次折半，反复利用复合求积公式，直到满足精度要求为止。 ; 2.变步长复合辛普森公式 逐次将区间［a，b］分成21，22，…，2m等份，并按复合抛物线公式逐次计算积分得到S1，S2，…，Sm，而; 2.变步长复合辛普森公式; 2.变步长复合辛普森公式;4.自动选择步长的算法步骤;§3 龙贝格积分方法;梯形公式,Simpson公式,Cotes公式的代数精度分别为;1. 复合梯形公式的递推化;各节点为;--------(3);(3)式称为递推的梯形公式;2. 龙贝格算法;即;复合辛普森公式;也就是，相邻的两个梯形值的外推，;同理，由复合辛普森公式的余项;一般的，若记;设以;T－数表; ;3. 算法设计;引言 求积公式 (1) 当求积系数 、求积节点 都可以自由选取时，其代数精确度最高可以达到多少次? 下面的引理可以回答上述问题。 ;