0

tdondur x 坐标，y 坐标 s。但我做不到。如果你能帮忙。谢谢大家。现在只绘制 tdondur 值。如何绘制 s 的值。

from math import e, exp,pow 
import pylab
from gaussxw import gaussxwab

Kb=1.3806505*10**-23
V=10**-3
p=6.022*10**28
Qd=428.0


def f(x):
    isisigasi=9*V*p*Kb*(T/Qd)**3
    return(isisigasi*x**4*(exp(x))/(exp(2)-1)**2)

tut=[] 

Tdondur=range(5,501)
tut=[None for T in Tdondur] 

for tsirala, T in enumerate (Tdondur):
    N=50
    a=0
    b=Qd/T
    x,w=gaussxwab(N,a,b)
    s=0.0
    for k in range (N):
        s+=w[k]*f(x[k])    
    print(T,"'in Isı Sığası : ",s)
    print("----------------------------------------------------------")


pylab.plot(Tdondur)
pylab.title('Isı Sığasının Sıcaklığa Göre Değişimi')
pylab.xlabel('Sıcaklık (T)')
pylab.ylabel('Isı Sığası')
pylab.gcf().savefig('isisigasi.png')
pylab.show()
4

1 回答 1

1

这里的问题是您在s循环的每次迭代中计算一个新值，但您永远不会保存您计算的内容。s相反，在进入循环之前创建一个数组，并s[k]在循环内设置。

from math import e, exp,pow 
import pylab
from gaussxw import gaussxwab

Kb=1.3806505*10**-23
V=10**-3
p=6.022*10**28
Qd=428.0


def f(x):
    isisigasi=9*V*p*Kb*(T/Qd)**3
    return(isisigasi*x**4*(exp(x))/(exp(2)-1)**2)

tut=[] 

Tdondur=range(5,501)
tut=[None for T in Tdondur] 

s = [0. for i in Tdondur]  # initialize s before entering the loop
for tsirala, T in enumerate (Tdondur):
    N=50
    a=0
    b=Qd/T
    x,w=gaussxwab(N,a,b)
    for k in range (N):
        s[k]+=w[k]*f(x[k])    # work on the values within the loop

pylab.plot(Tdondur, s)
#s_np = pylab.array(s)
#pylab.plot(Tdondur, pylab.log10(s_np))
pylab.title('title')
pylab.xlabel('T')
pylab.ylabel('Y')
#pylab.gcf().savefig('')
pylab.show()

结果如下（注意 x 轴被放大了）：

在此处输入图像描述

通常对于这样的数据，您会想要绘制log数据的 ，为此，请取消注释之后的两行：

s_np = pylab.array(s)
pylab.plot(pylab.log10(s_np))

（请注意，上面的行使用的是 numpy，它是通过 pylab 导入的，并且比此类问题的列表推导更容易。）

在此处输入图像描述

于 2013-06-10T16:18:00.597 回答