数值分析课件教学模板.pptVIP

*************************************辛普森法则1基本辛普森法则辛普森法则使用二次多项式（抛物线）近似区间上的被积函数。基本辛普森法则的公式为：∫??f(x)dx≈(b-a)/6[f(a)+4f((a+b)/2)+f(b)]，相当于用抛物线下的面积近似曲线下的面积。2精度与误差基本辛普森法则对三次或更低次多项式的积分是精确的。对一般函数，其截断误差为O(h?)，优于梯形法则。误差项与区间宽度的四次方和函数的四阶导数有关。复合辛普森法则的截断误差为O(h?)。3复合辛普森法则复合辛普森法则将积分区间均分为2n个子区间，每两个子区间应用一次基本辛普森法则。公式为：∫??f(x)dx≈h/3[f(a)+4∑?f(x????)+2∑?f(x??)+f(b)]，其中h=(b-a)/(2n)。n必须是偶数。龙贝格积分R(0,0)R(0,1)R(0,2)R(0,3)R(1,0)R(1,1)R(1,2)R(2,0)R(2,1)R(3,0)龙贝格积分是一种高效的数值积分方法，它结合了复合梯形法则和Richardson外推技术，能够显著提高积分精度。该方法构建一个三角形的积分表（龙贝格表），每行代表一个外推层次，通过逐步细化积分区间和外推操作来提高精度。龙贝格表的第一列R(i,0)是使用2^i个子区间的复合梯形法则结果。后续列通过递推公式计算：R(i,j)=R(i,j-1)+(R(i,j-1)-R(i-1,j-1))/(4^j-1)。这一过程相当于消除误差的低阶项，可以证明R(i,j)的截断误差为O(h^(2j+2))，精度随着j的增加而快速提高。高斯求积法积分点选择高斯求积法的关键在于精心选择积分点和权重，使得n点公式能够精确积分最高阶次为2n-1的多项式，这是同样使用n个函数求值的方法所能达到的最高精度。积分点是特定多项式（勒让德多项式）的零点，而非等距分布。公式与权重高斯-勒让德求积法表示为：∫??1f(x)dx≈∑?w?f(x?)，其中x?是勒让德多项式P?(x)的零点，w?是对应的权重。权重可通过解线性方程组或使用解析公式计算。积分区间[a,b]可通过变量替换转换为[-1,1]。变体与应用高斯求积法有多种变体，适用于不同的权函数和积分区间：高斯-切比雪夫求积法适用于无穷区间积分；高斯-埃尔米特适用于带指数权函数的积分；高斯-拉盖尔适用于半无穷区间积分。这些变体在量子力学计算、统计分析等领域有广泛应用。数值微分概述数值微分的必要性数值微分在函数解析表达式未知、仅有离散数据点，或者导数解析计算过于复杂时特别有用。它是数值模拟、优化算法、微分方程求解等领域的基础工具。然而，数值微分本质上是病态问题，对数据噪声极为敏感。基于差分的近似方法数值微分的基本思路是用差分近似微分。常见的差分格式包括前向差分、后向差分和中心差分。这些方法可通过泰勒展开推导，并分析其截断误差阶数。中心差分通常具有更高的精度，被广泛采用。高阶导数与高精度方法对于高阶导数，可以通过反复应用一阶差分公式获得，但这会导致误差累积。更好的方法是直接使用专门的高阶导数差分公式。提高精度的方法包括使用多点差分公式和Richardson外推法，后者可显著提高数值微分的精度。差分公式1前向差分前向差分公式近似一阶导数：f(x)≈[f(x+h)-f(x)]/h。它的截断误差为O(h)，是一阶精度。前向差分适用于计算区间起点的导数，或者当前向的函数值容易获取时。它的计算简单，但精度相对较低。2后向差分后向差分公式为：f(x)≈[f(x)-f(x-h)]/h，同样具有O(h)的截断误差。后向差分适用于计算区间终点的导数，或者当后向的函数值容易获取时。它的特性与前向差分类似，是一阶精度。3中心差分中心差分公式为：f(x)≈[f(x+h)-f(x-h)]/(2h)，具有O(h2)的截断误差，是二阶精度。中心差分通常比前向和后向差分更精确，是实际计算中常用的选择。它的缺点是需要计算两个函数值。4高阶差分高阶导数的差分公式可以通过重复应用一阶差分获得，也可以直接使用多点公式。例如，二阶导数的中心差分公式为：f(x)≈[f(x+h)-2f(x)+f(x-h)]/h2，具有O(h2)的截断误差。Richardson外推法精度提升原理Richardson外推法利用误差的渐近行为提高数值微分的精度。对于具有偶数阶截断误差的方法（如中心差分），其误差可以表示为：E(h)=c?h2+c?h?+...。通过组合不同步长的近似值，可以消除误差的低阶项。外推公式