6、滤波器组优化与MDCT滤波器组实现

滤波器组优化与MDCT滤波器组实现

1. 滤波器组基础与优化

1.1 折叠矩阵与延迟矩阵乘法

首先，我们有一个折叠矩阵 Fa 示例：

Fa[:,:,0]
Out:
array([[ 0.        , 0.19509032, 0.98078528, 0.        ],
       [ 0.55557023, 0.        , 0.        , 0.83146961],
       [ 0.83146961, 0.        , -0.        , -0.55557023],
       [ 0.        , 0.98078528, -0.19509032, -0.        ]])

将这个折叠矩阵 Fa 与延迟矩阵 D(z) 相乘，操作如下：

from Dmatrix import *
from polmatmult import *
D = Dmatrix(N)
Faz = polmatmult(Fa, D)

相乘后得到的 Faz 不同维度的结果如下：

Faz[:,:,0]
Out:
array([[ 0.        , 0.        , 0.98078528, 0.        ],