计算方法复习重点1.pptxVIP

引论一、误差、有效数字定义1 绝对误差，简称误差：误差限：相对误差：相对误差限：引论定义2 例1 42.195, 0.0375551, 8.00033, 2.71828,按四舍五入写出上述各数具有四位有效数字的近似数.，引论定理１ ，例3 解：设取n位有效数字，相对误差限?*r=引论二、数值运算的误差估计引论引论三、病态问题与条件数 引论四、算法的数值稳定性考虑初始数据误差在计算中的传播问题. 引论五、避免误差危害的若干原则1.避免‘大数’除以‘小数’ 绝对值太小的数不宜做除数这里分子的误差被扩大104倍,再如若将分母变为0.0011,即分母只有0.0001的变化时,计算结果却有了很大变化 引论2、两个相近的数相减，会严重损失有效数字 例如x =1958.75，y =1958.32都具有五位 有效数字，但x-y=0.43只有两位有效数字 通常采用的方法是改变计算公式,例如当 很接近时,由于用右端代替左端公式计算,有效数字就不会损失 引论当x很大时可作相应的变换 则用右端来代替左端。 引论引论３、防止‘大数’吃‘小数’引论六、简化计算步骤，减少运算次数 例：x255=xx2x4x8x16x32x64x128 原先要做254次乘法现只需14次即可 例： 如计算多项式 p(x)=anxn? an-1xn-1 ? … ? a1x ? a0 的值 若直接计算akxk,再逐项相加，一共要做 n+(n-1)+…+2+1=n(n+1)/2次乘法和n次加法 引论如果将前n项提出x，则有 p(x)=（anxn-1? an-1xn-2 ? … ? a1 ）x ? a0 =((anxn-2?an-1xn-3?…? a2)x?a1）x ? a0 =(…(anx ? an-1)x?…?a2)x ? a1)x ? a0写成递推公式 于是,这种多项式求值的算法称为秦九韶算法,只做n次乘法和n次加法,程序实现简单 插值法问题的提出函数解析式未知,通过实验观测得到的一组数据, 即在 某个区间[a, b]上给出一系列点的函数值 yi= f(xi)或者给出函数表xx0x1x2……xnyy0y1y2……yny=p(x)y=f(x)定理1 n次代数插值问题的解是存在且惟一的 单项式：拉格朗日插值：a?b 且依赖于x对于线性插值，其误差为对于抛物插值（二次插值），其误差为 例 已知f(x)的观测数据 x 1 2 3 4f(x) 0 -5 -6 3构造插值多项式 解: 四个点可以构造三次插值多项式, 将数据 代入插值公式，有例 已知=100, =121, 用线性插值估计 在x=115时的截断误差解: 由插值余项公式知 因为 f[x1,x2]- f[x0,x1]x2 – x0均差与牛顿插值：xif[xi]f[xi,xi+1]f[xi,xi+1,xi+2]f[xi,xi+1,xi+2]x0f(x0)x1f(x1)f[x0,x1]x2f(x2)f[x1,x2]f[x0,x1,x2]x3f(x3)f[x2,x3 ]f[x1,x2,x3]f[x0,x1,x2 ,x3]………例 2.12 已知 x = 1, 4, 9 的平方根值，求解：xif[xi]f[xi,xi+1]f[xi,xi+1,xi+2]114293N(x)=f(x0)+ (x- x0) f[x1,x0]+ (x- x0) (x-x1) f[x1,x0,x2]N2(7)=1+(7-1)*0.33333+ (7-1)*(7-4)*(-0.01667)= 2.69992例2.14已知 f(x) = x7+ x4+ 3x+ 1 求 f [20, 21, … 27 ] 及 f [20, 21, … 27, 28 ]解: 由差商与导数之间的关系埃尔米特插值：例6P36例： 已知函数 y= f(x) 的数据如下表所示, 求次数 不超过三次的Hermite的插值多项式H3(x)使 H3(xi) = yi (i=0,1,2) H′3(xi) = y′i 分段线性插值三次样条插值：函数逼近正交多项式：勒让德多项式 切比雪夫多项式定理6， 例3最佳平方逼近用正交函数族求最佳平方逼近最小二乘法：数值积分与数值微分定义 （代数精度） 设求积公式（4.1）对于一 切次数小于等于m的多项式(或)是准确的，而对于次数为m+1的多项式是不准确的，则称该求积公式具有m次代数精度（简称代数精度） 定理 n+1个节点的求积公式 为插值型求积公式的充要条件是公式 至少具有n次代数精度。 牛顿-柯特斯公式：n1 1/2 1/2 2 1/6 4/6 1