python求極值點(波峰波谷)

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

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 ------

 


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

更多相關文章
  • 下麵的代碼主要用於使用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的哪個節點來處理,可以使用網關來進行路由轉發( ...
一周排行
  • 《ASP.NET MVC 企業級實戰》 [作者] (中) 鄒瓊俊[出版] 清華大學出版社[版次] 2017年04月 第1版[印次] 2019年08月 第6次 印刷[定價] 89.00元 【第01章】 (P021) 只有在 Lambda 有一個輸入參數時,括弧才是可選的,否則括弧是必需的。 使用空括弧 ...
  • 上一篇(https://www.cnblogs.com/meowv/p/12971041.html)使用HtmlAgilityPack抓取壁紙數據成功將圖片存入資料庫,本篇繼續來完成一個全網各大平臺的熱點新聞數據的抓取。 同樣的,可以先預覽一下我個人博客中的成品:https://meowv.com/ ...
  • 前言 請了一天假後回公司,同事跟我說使用Newtonsoft.json序列化TreeView對象的時候出現報錯; 啊!什麼?這個類庫不是能夠序列化所有東西嗎?真的很懵逼,也是我第一次使用這個類庫出現問題! 問題異常 異常信息 : Newtonsoft.Json.JsonSerializationEx ...
  • 簡單瞭解下麵詞語的意思 節點:二叉樹中每個元素都稱為節點 葉子節點(簡稱:葉子):度為0的節點,葉子節點就是樹中最底段的節點,葉子節點沒有子節點,也叫終端結點 分枝節點:度不為0的結點 節點的度:二叉樹的度代表某個節點的孩子或者說直接後繼的個數,簡單說就是一個節點擁有的子樹數 樹的度: 樹中最大的結 ...
  • C# 中的LINQ 提供了兩種操作方式,查詢表達式和查詢操作符,所有的查詢表達式都有對應的查操作符類替代,查詢表達式有點“類” SQL,在代碼中寫SQL,總覺得不夠“優雅”,使用查詢操作符就顯得“優雅”很多, 本系列就來對所有的LINQ 標準操作符進行一個全面的總結,這些操作符和我上篇文章總結的Rx ...
  • 在Startup ConfigureServices 註冊本地化所需要的服務AddLocalization和 Configure<RequestLocalizationOptions> public void ConfigureServices(IServiceCollection services ...
  • 為什麼需要持久化,以及Redis持久化的RDB方式在這篇文章講的已經很透徹了,足以弔打面試官了。而且此篇內容需要RDB文章的內容支持,所以建議先看下:看完這篇還不懂Redis的RDB持久化,你們來打我! 一、什麼是AOF 它也是Redis持久化的重要手段之一,aof->Append Only Fil ...
  • 先上圖: @IT程式猿 微博網友評論: @迢書:前同事的,親眼見過 @AvenGeeker:Bug 404 @科技州:這是要逼死強迫症 @小島一瞥:哈哈哈哈哈我老家的車 最後小編整理了一套技術資料不僅能精準消除技術盲點、累計面試經驗,更可以攻剋JVM、Spring、分散式、微服務等技術難題。 海量電 ...
  • 概括來說,分三步: 1,首先找到是哪個進程的CPU占有率飆到了100%。 2,根據進程號pid,定位到是哪個線程,找到對應線程的tid。 3,導出對應線程的dump日誌文件,分析日誌文件定位具體代碼。 要解決這個問題,你應該具備以下技能: 1,linux的top命令。 2,jvm監控工具jps。 3 ...
  • 寫在最後 程式員為何害怕【別人的代碼】呢?這讓我想起一個段子。 寫這段代碼時 只有上帝和我知道他是幹嘛的 現在 只有上帝知道了 別人的代碼,似乎總意味著冗長、晦澀、凌亂,給人一種不想靠近的感覺。搞笑的是,對於一些程式員而言,即使是自己的代碼,在一段時間之後自己再拿來看,也成了【別人的代碼】... 作 ...