[数值分析1.ppt

第三章 数值积分与数值微分 §3.1 Newton-Cotes求积公式 总结 3.1.3 Newton-Cotes公式的误差分析 3.1.2 Newton-Cotes求积公式 3.1.1 插值型求积法 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 第 三 章 数 值 积 分 和 数 值 微 分 学习目标： 理解求积公式及代数精度概念，掌握确定求积公式的代数精度的方法，掌握 Newton-Cotes 求积公式、Romberg算法及Gauss求积公式的构造技术、特点及余项形式。掌握复化梯形求积公式、复化Simpson求积公式的构造技术及余项形式。了解上述求积公式的适用类型并会熟练使用这些公式做数值积分。了解数值微分法及 Richardson 加速技术，了解Newton-Cotes求积公式、Gauss 求积公式的稳定性问题。 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 一、数值求积的基本思想 积分 只要找到被积函数 f (x)原函数F(x)，便有 牛顿—莱布尼兹(Newton—Leibniz)公式 实际困难：大量的被积函数（ , sin x2 等）, 找不到用初等函数表示的原函数；另外, f (x)是（测量或数值计算出的）一张数据表时，牛顿—莱布尼兹公式也不能直接运用。 前言 例如，一块铝合金的横断面为正弦波，要求原材料铝合金板的长度。也就是f (x)=sinx 从x=0到x=b的曲线弧长L，可用积分表示为 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 这是一个椭圆积分计算问题。 卫星轨道的计算也一个椭圆积分计算问题。找不到被积函数的原函数。然而，用数值分析中的数值积分方法计算，并不是很难的计算问题。 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 　　本章讨论常用的数值求积公式及它们的误差估计和代数精度，而对数值微分只作简单介绍。 积分中值定理：在[a, b]内存在一点 ，有 f(?) 成立。 ? 就是说, 底为b-a 而高为f(?)的矩形面积恰等于所求曲边梯形的面积 . 问题 在于点?的具体位置一般是不知道的，因而 难以准确算出 f(?)的值．我们将f (?)称为区间[a, b]上的平均高度．这样,只要对平均高度f(?)提供一种算法，相应地便获得一种数值求积方法． Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 如果用两端点的“高度”f(a)与f(b)的算术平均作为平均高度f (?) 的近似值，这样导出的求积公式 : 便是我们所熟悉的梯形公式 . 而如果改用区间中点 的“高度”f (c)近似地取代平均高度f (?)，则又可导出所谓中矩形公式(今后简称矩形公式)： (1.1) (1.2) Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 3.1 Newton-Cotes求积分式 3.1.1 插值型求积法 近似计算 思路 利用插值多项式 则积分易算。 关键是f(x) Evaluation on