数学建模社区-数学中国
标题:
MATLAB ecel的错误
[打印本页]
作者:
yn120
时间:
2009-11-16 12:18
标题:
MATLAB ecel的错误
在matlab中求下列11阶下三角方阵的行列式出现BUG:正确应该为0,而结果为 -246.不知是什么原因?
' O) a1 h1 s- H! u3 h
b=[1 0 0 0 0 0 0 0 0 0 0
* t* L" m5 S. u! m2 }" O7 e/ w
-70 1 0 0 0 0 0 0 0 0 0
+ j( D( M+ {& n
-66 -70 1 0 0 0 0 0 0 0 0
" [' z- O/ q$ v2 W$ S; W5 k
-63 -67 -72 1 0 0 0 0 0 0 0
7 R" C2 V5 U( v& [ r0 L6 V0 m
-59 -61 -64 -66 1 0 0 0 0 0 0
9 b7 D( S6 ~" F" z
-59 -59 -61 -62 -69 1 0 0 0 0 0
! z# a \" r& [# n7 B* i
-59 -59 -59 -60 -66 -71 1 0 0 0 0
/ j# I! D. b K( M
-61 -59 -59 -59 -63 -67 -71 1 0 0 0
: \( _4 D" D$ j j6 x
-64 9 -59 -59 -61 -64 -67 -71 1 0 0
- J+ p# u4 @- Y3 e8 s& R8 K
-71 -66 -62 -60 -59 -60 -62 -64 -68 0 0
( c% |6 Y/ w/ E6 g1 d
0 -69 -65 0 -59 -59 -60 -62 -65 1231 1]
; w1 b& W, r6 v
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