Python如何找局部极值_scipy.signal.argrelextrema找波峰波谷

老宇同学_5665

老宇同学_5665

2026-03-31

1026人浏览

原创

scipy.signal.argrelextrema 默认要求严格大于/小于邻点,故平台段、噪声或高采样率下易漏峰;需先平滑、确保一维输入、显式指定axis和比较函数,并用order=3起调参;波峰波谷可统一用argrelextrema(-data)推导;索引须映射至原始x轴,注意边界截断。

python如何找局部极值_scipy.signal.argrelextrema找波峰波谷

scipy.signal.argrelextrema 为什么找不到明显的波峰?

默认只认「严格大于/小于邻点」的极值,如果信号里有平台段(连续相等值)、噪声毛刺或采样太密,argrelextrema 直接跳过——它不模糊判断,也不自动平滑。

实操建议:

  • 先用 scipy.signal.savgol_filter 或 np.convolve 做轻量平滑,尤其当原始数据带高频抖动时
  • 检查输入数组是否为一维 ndarray;传入 list 或二维数组会静默失败或报 ValueError: object too deep for desired array
  • 确认比较轴是 axis=0(默认),多维数组需显式指定,否则结果维度错乱

order 参数到底控制什么?不是“搜索窗口大小”

order 指的是在目标点左右各取 order 个点做比较,总跨度是 2 * order + 1。它不控制“局部范围”,而是硬性比较邻域内所有点——哪怕中间有个更陡的峰,只要没进这个固定窗口,就看不见。

常见错误现象:

  • 设 order=1 却想抓宽峰 → 峰顶稍平就漏掉
  • 设 order=50 处理长周期信号 → 计算变慢,且容易把缓坡误判为极值
  • 跨峰谷混合场景下,同一个 order 值对峰和谷敏感度不同(因符号差异)

建议:从 order=3 起试,再按信号特征微调;高频采样可略降,低频拉高。

怎么同时拿到波峰和波谷?别重复调用两次

直接复用同一组索引,靠符号翻转判断类型:对原数组取负后找极大值,等价于找原数组极小值。比两次调用 argrelextrema 更稳,避免因浮点误差或边界处理导致索引不一致。

macOS UI 自动化技能包 —— 用 Python 控制你的 Mac
macOS UI 自动化技能包 —— 用 Python 控制你的 Mac

让 macOS 自动化像说话一样简单 —— 用代码控制你的 Mac,无需学习 AppleScript 系统要求 macOS 14+ (Sonoma) Python 3.9+ 需要开启:辅助功能 + 屏幕录制权限

下载

示例逻辑:

import numpy as np
from scipy.signal import argrelextrema
<p>peaks = argrelextrema(data, np.greater, order=3)[0]
troughs = argrelextrema(-data, np.greater, order=3)[0]</p>

注意:np.greater 是必须显式传入的比较函数,不能省略;换成 np.less 反而会出反直觉结果。

argrelextrema 返回的索引怎么映射回原始 x 轴?

它只返回 y 值数组中的位置索引,不关心你原始 x 是时间、频率还是样本点编号。最容易被忽略的是:如果你的 x 轴非等距(比如 log 频率轴),这些索引不能直接当横坐标用。

正确做法:

  • 确保 x 和 y 长度一致,用 x[peaks]、x[troughs] 显式索引
  • 若 x 是 pd.Series,先转 np.array 再索引,否则可能触发 KeyError 或隐式对齐错误
  • 极值点靠近数组首尾时,order 过大会导致部分索引缺失(边界被截断),需手动补或提前裁剪

边界处理没标准答案,得看你的业务是否允许首尾半个周期参与判断。

Python免费学习笔记(深入):立即使用
在学习笔记中,你将探索 Python 的核心概念和高级技巧!

相关文章

PHP速学视频免费教程(入门到精通)
PHP速学视频免费教程(入门到精通)

PHP怎么学习?PHP怎么入门?PHP在哪学?PHP怎么学才快?不用担心,这里为大家提供了PHP速学教程(入门到精通),有需要的小伙伴保存下载就能学习啦!

下载

相关标签:

python

本站声明:本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系admin@php.cn

相关专题

更多
python打包成可执行文件
python打包成可执行文件

本专题为大家带来python打包成可执行文件相关的文章,大家可以免费的下载体验。

2023.07.20

1831

4

python能做什么
python能做什么

python能做的有:可用于开发基于控制台的应用程序、多媒体部分开发、用于开发基于Web的应用程序、使用python处理数据、系统编程等等。本专题为大家提供python相关的各种文章、以及下载和课程。

2023.07.25

4884

7

format在python中的用法
format在python中的用法

Python中的format是一种字符串格式化方法,用于将变量或值插入到字符串中的占位符位置。通过format方法,我们可以动态地构建字符串,使其包含不同值。php中文网给大家带来了相关的教程以及文章,欢迎大家前来阅读学习。

2023.07.31

1789

3

python教程
python教程

Python已成为一门网红语言,即使是在非编程开发者当中,也掀起了一股学习的热潮。本专题为大家带来python教程的相关文章,大家可以免费体验学习。

2023.08.03

29217

23

python环境变量的配置
python环境变量的配置

Python是一种流行的编程语言,被广泛用于软件开发、数据分析和科学计算等领域。在安装Python之后,我们需要配置环境变量,以便在任何位置都能够访问Python的可执行文件。php中文网给大家带来了相关的教程以及文章,欢迎大家前来学习阅读。

2023.08.04

3487

5

python eval
python eval

eval函数是Python中一个非常强大的函数,它可以将字符串作为Python代码进行执行,实现动态编程的效果。然而,由于其潜在的安全风险和性能问题,需要谨慎使用。php中文网给大家带来了相关的教程以及文章,欢迎大家前来学习阅读。

2023.08.04

3507

5

scratch和python区别
scratch和python区别

scratch和python的区别:1、scratch是一种专为初学者设计的图形化编程语言,python是一种文本编程语言;2、scratch使用的是基于积木的编程语法,python采用更加传统的文本编程语法等等。本专题为大家提供scratch和python相关的文章、下载、课程内容,供大家免费下载体验。

2023.08.11

1223

5

python合并两个列表
python合并两个列表

Python是一种强大的编程语言,具有许多方便的功能和工具。在Python中,有多种方法可以合并两个列表。php中文网给大家带来了相关的教程以及文章,欢迎大家前来学习阅读。

2023.08.10

616

4

python是前端还是后端
python是前端还是后端

Python属于前端也属于后端,其灵活性和丰富的生态系统使得开发人员能够在不同的领域中灵活运用。本专题为大家提供python相关的文章、下载、课程内容,供大家免费下载体验。

2023.08.11

2723

5

热门下载

更多
网站特效
/
网站源码
/
网站素材
/
前端模板

精品课程

更多
相关推荐
/
热门推荐
/
最新课程
Uvicorn官方安装与入门手册
Uvicorn官方安装与入门手册

共0课时 | 0人学习

Python uv官方项目管理指南
Python uv官方项目管理指南

共0课时 | 0人学习

Python uv官方安装与入门文档
Python uv官方安装与入门文档

共0课时 | 0人学习