实验一误差分析+刘子强+201930180234.docxVIP

PAGE 3 PAGE 数值分析实验报告 实验名称 误差传播与算法稳定性 实验时间 2021年 10月 9日 姓名 刘子强 班级 数学192 学号 201930180234 成绩 一、 实验目的 1．理解数值计算稳定性的概念。 2．了解数值计算方法的必要性。 3．体会数值计算的收敛性与收敛速度。 二、 实验内容 计算， 三、算法描述 由 ，知 则 可得递推关系 1． ， 2． ， 下面分别以1，2递推关系求解 算法1 ， 当时，㏑㏑1.2 , 递推公式为 = 1 \* GB2 （1） 算法2 ， 当时， 则 即 取递推初值 递推公式为 = 2 \* GB2 （2） 取递推公式 = 1 \* GB2 （1）中的初值 ，得 取递推公式 = 2 \* GB2 （2）中的初值 ，得 四、程序流程图 算法1 I0=0.182322 for n=1:20 I0=1/n-5*I0; disp([I,num2str(n),=,num2str(I0)]) end 算法2 I20=0.008730 for n=20:-1:1 I20=1/5*(1/n-I20); disp([I,num2str(n),=,num2str(I20)]) end 五、实验结果 计算结果如下： n 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 0.182322 0.088390 0.058050 0.043083 0.034583 0.027083 0.031250 -0.013393 0.191960 -0.848710 4.343600 -21.626800 108.217600 -541.011000 2705.126300 -13525.564600 67627.885500 -338139.368500 1690696.898300 -8453484.438700243500 0.182322 0.182320 0.088392 0.058039 0.043139 0.034306 0.028468 0.024325 0.021233 0.018837 0.016926 0.015368 0.014071 0.012977 0.012040 0.011229 0.010520 0.009898 0.009336 0.008876 0.008254 六、实验结果分析 由递推公式 = 1 \* GB2 （1）知当时，应当为精确解，递推公式的每一步都没有误差的取舍，但计算结果＞ ，出现负值，且之后随着n的增大，误差急剧增大。由此看出，当较大时，用递推公式 = 1 \* GB2 （1）中的近似是不正确的。 主要原因是初值不是精确值，设有误差，由递推公式 = 1 \* GB2 （1）知 则有 误差随的增大而迅速增加，增加到的倍。由此可见，递推公式计算的误差不仅取决于初值的误差，公式的精确性，还依赖于误差的传递即递推计算的稳定性。 由递推公式 = 2 \* GB2 （2）知 ，为估计值，并不精确，有，而由 得 误差随递推公式逐步缩小。 综上所述，在递推计算中，数值计算方法是非常重要的，误差估计、误差传播及递推计算的稳定性都会直接影响递推结果。 教 师 评 语 指导教师： 年 月 日 数值分析上机实验原始记录 实验名称：误差传播与算法稳定性 实验时间： 2021 年 10 月 9 日 姓名： 刘子强 学号： 201930180234 班级： 数学192 format long shiyan1 I0 = 0.182322000000000 I1=0.08839 I2=0.05805 I3=0.043083 I4=0.034583 I5=0.027083 I6=0.03125 I7=-0.013393 I8