| 查看: 1874 | 回复: 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 |
回复此楼» 猜你喜欢
时间戳变了，能看出什么问题？
已经有3人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有15人回复
-
filecode=后面第一个是大写字母
已经有9人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有80人回复
-
感觉是下周放榜了
已经有6人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有4人回复
-
今天系统多次维护，明天很可能放榜！
已经有9人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 【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人回复
4楼2017-02-17 16:49:39
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68166.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3292.3小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
2楼2017-02-17 16:33:03
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
5楼2017-02-17 16:51:15
25