import numpy as np
import networkx as nx
import pylab as plt
from acipy.sparse.linalg import eigs
L=[(1,2),(2,3),(3,4),(3,5),(3,6),
(4,1),(5,6),(6,1)]
G=nx.DiGraph()
G.add_nodes_from(range(1,7))
G.add_edges_from(L)
B=np.array(np.to_numpy_matrix(G))
plt.rc(‘font’,size=16); pos=nx.shell_layout(G)
nx.draw(G,pos,node_size=280,font_weight=‘bold’
node_color=‘r’,with_labels=True)
plt.savefig(“figure X.png”)
A=B/np.tile(B.sum(axis=1,keepdims=True),(1,B,shape[1]))
A=0.15/B.shape+0.85*A
print(“A=”,A)
W,V=eigs(A.T,1); V=V.real
V=V.flatten();
V=V/V.sum(),print(“V=”,V); plt.figure(2)
plt.bar(range(1,B.shape=[0]+1),V, width=0.6, color=‘b’)
plt.savefig(“figure x_2”);plt.show()
2021-06-02
最新推荐文章于 2022-06-16 19:04:49 发布