matlab2010a的bug：* 向量内积的错误

 

很偶然地发现了一个matlab2010a的一个bug：在某种非常特殊的情况下，matlab2010a的向量内积，即*指令会产生错误的结果。

下面陈述具体现象：

1.程序源代码

（注：运行该程序时需要把den1.mat和kern1.mat这两个matlab数据文件与程序放在同一个文件夹中）

clc;
clear;
load('den1.mat');
load('kern1.mat');
x = kern;
y = den(:);
z = 0;
for i = 1:length(y)
    z = z + x(1,i)*y(i,1);
end
z1 = x*y;
z2 = sum(y);

2.程序相关说明

（1）den1.mat是matlab数据文件（可从此处下载http://ishare.iask.sina.com.cn/f/22182234.html），其中数据是1*36维行向量，如下(其中的 - 是负号)：

[1 -34.9104720617963 591.960057485552 -6494.90695093987 51826.3768981548 -320501.857479884 1598411.70513468 -6605060.33215733 23058608.5581138 -68998984.4665517 178938782.362431 -405639572.104363 809205656.532103 -1428012719.05099 2238285526.53077 -3125592609.96192 3897008179.02336 -4344352067.53865 4333253669.23331 -3867217670.62412 3085871743.18171 -2198564376.37758 1395513287.06250 -786754091.156305 392372548.725918 -172203104.526507 66062863.9131647 -21964736.3140818 6259618.77763060 -1507085.68697902 300647.808748368 -48367.8227891117 6030.54857993965 -546.832712826108 32.0845535965159 -0.914362645545420]

（2）kern1.mat也是matlab数据文件（可从此处下载http://ishare.iask.sina.com.cn/f/22182235.html），其中数据是1*36维行向量，如下：

[1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1]
（3）den和kern分别是来自于den1.mat和kern1.mat的1*36维行向量，其中kern的所有元素都为1。

（4）很显然这个程序的作用是用三种程序不同的程序方法计算den中所有元素的和，即：

     z1使用向量内积计算，z2直接求元素和，z是用循环的方法求向量的内积。

3.bug现象描述

（1）在matlab7（Version 7.0.0.19920(R14)）下运行该程序所得z=z1=z2=5.3111e-008的完全相同，这符合我们的数学常识。

     但是令人奇怪的是，在matlab2010a中出现了奇怪的结果。

（2）在matlab210a中，z=z2=5.311088491222193e-08，z1=0.  

     是不是很奇怪！难道这是matlab2010a的bug！

（3）matlab210b与matlab210a出现同样问题。