算法的误差与稳定性.docVIP

算法的误差与稳定性 实验名称： 实验一 算法的误差与稳定性 指导教师： 数值分析实验组 实验时数： 2 实验设备：安装了Matlab 、C ++、VF 软件的计算机 实验日期：年月 实验地点： 第五教学楼北802或902 实验目的： 掌握舍入误差的概念，理解数值稳定性。 实验准备： 1. 在开始本实验之前，请回顾教科书的相关内容； 2. 需要一台准备安装Windows XP Professional操作系统和装有数学软件的计算机。 实验内容及要求 B 题 舍入误差在数值计算中是一个很重要的概念，在实际计算中，如果选用了不同的算法，由于舍入误差的影响，将会得到截然不同的结果，因此，选取数值稳定性的算法，在数值计算中是十分重要的。 对n =0,1,2, ,20计算定积分 y n =?x n e x -1d x 01 分别采用下面两个递推公式进行计算，并比较实验结果分析出哪个算法是稳定，并给出具体原因。 递推公式（1）y n =1-ny n -1(n =1,2, ,20) ； 递推公式（2）y n -1=1-y n (n =20,19, ,1) 。 n 说明：实验过程应包括对问题的简要分析、求解方法、求解步骤、程序及其必要的图表等内容。 实验过程： 本实验所选题为B 题 实验分析： x -1dx =1-e 递推公式为 方案1 y n =1-ny n (n =1,2,3....20) 当=0时y 0=?10e -1 ??y n =1-ny n -1(n =1, 2, , 20) ?-1??y 0=1-e 方案2 y n -1= ，则?1 0x n 1-y n (n =20,19, ,1) 当0 11111取递推初值y ≈[+]=(+1) 202e (20+1) 20+142e 1-y n ?y =(n =20,19, ,1) ??n -1n 递推公式为 ? 11?y =(+1) 20?42e ? n 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 Y1 Y2 0.63212055882856 0.63212055882856 0.36787944117144 0.36787944117144 0.26424111765712 0.26424111765712 0.20727664702865 0.20727664702865 0.17089341188538 0.17089341188538 0.14553294057308 0.14553294057308 0.12680235656152 0.12680235656153 0.11238350406936 0.11238350406930 0.10093196744509 0.10093196744560 0.09161229299417 0.09161229298962 0.08387707005829 0.08387707010385 0.07735222935878 0.07735222885769 0.07177324769464 0.07177325370770 0.06694777996972 0.06694770179987 0.06273108042387 0.06273217480180 0.05903379364190 0.05901737797298 0.05545930172957 0.05572195243238 0.05719187059731 0.05272680864948 19 20 实验结果分析： -0.02945367075154 0.05091744430931 1.55961974427919 0.03256855812313 由递推公式（1）知当y1(1)=1-exp(-1)时，yn 应当为精确解, 递推公式的每一步都没有误差的取舍，但计算结果y 18(*)=0.05719187059731y 17(*)=0.05545930172957 y19(*)出现负值。由此看出当n 较大时，用递推公式(*)中的y n (*) 近似yn 是不正确的。主要原因是初值y1=0.63212055882856不是精确值，设误差e (y 1(*))，由递推公式（*）知 则有 e (y n (*))=-ne (y n -1(*)), e (y n (*))=-10e (y n -1(*))=-n 2e (y 2(*))=(-10) n e (y 1(*)) 误差 e (y 1(*))的(-10) n 倍，由此可见，递推公式计算的误差不仅取决于初值的误差，公式的精确性，还依赖于误差的传递即递推计算的稳定性。 由递推公式（**）y 20≈0.032568558123130 由e (y n -1) =y n 为估计值，并不精