应用于导数的双线性变换(Tustin 方法)

信息处理 衍生物 双线性变换
2022-02-23 18:51:16

我希望我没有误解一些非常错误的东西,但是连续导数可以被认为是拉普拉斯空间中的传递函数,对吧?D=d/dtD(s)=s

因此,当我尝试使用双线性变换(Tustin 方法)对其进行离散化时,我很容易得到

D(z)=2T1z11+z1

当我将其应用于包含一个离散脉冲的序列时,响应以奈奎斯特频率振荡。附近的频谱是二次的,而不是像从导数中预期的编辑:后者只是由于舍入误差,因为低频幅度被奈奎斯特峰淹没ω=0iω

虽然我当然知道双线性变换在离散化各种过滤器方面的效果如何,但我很难理解,如果它对于像导数这样简单的东西似乎失败得如此惨烈,为什么它被认为是优越的,否则它可以是很容易用一阶有限差分表示:

D(z)=1z1T

甚至二阶(对称)有限差分

D(z)=1z22Tz1

我相信这一切都有一个非常简单的解释,但我看不到。

PS:更令人困惑的是:当我将双线性 D(z) 应用于阶跃函数时,结果(正确地)是一个峰值。因此,应用于脉冲的双线性 D(z) 的逆产生了阶跃函数,就像它必须的那样。那里发生了什么?

1个回答

如果双线性变换在尝试逼近导数时“失败得如此悲惨”,这是一个透视问题。首先,对于小近似值是二次的,这不是真的,它应该是线性的:ω

D(ejω)=2T1ejω1+ejω=2jTtan(ω2)jωTfor|ω|1

因此,对于小频率,离散时间系统表现为(离散时间)微分器。对于较大的频率,我们期望频率扭曲的影响,这正是这里发生的情况:无限连续时间频率范围映射到区间ω[0,π]

总之,双线性变换的行为与预期一致,但也很明显,与将双线性变换应用于连续时间传递函数相比,有更好的方法来逼近离散时间的导数。H(s)=s