| 查看: 1857 | 回复: 5 | |||
[求助]
matlab计算矩阵顺序问题和精度误差
|
|||
|
执行下面的代码，为什么y1和y2的结果不一样： format long m1 = [13 42;4 13]; m2 = [13 -42; -4 13]; x = [1.055;0]; y = m1*x; y1= m2*y y2 = m2*m1*x x-y1 x-y2 m1和m2矩阵互逆，通过matlab验证为单位阵，且其运算元素的精度并不长，但为什么分开计算（先算y，再算y1）和整体计算（y2）的结果不一样？ 试了几个PC机，都是如此。 y1 = 1.055000000000007e+00 0 y2 = 1.055000000000000e+00 0 ans = -6.883382752675971e-15 0 |
回复此楼» 猜你喜欢
小木虫没落了，除了祈祷帖子，几乎看不到有价值的帖子
已经有15人回复
-
面上的水很深
已经有18人回复
-
国社科通讯评审
已经有5人回复
-
5个%2F，非固定段22个字母
已经有9人回复
-
生命口C13面上会评时间
已经有11人回复
-
内心匮乏
已经有16人回复
-
生命口面上会评结束了吗大家有知道结果吗？现在还没有结果正常不？
已经有9人回复
-
投英文期刊如何查重
已经有6人回复
-
不要再数国自然申请书的 filecode 的分隔符个数了
已经有21人回复
-
伯胺上甲基求助
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 【MATLAB书籍159】《科学计算的介绍：使用MATLAB的矩阵-向量方法》(英文版)
已经有41人回复
- matlab计算数据的邻接矩阵
已经有0人回复
- 如何在matlab中计算对称矩阵的高精度本征函数？
已经有0人回复
- 如何将MATLAB中计算得到的矩阵保存到txt文档中
已经有5人回复
- MATLAB矩阵循环计算问题
已经有5人回复
- 我用matlab求矩阵的特性向量矩阵，怎么出错了？
已经有13人回复
- matlab简单程序问题（矩阵的计算）
已经有1人回复
- MATLAB计算的矩阵超出范围，比如1610*1610的矩阵，如何求解？
已经有5人回复
- matlab中std2函数计算矩阵的标准差时，分母是除以了n-1还是n？
已经有3人回复
- matlab处理矩阵为零的情况
已经有0人回复
- matlab 3维矩阵乘法求组
已经有30人回复
- Matlab矩阵计算
已经有3人回复
- matlab循环计算转为矩阵 30金币求助
已经有7人回复
- matlab计算
已经有2人回复
- 如何用Matlab画矩阵的散点图？
已经有3人回复
- 关于计算机内存和matlab产生矩阵能力的关系的困惑
已经有3人回复
- 用matlab计算刚性系数和柔性系数矩阵是否互逆？
已经有0人回复
- 一个关于用matlab计算矩阵的四阶矩的问题？
已经有0人回复
- 【求助】matlab 解矩阵方程 重金悬赏解决方案
已经有10人回复
yanglongren
新虫 (正式写手)
- 应助: 2 (幼儿园)
- 金币: 5308.7
- 红花: 2
- 帖子: 334
- 在线: 84.6小时
- 虫号: 3982702
- 注册: 2015-07-21
- 性别: GG
- 专业: 应用高分子化学与物理
【答案】应助回帖
cowlemon: 回帖置顶 2017-02-17 16:55:43
月只蓝: 应助指数+1 2017-02-18 17:14:14
月只蓝: 应助指数+1 2017-02-18 17:14:14
|
这是双精度浮点运算的误差吧，matlab中以64位二进制数表示，一旦有小数点的加和乘的运算，在转化为实数时难免会产生误差，若先y=m1*m2，y1=y*x.就不会存在这样的问题。这个误差并不是matlab的bug，C+等函数都有这个问题 发自小木虫Android客户端 |
3楼2017-02-17 16:45:50
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68161.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3292.2小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
2楼2017-02-17 16:33:03
4楼2017-02-17 16:49:39
5楼2017-02-17 16:51:15
匿名
用户注销 (职业作家)
Matlab专家
- 应助: 223 (大学生)
- 金币: 2891.3
- 散金: 5302
- 红花: 59
- 沙发: 10
- 帖子: 3308
- 在线: 945小时
- 虫号: 0
- 注册: 2012-05-01
- 性别: GG
- 专业: 力学
6楼2017-02-18 00:09:17
10