matlab如何做向量内积,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 -4344352