数值第七章数值微分与数值积分解答.pptVIP

上机作业 第260至261页 数值实验1题, 2题任选一题 ? 当区间[a, b] 2n等分时, 则有 在定理中以 h/2 代替 h 得 上式乘以4减去 T(h) 再除以3, 记之为 T1(h), 得 T1 (h) 逼近 I 的速度是 O ( h4 )阶, 效果比 T (h)好, 它不是别的, 就是Simpson序列. ? 类似地 上式乘以16减去 T1(h) 再除以15, 记之为 T2(h), 得 T2 (h) 逼近 I 的速度是 O ( h6 )阶, 效果比 T1 (h)好, 它不是别的, 就是Cotes公式序列. ? 对Cotes公式序列进行同样处理得到Romberg公式序列. Richardson外推加速方法 也称为Romberg求积算法 ? 收敛性说明: 如果 f (x) 充分光滑, 那么梯形公式序列, Simpson公式序列, Cotes公式序列, Romberg公式序列均收敛到所求的积分值. ? 对于 f (x)不充分光滑的函数也可用Romberg算法计算, 只是收敛慢一些. 也可以直接使用复化Simpson公式计算. 例 用Bomberg算法计算积分 解 在[0, 1]上仅是一次连续可微 用Romberg算法计算结果见下表 2 4 1 8 16 32 0.5 0.426777 0.407018 0.401812 0.400463 0.400118 0.402369 0.400432 0.400077 0.400014 0.400002 0.400302 0.400054 0.400009 0.400002 0.400050 0.400009 0.400002 §5 Gauss型求积公式 ? 基本思想 设计求积公式: 在节点数 n 固定时, 适当地选取求积节点 { xk }与求积系数{ Ak }, 使求积公式具有最高的代数精确度. 例 确定x1, x2, A1, A2, 使求积公式 具有最高次的代数精确度. x2 x1 1 ? 选取 (A1 , A2 , x1 , x2)使该求积公式对 f (x) = 1, x, x2, x3 时等号成立. ?1 ? 对 f = 1, x, x2, x3 积分精确成立 ? 四个方程四个未知数 梯形公式与Gauss求积公式的比较 ? 对1, x求积公式精确成立(1次代数精确度) ? 对1, x, x2, x3求积公式精确成立(3 次代数精确度) x3 x1 1 x2 ?1 ? 选取(A1, A2 , A3 , x1, x2 , x3) 使该求积公式对 f (x) = x0, x1, x2, x3, x4, x5 时等号成立. 区间[?1, 1]上的Gauss求积公式 ? 对 f = x0, x1, x2, x3, x4, x5 求积公式等号成立 syms A1 A2 A3 x1 x2 x3 eq1=A1+A2+A3-2 eq2=A1*x1+A2*x2+A3*x3 eq3=A1*x1^2+A2*x2^2+A3*x3^2-2/3 eq4=A1*x1^3+A2*x2^3+A3*x3^3 eq5=A1*x1^4+A2*x2^4+A3*x3^4-2/5 eq6=A1*x1^5+A2*x2^5+A3*x3^5 [A1,A2,A3,x1,x2,x3]=solve(eq1,eq2,eq3,eq4,eq5,eq6) Matlab求解 ? 选取(A1, A2, …, An , x1, x2 , …, xn)使该求积公式对f (x) = x0, x1, x2, …, x2n?1时等号成立. 区间[?1, 1]上的Gauss求积公式 ? 2n个方程2n个未知数, 非线性方程, 解是否存在唯一？若存在, 如何求解？ ? 对 f = x0, x1, x2, …, x2n?1求积公式等号成立 研究一般区间上的积分: 具有2n－1次代数精确度, 则称这组节点 {xk}为 Gauss 点, 上述公式称为Gauss型求积公式. 定义 如果一组节点x1, x2, …, xn?[a, b]能使求积公式 Gauss型求积公式 由上述 2n个方程确定全部 的2n个待定参数 xk , Ak ( k=1, …, n), 使求积公式至少具有 2n?1次代数精确度. 但上述方程组是非线性方程组, 求解十分困难. 一般利用 正交多项式来求出Gauss点与求积系数. Gauss型求积公式 ? 对 f = x0, x1, x2, …, x2n?1求积公式等号成立 定理 一求积公式的节点 x1, x2, …, xn是高斯点的充分必要条件是以这些节点为根(零点)的多项式 与任何次数不超过n－1的多项式P(x)正交 高斯求积基本定理 高斯求积公式的稳定性定理 定理 高斯求积公式的