很偶然地发现了一个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 -4344352