上一篇仿真从零极点的位置,对量化误差进行了形象的说明。本篇将用频率幅度响应来进行演示。
需要注意的是,scipy库中的signal模块下的tf2pzk函数,要求分母分子同阶:
a = np.array([1, 1.7, 0.745]) # denominator
b = np.array([0.05, 0, 0]) # numerator
为了计算频率响应,我们使用signal中的freqz方法,其用法可以使用signal.freqz?进行查看。
w1, h1 = signal.freqz(b, a) # compute the frequency response
上面是理想情况下的频率响应,如果对系数进行量化,引入误差后,会造成曲线的漂移,严重情况下甚