TA的每日心情 | 难过 2019-11-20 15:02 |
|---|
签到天数: 1 天 [LV.1]初来乍到
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
RT,我在求一个超大矩阵(10021*10021的稀疏图的邻接矩阵)的特征值和特征向量遇到如下问题:
* F9 C) N, q6 k& w, w' d9 L. J(1)特征值和特征向量计算时间过长,虽然可以运行出来,但是需要处理大约50分钟,这里使用matlab自带的eig函数运算,请问有无高效的求全部特征值所组成的对角阵和全部特征向量所组成的矩阵的运算方法?+ o; D6 k1 K' G+ u S+ @
9 w! \# Q" A- [5 m X8 P9 a: d( r
(2)在求得特征向量矩阵V以后,我将其求逆后与一个列向量X(10021*1)相乘,得到的是一个类奇异矩阵Y【Y=V^(-1)*X】,无法做到V与Y相乘来还原X【X=V*Y,算式应成立,但是Y奇异matlab无法运算】,请问我应该如何操作可以还原X矩阵?
6 F! y9 o; s. R, q- ]" i6 y# p1 r. ?7 c
0 [2 F1 Q& @6 ] t1 A" h8 k J7 y" S8 `) ]7 M
V^(-1)*X;" N" ~# C: a5 _2 a4 t" {# @& Y( n
2 X# j& V, C" f7 x |1 ~, J
警告: 矩阵为奇异工作精度。0 s- Y8 V: S) p/ x
, {5 h3 c- x3 \& I5 M2 B; c% c0 P! SIn matlab.internal.math.mpower.viaMtimes (line 35)
D( y2 r6 p+ A# f8 H$ f3 u+ @, W% W
1 j2 |: Y) B# m$ [- }5 H# K! U
9 J; n7 j5 Q. a, e2 F. Z
03.25更新:问题二V压根就是一个不可逆矩阵,描述有误,抱歉。/ F7 t8 Z8 d) m
使用pinv求逆时候svd内存不足,有方法解决吗?# P2 Q1 ]2 P: ~
* S8 z9 j t; T/ J, t* _! x
2 q- L( [5 q. ]; m' F4 l7 g2 ?
Y/ q4 \9 K0 W/ i: p7 U x; [1 ?
" q$ h& j" E4 o L4 Q' p; q7 s% r
|
|