很偶然地发现了一个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出现同样问题。
发布者：全栈程序员-站长，转载请注明出处：https://javaforall.net/189898.html原文链接：https://javaforall.net
