请帮我看看这个复数矩阵的行列式的值
矩阵M为一个8×8的复数方阵，我用det(M)求出的值为-.17202e28+.80487e-2*i
可是我用行列式的定义的方法求出的值却为
-.1124297554e25+.2345712643e25*i
很奇怪，两者的值差别很大。不知道究竟哪个值才是正确的？另外，我想问一下matlab求复数矩阵的行列式是怎么计算的？和实数一样的处理吗？
以下附矩阵M
Column 1
0 -1.388631098000000e+000i
1.373487414000000e+000 -2.045205351000000e-001i
0 +5.473233054000000e+000i
-5.413544840000000e+000 +8.061093798000000e-001i
0 -1.022401538000000e+002i
1.011251762000000e+002 -1.505814683000000e+001i
4.908771629000000e+008
7.229743031999999e+007 +4.855239136000000e+008i
Column 2
0 +1.388631098000000e+000i
1.373487414000000e+000 +2.045205351000000e-001i
0 -5.473233054000000e+000i
-5.413544840000000e+000 -8.061093798000000e-001i
0 +1.022401538000000e+002i
1.011251762000000e+002 +1.505814683000000e+001i
4.908771629000000e+008
7.229743031999999e+007 -4.855239136000000e+008i
Column 3
0 -8.516758816000000e-002i
-2.501301259000000e-002 -8.141171460000000e-002i
0 -8.580223028000000e-002i
-2.519940171000000e-002 -8.201836913000001e-002i
0 -4.040864364000000e-005i
-1.186768269000000e-005 -3.862663056000000e-005i
8.409365580000000e+005
8.038514244000000e+005 -2.469760758000000e+005i
Column 4
0 +8.516758816000000e-002i
-2.501301259000000e-002 +8.141171460000000e-002i
0 +8.580223028000000e-002i
-2.519940171000000e-002 +8.201836913000001e-002i
0 +4.040864364000000e-005i
-1.186768269000000e-005 +3.862663056000000e-005i
8.409365580000000e+005
8.038514244000000e+005 +2.469760758000000e+005i
Column 5
-3.428909785000000e+000
-2.104301496000000e-005
1.395799629000000e+001
8.565939121000000e-005
3.327269145000000e+002
2.041925241000000e-003
1.061039986000000e+008
6.511539145000000e+002
Column 6
-1.205975479000000e+000
-1.771000737000000e-002
4.751604739000000e+000
6.977833001000000e-002
-3.751887137000000e+001
-5.509726359000000e-001
-2.750086334000000e+008
-4.038560493000000e+006
Column 7
1.205975479000000e+000
8.212175327999999e+001
-4.751604739000000e+000
-3.235638857000000e+002
3.751887137000000e+001
2.554874084000000e+003
-2.750086334000000e+008
-1.872690742000000e+010
Column 8
3.428909785000000e+000
5.587327831000000e+005
-1.395799629000000e+001
-2.274422659000000e+006
-3.327269145000000e+002
-5.421706800000000e+007
1.061039986000000e+008
1.728939696000000e+013
以下是我用来求行列式A的值的程序：
%hanglieshijisuan
ZHENG1=M(1,1)*M(2,2)*M(3,3)*M(4,4)*M(5,5)*M(6,6)*M(7,7)*M(8,8)
ZHENG2=M(1,2)*M(2,3)*M(3,4)*M(4,5)*M(5,6)*M(6,7)*M(7,8)*M(8,1)
ZHENG3=M(1,3)*M(2,4)*M(3,5)*M(4,6)*M(5,7)*M(6,8)*M(7,1)*M(8,2)
ZHENG4=M(1,4)*M(2,5)*M(3,6)*M(4,7)*M(5,8)*M(6,1)*M(7,2)*M(8,3)
ZHENG5=M(1,5)*M(2,6)*M(3,7)*M(4,8)*M(5,1)*M(6,2)*M(7,3)*M(8,4)
ZHENG6=M(1,6)*M(2,7)*M(3,8)*M(4,1)*M(5,2)*M(6,3)*M(7,4)*M(8,5)
ZHENG7=M(1,7)*M(2,8)*M(3,1)*M(4,2)*M(5,3)*M(6,4)*M(7,5)*M(8,6)
ZHENG8=M(1,8)*M(2,1)*M(3,2)*M(4,3)*M(5,4)*M(6,5)*M(7,6)*M(8,7)
ZHENG=ZHENG1+ZHENG2+ZHENG3+ZHENG4+ZHENG5+ZHENG6+ZHENG7+ZHENG8
FU1=M(1,8)*M(2,7)*M(3,6)*M(4,5)*M(5,4)*M(6,3)*M(7,2)*M(8,1)
FU2=M(2,8)*M(3,7)*M(4,6)*M(5,5)*M(6,4)*M(7,3)*M(8,2)*M(1,1)
FU3=M(3,8)*M(4,7)*M(5,6)*M(6,5)*M(7,4)*M(8,3)*M(1,2)*M(2,1)
FU4=M(4,8)*M(5,7)*M(6,6)*M(7,5)*M(8,4)*M(1,3)*M(2,2)*M(3,1)
FU5=M(5,8)*M(6,7)*M(7,6)*M(8,5)*M(1,4)*M(2,3)*M(3,2)*M(4,1)
FU6=M(6,8)*M(7,7)*M(8,6)*M(1,5)*M(2,4)*M(3,3)*M(4,2)*M(5,1)
FU7=M(7,8)*M(8,7)*M(1,6)*M(2,5)*M(3,4)*M(4,3)*M(5,2)*M(6,1)
FU8=M(8,8)*M(1,7)*M(2,6)*M(3,5)*M(4,4)*M(5,3)*M(6,2)*M(7,1)
FU=FU1+FU2+FU3+FU4+FU5+FU6+FU7+FU8,
ZHI=vpa(ZHENG-FU,10)
简单的算例可是我用了一个简单的3×3的复数矩阵A来计算，det（A）的值又能和用行列是定义的方法计算的相同。
为什么8×8的那个却不对，奇怪？
以下是我的算例：
A=[2+i,1+i,3;
2,1,1-i;
4+i,1+2i,1-2i]
AZHENG1=A(1,1)*A(2,2)*A(3,3)
AZHENG2=A(1,2)*A(2,3)*A(3,1)
AZHENG3=A(1,3)*A(2,1)*A(3,2)
AZHENG=AZHENG1+AZHENG2+AZHENG3
AFU1=A(1,3)*A(2,2)*A(3,1）
AFU2=A(2,3)*A(3,2)*A(1,1)
AFU3=A(3,3)*A(1,2)*A(2,1)
AFU=AFU1+AFU2+AFU3
B2=AZHENG-AFU
B1=det（A）
此时B1=B2=-5.000000000000000e+000 +5.000000000000000e+000i
[ 本帖最后由 eight 于 2008-5-13 11:27 编辑 ]
页:
[1]