内插指的是增加数字信号采样率的方法。
步骤:
- 插入零值样本(即在所需要插值的地方插入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.plot