查找等效数字滤波器

Question

我正在尝试为简单的 RC 滤波器找到等效的数字滤波器。波特图不对齐，我不知道为什么。

============================================= =======

import numpy as np
import scipy.signal as sig
import matplotlib.pyplot as plt

# Analog filter
tau = 0.001 # param
lti = sig.lti([1], [tau, 1])

# Equivalent digital filter
fs = 20000 # param

T = 1 / fs
wd = 1 / tau
wa = 2 / T * np.tan(wd * T / 2)
dtau = 1 / wa
dnum, dden = sig.bilinear([1], [dtau, 1], fs)
dlti = sig.dlti(dnum, dden)

w, mag, phase = sig.bode(lti)
dw, dmag, dphase = sig.dbode(dlti)

plt.figure()

plt.subplot(211)
plt.semilogx(w, mag)    # Bode magnitude plot
plt.semilogx(dw, dmag)

plt.subplot(212)
plt.semilogx(w, phase)  # Bode phase plot
plt.semilogx(dw, dphase)

plt.show()

Answer 1

sig.bode returns第一个参数频率数组[rad/s].

sig.dbode returns代替频率数组[rad/time_unit].

数字信号的时间单位是采样，所以需要乘以采样率才能将rad/sample转换为rad/second:

plt.subplot(211)
plt.semilogx(w, mag)    # Bode magnitude plot
plt.semilogx(dw * fs, dmag)

plt.subplot(212)
plt.semilogx(w, phase)  # Bode phase plot
plt.semilogx(dw * fs, dphase)

查找等效数字滤波器

Find equivalent digital filter

signal-processing

scipy