数学建模社区-数学中国

标题: MATLAB ecel的错误 [打印本页]

作者: yn120    时间: 2009-11-16 12:18
标题: MATLAB ecel的错误
在matlab中求下列11阶下三角方阵的行列式出现BUG:正确应该为0,而结果为  -246.不知是什么原因?
) a: m' o6 Q/ {+ W+ X& m" T3 Kb=[1           0           0           0           0           0           0           0           0           0           0
- i2 t, ~+ W+ a1 D% l   -70           1           0           0           0           0           0           0           0           0           0
& |( ~! @, A5 }   -66         -70           1           0           0           0           0           0           0           0           0
! X' v1 Y: L/ v   -63         -67         -72           1           0           0           0           0           0           0           0
$ B" \/ N& b5 ~; M: f   -59         -61         -64         -66           1           0           0           0           0           0           0
2 e1 D* f; S& J4 n7 J( K' E) z8 F& N   -59         -59         -61         -62         -69           1           0           0           0           0           0
# p6 D0 [+ O$ w- z   -59         -59         -59         -60         -66         -71           1           0           0           0           0
9 ^. S: w9 j/ b0 c3 F2 P   -61         -59         -59         -59         -63         -67         -71           1           0           0           0
' V4 D0 a) \: S2 R   -64           9         -59         -59         -61         -64         -67         -71           1           0           0, a. w7 Y/ ~' V8 s
   -71         -66         -62         -60         -59         -60         -62         -64         -68           0           0
  _% k* c3 B; \6 ]    0         -69         -65           0         -59         -59         -60         -62         -65        1231           1]$ T% P: e5 a6 k
d=det(b)
作者: buct    时间: 2009-11-16 17:33
可能是病态矩阵 maltab不一定是对的 问问学数学的专业人士
作者: 厚积薄发    时间: 2009-11-19 18:00
我测试了一下,问题出在b(11,10)=1231上,当适当的改变这个数时候,det(A)可能等于0
作者: yn120    时间: 2009-11-19 22:43
看来,matlab软件,包含excel(也有同样错误):在数值计算中含有BUG,我不知道将来好能不能相信这些应用软件,是不是开发者为防止盗版故意留下BUG。这么简单的行列式运算都会出现问题。




欢迎光临 数学建模社区-数学中国 (http://www.madio.net/) Powered by Discuz! X2.5