数值分析4-06.ppt

上一页 下一页 湘潭大学数学与计算科学学院 * §6 奇异积分计算 一、分部积分与变量替换法 二、区间截去法 三、Gauss型求积法 四、乘积积分法 湘潭大学数学与计算科学学院 * 几种主要类型： 1、被积函数在积分上限或下限处的极限存在， 但函数值不存在。 例 2、积分上限是 ，或积分下限是 。 3、在积分限的两端存在积分奇异点； 例如 湘潭大学数学与计算科学学院 * 4.在积分上限和下限之间的某已知点处存在积分奇点。 5.在积分区间内某未知点处存在积分奇点。 标准数值求积方法不能直接应用的。 湘潭大学数学与计算科学学院 * 一、分部积分与变量替换法 基本思想 ： 采用变量替换法将无穷区间变为有限区间； 采用分部积分或变量替换法 消除或降低被积函数的奇异性 。 湘潭大学数学与计算科学学院 * 例1 积分 其中 在 上是光滑的， 作变量替换 得 积分区间变为有限区间. 在 处的性态为： 假设 湘潭大学数学与计算科学学院 * 则 积分 在 处为光滑的被积函数。 例2 积分 其中 在 上是一阶连续可微函数. 湘潭大学数学与计算科学学院 * 对I作分部积分，则 因此，I原有的奇异性得以消除. 湘潭大学数学与计算科学学院 * 二、区间截去法 基本思想 ： 若能够选取小数 (或大数R), 使得 或 上的积分值处在允许误差范围之内， 即 或 则可将正常积分 或 作为 （或 的近似. 湘潭大学数学与计算科学学院 * 例3 计算积分 上连续， 由于在[0,1]上， 所以 假如精度要求为 则可取 从而 湘潭大学数学与计算科学学院 * 例4 计算积分 因为 当 时有 所以对于 的精度要求 湘潭大学数学与计算科学学院 * 三、Gauss型求积法 Gauss型求积公式也可计算某些奇异积分 或无穷区间上的积分 例：Gauss-Chebyshev求积公式 其中 其本身就是对奇异积分 的一种算法. 湘潭大学数学与计算科学学院 * 四、乘积积分法 是接近奇异或奇异的函数， 是光滑函数. 考虑积分 其中 基本思想 ： 满足下列条件： 湘潭大学数学与计算科学学院 * 2、 非常容易计算。 由此可得 所以 因此， 的一个逼近序列。 是 湘潭大学数学与计算科学学院 * 注意： 为 上 的步长为 的分段多项式插值. 在实际应用中，常取 例5、计算积分 取 为分段线性插值多项式： 其中 湘潭大学数学与计算科学学院 * 由Lagrange插值余项定理知： 在 上二阶连续可微时， 当 有 因此 可见，当n充分大时， 可作为 的近似值. 湘潭大学数学与计算科学学院 * 其中 湘潭大学数学与计算科学学院 * 作变量替换 则 因此 可算。 上一页 下一页