这是我自己编的，我也是个初学者，希望对你有用，有些地方也可以改一下，不完善。 Plank1: # encoding: utf-8 import matplotlib.pyplot as plt from sympy import * import numpy as np k=1.38054e-23 h=6.626e-34 c=2.e8 sigma=5.6688e-008 b=0.0029 t1 = 1400 t2 = 1500 t3 = 1600 t4 = 1700 t5 = 1800 t6 = 1900 t7 = 2000 t8 = 6000 array 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()到此这篇python函数图像绘制程序（python绘制函数图像代码）的文章就 介绍到这了,更多相关内容请继续浏览下面的相关 推荐文章，希望大家都能在编程的领域有一番成就！版权声明：
本文来自互联网用户投稿，该文观点仅代表作者本人，不代表本站立场。本站仅提供信息存储空间服务，不拥有所有权，不承担相关法律责任。
如若内容造成侵权、违法违规、事实不符，请将相关资料发送至xkadmin@xkablog.com进行投诉反馈，一经查实，立即处理！
转载请注明出处，原文链接：https://www.xkablog.com/pythonbc/27493.html