x*z
发帖数: 381 | 1 我正在用matlab的lmisolver解线性矩阵不等式，碰到一个问题如下：
我需要求解A(X1,X2)<0, where X1 and X2 are unknown matrices,
所以我用lmisolver已经得到一个数值解X1_0和X2_0，为了验证结果，
我计算矩阵A(X1_0,X2_0)的eigenvalues,发现一些在-3E-6附近，
现在我老板的担心是，这些特征值这么小，可能是由于Matlab的计算误差使得
所有特征值都满足小于0的要求。所以希望我得到一些别的数值解使得A矩阵的
特征值离y轴更远点。
我想请教大家，我老板的担心是多余的吗？难道Matlab连1E-6这种计算精度都达不到？
谢谢。 | r****y
发帖数: 1437 | 2 the precision is in relative sense.
default matlab is double precision, so ~1e-14
for your case, take the ratio of the largest eigenvalue (in
magnitude)
and the smallest eigenvalue, if the ration is > 1e10, I would suggest to use
it with caution.
【在 x*z 的大作中提到】
: 我正在用matlab的lmisolver解线性矩阵不等式，碰到一个问题如下：
: 我需要求解A(X1,X2)<0, where X1 and X2 are unknown matrices,
: 所以我用lmisolver已经得到一个数值解X1_0和X2_0，为了验证结果，
: 我计算矩阵A(X1_0,X2_0)的eigenvalues,发现一些在-3E-6附近，
: 现在我老板的担心是，这些特征值这么小，可能是由于Matlab的计算误差使得
: 所有特征值都满足小于0的要求。所以希望我得到一些别的数值解使得A矩阵的
: 特征值离y轴更远点。
: 我想请教大家，我老板的担心是多余的吗？难道Matlab连1E-6这种计算精度都达不到？
: 谢谢。
| x*z
发帖数: 381 | 3 thanks!
use
【在 r****y 的大作中提到】
: the precision is in relative sense.
: default matlab is double precision, so ~1e-14
: for your case, take the ratio of the largest eigenvalue (in
: magnitude)
: and the smallest eigenvalue, if the ration is > 1e10, I would suggest to use
: it with caution.
| f**l
发帖数: 2041 | 4 matlab遵从的是IEEE标准. eps = 2^-52 .= 1e-16.
use
【在 r****y 的大作中提到】
: the precision is in relative sense.
: default matlab is double precision, so ~1e-14
: for your case, take the ratio of the largest eigenvalue (in
: magnitude)
: and the smallest eigenvalue, if the ration is > 1e10, I would suggest to use
: it with caution.
|
|
