数值计算第七讲.pptxVIP

数值分析与算法 (7) Numerical Analysis Algorithms 清华大学 计算机系 喻文健 Wenjian Yu 2 第七章 数值积分与数值微分 Wenjian Yu 3 数值积分的基本概念 Wenjian Yu 4 数值积分 目的与用途 经典问题: 算几何形体的面积、体积，力学中 物体的重心位置 例:铝制波纹瓦的长度问题 由一块平整的铝板压制而成. 若每个波纹的高度(自中心线)为1英寸, 周期为 2英寸, 做4英尺长波纹瓦需多长铝板? 第二类椭圆积分, 无法解析求出 ! Wenjian Yu 5 数值积分基本思想 . . . 积分系数积分节点 希望用较少的计算量得到较准确的结果 Wenjian Yu 6 插值型求积公式 中矩形公式 梯形公式 Wenjian Yu 7 积分余项与代数精度 反映了计算的截断误差 插值余项的积分 衡量求积公式准确度的另一个指标 注意: 对某些情况，代数精度并不是越高越好 Wenjian Yu 8 积分余项与代数精度 (至少0次代数精度) 它至少有n次代数精度 即插值型求积公式 Wenjian Yu 9 求积公式的收敛性与稳定性 (一系列求积公式的性质) 积分问题一般不太敏感 Wenjian Yu 10 求积公式的收敛性与稳定性 这是控制数值计算误差能达到的最佳情况 要尽量寻求稳定的求积公式 Wenjian Yu 11 牛顿-柯特斯公式 Wenjian Yu 12 Newton-Cotes公式 就是n阶牛顿-柯特斯公式 n=1, 1/2, 1/2 n=2, 1/6, 2/3, 1/6 n=4, 7/90, 16/45, 2/15, 16/45, 7/90 Cotes系数表 一系列求积公式 便于使用 Wenjian Yu 13 Newton-Cotes公式 n=1, 1/2, 1/2 n=2, 1/6, 2/3, 1/6 n=4, 7/90, 16/45, 2/15, 16/45, 7/90 n=8, Cotes系数表 思考题 梯形公式 Simpson公式 Cotes公式 中矩形公式可看成是n=0时的特例. 例题(板书) Wenjian Yu 14 Newton-Cotes公式 关键看积分: (n阶公式至少有n次代数精度) 一般不用n=3对应的N-C公式 Wenjian Yu 15 低阶N-C公式的积分余项 不保号, 无法用积分中值定理 2 详细过程自己看书 Wenjian Yu 16 稳定性、收敛性 n=8, 实际只使用n8的偶数阶N-C公式, 也看出代数精度不是越高越好 Wenjian Yu 17 复合求积公式 Wenjian Yu 18 复合求积公式 (composite quadrature) 积分误差: n增大,误差减小 仍是机械求积公式 Wenjian Yu 19 复合求积公式 2阶准确度 Wenjian Yu 20 复合求积公式 与复合梯形公式对比, 例7.4 Wenjian Yu 21 复合求积公式 步长折半的复合求积公式计算 积分余项公式包含被积函数的高阶导数, 很难应用. 常常动态地确定步长h 常用的动态减小步长策略是: 步长折半, 利用已算出结果 复合梯形公式的情况 递推化的复合梯形公式: (逐渐减小, 直到满足精度要求) 只需再计算新增节点的函数值 Wenjian Yu 22 复合求积公式 步长折半的复合求积公式计算 复合Simpson公式的情况 很少使用中矩形公式的原因 与梯形公式有相同的代数精度/准确度, 计算量更小 可类似构造复合中矩形公式, 但在步长折半时, 无法重用以前的结果 Wenjian Yu 23 Remberg积分算法 Wenjian Yu 24 复合梯形公式的余项展开式 (可与Richardson外推结合) Th7.5 所有小区间的积分求和: Wenjian Yu 25 复合梯形公式的余项展开式 Th7.5 所有小区间乘h/2求和: Wenjian Yu 26 Richardson Extrapolation (“0”代表未经外推的原始公式) Wenjian Yu 27 Richardson Extrapolation Wenjian Yu 28 Romberg算法 龙贝格算法 列三角形表格, 按行依次计算 计算公式 可证明: 具有2k+1次代数精度 (类似高阶差商的计算) Wenjian Yu 29 Romberg算法 做等距的节点分布, 充分利用其上函数值得到最高准确度的结果 要求被积函数充分光滑! Wenjian Yu 30 Romberg算法 Wenjian Yu 31 Romberg算法 Wenjian Yu 32