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

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

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

MatplotlibはPythonのおよび、NumPy用のグラフ描画ライブラリです。多くの場合、IPythonと連携して使われます。

NumPy

NumPyはPythonのプログラミング言語の科学的と数学的なコンピューティングに関する拡張モジュールです。

Python 3.x

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

pandas

Pandasは、PythonでRにおけるデータフレームに似た型を持たせることができるライブラリです。 行列計算の負担が大幅に軽減されるため、Rで行っていた集計作業をPythonでも比較的簡単に行えます。 データ構造を変更したりデータ分析したりするときにも便利です。

Q&A

1回答

8504閲覧

外部データをFourier変換することで周波数成分を可視化したい

55cm1m

総合スコア6

Matplotlib

MatplotlibはPythonのおよび、NumPy用のグラフ描画ライブラリです。多くの場合、IPythonと連携して使われます。

NumPy

NumPyはPythonのプログラミング言語の科学的と数学的なコンピューティングに関する拡張モジュールです。

Python 3.x

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

pandas

Pandasは、PythonでRにおけるデータフレームに似た型を持たせることができるライブラリです。 行列計算の負担が大幅に軽減されるため、Rで行っていた集計作業をPythonでも比較的簡単に行えます。 データ構造を変更したりデータ分析したりするときにも便利です。

0グッド

3クリップ

投稿2018/09/05 04:15

編集2022/01/12 10:55

前提・実現したいこと

ここに質問の内容を詳しく書いてください。
センサーから得られた文字列の外部データの値を読み込んでプロットしようとしています。
そのデータをフーリエ変換することで、周波数成分を可視化したいのですが、エラーが出てしまいました。

どのようにしたらこのwarningを解消して、周波数成分を可視化することができるでしょうか。
初歩的な質問かもしれませんが、よろしくお願いいたします。

i am trying a Matlab code in Python
i want to see FREQUENCY COMPONENT.
so i try to plot by reading data and do Fourier trans form.
but my code gives a warning

発生している問題・エラーメッセージ

C:\Users******\Anaconda3\lib\site-packages\numpy\core\numeric.py:492: ComplexWarning: Casting complex values to real discards the imaginary part return array(a, dtype, copy=False, order=order)

該当のソースコード

python

1import pandas as pd 2import numpy as np 3import matplotlib.pyplot as plt 4 5if __name__ == '__main__': 6 df = pd.read_csv("sleep.csv", names=('time','values','-400')) 7 8 x = df["time"] 9 y = df["-400"] 10 z = df["values"] 11 12 f = np.fft.fft(y) 13 f_abs = np.abs(f) 14 print(f_abs) 15 16 plt.subplot(211) 17 plt.plot(y) 18 19 plt.subplot(212) 20 plt.plot(f_abs) 21 22 plt.show()

試したこと

ここに問題に対して試したことを記載してください。

補足情報(FW/ツールのバージョンなど)

人物の転倒検出のために、
横になっている、座っている、転倒するなどの区別をできるような
周波数成分を知りたいです。

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

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

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

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

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

guest

回答1

0

フーリエ変換した結果は、周波数ごとのベクトル(大きさと向き(位相)))、もしくは複素数(実数部と虚数部なら成る)になります。

複素数は実数部と虚数部を持ちますが、それを実数にキャストすると虚数部の情報が失われます。
それを警告しているのが"Casting complex values to real discards the imaginary part"という一文です。

どのような周波数成分を見たい(得たい)のでしょうか?
周波数ごとの大きさ(ベクトルの長さ)だけ判れば良いのですか?

投稿2018/09/05 04:51

coco_bauer

総合スコア6915

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

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

55cm1m

2018/09/05 05:22

回答ありがとうございます。 worningの意味はそういうことだったんですね。納得しました。 今の段階でどのような周波数成分を知ることが最適なのかがわかっていないので、 はっきりとはお答えできないのですが(すみません)、、 まずは、周波数ごとの大きさを知りたいです。 追加の質問になってしまうのですが、 他にどのような周波数成分を知ることができるのでしょうか?
guest

あなたの回答

tips

太字

斜体

打ち消し線

見出し

引用テキストの挿入

コードの挿入

リンクの挿入

リストの挿入

番号リストの挿入

表の挿入

水平線の挿入

プレビュー

まだベストアンサーが選ばれていません

会員登録して回答してみよう

アカウントをお持ちの方は

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

ただいまの回答率
85.50%

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

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

質問する

関連した質問