第3章节 数值积分1.pptVIP

第3章 数值积分 引言 引言 引言 引言 引言 由定积分定义 引言 3.1 Newton-Cotes求积公式 由Lagrange插值,任何一的函数 都可以近似的表示成 其中 为简便起见,取节点为等分 现在关键是求 以此类推得Cotes系数表: Newton Cotes积分公式 常用的几个积分公式 梯形公式(n=1) Simpson公式(n=2) Newton公式(n=3) Cotes公式(n=4) 例题 3.1.2 Newton-Cotes公式截断误差及代数精度 几个常用的求积公式的代数精度 1.T 公式的代数精度 2. S–公式的代数精度 因此S-公式具有三次代数精度。 同理可得N-公式具有三次代数精度，C-公式具有五次代数精度。 3.2 复化求积公式 将区间[a,b]适当分割成若干个字区间，对每个子区间使用求积公式，构成所谓的复化求积公式，这是提高积分精度的一个常用的方法。 3.2.1定步长复化求积公式 1.复化梯形求积公式 一般地将[a,b]区间n等分，则 所以 而 定步长复化梯形求积公式算法 2.复化Simpson公式 类似于梯形公式： 定步长复化Simpson求积公式算法 例题 3.2.2变步长求积公式 变步长梯形求积公式 算法5.2.1变步长梯形求积算法 变步长Simpson求积公式 变步长Sinmpson求积算法 例题 f=inline((1+x).^-1); quad(f,0,1) ans = 0.6931 3 Romberg求积公式 3.1外推法基本思想 以较小的计算量为代价,达到提高数值结果的精度是外推法的中心思想. §2.5 Romberg求积算法 一.对低精度公式经过组合构造高精度公式 用T(2n)做为积分Ｉ的近似值时，其误差大约为 若我们不是把误差舍掉，而是把误差补偿到T(2n)后，我们会得到精度更好的求积公式－Simpson公式： 由（1），（2），（3）组成的方法称为Romberg算法。 序列{T(n)},{S(n)},{C(n)}和{R(n)}分别称为梯形序列，Simpson序列，Cotes序列和Romberg序列。 上述用若干个积分近似值推算出更为精确的积分近似值的方法，称为外推算法。得到Romberg序列后还可以继续外推，得到新的求积序列，称为Richardson外推算法。但由于在新的求积序列中，其线性组合的系数分别为： 因此，新的求积序列与前一个序列结果相差不大，故通常外推到Romberg序列为止。 可以证明，梯形序列，Simpson序列，Cotes序列和Romberg序列均收敛到积分值，且每次外推可使误差阶提高二阶。 Romberg求积算法 function [ s,n ] = rombg2( a,b,eps1) f=inline(4/(1+x^2)); a=0; b=1; eps1=0.000005; s=10; s0=0; k=2; t=[]; t(1,1)=(b-a)*(f(b)+f(a))/2; while abs(s-s0)eps1 h=(b-a)/2^(k-1); w=0; if h~=0 for i=1:2^(k-1)-1 w=w+f(a+i*h); end t(k,1)=h*(f(a)/2+w+f(b)/2); for l=2:k for i=1:(k-l+1) t(i,l)=(4^(l-1)*t(i+1,l-1)-t(i,l-1))/(4^(l-1)-1); end end s=t(1,k); s0=t(1,k-1); k=k+1; n=k; else s=s0 n=-k end end 两个二阶精度的梯形公式组合起来就得到一个四阶精度的Simpson公式，截断误差提高了二阶，起到了加速收敛的作用。 类似的，可以证明： 解：[s,n]=rombg2(0,1,5e-6) s = 3.1416 n = 7 * *