| 查看: 3050 | 回复: 15 | |||
| 当前主题已经存档。 | |||
[交流]
【求助】matlab里两个相等的数相减为什么不等于0呢
|
|||
|
自己编了一个小程序，调试时一个本应该为零的结果总是不等于零。认真检查中间的过程，没有发现错误。最后把每一个数拆开发现，运算中两个表达式计算出来的结果完全相等，但是它俩的差值却不等于零。请问这是什么原因造成的。这两个数都是通过矩阵中的元素运算得出来的。另：我将俩数相除，结果等于1。 MATLAB接触不多，还请各位帮忙解答一下，谢谢！ |
回复此楼» 猜你喜欢
40+中年人失业后，还能找到好工作吗？
已经有11人回复
-
用NBS溴化时溶剂除了CCl4外还有什么
已经有6人回复
-
申请科研项目纯靠实力能成功吗？
已经有11人回复
-
极限拉扯，如何抉择
已经有15人回复
-
二审4个多月了还是Awaiting Reviewer Scores
已经有4人回复
-
OpenAI又攻克数学难题了，离AI做有机合成还有多远？
已经有13人回复
-
课题组会帮你通过考核么？
已经有7人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有9人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
nono2009
超级版主 (文学泰斗)
No gains, no pains.
-
专家经验: +21105 - 仿真EPI: 1
- 应助: 28684 (院士)
- 贵宾: 513.911
- 金币: 2555230
- 散金: 27828
- 红花: 2148
- 沙发: 66666
- 帖子: 1602255
- 在线: 65200.9小时
- 虫号: 827383
- 注册: 2009-08-13
- 性别: GG
- 专业: 工程热物理与能源利用
- 管辖: 科研家筹备委员会
2楼2009-09-07 11:00:38
sitonyl
铁杆木虫 (文坛精英)
- 应助: 33 (小学生)
- 金币: 12935.9
- 红花: 3
- 帖子: 10448
- 在线: 338.8小时
- 虫号: 722261
- 注册: 2009-03-14
- 性别: GG
- 专业: 导航、制导与传感技术
3楼2009-09-07 11:13:46
|
程序比较长，数据也很多，截其中一部分吧，具体数据就不写出来了 就是很简单的矩阵运算，求解一个四元方程，得到的解回带回去就不等于0. b(1,1)=-U1(1,5); b(2,1)=-U12(1,5); b(3,1)=-U123(1,5); b(4,1)=-U123(3,5); A(1,1)=U1(1,2); A(2,1)=U12(1,2); A(3,1)=U123(1,2); A(4,1)=U123(3,2); A(1,2)=U1(1,4); A(2,2)=U12(1,4); A(3,2)=U123(1,4); A(4,2)=U123(3,4); A(2,3)=U2(1,4); A(3,3)=U3(1,1)*U2(1,4)+U3(1,2)*U2(2,4)+U3(1,3)*U2(3,4)+U3(1,4)*U2(4,4); A(4,3)=U3(3,1)*U2(1,4)+U3(3,2)*U2(2,4)+U3(3,3)*U2(3,4)+U3(3,4)*U2(4,4); A(3,4)=U3(1,4); A(4,4)=U3(3,4); X1=A\b; XL(1,1)=A(3,1)*X1(1,1)+A(3,2)*X1(2,1)+A(3,3)*X1(3,1)+A(3,4)*X1(4,1)-b(3,1) 涉及到的就是这一小部分，都是很简单的运算。 XL(1,1)应该等于0，但结果却不是。 单独计算了A(3,1)*X1(1,1)+A(3,2)*X1(2,1)+A(3,3)*X1(3,1)+A(3,4)*X1(4,1)和b(3,1)，两者结果相等，但是相减后就不等于0了。 情况就是这样的，非常简单的情况，但是却一直得不到正确结果。谢谢解答！ |
4楼2009-09-07 11:14:53
5楼2009-09-07 11:17:12
nono2009
超级版主 (文学泰斗)
No gains, no pains.
-
专家经验: +21105 - 仿真EPI: 1
- 应助: 28684 (院士)
- 贵宾: 513.911
- 金币: 2555230
- 散金: 27828
- 红花: 2148
- 沙发: 66666
- 帖子: 1602255
- 在线: 65200.9小时
- 虫号: 827383
- 注册: 2009-08-13
- 性别: GG
- 专业: 工程热物理与能源利用
- 管辖: 科研家筹备委员会
6楼2009-09-07 14:09:33
aaron1988
木虫 (正式写手)
- 应助: 15 (小学生)
- 金币: 2261.1
- 红花: 3
- 帖子: 938
- 在线: 135.2小时
- 虫号: 413419
- 注册: 2007-06-26
- 专业: 化工系统工程
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
wuguocheng(金币+2,VIP+0): 谢谢 9-8 21:59
小木虫(金币+0.5):给个红包，谢谢回帖交流
wuguocheng(金币+2,VIP+0): 谢谢 9-8 21:59
|
eps问题...(help eps能查到) 因为matlab用的是floating point.....默认的是eps=2.26e-16 也就是说实际上的1跟电脑储存的1的差距 你可以试下1+eps-1等于什么...................................(2.26e-16) 但是2+eps-2等于..................................................(0) 因为对2来说他的eps是eps(2)=2*eps>eps... 当我们计算a+b要是b是小于a的eps的话...我们可以做a+b=a来看(b=0) 就好象我们平时做计算的时候,要是a>>b的话...那么a+b约等于a |
7楼2009-09-08 21:34:23
snipher950
木虫 (正式写手)
- 应助: 2 (幼儿园)
- 金币: 3024.6
- 红花: 3
- 帖子: 797
- 在线: 22.5小时
- 虫号: 522854
- 注册: 2008-03-12
- 性别: MM
- 专业: 化工系统工程
8楼2009-09-08 21:42:19
9楼2009-09-09 09:22:45
10楼2009-10-23 09:23:41
200