这是我自己编的,我也是个初学者,希望对你有用,有些地方也可以改一下,不完善。
Plank1:
encoding: utf-8
x=list((x*0.01+1e-10)*1e-6 for x in range(0,500)) m=np.zeros(500)
def f(a):
for i in range(0,500): m[i]=2*np.pi*h*c2/(exp(h*c/(k*a*x[i]))-1)/x[i]5 return m
讯享网
g = [x[np.argmax(f(t1))], x[np.argmax(f(t2))], x[np.argmax(f(t3))], x[np.argmax(f(t4))], x[np.argmax(f(t5))], x[np.argmax(f(t6))] , x[np.argmax(f(t7))]] n = [max(f(t1)) , max(f(t2)) ,max(f(t3)) ,max(f(t4)) ,max(f(t5)) ,max(f(t6)) ,max(f(t7))]
print(np.argmax(f(t1)) , max(f(t1))) print(np.argmax(f(t2)) , max(f(t2))) print(np.argmax(f(t3)) , max(f(t3))) print(np.argmax(f(t4)) , max(f(t4))) print(np.argmax(f(t5)) , max(f(t5))) print(np.argmax(f(t6)) , max(f(t6))) print(np.argmax(f(t7)) , max(f(t7)))
plt.plot(x,f(t1),label=‘T=1400K’) plt.plot(x,f(t2),label=‘T=1500K’) plt.plot(x,f(t3),label=‘T=1600K’) plt.plot(x,f(t4),label=‘T=1700K’) plt.plot(x,f(t5),label=‘T=1800K’) plt.plot(x,f(t6),label=‘T=1900K’) plt.plot(x,f(t7),label=‘T=2000K’) plt.plot(g , n , color=‘red’, linewidth=2, linestyle=‘–’) plt.stem(g , n , linefmt=‘b:’ , markerfmt=‘C3.’, basefmt=‘r-’)
for a, b in zip(g, n):
讯享网plt.text(a, b, a, ha='left', va='bottom', fontsize=10) plt.title(u‘普朗克黑体单色辐射能力与波长的关系’) plt.xlabel(u‘波长λm’) plt.ylabel(u‘单色辐射能力M(λ,T)/(W▪m^)’)
plt.legend() plt.show()

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容,请联系我们,一经查实,本站将立刻删除。
如需转载请保留出处:https://51itzy.com/kjqy/164722.html