24小时热门版块排行榜

返回列表

当前主题已经存档。

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

guoguoxia

铜虫 (初入文坛)

应助: 0 (幼儿园)
金币: 278.9
帖子: 12
在线: 14.3小时
虫号: 570626
注册: 2008-06-07
专业: 土木工程

[交流] 【求助】matlab里两个相等的数相减为什么不等于0呢

自己编了一个小程序，调试时一个本应该为零的结果总是不等于零。认真检查中间的过程，没有发现错误。最后把每一个数拆开发现，运算中两个表达式计算出来的结果完全相等，但是它俩的差值却不等于零。请问这是什么原因造成的。这两个数都是通过矩阵中的元素运算得出来的。另：我将俩数相除，结果等于1。
MATLAB接触不多，还请各位帮忙解答一下，谢谢！

回复此楼

» 猜你喜欢

1楼 2009-09-07 10:41:45

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

guoguoxia

铜虫 (初入文坛)

应助: 0 (幼儿园)
金币: 278.9
帖子: 12
在线: 14.3小时
虫号: 570626
注册: 2008-06-07
专业: 土木工程

程序比较长，数据也很多，截其中一部分吧，具体数据就不写出来了
就是很简单的矩阵运算，求解一个四元方程，得到的解回带回去就不等于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了。

情况就是这样的，非常简单的情况，但是却一直得不到正确结果。谢谢解答！