数字信号内插方法的python实现

最新推荐文章于 2024-05-29 10:48:51 发布

原创

最新推荐文章于 2024-05-29 10:48:51 发布 · 2.5k 阅读

7 ·

CC 4.0 BY-SA版权

文章标签：

#python #数字信号处理 #内插

本文介绍了数字信号内插的基本步骤，包括插入零值样本、低通滤波以及为应对幅值变化进行的信号幅度调整。通过n倍内插，需要在滤波后将每个信号值乘以n以保持原始幅值。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

内插指的是增加数字信号采样率的方法。
步骤：

插入零值样本（即在所需要插值的地方插入0）

低通滤波器滤波
–来自：理查德G.莱昂斯, D.李・富盖尔, RichardG.Lyons,等. 数字信号处理精要[M]. 机械工业出版社, 2016.

但是我发现这种方法会导致信号的幅值改变。n倍内插相当于把幅值缩小了n倍。所以应该加入第3个步骤：n倍内插在滤波之后信号的幅值乘n。

import numpy as np
import math
import scipy.signal as signal
import pylab as pl
import matplotlib.pyplot as plt

sampling_rate=1000 # 采样率设置为1000Hz，内插之前的采样率为333Hz
t1=np.arange(0, 1.0, 1.0/sampling_rate) # 时间长度为1s
x1=np.sin(15.4*np.pi*t1)*8.9+np.sin(31*np.pi*t1)*35.9+np.sin(56*np.pi*t1)*29.3
plt.figure(figsize=(16,4))
plt.plot(t1,x1)
plt.