EDA365电子论坛网

标题: 行随机矩阵计算时谱半径大于1? [打印本页]

---

作者: domenica    时间: 2020-11-16 15:56
标题: 行随机矩阵计算时谱半径大于1?
用matlab计算行随机矩阵R的谱半径结果为什么大于1?新手求教
9 y1 c$ S3 ^" q$ J2 P$ O' P>> 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;
" X% j: D. O+ y- H# [* I    0 1/2 1/2 0 0 0 0 0 0 0;
    0 0 1/2 1/2 0 0 0 0 0 0;
& `' w$ b5 f  X( o    0 0 0 1/2 1/2 0 0 0 0 0;
# ?5 C( E2 U( U6 z' v* C2 |    0 0 0 0 1/2 1/2 0 0 0 0;
# R4 U! [" m' G. G( H  R    0 0 0 0 0 1/2 1/2 0 0 0;
% b' _! @) j+ L8 T    0 0 0 0 0 0 1/2 1/2 0 0;
" u+ V6 |' G6 d2 M! \# Q9 D& m    0 0 0 0 0 0 0 1/2 1/2 0;
2 ?0 Z- k5 g& v2 f; n+ f; @) \    1/3 0 0 0 0 0 0 0 1/3 1/3];
% M  p$ i* Z! d3 k$ h' o# E* N( o>> disp('计算谱半径')
& M' `! r# C1 F; q7 X1 o, K( w计算谱半径
>> max(abs(eig(R)))

  E$ S: I; z( u0 Q! z- ~ans =
) f8 M4 X) A1 S  i% o9 I3 v
   1.000000000000001
9 W( L& g* d7 k; |
2 ?3 Y# s) d3 r0 s3 u0 M

---

作者: yin123    时间: 2020-11-16 17:28
help哈哈哈

---

作者: regngfpcb    时间: 2020-11-16 17:28
支持楼上，哈哈哈

---

作者: mutougeda    时间: 2020-11-16 17:32
我觉得这应该是由于计算机中都是以二进制形式存储数据，导致产生了浮点误差，这是数值计算的固有问题

---

欢迎光临 EDA365电子论坛网 (https://bbs.eda365.com/)

Powered by Discuz! X3.2