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

其他略。

相关推荐
共享家95271 天前
搭建 AI 聊天机器人:”我的人生我做主“
前端·javascript·css·python·pycharm·html·状态模式
Hgfdsaqwr1 天前
Python在2024年的主要趋势与发展方向
jvm·数据库·python
一晌小贪欢1 天前
Python 测试利器:使用 pytest 高效编写和管理单元测试
python·单元测试·pytest·python3·python测试
小文数模1 天前
2026年美赛数学建模C题完整参考论文(含模型和代码)
python·数学建模·matlab
ghie90901 天前
使用经验模态分解(EMD)处理振动信号,并结合样本熵进行特征提取
信号处理
Halo_tjn1 天前
基于封装的专项 知识点
java·前端·python·算法
Hgfdsaqwr1 天前
掌握Python魔法方法(Magic Methods)
jvm·数据库·python
weixin_395448911 天前
export_onnx.py_0130
pytorch·python·深度学习
s1hiyu1 天前
使用Scrapy框架构建分布式爬虫
jvm·数据库·python
2301_763472461 天前
使用Seaborn绘制统计图形:更美更简单
jvm·数据库·python