绘制正余弦函数
import matplotlib.pylab as plt
import numpy as np
x=np.linspace(-np.pi,2*np.pi,200)
y=np.sin(x)
z=np.cos(x)
plt.plot(x,y,'r')
plt.plot(x,z,'b--',linewidth=2.5)
plt.show()
XRD 谱图绘制
import matplotlib.pyplot as plt
import numpy as np
with open('XRD_AFO.csv','r',encoding='UTF-8') as f:
x,y=[],[]
for i in f.readlines()[1:]:
l=i.strip().split(',')
x.append(eval(l[0]))
y.append(eval(l[1]))
plt.plot(x,y)
plt.xlabel('Position(2-Theta)')
plt.ylabel('Intensity')
plt.show()
分态密度图
import matplotlib.pyplot as plt
import numpy as np
with open('DOSofBaTiO3PDOSTitle.csv', 'r', encoding='utf-8') as f:
ls=[]
for i in f.readlines()[1:]:
ls.append(i.strip().split(','))
m,n=map(int,input().split())
if m>n:
m,n=n,m
x,y=[],[]
for i in range(0,7,2):
x=[eval(j[i]) for j in ls if m<=eval(j[i])<=n and eval(j[i])!=0]
y=[eval(j[i+1]) for j in ls if m<=eval(j[i])<=n and eval(j[i])!=0]
plt.plot(x,y)
