找回密码
 注册
关于网站域名变更的通知
查看: 387|回复: 3
打印 上一主题 下一主题

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

[复制链接]

该用户从未签到

跳转到指定楼层
1#
发表于 2020-11-16 15:56 | 只看该作者 回帖奖励 |倒序浏览 |阅读模式

EDA365欢迎您登录!

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

x
用matlab计算行随机矩阵R的谱半径结果为什么大于1?新手求教; D7 y/ @7 ?0 i# p4 I+ }  g7 g/ r
>> R=[1/3 0 0 0 1/3 0 0 0 0 1/3;& s2 I6 I# Z4 Y7 {/ ~( A
    1/2 1/2 0 0 0 0 0 0 0 0;
: G+ q6 M; S" y4 [' V3 p1 N    0 1/2 1/2 0 0 0 0 0 0 0;
: A) ]/ q# q3 n    0 0 1/2 1/2 0 0 0 0 0 0;* b/ R8 S& Q$ J6 r: v2 C4 b
    0 0 0 1/2 1/2 0 0 0 0 0;
0 z( A4 s5 A+ q+ W    0 0 0 0 1/2 1/2 0 0 0 0;8 E# Q  F; x( @' p( x1 n
    0 0 0 0 0 1/2 1/2 0 0 0;8 w- g7 Q1 o  H1 @8 _- y  a2 K( S
    0 0 0 0 0 0 1/2 1/2 0 0;% O: V4 P* @" B
    0 0 0 0 0 0 0 1/2 1/2 0;2 T9 _- ?( J  @1 E2 G8 R
    1/3 0 0 0 0 0 0 0 1/3 1/3];: v* c1 q1 w% z, y: O+ w
>> disp('计算谱半径')" \3 h& x* [/ s+ ~3 E
计算谱半径& d3 w% F1 k# o  o8 U
>> max(abs(eig(R)))  O8 j$ w; _1 y, p% a$ v% d; y
+ z" i, K5 i0 l5 R  C: W( c, M3 D
ans =8 o) j8 \+ v3 H* y: T9 u
$ {5 g/ n. f% E( U& f3 e
   1.000000000000001
- d7 s4 n6 y) Q2 o/ Q: D+ C, o$ t+ D! p$ {" i: ]
  • TA的每日心情

    2019-11-29 15:37
  • 签到天数: 1 天

    [LV.1]初来乍到

    2#
    发表于 2020-11-16 17:28 | 只看该作者
    help哈哈哈

    该用户从未签到

    3#
    发表于 2020-11-16 17:28 | 只看该作者
    支持楼上,哈哈哈
  • TA的每日心情

    2019-11-20 15:22
  • 签到天数: 2 天

    [LV.1]初来乍到

    4#
    发表于 2020-11-16 17:32 | 只看该作者
    我觉得这应该是由于计算机中都是以二进制形式存储数据,导致产生了浮点误差,这是数值计算的固有问题
    您需要登录后才可以回帖 登录 | 注册

    本版积分规则

    关闭

    推荐内容上一条 /1 下一条

    EDA365公众号

    关于我们|手机版|EDA365电子论坛网 ( 粤ICP备18020198号-1 )

    GMT+8, 2025-11-1 11:18 , Processed in 0.140625 second(s), 23 queries , Gzip On.

    深圳市墨知创新科技有限公司

    地址:深圳市南山区科技生态园2栋A座805 电话:19926409050

    快速回复 返回顶部 返回列表