数值计算方法”上机实验指导书——实验一 误差分析

实验要求：
**1、**选择充分小的 ess，反复进行上述实验，记录结果的变化并分析它们。如果扰动项的系数 ε 很小，我们自然感觉（1.1）和（1.2）的解应当相差很小。计算中你有什么出乎意料的发现？表明有些解关于如此的扰动敏感性如何？
令ess分别等于0.000000001（1e-8）、0.000000000001（1e-11）和0.000000000000001（1e-14），计算式（1.2）中的解，并使用MATLAB在坐标轴中将各点分别使用绿色、红色和蓝色连接起来。图示如下：

从图中可以看到，随着ess值变小，（1.2）的解的变化幅度也越来越小。
当ess=0.000000000000001（1e-14）时，如图中蓝线所示，此时扰动项的系数ε 很小，式（1.0）和（1.2）的解相差很小。
当解x<=6时，此时无论扰动项的系数ε 怎样变化，式（1.2）的解值不变与（1.1）相同，所有线段均呈直线。说明这些解的扰动敏感性较低。
代码如下：

ess = 0.000000001;
ve = zeros(1,21);
ve(2) = ess;
a=poly(1:20)+ve;
b=roots(poly(1:20)+ve);
ess = 0.000000000001;
ve = zeros(1,21);
ve(2) = ess;
c=poly(1:20)+ve;
d=roots(poly(1: