拉普拉斯矩阵——复杂链路预测

L=D-A

A是一个图的矩阵表示

A=[[0,1,0,0,1,0]
,[1,0,1,0,1,0]
,[0,1,0,1,0,0]
,[0,0,1,0,1,1]
,[1,1,0,1,0,0]
,[0,0,0,1,0,0]]
D是矩阵每个节点的度

[[2. 0. 0. 0. 0. 0.]
[0. 3. 0. 0. 0. 0.]
[0. 0. 2. 0. 0. 0.]
[0. 0. 0. 3. 0. 0.]
[0. 0. 0. 0. 3. 0.]
[0. 0. 0. 0. 0. 1.]]

L=D-A

[[ 2. -1. 0. 0. -1. 0.]
[-1. 3. -1. 0. -1. 0.]
[ 0. -1. 2. -1. 0. 0.]
[ 0. 0. -1. 3. -1. -1.]
[-1. -1. 0. -1. 3. 0.]
[ 0. 0. 0. -1. 0. 1.]]

python 求逆矩阵
Ln=np.linalg.pinv(L)

n(x,y)=M(Ixx+Iyy-2*Ixy)
Ixy为所求逆矩阵Ln对应的元素
M为常数可以忽略
n(x,y)为x->y 的平均通勤时间

ACT相似性为
Sxy_act=1/(Ixx+Iyy-2*Ixy)