python计算波峰波谷值的方法(极值点)


Posted in Python onFebruary 18, 2020

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

得到的图形是:

python计算波峰波谷值的方法(极值点)

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

python计算波峰波谷值的方法(极值点)

3. 在可导的情形下,可以求导来求极值点,同时得到极大值和极小值点:np.polyder

yyyd = np.polyder(p1,1) # 1表示一阶导
print(yyyd)

此时:yyyd.r 即可就得导数为0的点,可以与上述的极大值点对应比较

python计算波峰波谷值的方法(极值点)

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计算波峰波谷值的方法(极值点)

总结

以上所述是小编给大家介绍的python计算波峰波谷值的方法(极值点),希望对大家有所帮助,也非常感谢大家对三水点靠木网站的支持!

Python 相关文章推荐
python下如何让web元素的生成更简单的分析
Jul 17 Python
详细介绍Python函数中的默认参数
Mar 30 Python
在Python中使用matplotlib模块绘制数据图的示例
May 04 Python
python使用Apriori算法进行关联性解析
Dec 21 Python
Python+Pandas 获取数据库并加入DataFrame的实例
Jul 25 Python
Python字典对象实现原理详解
Jul 01 Python
python判断一个对象是否可迭代的例子
Jul 22 Python
Python ORM编程基础示例
Feb 02 Python
python3实现往mysql中插入datetime类型的数据
Mar 02 Python
自定义Django_rest_framework_jwt登陆错误返回的解决
Oct 18 Python
Python词云的正确实现方法实例
May 08 Python
Pytorch中TensorBoard及torchsummary的使用详解
May 12 Python
Python表达式的优先级详解
Feb 18 #Python
使用Tkinter制作信息提示框
Feb 18 #Python
Python中import导入不同目录的模块方法详解
Feb 18 #Python
python调用HEG工具批量处理MODIS数据的方法及注意事项
Feb 18 #Python
将pytorch转成longtensor的简单方法
Feb 18 #Python
python实现查找所有程序的安装信息
Feb 18 #Python
Python3中configparser模块读写ini文件并解析配置的用法详解
Feb 18 #Python
You might like
PHP 数组实例说明
2008/08/18 PHP
PHP中::、->、self、$this几种操作符的区别介绍
2013/04/24 PHP
fetchAll()与mysql_fetch_array()的区别详解
2013/06/05 PHP
php 计算两个时间相差的天数、小时数、分钟数、秒数详解及实例代码
2016/11/09 PHP
thinkPHP内置字符串截取函数用法详解
2016/11/15 PHP
jQuery温习篇 强大的JQuery选择器
2010/04/24 Javascript
基于jquery实现的可以编辑选择的下拉框的代码
2010/11/19 Javascript
jQuery滚动条插件nanoscroller使用指南
2015/04/21 Javascript
一起学写js Calender日历控件
2016/04/14 Javascript
把多个JavaScript函数绑定到onload事件处理函数上的方法
2016/09/04 Javascript
JS双击变input框批量修改内容
2016/12/12 Javascript
js实现的简练高效拖拽功能示例
2016/12/21 Javascript
详解webpack打包vue时提取css
2017/05/26 Javascript
小程序实现多选框功能
2018/10/30 Javascript
react高阶组件添加和删除props
2019/04/26 Javascript
微信小程序 如何获取网络状态
2019/07/26 Javascript
重置Redux的状态数据的方法实现
2019/11/18 Javascript
微信小程序通过websocket实时语音识别的实现代码
2020/08/19 Javascript
JSON stringify方法原理及实例解析
2020/10/23 Javascript
vue created钩子函数与mounted钩子函数的用法区别
2020/11/05 Javascript
Python中用format函数格式化字符串的用法
2015/04/08 Python
Python利用pandas计算多个CSV文件数据值的实例
2018/04/19 Python
Python3.5文件修改操作实例分析
2019/05/01 Python
如何通过50行Python代码获取公众号全部文章
2019/07/12 Python
Django-Model数据库操作(增删改查、连表结构)详解
2019/07/17 Python
如何通过命令行进入python
2020/07/06 Python
HTML5 层的叠加的实现
2020/07/07 HTML / CSS
德国健康生活方式网上商店:Landkaufhaus Mayer
2019/03/12 全球购物
小区门卫岗位职责
2013/12/31 职场文书
公司经理聘任书
2014/03/29 职场文书
4s店活动策划方案
2014/08/25 职场文书
2016年庆“七一”主题党日活动总结
2016/04/05 职场文书
python基于tkinter制作m3u8视频下载工具
2021/04/24 Python
python脚本框架webpy的url映射详解
2021/11/20 Python
Python实现抖音热搜定时爬取功能
2022/03/16 Python
青岛市的收音机研制与生产
2022/04/07 无线电