数学建模社区-数学中国
标题:
MATLAB ecel的错误
[打印本页]
作者:
yn120
时间:
2009-11-16 12:18
标题:
MATLAB ecel的错误
在matlab中求下列11阶下三角方阵的行列式出现BUG:正确应该为0,而结果为 -246.不知是什么原因?
* E8 T# U2 ?3 S/ O; p4 y: J7 o( N) a
b=[1 0 0 0 0 0 0 0 0 0 0
% y# D: t' V4 N' m0 D0 a" G4 U) a3 _
-70 1 0 0 0 0 0 0 0 0 0
0 P, ^7 z6 a5 f/ j* `1 M" x
-66 -70 1 0 0 0 0 0 0 0 0
! u6 j8 I1 @" N' M
-63 -67 -72 1 0 0 0 0 0 0 0
2 l v, }. R( @# F! t: X. F% R
-59 -61 -64 -66 1 0 0 0 0 0 0
0 u: P+ r/ y- v \
-59 -59 -61 -62 -69 1 0 0 0 0 0
3 P1 l4 P8 S- T* P6 ]0 ~; o
-59 -59 -59 -60 -66 -71 1 0 0 0 0
( m& }4 N, g' b0 o
-61 -59 -59 -59 -63 -67 -71 1 0 0 0
3 Z) n' f! \, y# \
-64 9 -59 -59 -61 -64 -67 -71 1 0 0
1 D2 z; V$ O6 v8 k) \# h
-71 -66 -62 -60 -59 -60 -62 -64 -68 0 0
& f: c: ]$ V. E/ d) n9 U2 {
0 -69 -65 0 -59 -59 -60 -62 -65 1231 1]
% ^$ w; f( j8 \& K6 u6 z
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