arXiv:1502.01588
画thermal dust的亮温谱
import matplotlib.pyplot as plt
import numpy as np
Td = 23 #K
kB = 1.3806505*10**(-23) #J/K
h = 6.6260693*10**(-34+9) #J·s,+9后单位变成J/GHz
v0 = 545 #GHz
betad = 1.55
Ad = 10
r = h/(kB*Td)
v = np.linspace(2,1200)
sd = Ad*(v/v0)**(betad+1)*((np.exp(r*v0)-1)/(np.exp(r*v)-1))
plt.loglog(v,sd,linewidth='8',color='r')
plt.show()
结果为: