matlab如何做向量内积,matlab2010a的bug:* 向量内积的错误解决方案「建议收藏」

matlab如何做向量内积,matlab2010a的bug:* 向量内积的错误解决方案「建议收藏」很偶然地发现了一个matlab2010a的一个bug:在某种非常特殊的情况下,matlab2010a的向量内积,即*指令会产生错误的结果。下面陈述具体现象:1.程序源代码(注:运行该程序时需要把den1.mat和kern1.mat这两个matlab数据文件与程序放在同一个文件夹中)clc;clear;load(‘den1.mat’);load(‘kern1.mat’);x=kern;y=d…

大家好,又见面了,我是你们的朋友全栈君。如果您正在找激活码,请点击查看最新教程,关注关注公众号 “全栈程序员社区” 获取激活教程,可能之前旧版本教程已经失效.最新Idea2022.1教程亲测有效,一键激活。

Jetbrains全系列IDE稳定放心使用

很偶然地发现了一个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出现同样问题。

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 举报,一经查实,本站将立刻删除。

发布者:全栈程序员-用户IM,转载请注明出处:https://javaforall.cn/189898.html原文链接:https://javaforall.cn

【正版授权,激活自己账号】: Jetbrains全家桶Ide使用,1年售后保障,每天仅需1毛

【官方授权 正版激活】: 官方授权 正版激活 支持Jetbrains家族下所有IDE 使用个人JB账号...

(0)


相关推荐

  • 常用的Map集合

    常用的Map集合Map集合Map集合常用方法增加删除判断获取Map集合中TreeSet的应用Map集合中HashSet的应用Map集合常用方法Map集合和Collection集合的区别Map集合是有Key和Value的,Collection集合是只有Value。Collection集合底层也是有Key和Value,只是隐藏起来。1.增加put(Kkey,Vvalue)2.删除…

  • python编写nc的思考

    0x00前言发现自己学习python已经有半个月了,也开发了自己的一些渗透的小脚本,但觉得还是不够,我个人觉得工具和脚本还有框架是个本质上的区别。脚本的话,不会考虑到其他的一些因素,例如报错和交互

    2021年12月11日
  • 408计算机考研

    408计算机考研为什么要考研?考研要先了动机,没有动机很难在枯燥且艰辛漫长的考研之路坚持到最后。有的人可能说,计算机薪资高,是钱让我充满了斗志!对此云泽只能说,真是庸俗且真实。可是你知道那掉头发的速度和你赚钱的速度是成正比的吗考研要先了动机,没有动机很难在枯燥且艰辛漫长的考研之路坚持到最后。有的人可能说,计算机薪资高,是钱让我充满了斗志!大家都知道,计算机本科就很好就业,所以我们面临的就是读研深造自己和就业的诱惑之间摇摆。我们目前主要就是数理,读研和就业的问题。读研能带给你什么东

  • 事务隔离级别与锁的对应关系_X1800S

    事务隔离级别与锁的对应关系_X1800S隔离级别org.springframework.transaction.annotation.Isolationpublic enum Isolation { DEFAULT(-1), READ_UNCOMMITTED(1), READ_COMMITTED(2), REPEATABLE_READ(4), SERIALIZABLE(8);}DEFAULT :这是默认值,表示使用底层数据库的默认隔离级别。对大部分数据库而言,通常这值就是: READ_CO..

  • JavaScript正则表达式简单教程「建议收藏」

    JavaScript正则表达式简单教程「建议收藏」1.常见的正则表达式符号?.匹配除换行符以外的任意字符\w匹配字母或数字或下划线或汉字\s匹配任意的空白符\d匹配数字\b匹配单词的开始和结束^匹配字符串的开始$匹配字符串的结束*重复零次或更多次+重复一次或更多次?重复零次或一次{n}重复n次{n,}重复n多次{n,m}重复n到m词\W匹配任意不是字母,数字,下划线,汉字的字符\S匹配任意不是空白符的字符\D匹配任意非数字的字符\B匹配不是单词开头或结尾的位置【^x】匹配除了x

  • putty 报server unexpectedly closed network connection错误

    putty 报server unexpectedly closed network connection错误由于IP变了,再用putty访问组内的服务器时,竟然总是出错,报serverunexpectedlyclosednetworkconnection错误。查了下防火墙设置,原来里面指定了allow的IP端,真的是晕死了,改好,一试,果然OK了

    2022年10月20日

发表回复

您的电子邮箱地址不会被公开。

关注全栈程序员社区公众号