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

domenica · 发表于 2020-11-16 15:56

EDA365欢迎您登录！

您需要登录才可以下载或查看，没有帐号？注册

x

用matlab计算行随机矩阵R的谱半径结果为什么大于1?新手求教
>> 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;
1/3 0 0 0 0 0 0 0 1/3 1/3];
>> disp('计算谱半径')
计算谱半径
>> max(abs(eig(R)))

ans =

1.000000000000001

yin123 · 发表于 2020-11-16 17:28

help哈哈哈

regngfpcb · 发表于 2020-11-16 17:28

支持楼上，哈哈哈

mutougeda · 发表于 2020-11-16 17:32

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

帐号		自动登录	找回密码
密码			注册

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

EDA365欢迎您登录！

浏览过的版块

推荐内容 /1