数值计算期末复习资料.docVIP

PAGE PAGE 7 引言 例： 利用秦九韶算法求多项式p(x)=x5-3x4+4x2-x+1,在x=3 时的值。（课本11页习题3） 解：p(x)=(x-3)x4+4x2-x+1 =((x-3)x2+4)x2-x+1 =(((x-3)x2+4)x-1)x+1 =((0*x2+4)x-1)x+1 =(4x-1)x+1 =(4*3-1)*3+1 =34 练习. x=3, f=2x5 – 5x4 – 4 x3 + 3x2 – 6 x + 7 解：f=(2x-5)x4 – 4 x3 + 3x2 – 6 x + 7 =((2x-5)x – 4 ) x3 + 3x2 – 6 x + 7 =(((2x-5)x – 4 ) x + 3)x2 – 6 x + 7 = ((((2x-5)x – 4 ) x + 3)x – 6) x +7 =((((2*3-5)*3 – 4 )* 3 + 3)*3 – 6) *3 +7 =-11 练习：用二分法求方程ex+10x-2=0在[0,1]内的近似根，要求误差不超过1/2*10-2.(课本11页，习题1） k ak bk xk f(ak)* f(xk) 0 0 1 0.5000 0 1 ? 0.5000 0.2500 0 2 ? 0.2500 0.1250 0 3 ? 0.1250 0.0625 0 4 0.0625 ? 0.0937 0 5 ? 0.0937 0.0781 0 6 0.0781 ? 0.0859 0 7 0.0859 ? 0.0898 ? 误差的来源： 误差限和有效数字： 误差限： 设以x代表x*的近似值，则 绝对误差为：|x-x*|，若有|x-x*|≤ε 则ε称为近似值x的绝对误差限，简称误差限，或称精度。 有效数字： 对于x*的近似值（规格化形式） x=±0.a1a2…an×10m (其中a1，a2，…，an是0-9之间的自然数） 如果误差|x-x*|≤1/2×10m-l,1 ≤l ≤n 则称近似值x有l位有效数字，或称x准确到第l位。 凡是由准确值经过四舍五入而得到的近似值，其绝对误差限等于该近似值末位的半个单位。 1、用四舍五入法取准确值的前n位作为近似值，则x*必有n位有效数字 2、有效数字相同的两个近似数，绝对误差限不一定相同。例如，设x1*=12345及x2*=0.12345均有五位有效数字，而它们的绝对误差限分别为0.5和0.000005 3、准确值被认为具有无穷位有效数字 相对误差与有效数字的联系： 设x,y分别是数x*,y*的近似值，考察x-y的相对误差 相对误差与有效数字： 由此可见：两个值相近的近似数相减，可能会造成有效数字的严重损失。 在实际计算时，需要加工计算公式，以避免这种情况发生。 8、已知e=2.71828…,试问其近似值x1=2.7,x2=2.71,x3=2.718各有几位有效数字？并给出它们的相对误差限。 解：x1=0.27×101,x2=0.271×101, x3=0.2718×101 |x1-e|=0.01828…≤ 1/2×101-2,有效数字2位， |x2-e|/=0.00828…≤1/2×101-2有效数字2位， |x3-e|=0.00028…≤1/2×101-4,有效数字4位 算法设计中，避免误差危害的要点： 避免除数绝对值远小于被除数绝对值的除法，否则可能会扩大舍入误差。 避免两相近数相减，否则会使有效数字严重损失。 尽可能防止大数“吃掉”小数，否则可能影响计算结果的可靠性。 简化计算步骤，将少运算次数。 选用数值稳定的算法 第一章 差值方法 例1： 解： 带入x0=100, 得 用p1(x)计算： 误差估计： 用p2(x)计算： 练习：课本54页1，求作f(x)=sin x在节点x0=0的5次泰勒多项式，并估计插值误差。 2、拉格朗日插值 求作n次多项式pn(x),使满足条件 pn(xi)=yi , i=0,1, … , n 这就是所谓拉格朗日插值 点xi（互不相同），称为插值节点。 1.2 拉格朗日插值公式 1、线性插值 问题3：求作一次式p1(x),使满足条件 p1(x0)=y0 , p1(x1)=y1 几何意义：y=p1(x)表示通过两点（x0,y0),(x1,y1)的直线，因此，一次插值又称线性插值 它的解p1(x)可表为点斜式： 解：已知x0=100,y0=10,x1=121,y1=11,x=115, 带入线性插值公式，得 插值基函数： 线性插值公式亦可表为对