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
投稿猜你喜欢
- 在学习asp过程中相信很多初学者对Sub与Function的用法有些疑惑，好像它们没什么区别都可以使用。呵呵，看了本文的介绍您就可以了解了S
安装laravel框架命令行cd进入指定目录下，执行composer create-project --prefer-dist larave- MongoDB是一个文档型数据库，是NOSQL家族中最重要的成员之一，以下代码封装了MongoDB的基本操作。MongoDBConfig.j
- uniapp页面跳转的几种方式一、uni.navigateTo定义：保留当前页面，跳转到应用内的某个页面，使用uni.navigateBac
- 下面是代码，如果看不懂，建议先把表格的一些<tr><td>的表格原理弄清楚了，就可以了代码如下:<table&
- 我想把本篇作为css基础教程的序曲，从今天开始翻译和整理我从设计网页以来所学到的css基础知识。本教程会分成N个部分，单独发表，其间可能会插
- 正题： 1.1 javascript的灵活性 面向对象对象的Javascript编程模式：1、可以保存状态 2、具有对象内部才能调用的方法
- 先 Create table 吧create table emp(id int not null primary key,name varc
- 1、问题现象：有个网站是通过session验证的，前端时间访问正常，但近期后台session老是失效很快，根本没法操作，我设置timeout
- 疫情肆虐，憋在家实在无聊，索性写点东西，于是就有了这个极极极极极简的音乐播放器。这个极极极简的音乐播放器类似于“阅后即焚”的软件，播放器可以
- 所谓“评论”不是必须得有文本框，用户录入提交数据才算，广义上的评论包括用户的任何“表态”，典型如打分。我接触电子商务时间并不长，最早应该是0
- 什么是deferdefer用来声明一个延迟函数，把这个函数放入到一个栈上， 当外部的包含方法return之前，返回参数到调用方法之前调用，也
- 1、创建存储过程 create or replace procedure test(var_name_1 in type,var_name_
- 现在jquery应用的越来越多, 有些同学在享受爽快淋漓coding时就将性能问题忽略了, 比如我. jquery虽
- 你可能在使用MySQL过程中，各种意外导致数据库表的损坏，而且这些数据往往是最新的数据，通常不可能在备份数据中找到。本章将讲述如何检测MyS
- 我们在切换选项卡的时候，如果使用的是ajax技术，会碰到如下情况：点击tab1选项，服务器发出一个Ajax请求获取该选项tab1的内容数据。
- PHP fprintf() 函数实例把一些文本写入到名为 "test.txt" 的文本文件：<?php $numb
- 阅读上一篇：WEB2.0网页制作标准教程（11）不用表格的菜单辛苦了好多天，我们努力学习使用XHTML+CSS来重新设计我们的网站。那么我们
- 方法一：在php中，抓取https的网站，提示如下的错误内容：Warning: file_get_contents() [function.
- 语法： text-overflow : clip | ellipsis 参数： clip : 不显示省略标记（...），而是简单的裁切 el
