b********e
发帖数: 109 | 1 我正在写一个大的计算程序，其中涉及到4*4矩阵b的行列式值计算。发现计算结果为0，
而实际上值很大（如下所示，在Matlab中单独计算的结果）。 计算Condition number发
现不小，但为什么程序中计算不出正确的结果，而单独算却可以呢？
b=[
-0.10502-0.052579i, -0.12305 + 0.069634i, 0 , 0 ;
0.00349+0.261i, 0.303+0.0028i, -2.98e+013 + 1.39e+013i, -1.39e+013
-2.98e+013i;
-0.3099+0.0257i, 0.023+0.263i,2.72e+013+5.85e+013i,-5.8558e+013 +
2.7219e+013i;
0 , 0 , -7.4743e+013 - 1.6081e+013i, 1.6081e+013 - 7.4743e+013i]
ans =
-1.8097e+010 +2.5265e+010i
ans calculated in the my compu | c*******e
发帖数: 8624 | 2 大才要出问题嘛.
【在 b********e 的大作中提到】
: 我正在写一个大的计算程序，其中涉及到4*4矩阵b的行列式值计算。发现计算结果为0，
: 而实际上值很大（如下所示，在Matlab中单独计算的结果）。 计算Condition number发
: 现不小，但为什么程序中计算不出正确的结果，而单独算却可以呢？
: b=[
: -0.10502-0.052579i, -0.12305 + 0.069634i, 0 , 0 ;
: 0.00349+0.261i, 0.303+0.0028i, -2.98e+013 + 1.39e+013i, -1.39e+013
: -2.98e+013i;
: -0.3099+0.0257i, 0.023+0.263i,2.72e+013+5.85e+013i,-5.8558e+013 +
: 2.7219e+013i;
: 0 , 0 , -7.4743e+013 - 1.6081e+013i, 1.6081e+013 - 7.4743e+013i]
| r****y
发帖数: 1437 | 3
condition should be less than 1e14, if you wanna a reliable calculation.
【在 b********e 的大作中提到】
: 我正在写一个大的计算程序，其中涉及到4*4矩阵b的行列式值计算。发现计算结果为0，
: 而实际上值很大（如下所示，在Matlab中单独计算的结果）。 计算Condition number发
: 现不小，但为什么程序中计算不出正确的结果，而单独算却可以呢？
: b=[
: -0.10502-0.052579i, -0.12305 + 0.069634i, 0 , 0 ;
: 0.00349+0.261i, 0.303+0.0028i, -2.98e+013 + 1.39e+013i, -1.39e+013
: -2.98e+013i;
: -0.3099+0.0257i, 0.023+0.263i,2.72e+013+5.85e+013i,-5.8558e+013 +
: 2.7219e+013i;
: 0 , 0 , -7.4743e+013 - 1.6081e+013i, 1.6081e+013 - 7.4743e+013i]
|
|
