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

domenica

用matlab计算行随机矩阵R的谱半径结果为什么大于1?新手求教
>> R=[1/3 0 0 0 1/3 0 0 0 0 1/3;
8 i4 z. f4 A: K0 a: C1 I    1/2 1/2 0 0 0 0 0 0 0 0;
. V+ k1 y/ s1 Y4 c3 W    0 1/2 1/2 0 0 0 0 0 0 0;
    0 0 1/2 1/2 0 0 0 0 0 0;
    0 0 0 1/2 1/2 0 0 0 0 0;
    0 0 0 0 1/2 1/2 0 0 0 0;
    0 0 0 0 0 1/2 1/2 0 0 0;
' ?# B, c  C/ t    0 0 0 0 0 0 1/2 1/2 0 0;
3 v$ o# L0 n! g    0 0 0 0 0 0 0 1/2 1/2 0;
2 R3 P' g1 p* ]+ @: M2 [    1/3 0 0 0 0 0 0 0 1/3 1/3];
  c$ c/ {+ P/ M) y" b% V5 y1 Y>> disp('计算谱半径')
- D3 P/ @: ~+ K. U- O计算谱半径
' p+ K$ D: k0 T1 e2 Q' i8 n>> max(abs(eig(R)))

ans =
+ H5 ~% t8 P, h
   1.000000000000001
% Q# h  X5 p2 G- z

mutougeda

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

regngfpcb

支持楼上，哈哈哈

yin123

help哈哈哈