数学建模社区-数学中国

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

作者: yn120    时间: 2009-11-16 12:18
标题: MATLAB ecel的错误
在matlab中求下列11阶下三角方阵的行列式出现BUG:正确应该为0,而结果为  -246.不知是什么原因?! \9 k  i3 A. a+ A1 C; a
b=[1           0           0           0           0           0           0           0           0           0           02 A3 k7 M  r/ d
   -70           1           0           0           0           0           0           0           0           0           0( x& j' {, `! K% c2 a; y/ j/ C
   -66         -70           1           0           0           0           0           0           0           0           0
& ^+ ^% E8 v5 B" `2 |, s   -63         -67         -72           1           0           0           0           0           0           0           0( ~/ |# `0 ]  V$ B: R. n8 Q5 Y: v
   -59         -61         -64         -66           1           0           0           0           0           0           0
2 a9 M4 F7 m, K, t   -59         -59         -61         -62         -69           1           0           0           0           0           09 t% ]2 [/ l! \, @; `
   -59         -59         -59         -60         -66         -71           1           0           0           0           0
; @. M5 `# A: v+ {' S   -61         -59         -59         -59         -63         -67         -71           1           0           0           0
$ }: \+ I( u' g! Q; U- y   -64           9         -59         -59         -61         -64         -67         -71           1           0           06 W, m6 K0 E: `$ l5 b2 S) A
   -71         -66         -62         -60         -59         -60         -62         -64         -68           0           0
1 v% k4 Y4 n" E    0         -69         -65           0         -59         -59         -60         -62         -65        1231           1]% J. {! s! L6 B! o$ x) j7 V% O4 q
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