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

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

[复制链接]

该用户从未签到

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

EDA365欢迎您登录!

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

x
用matlab计算行随机矩阵R的谱半径结果为什么大于1?新手求教
  C( O2 t% |: |" e( ?# V>> R=[1/3 0 0 0 1/3 0 0 0 0 1/3;
, c. Q9 R5 I5 I8 b    1/2 1/2 0 0 0 0 0 0 0 0;3 a( z; j3 h" b; Z" F$ C+ r, `7 a
    0 1/2 1/2 0 0 0 0 0 0 0;
9 Q% C7 B9 D0 m0 ?; g- ^. n    0 0 1/2 1/2 0 0 0 0 0 0;. I2 K# K6 m& A* @5 E3 {
    0 0 0 1/2 1/2 0 0 0 0 0;7 X) q0 i  u$ S+ b
    0 0 0 0 1/2 1/2 0 0 0 0;
0 N- z. \9 G0 y6 E" e$ w& \3 r2 o    0 0 0 0 0 1/2 1/2 0 0 0;
2 a5 |0 m, y& [: c( g& x    0 0 0 0 0 0 1/2 1/2 0 0;
- G3 C0 U/ f' Y* {2 ^  n) Y    0 0 0 0 0 0 0 1/2 1/2 0;6 G/ i4 }2 ^+ o$ }  @3 d# \
    1/3 0 0 0 0 0 0 0 1/3 1/3];, S, p" _$ R/ J) t% [; F1 G, @
>> disp('计算谱半径')( S  `) N% H3 K/ J1 [
计算谱半径$ o( D0 e. M* d' |% l# M) q
>> max(abs(eig(R)))2 v/ n3 u1 {# S# ~' V

( `: i9 j+ y2 G) K6 _4 E  F  `% U* ?ans =
' H: k- v$ U  ^+ Z) A
/ Q/ Q1 o% e$ [   1.000000000000001( B0 v* C* z3 d9 L! o, A  W$ C

0 h/ O/ `0 E2 h3 q" ~6 ~/ }* h
  • 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-6-21 10:23 , Processed in 0.078125 second(s), 23 queries , Gzip On.

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

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

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