第5章_数值积分和数值微分教材教学课件.pptVIP

计算定积分有微积分基本公式;数值积分/Numerical Integration /;1、解决第一个问题；节点xi 和系数Ai如何选取，即选取原则; 定义5.1 若求积公式;例 1 确定形如; 解得:A0=A2=1/3, A1=4/3.;例3 试确定参数A0, A1和x0, x1,使求积公式;5.1 插值型求积公式;例：对于[a, b]上1次插值，有; 为了简化计算,取等距节点xk=a+kh,(k=0,1,2,…,n,; 由于构造Newton-Cotes公式需要Cotes系数,将其列表如下:; （3）牛顿求积公式：;例3 求n =4的Newton-Cotes公式及误差.; 一般地 , Newton-Cotes公式的截断误差为;5. 3 复化求积公式;可见,复化梯形公式是收敛的。而且,要使|R[Tn]|? ,只要;? 复化 Simpson 公式：;如果记M4=; 例 已知函数; 例 利用复化梯形公式和复化Simpson公式分别计算上例中定积分,若使精度? =10-6,问各需取n为多少?;变步长求积方法;可见步长减半时;对于复化梯形公式;由此得;其中,xk=a+kh,k=0,1,2,…,n,h=(b-a)/n,;复化Simpson公式能加工成更高精度的公式吗? ;所以有;所以有; 用Tm(k)(m=1,2,3,4)分别表示把区间2k等分的复化梯形公式,复化Simpson公式, 复化Cotes公式和Romberg求积公式.;一般有：; 实际计算可按下表顺序进行 ; 解 按递推公式计算,结果如下;? 理查德森外推法 /* Richardson’s extrapolation */; 前面介绍的 n+1个节点的 Newton -Cotes求积公式， 其特征是节点是等距的。这种特点使得求积公式便于 构造，复化求积公式易于形成。但同时也限制了公式 的精度。 n是偶数时，代数精度为n+1， n是奇数时， 代数精度为n 。;构造具有2n+1次代数精度的求积公式;这样,高斯求积公式是 ; 定理;(2)求出pn(x)的n个零点x0 , x1 , … xn 即为Gsuss点. ;解 按 Schemite 正交化过程作出正交多项式: ;故两点Gauss公式为 ; 区间[-1,1]上权函数?(x)=1的Gauss型求积公式,称为Gauss-Legendre求积公式,其Gauss点为Legendre多项式的零点. ;例 用3点Gauss公式计算积分 ;由于; 区间[0,??)上权函数?(x)=e-x的Gauss型求积公式,称为Gauss-Laguerre求积公式,其Gauss点为Laguerre多项式的零点. ;Gauss-Laguerre求积公式为 ; 区间(-?,??)上权函数?(x)= 的Gauss型求积公式，;数值微分就是用离散方法近似地求出函数在某点的导数值. 按照Taylor展开原理可得 其中h步长。; 设Ln(x)是?(x)以a?x0x1…xn?b为节点的n次Lagrange插值多项式, 则取;如取n=1的线性插值L1(x)=[(x-x0)?(x1)-(x-x1)?(x0)]/h, (其中h=x1-x0)可得数值微分的二点公式:;可得数值微分的三点公式:;练习 题