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

domenica

用matlab计算行随机矩阵R的谱半径结果为什么大于1?新手求教
: v5 M6 N% H' Z9 ~>> R=[1/3 0 0 0 1/3 0 0 0 0 1/3;
    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 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;
4 W& G6 S& M- E" B& f    1/3 0 0 0 0 0 0 0 1/3 1/3];
>> disp('计算谱半径')
- ?# K$ Q7 H7 d# s2 k2 j) v计算谱半径
+ f+ f0 m/ B7 U! b+ W4 r>> max(abs(eig(R)))

. v+ ?: S5 `: @3 a: \+ w* C' ians =

   1.000000000000001
  P4 g4 _9 |' j9 v2 t! s. w8 K2 d. j

yin123

help哈哈哈

regngfpcb

支持楼上，哈哈哈

mutougeda

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