第四讲数值积分.ppt

第四章：数值积分与数值微分 1、数值微分在常微分方程离散化的应用 例 记 第四章：数值积分与数值微分 整理得 第四章：数值积分与数值微分 即有 其中 事实上，本方程的解析解为 第四章：数值积分与数值微分 解析解 近似解 h=(b-a)/10 近似解 h=(b-a)/20 误差 h=(b-a)/10 误差 h=(b-a)/20 u1 1.4552 1.4974 1.4769 0.0422 0.0217 u2 1.9033 1.9825 1.9441 0.0792 0.0408 u3 2.3384 2.4476 2.3946 0.1092 0.0562 u4 2.7555 2.8861 2.8227 0.1306 0.0672 u5 3.1510 3.2929 3.2240 0.1419 0.0730 u6 3.5229 3.6647 3.5959 0.1418 0.0730 u7 3.8715 4.0002 3.9377 0.1287 0.0662 u8 4.1988 4.3003 4.2510 0.1015 0.0522 u9 4.5094 4.5684 4.5398 0.0590 0.0304 第四章：数值积分与数值微分 N=10时 第四章：数值积分与数值微分 N=20时 第四章：数值积分与数值微分 N=20时 N=10时 如果已知方程 数值微分 整理得 第四章：数值积分与数值微分 其中 1、数值微分在偏微分方程离散化的应用 例： 离散得 第四章：数值积分与数值微分 扩散方程初边值问题 令 第四章：数值积分与数值微分 第四章：数值积分与数值微分 第四章：数值积分与数值微分 第四章：数值积分与数值微分 * 性 * 性 * 性 * 性 * 性 * 性 第四章：数值积分与数值微分 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算法 第四章：数值积分与数值微分 考察 可以得到两个结果 结果一： 区间二分后的误差是二分前后差值的三分之一。 第四章：数值积分与数值微分 结果二： 其中 二分前的步长 第四章：数值积分与数值微分 即 令 则 第四章：数值积分与数值微分 同理，由 直接验证得 又 第四章：数值积分与数值微分 记 称Romberg公式 解： 令 则令 因此 例：计算 对照值 （1）二分前的步长 （2）二分前的区间中值 二分点 第四章：数值积分与数值微分 （1）二分前的步长 （2）二分前的区间中值 第四章：数值积分与数值微分 　 0.9397933 0.9460831 0.0062898 　 0.9445135 0.9460831 0.0015696 　 0.9456909 0.9460831 0.0003922 　 0.9459850 0.9460831 0.0000981 　 0.9460596 0.9460831 0.0000235 　 0.9460769 0.9460831 0.0000062 　 0.9460815 0.9460831 0.0000016 　 0.9460827 0.9460831 0.0000004 　 0.9460830 0.9460831 0.0000001 　 0.9460831 0.9460831 0.0000000 二分次数 区间个数数 第四章：数值积分与数值微分 第四章：数值积分与数值微分 第四章：数值积分与数值微分 　 0.9397933 0.0062898 　 0.9445135 0.001