行随机矩阵计算时谱半径大于1?

domenica

用matlab计算行随机矩阵R的谱半径结果为什么大于1?新手求教
) V) a9 Z0 Q* s& `5 b9 P3 U/ B>> R=[1/3 0 0 0 1/3 0 0 0 0 1/3;
" k9 ~& i$ f2 }    1/2 1/2 0 0 0 0 0 0 0 0;
7 g  g9 B& C2 Y0 i    0 1/2 1/2 0 0 0 0 0 0 0;
, t# G1 E1 W" |7 \6 R' U    0 0 1/2 1/2 0 0 0 0 0 0;
    0 0 0 1/2 1/2 0 0 0 0 0;
) n0 G* G7 v) B& \% k7 B6 ~    0 0 0 0 1/2 1/2 0 0 0 0;
. L+ h  v! m8 a8 K    0 0 0 0 0 1/2 1/2 0 0 0;
    0 0 0 0 0 0 1/2 1/2 0 0;
/ B2 u- B& q' y, N! b    0 0 0 0 0 0 0 1/2 1/2 0;
    1/3 0 0 0 0 0 0 0 1/3 1/3];
4 u& c+ Z. u5 }6 }6 Y>> disp('计算谱半径')
. r) v3 z: j. C' Y* ^8 s$ T- f$ f计算谱半径
>> max(abs(eig(R)))
' Y0 Z7 h0 }. Z( c
ans =
0 E1 m8 W& |4 u- N
  P$ H+ m' K9 G0 r   1.000000000000001
- ]. t( y/ K7 |2 a
2 ]& ]' ]  w* `2 b" V, b8 }+ ~

yin123

help哈哈哈

regngfpcb

支持楼上，哈哈哈

mutougeda

我觉得这应该是由于计算机中都是以二进制形式存储数据，导致产生了浮点误差，这是数值计算的固有问题