数学建模社区-数学中国
标题:
MATLAB ecel的错误
[打印本页]
作者:
yn120
时间:
2009-11-16 12:18
标题:
MATLAB ecel的错误
在matlab中求下列11阶下三角方阵的行列式出现BUG:正确应该为0,而结果为 -246.不知是什么原因?
, d& I& L; Z4 l7 i
b=[1 0 0 0 0 0 0 0 0 0 0
3 B3 v- ~2 G- U
-70 1 0 0 0 0 0 0 0 0 0
0 B( c: P+ P1 V& J4 f3 D, P( p
-66 -70 1 0 0 0 0 0 0 0 0
* T1 e |/ {/ {
-63 -67 -72 1 0 0 0 0 0 0 0
[3 v1 z; \+ M6 D9 a7 E$ _& y
-59 -61 -64 -66 1 0 0 0 0 0 0
& Y- Q. @8 Z; O3 i; B
-59 -59 -61 -62 -69 1 0 0 0 0 0
) k, X# h0 J0 S; o9 c" S" T
-59 -59 -59 -60 -66 -71 1 0 0 0 0
I j3 p0 a% b1 g3 H" J
-61 -59 -59 -59 -63 -67 -71 1 0 0 0
& Y5 L! P: u' E7 c* i- H8 k
-64 9 -59 -59 -61 -64 -67 -71 1 0 0
! k0 T; D- H3 m' c
-71 -66 -62 -60 -59 -60 -62 -64 -68 0 0
0 `5 M9 P/ n6 I( p
0 -69 -65 0 -59 -59 -60 -62 -65 1231 1]
+ _& Z! V( ?9 X9 `2 F9 O( j- i0 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