matlab如何做向量内积,matlab2010a的bug：* 向量内积的错误解决方案

最新推荐文章于 2023-08-16 04:33:53 发布

转载

最新推荐文章于 2023-08-16 04:33:53 发布 · 1.1k 阅读

文章标签：

#matlab如何做向量内积

本文揭示了MATLAB 2010a的一个罕见bug，当进行特定条件下的向量内积计算时，*操作符会产生错误结果。通过对比不同版本MATLAB的运行结果，作者发现2010a和2010b在计算向量内积时返回了0，而正确结果应与直接求和一致。提供了一个示例程序，详细介绍了问题的现象，并附带了相关数据文件链接。

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