第9章_数值积分与数值微分精品.ppt

结束 * 如认为 则有 于是有: 记 这样我们从{Cn}序列又推出了收敛更快的{Rn}序列. {Rn}序列也称为龙贝格序列. 这样我们从收敛较慢的{Tn}序列只用了一些四则运算,便推出了收敛更快的{Sn}序列, {Cn}序列和{Rn}序列. 这个过程还可继续下去,但已意义不大.我们常将这四个序列排成如下的三角形数表(表7-2) 结束 * T1 T2 S1 T4 S2 C1 T8 S4 C2 R1 T16 S8 C4 R2 ﹕ ﹕ ﹕ ﹕ 表7-2 该表四个序列都是收敛的. 结束 * 例5 利用龙贝格方法计算 解:计算结果列如下表: i 2i T序列 S序列 C序列 R序列 0 1 3.00000 1 2 3.10000 3.13333 2 4 3.13118 3.14157 3.14212 3 8 3.13899 3.14159 3.14159 3.14159 4 16 3.14094 3.14159 3.14159 3.14159 这一结果与I=π相比较已有较好的精度. 算法公式： i=2,3,… 其中： 算法： 例5 利用龙贝格方法计算 程序：先建立函数文件，文件名, ff.m function y=ff(x) y=1/(1+x) 主程序：Romberg1.m a=0,b=1.5,m=5; h=b-a; T=zeros(m); T(1,1)=h*(ff(a)+ff(b))/2; for i=2:m for k=1:2^(i-2) ss(k)=ff(a+(k-0.5)*h); end T(i,1)=(T(i-1,1)+h*sum(ss))/2; for j=2:i T(i,j)=(4^(j-1)*T(i,j-1)-T(i-1,j-1))/(4^(j-1)-1); end if abs(T(i,j)-T(i,j-1))tol break end h=h/2; end T 结束 * §7.5 高斯型求积公式 由前面的讨论已经知道，以a=x0x1…xn=b为节点的N-C求积公式的代数精度一般为n或n+1,这时节点简单地按照闭式等距的方式确定。对一个求积公式而言，如果不固定节点的位置，在节点数目不变的情况下，代数精度能否提高，最多能达到多少?高斯型求积公式讨论的就是最高代数精度的求积公式.先看一个简单的例子，考虑两个节点的求积公式. 7.5.1 最高代数精度的求积公式 这里积分区间选为[-1,1] 不失一般性，因为对区间[a,b]，总可用变量替换 将化为[-1,1]上的积分.从而用（7.35）式得到解决。 结束 * 现在不固定节点x0,x1的位置为区间端点-1,1，而允许取在(-1,1)内，即选取x0,x1及A0，A1，使(7.35)具有最高的代数精度.为此，分别取 代入(7.35)式，得方程组 取前四个方程构成此方程组是因为只有4个待定数x0,x1，A0, A1，即代数精度m=3. 方程组(7.36)有一组解为 即求积公式 具有的代数精度为3 结束 * 从上面这一简单的例子可以看到，节点数目不变的情况下，求积公式的代数精度是可以提高的.下面就对一般问题进行讨论，即当节点数目为n+1时，求积公式的代数精度最高能达多少，怎样才能达到这一最高的代数精度. 设节点为x0,x1,…,xn，求积公式为： 共2n+2未知参数,可列2n+2个方程,由此可以推设在n+1个节点上的求积公式，其代数精度最多为2n+1，即 (7.38)的代数精度可以达到2n+1.通常确定节点xk与求积系数Ak不是通过解非线性方程组，而是利用正交多项式的性质来求得的.下面讨论更一般情况的高斯型求积公式 结束 * 其中ρ(x)≥0为权函数，为使此公式对f(x)为不超过2n+1次的多项式时能精确成立，记: 用ω(x)去除f(x)，则可表示成 其中，q(x)为商，r(x)为余，q(x),r(x)均为不超过n次的多项式，于是有 如果对任何不超过n次的多项式q(x)都有 结束 * 又因为 则 即求积公式(7.40)对不超过2n+1次的多项式能精确成立,而满足条件(7.43)只需q(x)与ω(x)为区间［a,b］上关于权函数ρ(x)正交，若取节点xk(k=0,1,2,…,n)恰为此正交多项式系中n+1次多项式的n+1个零点，而由正交多项式的性质可知这些根均为实根，无重根，且全部分布在(a,b)内,这样，对于给定的权函数总能构造出关于此权函数在［a,b］区间上的正交多项