数学建模社区-数学中国
标题:
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 0
2 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 0
9 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 0
6 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