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

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

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

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

Q&A

解決済

1回答

669閲覧

FFT変換の値の表示について

hirohiro1999

総合スコア5

Python

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

0グッド

0クリップ

投稿2021/11/08 06:56

編集2021/11/09 04:57

・やってること
現在、単純なフーリエ変換をしていますが、振幅スペクトルであるy軸の最大値が以下の図のように最大1となってしまいます。
この振幅が100などの数値になるはずなのですが、そのように変換する方法がわかりません。
どのようにコードを変換すると数値にすることができますか?
わかる方がいらっしゃいましたら教えていただきたいです。
イメージ説明

・使っているコード

Python

1#! /usr/bin/env python 2# -*- coding: utf-8 -*- 3from pydub import AudioSegment as AS 4import ffmpeg 5import numpy as np 6import matplotlib.pyplot as plt 7import matplotlib as mpl 8import librosa 9import wave 10import math 11 12 13file = "2_cut5_bp" 14sound = AS.from_file(file + ".wav", "wav") 15data = np.array(sound.get_array_of_samples()) 16spec = np.fft.fft(data) #2次元配列(実部,虚部) 17freq = np.fft.fftfreq(data.shape[0], 1.0/sound.frame_rate) 18spec = spec[:int(spec.shape[0]/2 + 1)] #周波数がマイナスになるスペクトル要素の削除 19freq = freq[:int(freq.shape[0]/2 + 1)] #周波数がマイナスになる周波数要素の削除 20g = wf.readframes(wf.getnframes()) 21g = np.frombuffer(g, dtype= "int16")/32768.0 # -1~1に正規化 22max_spec=max(np.abs(spec)) #最大音圧を取得(音圧を正規化するために使用) 23 24fig = plt.figure() 25 26plt.plot(freq, np.abs(spec)/max_spec) 27plt.grid() 28plt.xlim([0,20000]) #グラフに出力する周波数の範囲[Hz] 29plt.xticks([2000, 2500, 3000, 3500, 4000, 4500, 5000]) 30plt.xlabel("周波数[Hz]",fontname="MS Gothic") 31plt.ylabel("振幅スペクトル",fontname="MS Gothic") 32 33#plt.yscale("log") 34output ="C:\Users\hiro2\Desktop\sound\kara" 35fig.savefig(file + ".png") 36plt.close()

追記-max_specを省いた場合
イメージ説明
振幅スペクトルが変わらずに左上に1e8と表示されてしまいます。

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

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

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

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

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

guest

回答1

0

ベストアンサー

python

1plt.plot(freq, np.abs(spec)/max_spec) # ここでmax_specで割ってる

投稿2021/11/08 07:13

ozwk

総合スコア13528

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

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

hirohiro1999

2021/11/09 04:55

なるほど! しかしそこでmax_specで割る処理を省いても表では反映されずに追記の図のようになってしまいます。 ここを数値にする方法を教えていただきたいです。。。
ozwk

2021/11/09 04:59

1e8は10の8乗という意味です
hirohiro1999

2021/11/09 05:09

左にプロットされている0.2や0.4などの数字に10の8乗をすればちゃんとした数値になるということでしょうか?
ozwk

2021/11/09 05:14

はい
hirohiro1999

2021/11/09 05:18

つまり例えば1.0×10の8乗の場合大きすぎて表示できないから1e8と表示されているということでしょうか?
ozwk

2021/11/09 06:38

はい、"100000000"だと長いから"1e8"と表記してます。 Excelなどでも出てくる割と一般的な表記だと思います。
hirohiro1999

2021/11/09 06:43

なるほど! ありがとうございます!!!
guest

あなたの回答

tips

太字

斜体

打ち消し線

見出し

引用テキストの挿入

コードの挿入

リンクの挿入

リストの挿入

番号リストの挿入

表の挿入

水平線の挿入

プレビュー

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

ただいまの回答率
85.48%

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

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

質問する

関連した質問