Pythonでのfor文の使い方についての質問です。
以下のプログラムについて、[波形を12分割し~]の部分をfor文を用いて12個目まで計算し、それぞれの結果をグラフに示し、さらに各グラフの最も価の高い部分をつなげたグラフを作りたいと考えています。
#coding:utf-8 import numpy as np import pylab import wave def wavread(filename): wf = wave.open(filename, "r") fs = wf.getframerate() x = wf.readframes(wf.getnframes()) x = np.frombuffer(x, dtype="int16") / 32768.0 # (-1, 1)に正規化 wf.close() return x, float(fs) if __name__ == "__main__": # 波形を表示 wav, fs = wavread("sinwave_440hz.wav") t = np.arange(0.0, len(wav) / fs, 1/fs) pylab.plot(t * 1000, wav) pylab.xlabel("time [ms]") pylab.ylabel("amplitude") pylab.show() #波形を12分割し、1つ目から2つ目について検討する cuttime = len(wav)/12 # 波形を12分割 wavdata = wav[int(1*cuttime): int(2*cuttime)] time = t[int(1*cuttime): int(2*cuttime)] pylab.subplot(211) pylab.plot(time * 1000, wavdata) pylab.ylabel("amplitude") # ハニング窓をかける hanningWindow = np.hanning(len(wavdata)) wavdata = wavdata * hanningWindow pylab.subplot(212) pylab.plot(time * 1000, wavdata) pylab.xlabel("time [ms]") pylab.ylabel("amplitude") pylab.show() # 切り出した波形のスペクトルを求める n = 2048 # FFTのサンプル数 # 離散フーリエ変換 dft = np.fft.fft(wavdata, n) # 振幅スペクトル Adft = np.abs(dft) # パワースペクトル Pdft = np.abs(dft) ** 2 # 周波数スケール fscale = np.fft.fftfreq(n, d=1.0 / fs) # プロット pylab.subplot(211) pylab.plot(fscale[0:int(n / 2)], Adft[0:int(n / 2)]) pylab.xlabel("frequency [Hz]") pylab.ylabel("amplitude spectrum") pylab.xlim(0, 1000) pylab.subplot(212) pylab.plot(fscale[0:int(n / 2)], Pdft[0:int(n / 2)]) pylab.xlabel("frequency [Hz]") pylab.ylabel("power spectrum") pylab.xlim(0, 1000) pylab.show()
まず、for文を用いて次のようにしたところ、エラー(後記)が発生しました。
for k in(12): wavdata = wav[int(k*cuttime): int((k+1)*cuttime)] time = t[int(k*cuttime): int((k+1)*cuttime)] pylab.subplot(211) pylab.plot(time * 1000, wavdata) pylab.ylabel("amplitude") # ハニング窓をかける hanningWindow = np.hanning(len(wavdata)) wavdata = wavdata * hanningWindow pylab.subplot(212) pylab.plot(time * 1000, wavdata) pylab.xlabel("time [ms]") pylab.ylabel("amplitude") pylab.show() # 切り出した音声のスペクトルを求める n = 2048 # FFTのサンプル数 # 離散フーリエ変換 dft = np.fft.fft(wavdata, n) # 振幅スペクトル Adft = np.abs(dft) # パワースペクトル Pdft = np.abs(dft) ** 2 # 周波数スケール fscale = np.fft.fftfreq(n, d=1.0 / fs) # プロット pylab.subplot(211) pylab.plot(fscale[0:int(n / 2)], Adft[0:int(n / 2)]) pylab.xlabel("frequency [Hz]") pylab.ylabel("amplitude spectrum") pylab.xlim(0, 1000) pylab.subplot(212) pylab.plot(fscale[0:int(n / 2)], Pdft[0:int(n / 2)]) pylab.xlabel("frequency [Hz]") pylab.ylabel("power spectrum") pylab.xlim(0, 1000) pylab.show()
Traceback (most recent call last): File "C:/Users/KURATAMA/Desktop/卒研/卒研/sample/a4.py", line 27, in <module> for k in(12): TypeError: 'int' object is not iterable
for文の次の行にあるintが問題なのかと思って削除してみましたが、結果は変わりませんでした。
for文の使い方に問題があるのか、そもそもfor文ではできないことをやろうとしているのかなど、解決方法を教えていただけると幸いです。
回答1件
あなたの回答
tips
プレビュー
バッドをするには、ログインかつ
こちらの条件を満たす必要があります。
2020/01/11 00:07