[数字信号处理][Python] numpy.gradient()函数的算法实现

先看实例

python 复制代码
import numpy as np
signal = [3,2,1,3,8,10]
grad = np.gradient(signal)
print(grad)

输出结果是
[-1. -1. 0.5 3.5 3.5 2. ]

这个结果是怎么来的呢?

np.gradient 计算信号的数值梯度,也就是信号值的变化率。它使用中心差分法来计算中间点的梯度,并使用前向差分法和后向差分法来计算边界点的梯度。

一、边界点:
g r a d [ 0 ] = s i g n a l [ 1 ] − s i g n a l [ 0 ] = 2 − 3 = − 1 grad[0] = signal[1] - signal[0] = 2 - 3 = -1 grad[0]=signal[1]−signal[0]=2−3=−1
g r a d [ 5 ] = s i g n a l [ 5 ] − s i g n a l [ 4 ] = 10 − 8 = 2 grad[5] = signal[5] - signal[4] = 10-8 = 2 grad[5]=signal[5]−signal[4]=10−8=2

二、中间点:
g r a d [ i ] = s i g n a l [ i + 1 ] − s i g n a l [ i − 1 ] 2 grad[i] =\frac{signal[i+1]-signal[i-1]}{2} grad[i]=2signal[i+1]−signal[i−1],由此可知
g r a d [ 1 ] = s i g n a l [ 2 ] − s i g n a l [ 0 ] 2 = − 1 grad[1] =\frac{signal[2]-signal[0]}{2} = -1 grad[1]=2signal[2]−signal[0]=−1

其他略。

相关推荐
习习.y8 分钟前
关于python中的面向对象
开发语言·python
hmbbcsm17 分钟前
练习python题目小记(六)
开发语言·python
wow_DG26 分钟前
【Python✨】VS Code 秒开 Python 类型检查:一招 mypy + settings.json 让你的 Bug 原地现形!
python·json·bug
Aspect of twilight43 分钟前
LeetCode华为大模型岗刷题
python·leetcode·华为·力扣·算法题
空影星1 小时前
高效追踪电脑使用时间,Tockler助你优化时间管理
python·django·flask
LiLiYuan.1 小时前
【Lombok库常用注解】
java·开发语言·python
不去幼儿园2 小时前
【启发式算法】灰狼优化算法(Grey Wolf Optimizer, GWO)详细介绍(Python)
人工智能·python·算法·机器学习·启发式算法
二川bro2 小时前
数据可视化进阶:Python动态图表制作实战
开发语言·python·信息可视化
青青子衿_212 小时前
TikTok爬取——视频、元数据、一级评论
爬虫·python·selenium
徒慕风流2 小时前
GeoSight:基于 Open3D 与 PySide6 的参数化 3D 模型处理与实时点云监控工具
计算机视觉·3d·信号处理