第七章数值微分与数值积分初次修改稿.pptVIP

2021/2/11 * 上机作业 第207页 26题,27题任选一题 * * 2021/2/11 * 停机准则: 梯形值序列 Simpson序列 Cotes序列 Romberg序列 ? Romberg求积公式 2021/2/11 * 例 计算 =0.9207355 =0.9397933 =0.9445135 =0.9456909 解 先求梯形值序列 2021/2/11 * 2 4 1 8 0.9207355 0.9397933 0.9445135 0.9456909 0.9461459 0.9560869 0.9460833 0.9460830 0.9460831 0. ? 利用只有两三位有效数字的T1, … ,T8 经过三次外推得到7位有效数字. 可见加速的效果十分显著. 用Romberg算法计算如下 2021/2/11 * ? 理论依据: 复化梯形公式的余项展开. 记 ?定理 设 则 其中系数 ?k ( k=0, 1, …)是与 h 无关的常数. ? T (h) 逼近 I 的速度是 O ( h2 )阶. 2021/2/11 * ? 当区间[a, b] 2n等分时, 则有 在定理中以 h/2 代替 h 得 上式乘以4减去 T(h) 再除以3, 记之为 T1(h), 得 T1 (h) 逼近 I 的速度是 O ( h4 )阶, 效果比 T (h)好, 它不是别的, 就是Simpson序列. 2021/2/11 * ? 类似地 上式乘以16减去 T1(h) 再除以15, 记之为 T2(h), 得 T2 (h) 逼近 I 的速度是 O ( h6 )阶, 效果比 T1 (h)好, 它不是别的, 就是Cotes公式序列. 2021/2/11 * ? 对Cotes公式序列进行同样处理得到Romberg公式序列. Richardson外推加速方法 也称为Romberg求积算法 ? 收敛性说明: 如果 f (x) 充分光滑, 那么梯形公式序列, Simpson公式序列, Cotes公式序列, Romberg公式序列均收敛到所求的积分值. ? 对于 f (x)不充分光滑的函数也可用Romberg算法计算, 只是收敛慢一些. 也可以直接使用复化Simpson公式计算. 2021/2/11 * 例 用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 2021/2/11 * §5 Gauss型求积公式 ? 基本思想 设计求积公式: 在节点数 n 固定时, 适当地选取求积节点 { xk }与求积系数{ Ak }, 使求积公式具有最高的代数精确度. 2021/2/11 * 例 确定x1, x2, A1, A2, 使求积公式 具有最高次的代数精确度. x2 x1 1 ? 选取 (A1 , A2 , x1 , x2)使该求积公式对 f (x) = 1, x, x2, x3 时等号成立. ?1 2021/2/11 * ? 对 f = 1, x, x2, x3 积分精确成立 ? 四个方程四个未知数 2021/2/11 * 梯形公式与Gauss求积公式的比较 ? 对1, x求积公式精确成立(1次代数精确度) ? 对1, x, x2, x3求积公式精确成立(3 次代数精确度) 2021/2/11 * x3 x1 1 x2 ?1 ? 选取(A1, A2 , A3 , x1, x2 , x3) 使该求积公式对 f (x) = x0, x1, x2, x3, x4, x5 时等号成立. 区间[?1, 1]上的Gauss求积公式 2021/2/11 * ? 对 f = x0, x1, x2, x3, x4, x5 求积公式等号成立 2021/2/11 * 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求解 2021/2/11