数值积分分和数值微分.pptVIP

第五章数值积分与数值微分 5.1引言 5.2 Newton- Cotes公式 5,3复化求积公式 5.4龙贝格求积公式 5.5高斯型积分 56数值微分 2013326数计学《数计算》课程建设生QAB率 合*  【本章重点】 1.求积公式代数精确度定义,应用此定义建立求积公式 2.梯形公式, Simpson/公式及它们的复合公式及余项表达式 和误差估计。 3.外推原理及 Romberg求积公式。 4. Gauss型求积公式及求积节点(即 Gauss点)的充要条件, Gauss-- Legendre求积公式与 Gauss- chebysher求积公式 5.求积公式收敛性与稳定性概念和基本结论。 2013326数计学《数计算》课程建设生QAB率 合*  课前思考】 1.什么是代数精确度?如何确定和验证求积公式的代数精确 度和验证求积公式的代数精确度次数?如何利用代数精确 度概念确定求积公式系数与节点? 2.写出梯形公式、 Simpson公式及复合梯形公式及复合 Simpson 公式和它们的截断误差。复合梯形公式和复合 Simpson公式的 误差是步长h的几阶小量?若给定误差,如何根据复合求积 公式的误差估计求积区间等分数 3.什么是 Gauss型求积公式?什么是 Gauss点?如何求 Gauss 点? 2013326数计学《数计算》课程建设生QAB率 合*  §1引言 一、数值求积的基本思想 计算定积分的方法: 1、牛抓菜布尼兹公式了f(=Fa)-F( 但是求函数fx)的原函数F(x)不一定比计算积分容易, 例如函数,sinx2找不到用初等函数表示的原函数 另外若给出的函数fx)是数据表,也不好求函数的积分。 2013326数计学《数计算》课程建设生QAB率 合*  2、积分中值定理『f(xk=(b-a)(5) 但是点ξ的具体位置一般不知道,故难以准确算出∫f的值。 有两种近似方法 种是梯形法,用 f(a)+f(b) 代替∫( 得到[r(xlk=(b-a)(a)+fb 梯形公式 另一种是矩形法, 一(ok=b-00—→在矩形公式 atb f(rdx=(b-a)f( 中矩形公式 F(三 f(xdx 今右矩形简称矩形公式 2013326数计学数自计算课程建设组QAB DESk  义:公式(k=A/() 叫做数值求积公式(机械求积) 其中x称为求积节点, A称为求积系数,亦称伴随节点x的权。 2013326数计学《数计算》课程建设生QAB率 合*  代数精度的概念 定义若某个求积公式对于次数不超过m的多项式均能准确成立, 但对于m1次的多项式就不能准确成立,则称此求积公式的代数 精度为m。 般要使求积公式「f(x)dx≈∑A(x)具有m次 精度,只要令它使∫(x)=1,x,x2,…,都能准确成立,即: A-a Cb 大-k (em+ a+1 2+1  例1考察「f(x)d≈(1)+2f()+f0有几次代数精度。 解:逐次检查公式是否精确成立 当f(x)=1:∫1x=2=(+2+1 当fx)=:∫,(hx=0=(+1+2×0+ 当fx)=x2:如=≠(-1)2+2x02+1=1 故:代数精度=1 2013326数计学《数计算》课程建设生QAB率 合*  三、插值型的求积公式 设给定一组节点:4≤xx2…x≤b 且已知函数在这些节点处的函数值f(x)(i-0,1…,n,由第二 章知可以作插值函数Ln(x,由于L1(x)为多项式,所以其积 分可以很容易求得: 1()=L(d=∑f(x)(xx=∑可广1(x)df(x) 记A=(x)则上式=∑A F foda 定义称公式:1()=f(x)d∑Af(x)=L( 为插值型的求积公式 其中A=1()(k=0,1,…,m)为求积系数 其余项为:R「门-=1-1,=(2m( 在合4s  定理含有n+1个节点的插值型求积公式2()=∑4f(x) k=0 至少具有n次代数精度的充分必要条件是它是插值型的。 证明 显然对于次数不超过n次的多项式,R印等于0,所以说明, 插值型求积至少是n次代数精度的 假设∫f(=∑A()至少具有次代数精度, 所以用插值基函数l(x)当作(x)代入,上式精确成立,即 4(x)d==∑4(x)=A 所以1=∑Af(x)为插值型的求积 k=0 2013326数计学《数计算》课程建设生QAB率 合*