使用python怎么計(jì)算波峰波谷值

使用python怎么計(jì)算波峰波谷值？針對這個(gè)問題，這篇文章詳細(xì)介紹了相對應(yīng)的分析和解答，希望可以幫助更多想解決這個(gè)問題的小伙伴找到更簡單易行的方法。

1. 首先可先選擇一個(gè)函數(shù)或者擬合一個(gè)函數(shù)，這里選擇擬合數(shù)據(jù)：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次多項(xiàng)式擬合
p1 = np.poly1d(z1) #多項(xiàng)式系數(shù)
print(p1) # 在屏幕上打印擬合多項(xiàng)式
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表極大值點(diǎn)的距離至少大于等于10個(gè)水平單位
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. 在可導(dǎo)的情形下，可以求導(dǎo)來求極值點(diǎn)，同時(shí)得到極大值和極小值點(diǎn)：np.polyder

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

此時(shí)：yyyd.r 即可就得導(dǎo)數(shù)為0的點(diǎn)，可以與上述的極大值點(diǎn)對應(yīng)比較

4. 直接函數(shù)分別求極大值和極小值：signal.argrelextrema 函數(shù)

print(yvals[signal.argrelextrema(yvals, np.greater)]) #極大值的y軸, yvals為要求極值的序列
print(signal.argrelextrema(yvals, np.greater)) #極大值的x軸
peak_ind = signal.argrelextrema(yvals,np.greater)[0] #極大值點(diǎn)，改為np.less即可得到極小值點(diǎn)
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) #極大值點(diǎn)
plt.plot(signal.argrelextrema(yvals,np.less)[0],yvals[signal.argrelextrema(yvals, np.less)],'+', markersize=10) #極小值點(diǎn)
plt.show()

Python主要用來做什么

Python主要應(yīng)用于：1、Web開發(fā)；2、數(shù)據(jù)科學(xué)研究；3、網(wǎng)絡(luò)爬蟲；4、嵌入式應(yīng)用開發(fā)；5、游戲開發(fā)；6、桌面應(yīng)用開發(fā)。

關(guān)于使用python怎么計(jì)算波峰波谷值問題的解答就分享到這里了，希望以上內(nèi)容可以對大家有一定的幫助，如果你還有很多疑惑沒有解開，可以關(guān)注億速云行業(yè)資訊頻道了解更多相關(guān)知識。