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

实验一：误差传播与算法稳定性 一：实验内容 考虑一个简单由积分定义的序列： 显然当n=1时，。而对于时，利用分步积分易得： 另一方面，我们有。 由以上递推关系，我们可以得到计算序列的两种方法。 (Ⅰ) ， (Ⅱ) , 二：实验要求及实验结果 分别用算法(Ⅰ)、(Ⅱ)计算，并且在计算机中分别采用5位、6位和7位有效数字，请判断哪种算法能给出更精确的结果。 实验过程： ⅰ）编写MATLAB程序如下： a= input (请输入有效位数a:); %设定有效数字位数 syms n In In=vpa((exp(-1)),a) %vpa设定结果有效数字 for n=2:10; In=vpa((1-n*In),a) %循环计算 End 运行文件，输入有效数字a分别为5位、6位和7位，得到运算结果如下表格所示： 有效位数 5位 6位 7位 I1 .36788 .367879 .3678794 I2 .26424 .264242 .2642412 I3 .20728 .207274 .2072764 I4 .17088 .170904 .1708944 I5 .14560 .145480 .1455280 I6 .12640 .127120 .1268320 I7 .11520 .110160 .1121760 I8 .7840e-1 .118720 .1025920 I9 .29440 -.68480e-1 .766720e-1 I10 -1.9440 1.68480 .2332800 ⅱ）编写MATLAB程序如下： function In=NO1B b= input (请输入有效位数b:); syms n En En=vpa(0,b) for n=10:-1:2; En=vpa(((1-En)/n),b) End 运行文件，输入有效数字a分别为5位、6位和7位，得到运算结果如下表格所示： 有效位数 5位 6位 7位 E11 0. 0. 0. E10 .10000 .100000 .1000000 E9 .10000 .100000 .1000000 E8 .11250 .112500 .1125000 E7 .12679 .126786 .1267857 E6 .14554 .145536 .1455357 E5 .17089 .170893 .1708929 E4 .20728 .207277 .2072768 E3 .26424 .264241 .2642411 E2 .36788 .367880 .3678794 由以上两种算法所得到的数据可知，对算法，从开始，结果变得无规律，各个有效位数计算结果都不一样，这是因为随着计算的n增大，误差会越来越大。而对, ，5位、6位和7位结果相近，随着有效数字位数的增加，结果越来越精确。 两种算法的优劣，与你第一感觉是否吻合。请从理论上证明你的实验得出的结果，解释实验得到的结果，算法(Ⅰ)中的计算误差为，由递推计算的误差为；算法(Ⅱ)中的计算误差为，由向前递推计算（）的误差为。如果在上述两算法中都假定后面的计算不再引入其它误差，试给出与的关系和与的关系。 ⅰ）。 ⅱ）因为由此类推，对nN有：。 算法(Ⅰ)中的会很小，当n增大时，的变化趋势如何？算法(Ⅱ)中通常相对较大，当n减小时，误差又是如何传播的？也就是说比较一下上述两个算法，当某一步产生误差后，该误差对后面的影响是衰减还是扩张的。 ⅰ）算法(Ⅰ)中，。当很小时，随着n的增大，以阶乘为系数迅速增大。所以当某一步产生误差后，该误差最后面的影响是扩张的。 ⅱ）算法(Ⅱ)中，虽然开始很大，但是随着n的增大，以阶乘为除数迅速减小。所以当某一步产生误差后，该误差最后面的影响是衰减的。 通过理论分析与计算实验，针对(Ⅰ)和(Ⅱ)的稳定性给出你的结论。 ⅰ）算法(Ⅰ)中，。当某一步产生误差后，该误差最后面的影响是扩张的，所以该算法是不稳定的。 ⅱ）算法(Ⅱ)中，虽当某一步产生误差后，该误差最后面的影响是衰减的，所以该算法是稳定的。