python求極值點(波峰波谷)

来源:https://www.cnblogs.com/qi-yuan-008/archive/2020/02/17/12323710.html
-Advertisement-
Play Games

python求極值點主要用到scipy庫。 1. 首先可先選擇一個函數或者擬合一個函數,這裡選擇擬合數據:np.polyfit import pandas as pd import matplotlib.pyplot as plt import numpy as np from scipy impo ...


python求極值點主要用到scipy庫。

1. 首先可先選擇一個函數或者擬合一個函數,這裡選擇擬合數據:np.polyfit

import pandas as pd
import matplotlib.pyplot as plt
import numpy as np
from scipy import signal   #濾波等

xxx = np.arange(0, 1000)
yyy = np.sin(xxx*np.pi/180)

z1 = np.polyfit(xxx, yyy, 7) # 用7次多項式擬合
p1 = np.poly1d(z1) #多項式繫數
print(p1) # 在屏幕上列印擬合多項式
yvals=p1(xxx) 

plt.plot(xxx, yyy, '*',label='original values')
plt.plot(xxx, yvals, 'r',label='polyfit values')
plt.xlabel('x axis')
plt.ylabel('y axis')
plt.legend(loc=4)
plt.title('polyfitting')
plt.show()

得到的圖形是:

 

2. 求波峰值,也就是極大值,得到:signal.find_peaks

# 極值
num_peak_3 = signal.find_peaks(yvals, distance=10) #distance表極大值點的距離至少大於等於10個水平單位
print(num_peak_3[0])
print('the number of peaks is ' + str(len(num_peak_3[0])))
plt.plot(xxx, yyy, '*',label='original values')
plt.plot(xxx, yvals, 'r',label='polyfit values')
plt.xlabel('x axis')
plt.ylabel('y axis')
plt.legend(loc=4)
plt.title('polyfitting')
for ii in range(len(num_peak_3[0])):
    plt.plot(num_peak_3[0][ii], yvals[num_peak_3[0][ii]],'*',markersize=10)
plt.show()

 

3. 在可導的情形下,可以求導來求極值點,同時得到極大值和極小值點:np.polyder

yyyd = np.polyder(p1,1) # 1表示一階導
print(yyyd)

此時:yyyd.r 即可就得導數為0的點,可以與上述的極大值點對應比較

 

4. 直接函數分別求極大值和極小值:signal.argrelextrema 函數

print(yvals[signal.argrelextrema(yvals, np.greater)]) #極大值的y軸, yvals為要求極值的序列
print(signal.argrelextrema(yvals, np.greater))  #極大值的x軸
peak_ind = signal.argrelextrema(yvals,np.greater)[0] #極大值點,改為np.less即可得到極小值點
plt.plot(xxx, yyy, '*',label='original values')
plt.plot(xxx, yvals, 'r',label='polyfit values')
plt.xlabel('x axis')
plt.ylabel('y axis')
plt.legend(loc=4)
plt.title('polyfitting')
plt.plot(signal.argrelextrema(yvals,np.greater)[0],yvals[signal.argrelextrema(yvals, np.greater)],'o', markersize=10)  #極大值點
plt.plot(signal.argrelextrema(yvals,np.less)[0],yvals[signal.argrelextrema(yvals, np.less)],'+', markersize=10)  #極小值點
plt.show()

## ----- end ------

 


您的分享是我們最大的動力!

-Advertisement-
Play Games
更多相關文章
  • 下麵的代碼主要用於使用python語言調用NASA官方的MODIS處理工具HEG進行投影坐標轉換與重採樣批量處理 主要參考 1. HEG的用戶手冊:https://newsroom.gsfc.nasa.gov/sdptoolkit/HEG/HEG215/EED2 TP 030_Rev01_HEG_U ...
  • 內部類的訪問規則 + 內部類可以直接訪問外部類中的成員,包括私有成員。 因為內部類中持有了一個外部類的引用,格式為:外部類名.this + 外部類要訪問內部類,必須要建立內部對象。 運行結果: 訪問 1. 當內部類定義在外部類的成員位置上,而且非私有,可以在外部其他類中直接建立內部類對象。 格式: ...
  • 消費者和消費者組 如何創建消費者 如何消費消息 消費者配置 提交和偏移量 再均衡 結束消費 ...
  • 一、正面回答 潛移默化中改變 二、日常成長 關心造輪子和用輪子,對於程式員而言,任何一件需要重覆做兩次的事情,就可以考慮通過寫程式把它自動化處理,所以程式員們不太能夠忍受重覆性高的工作。以Java為例,先後有ssh、ssm、spring boot、spring cloud、dubbo等,相關的一些任 ...
  • 1、什麼是樹? 什麼是樹?就是又數據組成的一種數據結構,按照平面展示出來就像一顆倒過來的樹。 2、樹的基本概念 1、樹的節點稱呼:節點、根節點(只有一個,就是第一個)、父節點、子節點、兄弟節點(兄弟節點只能是從同一個父節點的才能稱之為兄弟節點)。 2、一顆樹可以沒有任何節點,稱之為空樹。 3、一棵樹 ...
  • user 資料庫用戶名,用於連接資料庫……轉載自https://blog.csdn.net/wfanking/article/details/95504879 ...
  • JavaWeb Java Web 1、基本概念 1.1、前言 web開發: web,網頁的意思 , www.baidu.com 靜態web html,css 提供給所有人看的數據始終不會發生變化! 動態web 淘寶,幾乎是所有的網站; 提供給所有人看的數據始終會發生變化,每個人在不同的時間,不同的地 ...
  • 網關的概念 服務A、B都是暴露出來,供外部直接調用的, 有時候需要對請求進行過濾、校驗,比如檢驗用戶是否已登陸,可以寫在暴露出來的每個服務中,但要在多個服務中寫相同的代碼,太繁瑣,可以提出來,放在網關中。 如果A、B進行集群,需要負載均衡來確定使用A|B的哪個節點來處理,可以使用網關來進行路由轉發( ...
一周排行
    -Advertisement-
    Play Games
  • Timer是什麼 Timer 是一種用於創建定期粒度行為的機制。 與標準的 .NET System.Threading.Timer 類相似,Orleans 的 Timer 允許在一段時間後執行特定的操作,或者在特定的時間間隔內重覆執行操作。 它在分散式系統中具有重要作用,特別是在處理需要周期性執行的 ...
  • 前言 相信很多做WPF開發的小伙伴都遇到過表格類的需求,雖然現有的Grid控制項也能實現,但是使用起來的體驗感並不好,比如要實現一個Excel中的表格效果,估計你能想到的第一個方法就是套Border控制項,用這種方法你需要控制每個Border的邊框,並且在一堆Bordr中找到Grid.Row,Grid. ...
  • .NET C#程式啟動閃退,目錄導致的問題 這是第2次踩這個坑了,很小的編程細節,容易忽略,所以寫個博客,分享給大家。 1.第一次坑:是windows 系統把程式運行成服務,找不到配置文件,原因是以服務運行它的工作目錄是在C:\Windows\System32 2.本次坑:WPF桌面程式通過註冊表設 ...
  • 在分散式系統中,數據的持久化是至關重要的一環。 Orleans 7 引入了強大的持久化功能,使得在分散式環境下管理數據變得更加輕鬆和可靠。 本文將介紹什麼是 Orleans 7 的持久化,如何設置它以及相應的代碼示例。 什麼是 Orleans 7 的持久化? Orleans 7 的持久化是指將 Or ...
  • 前言 .NET Feature Management 是一個用於管理應用程式功能的庫,它可以幫助開發人員在應用程式中輕鬆地添加、移除和管理功能。使用 Feature Management,開發人員可以根據不同用戶、環境或其他條件來動態地控制應用程式中的功能。這使得開發人員可以更靈活地管理應用程式的功 ...
  • 在 WPF 應用程式中,拖放操作是實現用戶交互的重要組成部分。通過拖放操作,用戶可以輕鬆地將數據從一個位置移動到另一個位置,或者將控制項從一個容器移動到另一個容器。然而,WPF 中預設的拖放操作可能並不是那麼好用。為瞭解決這個問題,我們可以自定義一個 Panel 來實現更簡單的拖拽操作。 自定義 Pa ...
  • 在實際使用中,由於涉及到不同編程語言之間互相調用,導致C++ 中的OpenCV與C#中的OpenCvSharp 圖像數據在不同編程語言之間難以有效傳遞。在本文中我們將結合OpenCvSharp源碼實現原理,探究兩種數據之間的通信方式。 ...
  • 一、前言 這是一篇搭建許可權管理系統的系列文章。 隨著網路的發展,信息安全對應任何企業來說都越發的重要,而本系列文章將和大家一起一步一步搭建一個全新的許可權管理系統。 說明:由於搭建一個全新的項目過於繁瑣,所有作者將挑選核心代碼和核心思路進行分享。 二、技術選擇 三、開始設計 1、自主搭建vue前端和. ...
  • Csharper中的表達式樹 這節課來瞭解一下表示式樹是什麼? 在C#中,表達式樹是一種數據結構,它可以表示一些代碼塊,如Lambda表達式或查詢表達式。表達式樹使你能夠查看和操作數據,就像你可以查看和操作代碼一樣。它們通常用於創建動態查詢和解析表達式。 一、認識表達式樹 為什麼要這樣說?它和委托有 ...
  • 在使用Django等框架來操作MySQL時,實際上底層還是通過Python來操作的,首先需要安裝一個驅動程式,在Python3中,驅動程式有多種選擇,比如有pymysql以及mysqlclient等。使用pip命令安裝mysqlclient失敗應如何解決? 安裝的python版本說明 機器同時安裝了 ...