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出现同样问题。

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

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

(0)
全栈程序员-站长的头像全栈程序员-站长


相关推荐

  • squid反向代理

    squid反向代理

    2021年5月28日
    106
  • MongoDB(四)—-MongoDB的文档操作

    MongoDB(四)—-MongoDB的文档操作

    2020年11月12日
    163
  • Java两整数相除向上取整

    Java两整数相除向上取整前言:Java中两个整数相除,如果不能整除,默认是向下取整的。例如:11除以3的结果是3。然而,某些情况下(eg.把11个糖果,每3个分一堆,不足三个也分成一堆,可以分几堆?),我们需要向上取整,这样的情况该如果处理呢?方式一:添加三目运算符逻辑代码x/y+(x%y!=0?1:0);这种方法逻辑上很简单,如果x可以整除y,就将x/y的结果加0,不能整除y就将x/y的结果加1。方式二:使用ceil函数(int)Math.ceil((double.

    2022年6月21日
    107
  • 用EasyBoot轻松做启动光盘

    用EasyBoot轻松做启动光盘原文转自:电脑技术资料园BrianLiu学习之园原版系统安装盘的缺憾——不管是Windows98还是WindowsNT/2000/XP,仅能实现单一系统的初始安装,缺少调试维护、系统恢复、DOS杀毒等工具。虽然市面上出现了N合1光盘,但一般体积庞大,且无法满足自己的需要。  用EasyBoot刻盘正好可以解决这个问题。EasyBoot是一款集成…

    2022年7月26日
    1
  • 爬虫基本知识,如何发起请求,进行分析

    爬虫基本知识,如何发起请求,进行分析爬虫基础知识爬虫一个实战性很强的内容,下面是一些知识点,方便日后复习,具体还要去案例看看,随机应变。这是我的github爬虫仓库,欢迎大家clone进行学习和体验。一.网络爬虫概述定义网络蜘蛛(spider)、网络机器人(robot),抓取网络数据的程序其实就是用Python程序模仿人点击浏览器并访问网站,而且模仿的越像越好,让Web站点无法发现你不是人爬取数据的目的1、公司项目测试数据2、公司业务部门及其他部门所需数据3、数据分析企业获取数据方式1、公司自有数据2、第三方

    2022年10月3日
    0
  • auto.js淘宝秒杀脚本_京东秒杀脚本

    auto.js淘宝秒杀脚本_京东秒杀脚本AUTO.JS脚本实现小米、淘宝、京东抢购代码新的改变功能快捷键合理的创建标题,有助于目录的生成如何改变文本的样式插入链接与图片如何插入一段漂亮的代码片生成一个适合你的列表创建一个表格设定内容居中、居左、居右SmartyPants创建一个自定义列表如何创建一个注脚注释也是必不可少的KaTeX数学公式新的甘特图功能,丰富你的文章UML图表FLowchart流程图导出与导入导出导入代码你好!这是你第一次使用Markdown编辑器所展示的欢迎页。如果你想学习如何使用Markdown编辑器,可以仔细

    2022年5月3日
    294

发表回复

您的邮箱地址不会被公开。 必填项已用 * 标注

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