python中的for循环绘图

Question

我试图在同一个图中用 python 绘制 sinc 函数，这基本上是一个 OFDM 载波信号，将在第二个图中总结。

你能告诉我有什么问题吗？这是代码sn-p。

NoOfCarriers = 11 
interval = math.pi/50
f = np.arange((-5*math.pi),(5*math.pi),interval)  
fnoiseMax = 0.3
iMin = -(NoOfCarriers-1)//2
iMax =  (NoOfCarriers-1)//2
csum = np.zeros(len(f))
fList = [];cList = []

ax = plt.subplot(111)

for i in range(iMin,iMax):
    print("i = ", i)
    fnoise = fnoiseMax*(np.random.uniform(-1,1))
    fshift = (i * (1//math.pi) * math.pi) + fnoise
    c = np.sinc(f - fshift)
    csum = csum + c[i]
    fList = [fList,fshift]
    cList = [cList,max(c)]
    ax.plot(f, c)
    plt.grid(True)
    plt.show()

这是我得到的：

这是我所期望的：

• 我不知道如何在python中添加stem函数。干函数stem((i * (1/pi) * pi) + fnoise,1)的基本数学逻辑

从循环中取出 plt.show 后更新了绘图

Answer 1

尝试将plt.show 置于循环之外，如下所示：

import math
import numpy as np
import matplotlib.pyplot as plt

NoOfCarriers = 11
interval = math.pi/50
f = np.arange((-5*math.pi),(5*math.pi),interval)
fnoiseMax = 0.3
iMin = -(NoOfCarriers-1)//2
iMax =  (NoOfCarriers-1)//2
csum = np.zeros(len(f))
fList = [];cList = []

ax = plt.subplot(111)

for i in range(iMin,iMax):
    print("i = ", i)
    fnoise = fnoiseMax*(np.random.uniform(-1,1))
    fshift = (i * (1//math.pi) * math.pi) + fnoise
    c = np.sinc(f - fshift)
    csum = csum + c[i]
    fList = [fList,fshift]
    cList = [cList,max(c)]
    ax.plot(f, c)

plt.grid(True)
plt.show()

当并非所有地块都在ax 上时，问题是触发show 函数（阻塞）。