动画剧情更新

Question

我正在尝试查找正弦曲线的最大值并将其显示在一个子图中。并尝试将其更新为动画。但是最大值的子图给出所有的零值。当我打印数组时，它不是零。我认为它不更新y值。我无法弄清楚原因。任何帮助将不胜感激。

我将把我的代码，它是可执行的：

from pylab import * 
import time 

ion() 
fs = 1e6 
Ts = 1/fs 
SNR=10 
sinfreq=2*pi*1e5 
pack= 512 
t = Ts*arange(0,pack) 
f = fs*(arange(0,pack)-pack/2)/pack 
max_y = zeros (len(t)) 
y=sin(sinfreq*t)   
y=y+randn(size(y))/sqrt(10^(SNR/10)*2) 
subplot(211)  
line1, = plot(y) 

subplot(212) 
line2, = plot(max_y) 
for i1 in arange(1,1000): 
    y=sin(sinfreq*t) 
    y=y+randn(size(y))/sqrt(10^(SNR/10)*2) 

    line1.set_ydata(y) 

    mk=0 
    for mk in range(0,len(y)): 
     if y[mk] > max_y[mk]: 
      max_y[mk] = y[mk] 
    print max_y 
    line2.set_ydata(max_y) 
    draw()       
    waitforbuttonpress(timeout=0.5) 

Answer 1

好的，我偶然发现了解决方案，我不知道原因，但我先更新max_y图的任何其他值，然后进行实际更新，然后该图显示我的更改。除此之外，它没有显示。我尝试使用不同的for循环来绘制一个更新，但对于我的循环它想更新两次。

我也添加set_ylim来看看更好的限制。我把星星放在我改变的地方。我也把新的代码。我希望它能帮助那些有同样问题的人。

from pylab import * 
import time 

ion() 
fs = 1e6 
Ts = 1/fs 
SNR=10 
sinfreq=2*pi*1e5 
pack= 512 
t = Ts*arange(0,pack) 
f = fs*(arange(0,pack)-pack/2)/pack 
max_y = zeros (len(t)) 
y=sin(sinfreq*t)   
y=y+randn(size(y))/sqrt(10^(SNR/10)*2) 
subplot(211)  
line1, = plot(y) 

sub2=subplot(212) #**** 
line2, = plot(max_y) 
for i1 in arange(1,1000): 
    y=sin(sinfreq*t) 
    y=y+randn(size(y))/sqrt(10^(SNR/10)*2) 

    line1.set_ydata(y) 

    mk=0 
    for mk in range(0,len(y)): 
     if y[mk] > max_y[mk]: 
      max_y[mk] = y[mk] 
    #print max_y 
    line2.set_ydata(zeros(len(max_y)))#****  
    line2.set_ydata(max_y) 
    sub2.set_ylim(min(max_y),max(max_y)) #**** 
    draw()       
    waitforbuttonpress(timeout=0.5) 

Answer 2

你忘了在第二个for循环的缩进。 这实际上应该给你一个IndentationError，所以我不明白你怎么可以说程序是可执行的（实际上，我编辑了你的条目以删除剩余的“在此输入代码”语句;如果你已经检查并且复制粘贴你的条目，你可能会发现两个错误）。

但是，你肯定不希望只是

max_y[:] = max(y) 

的而不是for循环？