时间:2021-05-22
python求极值点主要用到scipy库。
1. 首先可先选择一个函数或者拟合一个函数,这里选择拟合数据:np.polyfit
import pandas as pdimport matplotlib.pyplot as pltimport numpy as npfrom 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()总结
以上所述是小编给大家介绍的python计算波峰波谷值的方法(极值点),希望对大家有所帮助,也非常感谢大家对网站的支持!
声明:本页内容来源网络,仅供用户参考;我单位不保证亦不表示资料全面及准确无误,也不保证亦不表示这些资料为最新信息,如因任何原因,本网内容或者用户因倚赖本网内容造成任何损失或损害,我单位将不会负任何法律责任。如涉及版权问题,请提交至online#300.cn邮箱联系删除。
波峰绘制需要一点技巧,利用基本图形来绘制即可,今天我们就来看看使用ppt绘制波峰图形的教程。软件名称:office2019专业增强版激活版v19.0(附离线安装
本文实例讲述了Python实现计算圆周率π的值到任意位的方法。分享给大家供大家参考,具体如下:一、需求分析输入想要计算到小数点后的位数,计算圆周率π的值。二、算
本文实例讲述了python计算对角线有理函数插值的方法。分享给大家供大家参考。具体实现方法如下:'''p=rational(xData,yData,x)Eval
本文实例讲述了python实现两个dict合并与计算操作。分享给大家供大家参考,具体如下:用pythonic的方法,将两个dict合并,并进行计算.如果key值
本文实例讲述了PHP排序二叉树基本功能实现方法。分享给大家供大家参考,具体如下:这里演示了排序二叉树节点的插入,中序遍历,极值的查找和特定值的查找的功能.基本没