def node_sst_test(vno):
NP=50
#计算压力灵敏度(2维数组,即数组的数组)
Dp=node_sst(vno,nCount=10,accuracy=0.01)
#print Dp
np.savetxt('Dp.txt',Dp)
whiten_Dp = whiten(np.abs(Dp))
#print whiten_Dp
np.savetxt('whiten_Dp.txt',whiten_Dp)
#分类后计算聚类重心的索引位置
print '--------压力灵敏度--------------'
centroids=fca(whiten_Dp,n)
print n,'-->',[vno.nodes[i+1].id for i in centroids]
#---------相对压力灵敏度--------------------#
ZDp=sst(Dp)
np.savetxt('ZDp.txt',ZDp)
#数据标准化(白化)
whiten_ZDp = whiten(ZDp)
np.savetxt('whiten_ZDp.txt',whiten_ZDp)
print '--------whiten yali--------------'
centroids=fca(whiten_ZDp,NP)
print n,'-->',[vno.nodes[i+1].id for i in centroids]
def edge_sst_test(vno):
NQ=50
#计算风量灵敏度(2维数组,即数组的数组)
Dq=edge_sst(vno,nCount=10,accuracy=0.01)
#print Dq
np.savetxt('Dq.txt',Dq)
#数据标准化(白化)
whiten_Dq = whiten(np.abs(Dq))
#print whiten_Dq
np.savetxt('whiten_Dq.txt',whiten_Dq)
#分类后计算聚类重心的索引位置
print '------------风量灵敏度-----------'
centroids=fca(whiten_Dq,NQ)
print n,'-->',[vno.edges[i].id for i in centroids]
#----------相对风量灵敏度------------#
ZDq=sst(Dq)
np.savetxt('ZDq.txt',ZDq)
#数据标准化(白化)
whiten_ZDq = whiten(ZDq)
np.savetxt('whiten_ZDq.txt',whiten_ZDq)
print '--------标准化后的相对风量灵敏度--------------'
centroids=fca(whiten_ZDq,NQ)
print n,'-->',[vno.edges[i].id for i in centroids]
def cvno_test():
print("[vno模块]测试...")
# 初始化vno数据(测试用)
vno = vno_data()
#initVNO(vno)
#initVNO2(vno)
#initVNO3(vno)
initVNO4(vno)
# 解算通风网络
ret = cvno(vno)
# 打印解算结果
if not ret:
print("解算失败")
else:
# 输出解算风量
PrintVN(vno)
# 计算风量灵敏度(2维数组,即数组的数组)
Dq=edge_sst(vno,nCount=10,accuracy=0.01)
#print Dq
np.savetxt('Dq.txt',Dq)
#数据标准化(白化)
whiten_Dq = whiten(np.abs(Dq))
#print whiten_Dq
np.savetxt('whiten_Dq.txt',whiten_Dq)
#分类后计算聚类重心的索引位置
print '------------fengliang-----------'
for n in range(1,9):
centroids=fca(whiten_Dq,n)
print n,'-->',[vno.edges[i].id for i in centroids]
# 计算压力灵敏度(2维数组,即数组的数组)
Dp=node_sst(vno,nCount=50,accuracy=0.01)
#print Dp
np.savetxt('Dp.txt',Dp)
whiten_Dp = whiten(np.abs(Dp))
#print whiten_Dp
np.savetxt('whiten_Dp.txt',whiten_Dp)
#分类后计算聚类重心的索引位置
print '--------yali--------------'
for n in range(1,9):
centroids=fca(whiten_Dp,n)
#print centroids
print n,'-->',[vno.nodes[i+1].id for i in centroids]
#np.savetxt('centroids_Dq.txt',[vno.edges[i].id for i in centroids])
#---------------------------------------#
ZDq=sst(Dq)
print ZDq
np.savetxt('ZDq.txt',ZDq)
#数据标准化(白化)
whiten_ZDq = whiten(ZDq)
#print whiten_ZDq
np.savetxt('whiten_ZDq.txt',whiten_ZDq)
print '--------whiten fengliang--------------'
for n in range(1,9):
centroids=fca(whiten_ZDq,n)
print n,'-->',[vno.edges[i].id for i in centroids]
#------------------
ZDp=sst(Dp)
np.savetxt('ZDp.txt',ZDp)
#数据标准化(白化)
whiten_ZDp = whiten(ZDp)
print whiten_ZDp
np.savetxt('whiten_ZDp.txt',whiten_ZDp)
print '--------whiten yali--------------'
for n in range(1,9):
centroids=fca(whiten_ZDp,n)
print centroids
print n,'-->',[vno.nodes[i+1].id for i in centroids]
