计算分子动力学模拟轨迹contact-map概率，并绘制热图

1.将矩阵文件处理成一个文件

mport numpy as np
import matplotlib.pyplot as plt
from matplotlib import cm
import seaborn as sns

2.将多个矩阵读入到一个列表sunmata中，且列表的每一个值为矩阵每一行

b = open('55.txt','w')    # 对文件进行写操作
sumdata = [] #建立列表
squ=[]
with open('44.txt','r')as f:  #文件读取操作
   data=f.readlines()     #逐行读取，且data是列表，列表的每一个元素是读取的每一行包括换行
sumdata+=[line.strip()for line in data]  #删除空格符，把列表元素存入sumdata列表

3.把sumdata列表中的元素以0/1写入squ列表中

 for i in range ( 0,len(sumdata) ):  
   for j in range ( 4 ):
       if sumdata[i][j] == str('n'):        #判断元素是否为n，是记为0 ，否记为1
           squ.append(0)   
       else:
           squ.append ( 1 )

4.把suq列表中元素按一个矩阵为单位写入v列表

q=[squ[x:x+4] for x in range(0