2025年python函数图像绘制程序(python绘制函数图像代码)

python函数图像绘制程序(python绘制函数图像代码)这是我自己编的 我也是个初学者 希望对你有用 有些地方也可以改一下 不完善 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

大家好,我是讯享网,很高兴认识大家。

 这是我自己编的,我也是个初学者,希望对你有用,有些地方也可以改一下,不完善。

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()


小讯
上一篇 2025-05-26 14:15
下一篇 2025-06-03 10:41

相关推荐

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