三种复化求积分算法的精度分析.docxVIP

【摘 要】分别利用复化梯形公式、复化 simpson 公式和复化 gauss-legendre i 型公式对定 积分进行运算，得到近似数值解，并对各算法的精度和计算复杂度进行了比较与分析。数值 举例结果表明，三种复化求积分算法的运算结果均在绝对误差限ε =5e-8 内，并且在相同的 精度下，复化 gauss-legendre i 型公式的步长和计算量最小。 【关键词】复化梯形公式；复化  simpson  公式；  gauss-legendre  公式 引言 数值积分是计算数学的基本内容，在工程技术和科学计算中起着十分重要的作用，当积 分的精确值不能不能求出时，数值积分就变得越来越重要。通常数值积分的计算常利用机械积分来实现，其基本思想为： 1） 理论模型 2.1 复化梯形求积公式 将区间 [a ，b] 划分成 n 等分，分点 xk=a+kh（， k=1，2，3?n），在每个子区间 [xk ，xk+1] k=1， 2， 3 ?n-1 ）上采用梯形式，则得到 2） 记 3） 上式（ 3）为复化梯形公式，其余项可由式 ，（ a≤η≤ b） （ 4） 得 ，η k∈ [xk ， xk+1] 由于 f （ x）∈ c2[a ， b]  （ 5） 且 ，（ 0≤ k≤ n-1 ） （ 6） 所以 ?∈（ a， b），使 （7） 于是复化梯形公式余项为 （8） 2.2 复化 simpson 求积公式 将区间 [a ，b] 划分为 n 等分，在每个子区间 [xk ， xk+1] 上采用 simpson 式，若记，则得 （9） 记 （10） 上式（ 10）为复化 simpson 求积公式，其余项可由式 ，（ a≤η≤ b） （ 11） 得 ，η k∈ [xk ， xk+1] （ 12） 于是当 f （x）∈ c4[a ， b] 时，与复化梯形公式相似有 ，η∈ [a ，b] （ 13） 2.3 复化 gauss-legendre i 型求积公式 gauss 型求积公式是具有最高代数精度的插值求积公式。通过适当选取求积公式（ 1）的 节点ε =5e-8 和求积系数 ak≥ 0 和 xk ∈ [a ， b] （ k=1， 2， 3?n），可使其代数精度达到最高 的 2n+1 次。利用特殊区间 [-1 ，1] 上 n+1 次 legendre 正交多项式的根作为节点， 我们可以建 gauss-legendre 型求积公式。 将区间 [a ，b] 划分成 n 等分，分点 xk=a+kh（，k=1，2，3 ? n），在每个子区间 [xk ，xk+1] （ k=1，2，3? n-1 ）上采用 2 点 gauss-legendre i 型求积公式 （14） 在[a ， b] 区间上的复化积分公式为 （15） 上式（ 15）称为复化 gauss-legendre i 型求积公式。 于是当 f （x）∈ c4[a ， b] ，时，复化 gauss-legendre i 型求积公式的余项表达式为 ，（ a≤η≤ b） （ 16） 3 数值举例 先考察下面等式（ 17）右边定积分的近似值 17） 分别用复化梯形公式、 复化 simpson 公式和复化 gauss-legendre i 型公式做运算， 求出其在绝对误差限为ε =5e-8 内的近似数值解。 假定 18） 因此 ， （ 19）所以 ， （ 20） 对于复化梯形公式有 21） 所以 n≥ 1791.6 （ 22） 因此取步长 n=1792 （ 23） 对于复化 simpson 求积公式有 24） 所以 n≥ 20.1 （25） 因此取步长 n=21 （ 26） 对于复化 gauss-legendre i 型求积公式有 27） 所以 n≥ 18.2 （28） 因此取步长 n=19 （ 29） 同理也可以考察等式 和 （ 30） 右端定积分的近似数值值，具体结果见表 1。 1 三种复化算法步长的事前估函数 复化梯形 求积公式 复化 simpson 求积公式 复化 gauss-legendre i 型求积公式 1792 21 19 2457 14 12 7019 24 22 表 2 三种复化算法的计算结果函数 复化梯形 求积公式 复化 simpson 求积公式 复化 gauss-legendre i 型求积公式 -0.405465126309431 -0.405465118046333 -0.405465098225125 1.820478483584408 1.820478477218769 1.820478423657262 7.389056127230221 7.389056126214707 7.389056073169591 3 三种复化算法的精度分析函数 复化梯形