找到等效的数字滤波器

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在第一个参数频率数组[rad / s ] 中返回。

sig.dbode返回频率数组[rad / time_unit ] 。

数字信号的时间单位是样本，因此您需要乘以采样率将rad / sample转换为rad /秒：

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)