matlab 单双精度引发的血案

今天在测试一个算法时遇到了少见的问题，觉得值得记录，总结起来，可以用下面几行代码来表示：

看到了么，这个数字在加减去1后都是不变的，然而，当bbb不是单精度时，答案就是正确的：

当然，在交互行中直接赋值的变量默认是double类型。这么写只是为了强调数据的格式。

matlab在用双精度减（或加）单精度类型时，返回的结果竟然存储格式是单精度的，而单精度的实际有效位数为24位，也就是能够表示的范围最大为：
2^24-1 = 16777215

而我的结果明显已经超过了这个范围，所以才会有这么怪异的输出。