|
最近开始使用matlab,在写代码过程中遇到了除数为0的情况，在分母位置加了eps就没有报错了，所以查了下eps函数的用法，在这里分享一下，也当作是自己的学习记录！
首先matlab中eps是一个函数，可以返回某一个数N的最小浮点数精度，形式例如eps(N)。eps(a)是|a|与大于|a|的最小的浮点数之间的距离，距离越小表示精度越高。 默认a=1，即eps = eps(1); 我们在matlab中敲入eps和eps(1)可以发现结果是一样的。eps = eps(1) = 2.2204e-16。eps为系统运算时计算机允许取到的最小值。例如对于函数y=cosx/x。由于在编程时分母可能出现为0，所以编程时要写成：y=cosx/(x+eps)。 下来解释下eps: 我们知道浮点数其实是离散的，有限的，而且间隔是不均匀的。我们可以说一个数旁边的数是什么，而它们之间的距离就反应了其精度。越靠近0，数和数之间就越密集，精度就越高。  我画了一个简图，黑色实心点代表的是数轴的数，例如1，红圈与1的距离代表的就是eps(1)，那么在1+eps(1)/2到1之间的所有数，都被认作1，MATLAB无法识别介于1和1+eps(1)/2之间的数。
例一，eps(1)   可以发现eps=eps(1)=2.2204e-16,由于matlab无法区分1到1+eps(1)/2之间的数，所以就将中间的这些数四舍五入到1或者1+eps(1),比如上图中，1+0.3eps(1)就被舍入到1，然后大于0.5时，1+0.7eps(1)就被近似为1+eps(1)。 例二，eps(0)  我们可以看到eps(0)比eps(1)小很多，别的精度和1一样，在0到eps(0)*0.5之间的所有数都是0，eps(0）*0.5以上到eps(0)之间的所有数都被近似为eps(0)。 例三，eps(2)以及eps(N)  我们可以发现最小精度eps(N)随着数量级增大，eps也在逐渐增大，这在计算的过程中，都要注意。防止大数吃小数。
|