質問をすることでしか得られない、回答やアドバイスがある。

15分調べてもわからないことは、質問しよう!

新規登録して質問してみよう
ただいま回答率
85.50%
for

for文は、様々なプログラミング言語で使われている制御構造です。for文に定義している条件から外れるまで、for文内の命令文を繰り返し実行します。

Python 3.x

Python 3はPythonプログラミング言語の最新バージョンであり、2008年12月3日にリリースされました。

Q&A

解決済

1回答

1771閲覧

Pythonでfor文を使った解析をしたい

tamaki6485

総合スコア6

for

for文は、様々なプログラミング言語で使われている制御構造です。for文に定義している条件から外れるまで、for文内の命令文を繰り返し実行します。

Python 3.x

Python 3はPythonプログラミング言語の最新バージョンであり、2008年12月3日にリリースされました。

0グッド

0クリップ

投稿2020/01/09 05:18

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文ではできないことをやろうとしているのかなど、解決方法を教えていただけると幸いです。

気になる質問をクリップする

クリップした質問は、後からいつでもMYページで確認できます。

またクリップした質問に回答があった際、通知やメールを受け取ることができます。

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

guest

回答1

0

ベストアンサー

kを0から11まで回したい場合は

python

1for k in range(12):

としてください。

投稿2020/01/09 05:20

jeanbiego

総合スコア3966

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

tamaki6485

2020/01/11 00:07

超イージーミスだったとは・・・。直したのち、intをつける部分を変更したら実行できました。ご協力ありがとうございました。
guest

あなたの回答

tips

太字

斜体

打ち消し線

見出し

引用テキストの挿入

コードの挿入

リンクの挿入

リストの挿入

番号リストの挿入

表の挿入

水平線の挿入

プレビュー

15分調べてもわからないことは
teratailで質問しよう!

ただいまの回答率
85.50%

質問をまとめることで
思考を整理して素早く解決

テンプレート機能で
簡単に質問をまとめる

質問する

関連した質問