数值计算方法　全套课件.ppt

6.2 牛顿-柯特斯求积公式 6.2.1 公式的导出 6.2.2 牛顿-柯特斯公式的代数精度 6.2.3 低阶求积公式的余项 6.2.4 复化求积法 6.2 牛顿-柯特斯求积公式6.2.1 公式的导出 2 柯特斯系数的求取 n ? ? ? ? ? ? ? ? 1 1/2 1/2 ? ? ? ? ? ? ? 2 1/6 4/6 1/6 ? ? ? ? ? ? 3 1/8 3/8 3/8 1/8 ? ? ? ? ? 4 7/90 16/45 2/15 16/45 7/90 ? ? ? ? 5 19/288 25/96 25/144 25/144 25/96 19/288 ? ? ? 6 41/840 9/35 9/280 34/105 9/280 9/35 41/840 ? ? 7 751/17280 3577/17280 1323/17280 2989/17280 2989/17280 1323/17280 3577/17280 751/17280 ? 8 989/ 28350 5888/ 28350 -928/ 28350 10496/ 28350 -4540/ 28350 10496/ 28350 -928/ 28350 5888/ 28350 989/ 28350 柯特斯求积系数表： 例如：n=1时，有 n=2时，有 柯特斯系数的性质 (2) 系数有对称性。 (3) 当n≥8时开始出现负值的柯特斯系数。 (1) 取f(x)≡1，则 f(n+1)(x)≡0, Rn(f)≡0，于是 梯形公式 当n=1时, 有 相当于用直线P(x)代替f(x)计算积分。 3 常用的低阶牛顿-柯特斯公式 5.3 可线性性化模型的最小二乘拟合 例 已知观测数据(1,–5),(2,0),(4,5),(5,6) ，试用最小二乘法求形如 的经验公式。 解 作超定方程组 即 法方程组为 求得 a=1.537650114 b= －6.432976311 所求经验公式为 5.5 多项式拟合 1 直线拟合 作超定方程组 n 记号?指对i从1到n取和 法方程组 2 二次拟合、抛物拟合 作超定方程组 有 法方程组 3 一般情形 超定方程组的系数矩阵 例 给定函数y=f(x)的实例数据表。 试用最小二乘法求二次拟合多项式。 解 设二次拟合多项式 写出其正则方程组 x 1 2 3 4 6 7 8 y 2 3 6 7 5 3 2 将计算结果代入正则方程组 解得 a0=- 1。3185 , a1= 3.4321 , a2=-0.3864 二次拟合曲线 第6章 数值积分和数值微分  本章的问题： 计算定积分∫abf(x)dx的近似值。 必要性: 如果f(x)的原函数是F(x)，则 等. 实际问题中常有些被积函数没有表达式,只是通过观测得到一些离散的数据点, 这样的定积分也只能用数值方法近似计算. （牛顿-莱布尼兹公式） 但有些定积分的被积函数的原函数不能用初等函数明显表示, 牛顿-莱布尼兹公式不能用.如 第6 章 数值积分和数值微分 6.1 数值积分概述 6.2 牛顿-柯特斯公式 6.3 变步长求积和龙贝格算法 6.4 高斯型求积公式 6.5 数值微分 6.1.2 代数精度 代数精度与节点数的关系 6.1.3 插值求积公式 6.1.4 构造插值求积公式的步骤 用待定系数法构造插值求积公式 按差商的定义 4.3.2牛顿插值公式 1。牛顿插值公式的建立 牛顿插值多项式 (f(x)的前n+1项） 牛顿插值余项 (f(x)的最后一项） 牛顿插值多项式的构成 零阶差商 一阶差商 二阶差商 乘积因子 1 ……. 2。牛顿插值的特点 (1) P(x)次数不超过n次，项数不超过n+1项。各项系数是差商表上对角线的各阶差商值。 (2) P(x)满足插值条件，在节点上f(xi )=P(xi ) . (3) 增加一个节点，只需增加一项。 n次牛顿插值多项式 计算牛顿插值多项式的步骤 (1)