数字信号处理--双线性变换法

双线性变换使用双线性映射将 s-平面映射到 z-平面，这样可以避免频率混叠。通常被称为Tustin方法。

具体步骤：

1. 从模拟域传递函数 H(s)开始。
2. 用s = ( 2 / T ) * ( ( 1 - ( z^-1 ) ) / ( 1 + ( z^-1 ) )替换s。
3. 简化表达式，得到离散化后的传递函数 H(z)。

传递函数的频率相应可以通过将z替换成e^(jωT)来计算。其中，ω是连续角频率（单位为弧度/秒），T是采样周期。
传递函数H(z)替换之后，变成H(e^(jωT))。其中|H(e^(jωT))|是幅频响应。arctan(H(e^(jωT)))也就是H(e^(jωT))的相角为相位响应。
特别注意：网上部分资料用z=e^(jω)，此时ω=2πf⋅Ts​=2πf/Fs代表数字频率（单位为弧度/样本）。万不可混淆。

频响计算例1

题目

计算传递函数H(z)=1/(1 - 0.3 * z^（-1) )。采样频率Fs=8000Hz。求1000Hz频响和相位。