数学建模社区-数学中国
标题:
MATLAB ecel的错误
[打印本页]
作者:
yn120
时间:
2009-11-16 12:18
标题:
MATLAB ecel的错误
在matlab中求下列11阶下三角方阵的行列式出现BUG:正确应该为0,而结果为 -246.不知是什么原因?
4 J! N( U+ D4 o
b=[1 0 0 0 0 0 0 0 0 0 0
7 X/ u/ ?+ A$ P1 D0 [- E/ f
-70 1 0 0 0 0 0 0 0 0 0
$ g3 {5 K; Q! A' x2 D) ^
-66 -70 1 0 0 0 0 0 0 0 0
; _: k3 f6 L2 z4 J( d( W3 F% n; T( w
-63 -67 -72 1 0 0 0 0 0 0 0
) {& N7 M H) @: S# |4 ]( y3 R. E
-59 -61 -64 -66 1 0 0 0 0 0 0
) _) m2 o. Q) s$ e r8 q# `
-59 -59 -61 -62 -69 1 0 0 0 0 0
( l, f9 D& W; O# y; f( _$ X
-59 -59 -59 -60 -66 -71 1 0 0 0 0
# N1 l: y5 J ^. ^! w- b2 ]
-61 -59 -59 -59 -63 -67 -71 1 0 0 0
7 s! a+ s7 U' O
-64 9 -59 -59 -61 -64 -67 -71 1 0 0
% p- g$ ?3 K9 ?3 k* D! y4 H
-71 -66 -62 -60 -59 -60 -62 -64 -68 0 0
5 C( s" f: K$ C( g. S: d# t. n" I" Y
0 -69 -65 0 -59 -59 -60 -62 -65 1231 1]
8 S) V2 _, O* ?' ]
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