8、数字滤波器与算法详解

数字滤波器与算法详解

1. 并行形式的状态空间实现

考虑传递函数，其可写为：
[H(z) = c_0 + \sum_{i = 1}^{N} \frac{c_i}{z - z_{p_i}}]
这里(z_p)是复极点。为简化，假设极点唯一，且部分分式展开形式如上式。选择状态变量为：
[V_i(z) = \frac{1}{z - z_{p_i}}X(z)]
由此可得：
[zV_i(z) = z_{p_i}X(z) + X(z)]
原传递函数可写成：
[Y(z) = c_0X(z) + \sum_{i = 1}^{N} c_iV_i(z)]
写成矩阵形式为：
[\begin{bmatrix}
v_1(nT + T) \
v_2(nT + T) \
\vdots \
v_N(nT + T)
\end{bmatrix} =
\begin{bmatrix}
z_{p_1} & 0 & \cdots & 0 \
0 & z_{p_2} & \cdots & 0 \
\vdots & \vdots & \ddots & \vdots \
0 & 0 & \cdots & z_{p_N}
\end{bmatrix}
\begin{bmatrix}
v_1(nT) \
v_2(nT) \
\vdots \
v_N(nT)
\end{bmatrix} +
\begin{bmatrix