Question

数值运算的误差（MATLAB实验）？

Answer 1

方法1、计算f(0.01)的真值。其代码

format long g %将输出格式设置为长固定十进制小数点格式

x=0.01;fx=(exp(x)-1-x)/x^2

format

方法2、使用六位有效数字计算

digitsOld = digits(7);

x=0.01;fx0=exp(x);

fx0=(vpa(fx0,7)-vpa(1+x,7))/ vpa(x^2,7);

fx1=vpa(fx0,6)

err=fx-fx1 %误差

%3、使用近似公式计算

x=0.01;fx2=1/2+vpa(x/6,7)+vpa(x^2/24,7);

fx2=vpa(fx2,6)

err=fx-fx2 %误差

从运算结果，我们可以看到方法2比方法3精确高