数值计算及其不确定性问题(简化)要点.ppt

变步长Sinmpson求积算法 例题 谢谢！ * X255=xx2x4x8x16x32x64x128 * * 要计算N+1个行列式，每一个需要n!(n-1)次乘法，计算N+1个行列式共需要(n+1) n!(n-1)次乘法,计算全部的xi需要(n+1) n!(n-1)次乘法，n次除法，总计算量为N=[(n2-1)n!+n]flop。 n=20,N= 9.707?1020flop,换算为时间，设在10亿（109)次/每秒的计算机上工作，则需9.707?1011秒，一天为24×3600 =86400秒，一年为86400 × 365秒， 9.707?1011秒÷86400 × 365秒=30781年（3万多年）。显然方法不可行；但若改用gauss消元法，只需3060flop,即使在一台小型机上，也可在一秒中内完成。 * 在109flop/秒(每秒10亿次）计算机上计算，T= 9.707?1011秒=3万多年。 * 1 计算机数系不同于我们熟悉的实数系 2 在数值计算中，经常将实数表示成浮点形式 * * * 以上错误,是因为计算机计算时做加减法要 “对阶”,“对阶”的结果使大数吃掉了小数.产生了误差.为了避免由于上述原因引起的计算结果严重失真,可以把某些算式改写成另一种等价的形式. 以上错误,是因为计算机计算时做加减法要 “对阶”,“对阶”的结果使大数吃掉了小数.产生了误差.为了避免由于上述原因引起的计算结果严重失真,可以把某些算式改写成另一种等价的形式. * Newton公式(n=3) Cotes公式(n=4) 例题 5.1.2 Newton-Cotes公式截断误差及代数精度 几个常用的求积公式的代数精度 1.T 公式的代数精度 2. S–公式的代数精度 因此S-公式具有三次代数精度。 同理可得N-公式具有三次代数精度，C-公式具有五次代数精度。 5.2 复化求积公式 将区间[a,b]适当分割成若干个字区间，对每个子区间使用求积公式，构成所谓的复化求积公式，这是提高积分精度的一个常用的方法。 5.2.1定步长复化求积公式 1.复化梯形求积公式 一般地将[a,b]区间n等分，则 所以 而 定步长复化梯形求积公式算法 2.复化Simpson公式 类似于梯形公式： 定步长复化Simpson求积公式算法 例题 5.2.2变步长求积公式 变步长梯形求积公式 算法5.2.1变步长梯形求积算法 变步长Simpson求积公式 §2 Newton’s Interpolation ? 牛顿插值 /* Newton’s Interpolation */ 1 2 … … … … n?1 1 + (x ? x0) ? 2 + … … + (x ? x0)…(x ? xn?1) ? n?1 Nn(x) Rn(x) ai = f [ x0, …, xi ] §2 Newton’s Interpolation 注：? 由唯一性可知 Nn(x) ? Ln(x)， 只是算法不同，故其余项也相同，即 ? 实际计算过程为 f (x0) f (x1) f (x2) … f (xn?1) f (xn) f [x0, x1] f [x1, x2] … … … … f [xn?1, xn] f [x0, x1 , x2] … … … … f [xn?2, xn?1, xn] f [x0, …, xn] f (xn+1) f [xn, xn+1] f [xn?1, xn, xn+1] f [x1, …, xn+1] f [x0, …, xn+1] Newton插值多项式的计算 表5-2 Newton插值多项式可按表5-2计算。 xi yi= f (xi) 一阶差商 二阶差商 … n阶差商 ? x0 y0 ? ? ? ? 1 x1 y1 f [x0,x1] ? ? ? x-x0 x2 y2 f [x1,x2] f [x0,x1,x2] ? ? ? x3 y3 f [x2,x3] f [x1,x2,x3] ? ? ? ? ? ? ? ? ? ? xn yn f [xn-1,xn] f [xn-2,xn-1,xn] ? f [x0,x1,…,xn] ? n次Newton插值多项式Nn(x)为表5-2中对角线上的差商值 与右端因子乘积的和。 Newton插值公式计算举例 例3 用Newton插值公式计算例1中的ln11.5。 [解] 如果仍取点x0=11, x1=12, x2=13,作抛物线插值， 按表5-2计算，结果如下： xi yi= lnxi 一阶差商 二阶差商 ? 11 2