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

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

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

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

if

if文とは様々なプログラミング言語で使用される制御構文の一種であり、条件によって処理の流れを制御します。

Python

Pythonは、コードの読みやすさが特徴的なプログラミング言語の1つです。 強い型付け、動的型付けに対応しており、後方互換性がないバージョン2系とバージョン3系が使用されています。 商用製品の開発にも無料で使用でき、OSだけでなく仮想環境にも対応。Unicodeによる文字列操作をサポートしているため、日本語処理も標準で可能です。

配列

配列は、各データの要素(値または変数)が連続的に並べられたデータ構造です。各配列は添え字(INDEX)で識別されています。

解決済

pythonでx=aの時のyの値を出力したい

tamaki6485
tamaki6485

総合スコア6

for

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

if

if文とは様々なプログラミング言語で使用される制御構文の一種であり、条件によって処理の流れを制御します。

Python

Pythonは、コードの読みやすさが特徴的なプログラミング言語の1つです。 強い型付け、動的型付けに対応しており、後方互換性がないバージョン2系とバージョン3系が使用されています。 商用製品の開発にも無料で使用でき、OSだけでなく仮想環境にも対応。Unicodeによる文字列操作をサポートしているため、日本語処理も標準で可能です。

配列

配列は、各データの要素(値または変数)が連続的に並べられたデータ構造です。各配列は添え字(INDEX)で識別されています。

1回答

0グッド

0クリップ

3543閲覧

投稿2020/01/11 14:40

pythonで、スペクトルの値(y軸)が最大値をとるときの周波数(x軸)を出力したいと考えています。
以下のプログラムで、それぞれのkについてスペクトルの最大値を出力することには成功しました。しかし、その時の周波数を求めるために書いた最後の行は意味をなしていないように思えます。また、mAdftが一回ごとに書き換えられてしまっており、配列の形になっていないと思っています。

#cording: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("eiyuu.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() # 母音の定常部分(中心部)のスペクトルを center = int(len(wav) / 2) # 中心のサンプル番号 cuttime = len(wav)/12 # 切り出す長さ [s] for k in range(12): wavdata = wav[(k*int(cuttime)): ((k+1)*int(cuttime))] time = t[(k*int(cuttime)): ((k+1)*int(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() mAdft=[max(Adft)] print(mAdft) print(fscale)

そこで、プログラム最終部分をif文を使って以下のように書き換えたところ、後述するようなエラーが起こってしまいました。

mAdft=[max(Adft)] print(mAdft) print(fscale) if Adft==mAdft: print(fscale)
if Adft==mAdft: ValueError: The truth value of an array with more than one element is ambiguous. Use a.any() or a.all()

おそらく、Adft,mAdft,fscaleがすべて配列であることが原因のエラーのようですが、解決策までは掴めませんでした。
繰り返しになりますが、for文のそれぞれの出力結果をすべて配列に入れる方法、および、スペクトルが最大値の時の周波数の値を出力する方法を教えていただきたいです。

以下のような質問にはグッドを送りましょう

  • 質問内容が明確
  • 自分も答えを知りたい
  • 質問者以外のユーザにも役立つ

グッドが多くついた質問は、TOPページの「注目」タブのフィードに表示されやすくなります。

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

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

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

下記のような質問は推奨されていません。

  • 間違っている
  • 質問になっていない投稿
  • スパムや攻撃的な表現を用いた投稿

適切な質問に修正を依頼しましょう。

回答1

0

ベストアンサー

AdftとmAdftの要素数合ってますか??

forループの中で毎回
mAdft=[max(Adft)]
を呼び出した場合、mAdftはkが最後の時のmax(Adft)が1つだけ入った配列になるかと思います。
ループの前で空のmAdftを定義してあげてから
np.append(mAdft,max(Adft))
としないといけないような気がします。

間違ってたらすみません

投稿2020/01/11 16:59

RyotaroIsoyama

総合スコア183

良いと思った回答にはグッドを送りましょう。
グッドが多くついた回答ほどページの上位に表示されるので、他の人が素晴らしい回答を見つけやすくなります。

下記のような回答は推奨されていません。

  • 間違っている回答
  • 質問の回答になっていない投稿
  • スパムや攻撃的な表現を用いた投稿

このような回答には修正を依頼しましょう。

回答へのコメント

RyotaroIsoyama

2020/01/11 17:07

あ、ちがう、そのときの周波数がほしいんですね。でしたら fscale[np.argmax(Adft)] で取れそうな気がします。 試してないですが、、
tamaki6485

2020/01/13 02:00

ありがとうございます。 mfscale=fscale[np.argmax(Adft)] print(mfscale) で実行できました。

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

ただいまの回答率
86.12%

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

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

質問する

関連した質問

同じタグがついた質問を見る

for

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

if

if文とは様々なプログラミング言語で使用される制御構文の一種であり、条件によって処理の流れを制御します。

Python

Pythonは、コードの読みやすさが特徴的なプログラミング言語の1つです。 強い型付け、動的型付けに対応しており、後方互換性がないバージョン2系とバージョン3系が使用されています。 商用製品の開発にも無料で使用でき、OSだけでなく仮想環境にも対応。Unicodeによる文字列操作をサポートしているため、日本語処理も標準で可能です。

配列

配列は、各データの要素(値または変数)が連続的に並べられたデータ構造です。各配列は添え字(INDEX)で識別されています。