MATLAB下邻接阵分解

在有关网络图的计算中,经常要计算邻接阵的分解,形如A=QXQT。因为A是实对陈阵,所以QT=Q-1。这样在计算A的power时就可以很方便计算那种超大矩阵了。

简单做了一下实验,首先用pajek打开.net文件,保存成mat的binary格式,然后用matlab的import data将其导入变成邻接阵。

分解采用的是EVD分解,即AV-1=VD,其中V是特征向量矩阵,D是三角阵,调用函数eig即可,如[v,d]=eig(A);

测试了一下发现matlab有个小问题,输入是整数,输出是浮点数,通过isequal函数比较矩阵会出现不等的情况,用round函数限制一下就好了。不过这个可能还有问题,以后再说。

posted @ 2013-04-21 11:16  NEU-Jimmy  阅读(243)  评论(0)    收藏  举报