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

三种复化求积分算法的精度分析 　　【摘 要】分别利用复化梯形公式、复化Simpson公式和复化Gauss-Legendre I型公式对定积分进行运算，得到近似数值解，并对各算法的精度和计算复杂度进行了比较与分析。数值举例结果表明，三种复化求积分算法的运算结果均在绝对误差限ε=5e-8内，并且在相同的精度下，复化Gauss-Legendre I型公式的步长和计算量最小。 　　【关键词】复化梯形公式；复化Simpson公式；Gauss-Legendre公式 　　1 引言 　　数值积分是计算数学的基本内容，在工程技术和科学计算中起着十分重要的作用，当积分的精确值不能不能求出时，数值积分就变得越来越重要。通常数值积分的计算常利用机械积分来实现，其基本思想为： 　　（1） 　　其中Ak≥0和xk∈[a，b]称为求积系数和求积节点，k=1，2，3 … n。Newton- Cot’s方法可以为数值积分提供思路，主要思路就是在节点确定时利用插值多项式的积分。对节点和系数都使用待定法使代数精度达到最高阶的一类方法是Gauss型公式，常见的有Gauss-Legendre公式、Gauss-Chebyshev公式等。本文利用复化梯形公式、复化Simpson公式和复化Gauss-Legendre I型公式来计算几例常见定积分的近似数值解，并且对每种算法的精确度、计算复杂度进行了比较与分析。 　　2 理论模型 　　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] （5） 　　由于 　　f（x）∈C2[a，b] 　　且 　　，（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 三种复化算法步长的事前估 　　函数 复化梯形 　　求积公式 复化Sim