实验一 误差传播与算法稳定性.docVIP

实验一 误差传播与算法稳定性 1.1 实验目的 体会稳定性在选择算法中的地位．误差扩张的算法是不稳定的，是我们所不期望的；误差衰竭的算法是稳定的．是我们努力寻求的，这是贯穿本课程的目标．通过上机计算，了解舍入误差所引起的数值不稳定性。 1.2 算法描述 概要：舍人误差在计算方法中是—个很重要的概念。在实际计算中，如果选用了不同的算法，由于舍入误差的影响．将会得到截然不同的结果。因此，选取稳定的算法，在实际计算中是十分重要的。 考虑一个简单的由积分定义的序列 ， (E.1) 利用分部积分易得 , 得递推公式 ， (E.2) 注意到 取 由利用(E.2)变形得到 (E.3) 计算方法：先估计一个,再反推要求的 ( )。 可取 （）。 1.3 实验内容 由递推关系(E.2)，可以得到计算(E.1)积分序列的两种算法．其一为(E.2)的直接应用，即 ， ， (E.4) 另外一种算法则是利用(E.2)变形得到的 ， ， (E.5) 1.4 实验要求 （1）分别用算法(E.4)、(E.5)并在计算中分别采用5位、6位和7位有效数字，请判断哪种算法能给出更精确的结果． （2）两种算法的优劣，与你的第一感觉是否吻合．请从理论上证明你实验得出的结果，解释实验的结果．设(E.4)中的计算误差为，由递推计算的误差为；算法(E.5)中的误差为，由向前递推计算的误差为．如果在上述两算法中都假定后面的计算不再引入其他误差，试给出与的关系和与的关系． （3）算法(E.4)中通常会很小，当增大时，的变化趋势如何?算法(E.5)中通常相对比较大，当减小时，误差又是如何传播的?也就是说比较一下上述两个算法，当某一步产生误差后．该误差对后面的影响是衰减还是扩张的． （4）通过理论分析与计算实验，针对算法(E.4)和(E.5)的稳定性，给出你的结论。 1.5 实验与分析 用算法(E.4)初始值采用7，8位有效数字 y(1)=vpa((1-1/exp(1)),7);for i=1:19 y(i+1)=1.0-(i+1)*y(i); end y y = 1.0e+017 * 0.000000000000000 -0.000000000000000 0.000000000000000 -0.000000000000000 0.000000000000000 -0.000000000000002 0.000000000000013 -0.000000000000107 0.000000000000959 -0.000000000009589 0.000000000105477 -0.000000001265719 0.000000016454352 -0.000000230360927 0.000003455413904 -0.000055286622459 0.000939872581808 -0.016917706472551 0.321436422978460 -6.428728459569207 y(1)=vpa((1-1/exp(1)),8);for i=1:19 y(i+1)=1.0-(i+1)*y(i); end format long y y = 1.0e+017 * 0.000000000000000 -0.000000000000000 0.000000000000000 -0.000000000000000 0.000000000000000 -0.000000000000002 0.000000000000013 -0.000000000000107 0.000000000000959 -0.000000000009589 0.000000000105477 -0.000000001265719 0.000000016454349 -0.000000230360892 0.000003455413381 -0.000055286614090