前提・実現したいこと
グラフの極値を抽出して点と点をつなぎたいです。
・ここに質問の内容を詳しく書いてください。
実験データの点と点を結んでグラフにしました。
次に取り組みたい作業として、グラフの極値を抽出して点と点をつなごうとしています。
極値の抽出として、Scipyやargrelmaxといったコードを使うといいとアドバイスをいただきましたが、元のコードにどのようにして付け加えるかわかりません。
発生している問題・エラーメッセージ
Traceback (most recent call last): File "FSP_f0all_plot_singleKID'.py", line 2, in <module> from scipy import signal File "C:\Users\User\anaconda\lib\site-packages\scipy\signal\__init__.py", line 307, in <module> from ._peak_finding import * AttributeError: module 'scipy.signal._peak_finding' has no attribute 'argrelextrema'
該当のソースコード
import numpy as np from scipy import signal import matplotlib.pyplot as plt ###------ Edit Here ----------------- filename = "C:/Users/User/OneDrive/デスクトップ/卒研/data1.txt" samplingrate = 2**5 ###------ Edit Here ----------------- PlotTitle = [] PlotTitle.append(filename.split(".")[1]) with open(filename) as f: data = np.loadtxt(f) up_points = [] up_time = [] up_f0 = [] down_points = [] down_time = [] down_f0 = [] for i in range(0,len(data)): up_points.append(data[i][0]) up_time.append(data[i][0]/samplingrate) up_f0.append(data[i][2]) down_points.append(data[i][1]) down_time.append(data[i][1]/samplingrate) down_f0.append(data[i][3]) y = up_f0 maxid = signal.argrelmax(y, order=1) minid = signal.argrelmin(y, order=1) plt.figure() plt.plot(up_time, up_f0,lw =1,label="up sweep") plt.plot(down_time, down_f0,lw = 1,label="down sweep") plt.grid(True) plt.title(PlotTitle[0]) plt.xlabel("time [s]") plt.ylabel("Resonance Frequency [MHz]") plt.plot(x[maxid],y[maxid],'ro',label='ピーク値') plt.plot(x[minid],y[minid],'bo',label='ピーク値(最小)') plt.legend(loc = "best") plt.show()
試したこと
補足情報(FW/ツールのバージョンなど)
実験データを折れ線グラフでつないだものです。
拡大したものを以下に貼り付けます。
このようなグラフの極大値、極小値を抽出して、点と点をつなげたいです。
完成図は以下のような感じになることが予想されます。
ここにより詳細な情報を記載してください。