第五篇数值积分.pptVIP

抛物型求积公式： 思考： 代数精度是否是越高越好？ 梯形公式的误差 定理证明 思考： 结论是什么？ 推广： 加权Gauss积分公式 Matlab 积分函数 Matlab 积分函数 求解思路： 用数值积分代替积分 q=quad(‘fun’,a,b,tol,trace,p1,p2,…) q=quad8(‘fun’,a,b,tol,trace,p1,p2,…) 参数‘fun’是被积函数,可以是表达式字符串、内联函数、M函数文件名，被积函数的自变量一般采用字母 x；a、b分别是积分的上、下限，都为确定的值；tol 是一二元向量，第一个元素控制相对误差，第二个元素控制绝对误差；trace若取非零值，将以动态图形展现积分的整个过程，若取零值，则不画图，其缺省值为0；p1、p2是向被积函数传递的参数。 在调用函数时，前三个参数是必须的，其余参数可缺省。 符号积分： int(f) — 对f表达式的缺省变量求积分 int(f,v) — 对f表达式的v变量求积分 int(f,v,a,b) — 对f表达式的v变量在（a,b) 区间求定积分 §5.3 积分方程的数值求解 怎么求解？ * * 第八章 数值积分 近似计算 思路 利用插值多项式 则积分易算。 问题的提出： ? 在[a, b]上取 a ? x0 x1 … xn ? b，做 f 的 n 次插值多项式 ，即得到 Ak 由 决定， 与 无关。 节点 f (x) 插值型积分公式 误差 例：对于[a, b]上1次插值，有 考察其代数精度。 f(x) a b f(a) f(b) 梯形公式 解：逐次检查公式是否精确成立 代入 P0 = 1： = 代入 P1 = x ： = 代入 P2 = x2 ： ? 代数精度 = 1 二次插值求积公式： Simpson公式 定 义 　　　若某个求积公式所对应的误差R[ f ]满足：R[ Pk ]=0 对任意 k ? n 阶的多项式成立，且 R[ Pn+1 ] ? 0 对某个 n+1 阶多项式成立，则称此求积公式的代数精度为 n 。 代数精度： 怎样验证代数精度： 注：形如 的求积公式至少有 n 次代数精度 ? 该公式为插值型（即： ） ? 辛普森(Simpson)求积公式的误差： 怎么办？ 复合求积： 高次插值有Runge 现象，故采用分段低次插值 ? 分段低次合成的 复合求积公式。 ? 复合梯形公式： 在每个 上用梯形公式： = Tn /*中值定理*/ 怎么办？ 4 4 4 4 4 = Sn 注：为方便编程，可采用另一记法：令 n’ = 2n 为偶数， 这时 ，有 ? 复合Simpson公式： 复化求积例： 复化求积例： 两种方法谁好？ 给定精度 ?，如何取 n ? 通常采取将区间不断对分的方法，即取 n = 2k 上例中2k ? 68 ? k = 7 注意到区间再次对分时 可用来判断迭代是否停止。 §5.2 高斯型积分 构造具有2n+1次代数精度的求积公式 将节点 x0 … xn 以及系数 A0 … An 都作为待定系数。令 f (x) = 1, x, x2, …, x2n+1 代入可求解，得到的公式具有2n+1 次代数精度。这样的节点称为Gauss 点，公式称为Gauss 型求积公式。 例：在两点数值积分公式中，如果积分点也作为未知量， 则有4个未知量,可以列出4个方程：（在[-1,1]为例） 可解出： 数值积分公式 具有3阶代数精 度，比梯形公式 1阶代数精度高 权函数 例：求 的 2 点 Gauss 公式。 解：设 ，应有 3 次代数精度。 ? + ? 1 0 1 1 0 0 ) ( ) ( ) ( x f A x f A dx x f x 代入 f (x) = 1, x, x2, x3 不是线性方程组，不易求解。 　　　 x0 … xn 为 Gauss 点 ? 与任意次数不大于n 的多项式 P(x) （带权）正交。 定理 求 Gauss