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

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

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

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

Q&A

0回答

2370閲覧

Pythonを用いてグラフの極大値を求めたい

tapio21

総合スコア1

Python

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

0グッド

0クリップ

投稿2021/12/14 05:25

編集2021/12/14 05:31

お世話になっております。
Pythonを使用して、極大値を算出したいです。
数値が記入されているExcelファイルにargrelmaxを使用して求めようと思いましたが、うまくいきません。
独学でプログラミングを始めたばかりですので、初歩的な質問で申し訳ございません。
やりたいこととしては、一つのグラフより、グラフの立ち上がりと立下りよりそれぞれを抜き出して保存していきたいです。

df1=df["water_lev"] print(df1) 0 -2.09 1 -2.09 2 -2.09 3 -2.08 4 -2.08 ... 8779 -2.33 8780 -2.35 8781 -2.37 8782 -2.39 8783 -2.40 Name: water_lev, Length: 8784, dtype: float64
import numpy as np from scipy import signal from scipy import fftpack from matplotlib import pyplot as plt df = pd.read_excel('test.xlsx') df1=df["water_lev"] # フーリエ変換して複素スペクトル、振幅、周波数を出力する関数 def fft(df1, dt): spec = fftpack.fft(df1) # フーリエスペクトル(複素数) freq = fftpack.fftfreq(len(spec), d=dt) # 周波数軸を生成 amp = np.sqrt(spec.real ** 2 + spec.imag ** 2) / (len(spec)/2) # 振幅を計算 return spec, amp, freq # 波形(x, y)からn個のピークを幅wで検出する関数(xは0から始まる仕様) def findpeaks(df1, n, w): index_all = list(signal.argrelmax(df1, order=w)) # scipyのピーク検出 index = [] # ピーク指標の空リスト peaks = [] # ピーク値の空リスト # n個分のピーク情報(指標、値)を格納 for i in range(n): index.append(index_all[0][i]) peaks.append(y[index_all[0][i]]) index = np.array(index) * x[1] # xの分解能x[1]をかけて指標を物理軸に変換 return index, peaks # フーリエ変換する spec, amp, freq = fft(df1, dt) # ピーク検出する index, peaks = findpeaks(freq, amp, n=100, w=6)

イメージ説明

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

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

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

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

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

guest

あなたの回答

tips

太字

斜体

打ち消し線

見出し

引用テキストの挿入

コードの挿入

リンクの挿入

リストの挿入

番号リストの挿入

表の挿入

水平線の挿入

プレビュー

まだ回答がついていません

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

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

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

ただいまの回答率
85.35%

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

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

質問する

関連した質問