[数学]第3章数值积分1.pptVIP

第3章 数值积分 引言 引言 引言 引言 引言 由定积分定义 引言 5.1 Newton-Cotes求积公式 由Lagrange插值,任何一的函数 都可以近似的表示成 其中 为简便起见,取节点为等分 现在关键是求 以此类推得Cotes系数表: Newton Cotes积分公式 常用的几个积分公式 梯形公式(n=1) Simpson公式(n=2) Newton公式(n=3) Cotes公式(n=4) 例题 3.1.2 Newton-Cotes公式截断误差及代数精度 几个常用的求积公式的代数精度 1.T 公式的代数精度 2. S–公式的代数精度 因此S-公式具有三次代数精度。 同理可得N-公式具有三次代数精度，C-公式具有五次代数精度。 3.2 复化求积公式 将区间[a,b]适当分割成若干个字区间，对每个子区间使用求积公式，构成所谓的复化求积公式，这是提高积分精度的一个常用的方法。 3.2.1定步长复化求积公式 1.复化梯形求积公式 一般地将[a,b]区间n等分，则 所以 而 定步长复化梯形求积公式算法 2.复化Simpson公式 类似于梯形公式： 定步长复化Simpson求积公式算法 例题 3.2.2变步长求积公式 变步长梯形求积公式 算法3.2.1变步长梯形求积算法 变步长Simpson求积公式 变步长Sinmpson求积算法 例题 3.3 Romberg求积公式 3.3.1外推法基本思想 以较小的计算量为代价,达到提高数值结果的精度是外推法的中心思想. 3.3.2 Romberg求积算法 Romberg求积算法 例 用四点复化梯形公式计算 解：四点复化梯形公式就是将区间[0，1]三等分得到 0 1 而梯形公式的结果为 例 试利用函数 的数据表分别用复化梯形公式、复化Simpson公式计算下列积分的近似值。 数据表 0.958811 1/2 0.8414710 1 0.9767267 3/8 0.8771926 7/8 0.9896158 1/4 0.9088517 3/4 0.9973978 1/8 0.9361556 5/8 1 0 也就是 解: 两种复化公式分别计算如下： 根据已知点的数据，需要用到九点复化梯形公式： 以上两种算法对区间采用不同等分，计算量大体一致，定积分精确到小数点后七位的值是0.9460831，Simpson公式精度要高一些。 对于复化抛物型公式： 在这里n=4 ,步长 例 利用复化梯形公式和复化Simpson公式分别求下列定积分 ，若要使精度 达到 ε=10-6 ，问各需将区间[0,1]多少等分？ 解 由于 从而 于是有 由复化梯形公式和复化Simpson公式的误差表示式 得到 根据上面的估计分别取 则只要 可分别解出 可见 满足同样的精度 要求 复化梯形公式 需将区 间 167 等分 复化抛 物 线 公式 只 需将区间 3 等分 * * * * 例 用梯形公式，Simpson公式和 Cotes 公式求积分 解:利用梯形公式 利用 Simpson 公式 得 利用Cotes公式得 而原积分为 相对而言，Cotes求积公式精度最高，梯形求积公式精度最低。 例 用梯形公式和Simpson公式 计算积分 解: 由梯形公式得 由Simpson公式得