Question

我正在尝试使用我从plotting multiple plots but whith offset ranges python学到的东西，但我似乎无法对我的Legendre绘图代码进行适当的调整：

import numpy as np
import pylab
from numpy.polynomial.legendre import leggauss, legval


def f(x):
    if 0 <= x <= 1:
        return 1
    if -1 <= x <= 0:
        return -1


f = np.vectorize(f)

deg = 1000
x, w = leggauss(deg)  #  len(x) == deg
L = np.polynomial.legendre.legval(x, np.identity(deg))
integral = (L * (f(x) * w)[None,:]).sum(axis=1)
xx = np.linspace(-1, 1, 500000)
csum = []


for N in [5, 15, 25, 51, 97]:
    c = (np.arange(1, N) + 0.5) * integral[1:N]
    clnsum = (c[:,None] * L[1:N,:]).sum(axis = 0)
    csum.append(clnsum)


fig = pylab.figure()
ax = fig.add_subplot(111)


for i in csum:
    ax.plot(x, csum[i])


pylab.xlim((-1, 1))
pylab.ylim((-1.25, 1.25))
pylab.plot([0, 1], [1, 1], 'k')
pylab.plot([-1, 0], [-1, -1], 'k')
pylab.show()

我正在使用csum为clnsum保留N = 5, 15, 25, 51, 97的每次迭代。然后我想绘制每个存储的clnsum，但我相信这就是问题发生的地方。

我相信

for i in csum:

是正确的设置，但ax.plot(x, csum[i])必须是绘制每次迭代的错误方法。至少，这是我所相信的，但也许整个设置是错误的或错误的。

如何为每个clnsum实现每个N的绘图？

Answer 1

for i in csum:
    ax.plot(x, csum[i])

这就是你的问题所在。我不是一个整数，它是一个数组。你可能意味着

for i in range(len(csum)):

您也可以

for y in csum:
    ax.plot(x, y)

Python：在同一窗口中绘制多个图

1 个答案: