ウェーブ パッケージのアニメーションを作成し、ムービーとして保存しようとしています。保存以外はすべて機能しています。私が間違っていることを教えてください。ラインに入ると、ani.save('MovWave.mp4')
彼は私にこう言います:
writer = writers.list()[0]
IndexError: list index out of range
もちろんググってみましたが、意味がわかりません。
更新:ffmpeg
コンソールで呼び出すことができるようになりました。ffmpeg バージョンが0.10.7-6:0.10.7-0jon1~precise
インストールされていると表示されます。コードを更新してプログラムを実行しましたが、次のエラーが発生しました。
Traceback (most recent call last):
ani.save('MovWave.mpeg', writer="ffmpeg")
writer.grab_frame()
dpi=self.dpi)
self.canvas.print_figure(*args, **kwargs)
self.figure.dpi = origDPI
self.dpi_scale_trans.clear().scale(dpi, dpi)
self._mtx = np.identity(3)
from numpy import eye
File "<frozen importlib._bootstrap>", line 1609, in _handle_fromlist
UnicodeDecodeError: 'utf-8' codec can't decode byte 0xff in position 0: invalid start byte
更新 2: Doctorlove が指摘したように、python 3.3 を使用するとバグがあるようです。現在、代わりに python 2.7 を使用しようとしています。現在、mpeg ファイルを作成しますが、再生できず、約 150 kB しかありません。
更新 3:わかりました。Win7 マシンでまったく同じコードを試してみましたが、Python 3.3 でも動作します。しかし、私は同じ問題を抱えています.以前にpython 2.7. 作成された mpeg ファイルは再生できず、数百 kB しかありません。
#! coding=utf-8
import numpy as np
import matplotlib.pyplot as plt
import matplotlib.animation as animation
import time
time.clock()
def FFT(x,y):
X = (x[-1]-x[0])/len(y)
f = np.linspace(-2*np.pi/X/2,2*np.pi/X/2,len(y))
F = np.fft.fftshift(np.fft.fft(y))/np.sqrt(len(y))
return(f,F)
def FUNCTION(k_0,dx,c,t):
y = np.exp(1j*k_0*(x-c*t))*np.exp(-((x-c*t)/(2*dx))**2 )*(2/np.pi/dx**2)**(1/4)
k,F = FFT((x-c*t),y)
return(x,y,k,F)
#Parameter
N = 1000
x = np.linspace(0,30,N)
k_0 = 5
dx = 1
c = 1
l = [k_0,c,dx]
fig = plt.figure("Moving Wavepackage and it's FFT")
sub1 = plt.subplot(211)
sub2 = plt.subplot(212)
sub2.set_xlim([-10,10])
sub1.set_title("Moving Wavepackage and it's FFT")
sub1.set_ylabel("$Re[\psi(x,t)]$")
sub1.set_xlabel("$t$")
sub2.set_ylabel("$Re[\psi(k_x,t)]$")
sub2.set_xlabel("$k_x$")
n = 50
t = np.linspace(0,30,n)
img = []
for i in range(n):
x,y,k,F = FUNCTION(k_0,dx,c,t[i])
img.append(plt.plot(x,np.real(y),color="red", axes=plt.subplot(211)))
img.append(plt.plot(k,np.real(F),color="red", axes=plt.subplot(212)))
ani = animation.ArtistAnimation(fig, img, interval=20, blit=True, repeat_delay=0)
ani.save('MovWave.mpeg', writer="ffmpeg")
print(time.clock())
plt.show()