[数字信号处理][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

其他略。

相关推荐
花好月圆春祺夏安1 小时前
基于odoo17的设计模式详解---装饰模式
数据库·python·设计模式
萧鼎2 小时前
深度探索 Py2neo:用 Python 玩转图数据库 Neo4j
数据库·python·neo4j
华子w9089258592 小时前
基于 Python Django 和 Spark 的电力能耗数据分析系统设计与实现7000字论文实现
python·spark·django
Rockson3 小时前
使用Ruby接入实时行情API教程
javascript·python
GG不是gg4 小时前
NumPy-核心函数np.matmul()深入解析
numpy
Tipriest_4 小时前
Python关键字梳理
python·关键字·keyword
im_AMBER5 小时前
学习日志05 python
python·学习
大虫小呓5 小时前
Python 处理 Excel 数据 pandas 和 openpyxl 哪家强?
python·pandas
哪 吒6 小时前
2025B卷 - 华为OD机试七日集训第5期 - 按算法分类,由易到难,循序渐进,玩转OD(Python/JS/C/C++)
python·算法·华为od·华为od机试·2025b卷
摸爬滚打李上进7 小时前
重生学AI第十六集:线性层nn.Linear
人工智能·pytorch·python·神经网络·机器学习