《数值分析》讲：数值积分.pptVIP

理学院 陈丽娟 《数值分析》第四讲 数值积分与微分 §4.1 引言 第四章：数值积分与数值微分 1、积分的概念 设 任取 做 则称 可积，极限值称为函数 在区间[a,b]上的 定积分，记为: Riemann积分 第四章：数值积分与数值微分 黎曼是世界数学史上最具独创精神的数学家之一,著作不多，却异常深刻，富于对概念的创造与想象，思想极其深邃难以理解。许多奠基性、创造性的工作，直接影响了19世纪以后的数学发展，在黎曼思想的影响下数学许多分支取得了辉煌成就。 ■ 黎曼几何、流形、微分流形、椭圆几何的创始人 爱因斯坦用黎曼几何将广义相对论几何化。黎曼几何是现代理论物理必备的数学基础。 第四章：数值积分与数值微分 ■ 完善微积分理论的出杰人物之一 微积分理论严谨性论证的杰出贡献者有：黎曼、波尔查诺、柯西、阿贝尔、狄利克莱、维尔斯特拉斯等等。柯西证明连续函数必定可积，黎曼指出可积函数不一定连续。黎曼推广了博里叶展开式成立的狄利克莱条件，即三角级数收敛的黎曼条件等等。 ■ 解析数论、与复变函数的里程碑 ■ 组合拓扑的开拓者 ■ 代数几何的奠基人 ■ 在数学物理、微分方程等领域贡献卓著 第四章：数值积分与数值微分 2、积分的计算 Riemann积分从定义上基本不可算 第四章：数值积分与数值微分 理论求解定积分基本看运气 第四章：数值积分与数值微分 3、数值积分的思想 第四章：数值积分与数值微分 所谓数值积分就是 4、代数精度 如果某求积公式对次数不超过m的多项式能够准确成立，对m+1次多项式不一定准确成立，则称该求积公式具有m次代数精度。 第四章：数值积分与数值微分 考查 第四章：数值积分与数值微分 即 第四章：数值积分与数值微分 如梯形公式 所以梯形公式具有1次代数精度 第四章：数值积分与数值微分 图示 第四章：数值积分与数值微分 5、插值型求积公式 第四章：数值积分与数值微分 定理 求积公式 充分性 具有n阶代数精度， 当且仅当该求积公式为插值型求积公式。 第四章：数值积分与数值微分 必要性 所以该求积公式为插值型。 §4.2 Newton-Cotes系数 第四章：数值积分与数值微分 1、Cotes系数 将区间[a,b] n 等分 第四章：数值积分与数值微分 第四章：数值积分与数值微分 第四章：数值积分与数值微分 第四章：数值积分与数值微分 第四章：数值积分与数值微分 由 得求积公式 第四章：数值积分与数值微分 第四章：数值积分与数值微分 此时，求积公式为 当 时可得 第四章：数值积分与数值微分 同理, 第四章：数值积分与数值微分 考虑Simpson公式 第四章：数值积分与数值微分 Simpson公式具有三次代数精度 第四章：数值积分与数值微分 至少具有n+1次代数精度。 2、数值积分的误差分析 以梯形公式误差为例 第四章：数值积分与数值微分 因为 保号且可积，由积分中值定理得 第四章：数值积分与数值微分 所以 第四章：数值积分与数值微分 同理 误差取决于区间[a,b]的长度。 第四章：数值积分与数值微分 3、复化求积法 第四章：数值积分与数值微分 （1）复化梯形公式 （2）复化Simpson公式 （3）复化Cotes公式 第四章：数值积分与数值微分 4、复化求积公式的误差 因为复化梯形公式为 第四章：数值积分与数值微分 同理，复化Simpson和复化Cotes公式的误差分别为 例 第四章：数值积分与数值微分 a=0; b=2.5; n=4; h=(b-a)/n; x=a:h:b; y=-x.^3+2*x.^2+2*x+1; s=0.0; for k=2:length(x)-1 s=s+y(k); end tn=h*(y(1)+y(length(x))+2*s)/2; 复化梯形程序 第四章：数值积分与数值微分 a=0;b=2.5;n=4;h=(b-a)/n; x=a:h:b; y=-x.^3+2*x.^2+2*x+1; s1=0.0; for k=2:length(x)-1 s1=s1+2*y(k); end s2=0.0; for i=1:length(x)-1 x0=(x(i+1)+x(i))/2; y0=-x0^3+2*x0^2+2*x0+1; s2=s2+4*y0; end tn=h*(y(1)+s1+s2+y(length(x)))/6; 第四章：数值积分与数值微分 复化Simpson程序 第四章：数值积分与数值微分 第四章：数值积分与数值微分 §4.3 Romberg算法 第四章：数值积分与数值微分 考察 可以得到两个结果 区间二分后的误差是二分前后差值的三分之一。 第四章：数值积分与数值微分 结果二