0

我是 python 和绘图的新手。我试图使用以下脚本生成一个情节。目标是为所有值绘制 Q 与 F 的关系图。

from pylab import *

n = 5
D = 13
B = 10

x = linspace(-6.5, 6.5, 1000)
y = 1/sqrt(2*pi)*exp(-(x)**2/2)

for i in range(1,n):
    F = sum(y*cos(2*pi*i*x/D)*exp(-i**2*B/(4*D**2)))
    print F

for j in range(1,n):
    Q = 2*pi*(j)/D
    print Q

plt.plot(Q,F,'rx')
plt.show()

当我运行脚本时,它只绘制一个数据点而不是全部。我敢肯定,我犯了一些愚蠢的错误。有人可以帮我吗?谢谢你。

4

2 回答 2

2

这是因为您将 F 和 Q 设置为每个循环上的值,而不是将值附加到数组的末尾。

from pylab import *

n = 5
D = 13
B = 10

x = linspace(-6.5, 6.5, 1000)
y = 1/sqrt(2*pi)*exp(-(x)**2/2)

F,Q = [],[]

for i in range(1,n):
    F.append(sum(y*cos(2*pi*i*x/D)*exp(-i**2*B/(4*D**2))))
    Q.append(2*pi*(j)/D)

plt.plot(Q,F,'rx')
plt.show()
于 2013-02-21T13:46:19.227 回答
1

您只是在循环中设置 Q、F 相等。

from pylab import *

n = 5
D = 13
B = 10

x = linspace(-6.5, 6.5, 1000)
y = 1/sqrt(2*pi)*exp(-(x)**2/2)

for i in range(1,n):
    F.append(sum(y*cos(2*pi*i*x/D)*exp(-i**2*B/(4*D**2))))
    print F
    Q.append(2*pi*(j)/D)
    print Q

plt.plot(Q,F, 'rx')    
plt.show()
于 2013-02-21T13:45:45.210 回答