python 寻找离散序列极值点的方法
作者：weijifen000 发布时间：2021-03-13 15:26:15
标签：python,离散序列,极值点
使用 scipy.signal 的 argrelextrema 函数(API)，简单方便
import numpy as np
import pylab as pl
import matplotlib.pyplot as plt
import scipy.signal as signal
x=np.array([
0, 6, 25, 20, 15, 8, 15, 6, 0, 6, 0, -5, -15, -3, 4, 10, 8, 13, 8, 10, 3,
1, 20, 7, 3, 0 ])
plt.figure(figsize=(16,4))
plt.plot(np.arange(len(x)),x)
print x[signal.argrelextrema(x, np.greater)]
print signal.argrelextrema(x, np.greater)
plt.plot(signal.argrelextrema(x,np.greater)[0],x[signal.argrelextrema(x, np.greater)],'o')
plt.plot(signal.argrelextrema(-x,np.greater)[0],x[signal.argrelextrema(-x, np.greater)],'+')
# plt.plot(peakutils.index(-x),x[peakutils.index(-x)],'*')
plt.show()
[25 15 6 10 13 10 20]
(array([ 2, 6, 9, 15, 17, 19, 22]),)
但是存在一个问题，在极值有左右相同点的时候无法识别，但是个人认为在实际的使用过程中极少会出现这种情况，所以可以忽略。
x=np.array([
0, 15, 15, 15, 15, 8, 15, 6, 0, 6, 0, -5, -15, -3, 4, 10, 8, 13, 8, 10, 3,
1, 20, 7, 3, 0 ])
plt.figure(figsize=(16,4))
plt.plot(np.arange(len(x)),x)
print x[signal.argrelextrema(x, np.greater)]
print signal.argrelextrema(x, np.greater)
plt.plot(signal.argrelextrema(x,np.greater)[0],x[signal.argrelextrema(x, np.greater)],'o')
plt.plot(signal.argrelextrema(x,np.less)[0],x[signal.argrelextrema(x, np.less)],'+')
plt.show()
[15 6 10 13 10 20]
(array([ 6, 9, 15, 17, 19, 22]),)
来源：https://blog.csdn.net/weijifen000/article/details/80070520
0
投稿猜你喜欢
- 比如：我们导入了某个客户的资料，我们知道此客户的姓名是ZhangShan，我们想知道，在我们的业务数据库(eg:NorthWind)中，有哪
- 前端开发部门在壮大，人员在增加，技术在进步。CSSer虽然不能代表前端开发团队，但CSSer是前端开发中非常重要的组成部分。一、文档写程序的
本文主要包括三大方面，大家仔细学习。1、导航栏中的表单导航栏中的表单不是使用 Bootstrap 表单 章节中所讲到的默认的 class，它- 目录1、typing介绍2、typing的作用3、常用类型3.1 代码示例4、typing模块的其他用法4.1 类型别名4.2 NewTyp
- RGB图像转灰度图RGB图像转换为灰度图时通常使用：进行转换，以下尝试通过其他对图像像素操作的方式将RGB图像转换为灰度图像。#includ
- 本文在上文的基础上重新实现支持多线程的服务器。以下为TCP客户端的程序代码：#!/usr/bin/env python3import sys
- 本文研究的主要是Python生成器及其应用，具体如下。一、定义可以理解为一种数据类型，这种数据类型自动实现了迭代器协议（其他的数据类型需要调
- 写在之前「装饰器」作为 Python 高级语言特性中的重要部分，是修改函数的一种超级便捷的方式，适当使用能够有效提高代码的可读性和可维护性，
- 使用celery在django项目中实现异步发送短信在项目的目录下创建celery_tasks用于保存celery异步任务。在celery_
- 本文介绍了python同时替换多个字符串方法示例，分享给大家，具体如下：import rewords = ''' &
- 本文所述实例为Python处理文本文件并生成指定格式文件的方法，具体实现功能代码如下所示：import osimport sysimport
- Python提供了一个内联模块buildin，该模块定义了一些软件开发中经常用到的函数，利用这些函数可以实现数据类型的转换、数据的计算、序列
- 事务安全 transaction事务 transaction 访问可能更新数据库中各种数据项的一个程序执行单元unit事务由事务开始（beg
- 一、简介简单记录一下存储过程的使用。存储过程是预编译SQL语句集合，也可以包含一些逻辑语句，而且当第一次调用存储过程时，被调用的存储过程会放
- 本文介绍我使用QQ得到服务器上回传的python代码的探索历程，面向的对象是对计算机网络有一定了解的读者。期待有兴趣的人和我一起探讨！需求来
- Python的3.0版本，常被称为Python 3000，或简称Py3k。相对于Python的早期版本，这是一个较大的升级。为了不带入过多的
- PHP PDO 事务与自动提交现在通过 PDO 连接上了，在开始进行查询前，必须先理解 PDO 是如何管理事务的。事务支持四大特性（ACID
- 这篇文章记录一个采样器都随机地从原始的数据集中抽样数据。抽样数据采用permutation。 生成任意一个下标重排，从而利用下标来提取dat
- 一、功能简述番茄钟即番茄工作法，番茄工作法是简单易行的时间管理工具，使用番茄工作法即一个番茄时间共30分钟，25分钟工作,5分钟休息；特点一
- 之前已经写过一篇关于Django外键的文章，但是当时并没有介绍如何根据外键对数据的操作，也就是如何通过主表查询子表或者通过子表查询主表的信息
