GEQ matlab实现流程

GEQ：

均衡器通常被音频工程师和消费者用来调整音频的频率响应。例如，它们可以用来补偿扬声器引入的偏差，或者为歌曲添加低音。它们本质上是一组滤波器，设计用于提供自定义的整体频率响应。 图形均衡器相比较参数均衡器更简单。
八倍频程、二个三分之一倍频程和三分之一八倍频程已经成为基于人耳行为的均衡器的常见带宽。IS0 266:1997（E）、ANSI S1.11-2004和IEC 61672-1:2013等标准定义了倍频程和分数倍频程滤波器的中心频率。这只留下一个参数需要调整：滤波器增益。图形均衡器在使用标准中心频率和公共带宽的同时提供对增益参数的控制。
matlab里audiotoolbox有两种geq实现，它们在组成滤波器的排列上有所不同：一个使用一组平行的倍频程或分数倍频程滤波器，另一个使用级联的双二次滤波器。两种实现方式中的中心频率均遵循ANSI S1.11-2004标准。
（https://www.docin.com/p-485047366.html）

带并行滤波器的图形均衡器

构造图形均衡器的一种方法是并行放置一组带通滤波器。每个滤波器的带宽是倍频程或分数倍频程，它们的中心频率被设置为一起覆盖[20-20k]Hz的音频频率范围。

参数EQOrder是单个均衡器频带的顺序，指定为正偶数。所有均衡器频带具有相同的顺序。每个带的滤波器设计都遵循octaveFilter design method。
————后补(octaveFilter design)
如果设置SampleRate属性，使奈奎斯特频率（SampleRate/2）小于ANSI S1.11-2004标准定义的最终带通边缘，则：最后的带通滤波器是其上带通边缘小于奈奎斯特频率的滤波器。最终滤波器被实现为通过designParamEQ函数设计的高通滤波器。
这个GEQ系统的传递函数是所有分支传递函数的总和，公式：

由于增益与滤波器设计无关，因此调整增益不会产生显著的计算成本。并行滤波器结构非常适合于并行硬件实现。带通滤波器的幅度响应应在其带宽之外的所有其他频率处接近于零，以避免滤波器之间的相互作用。虽然实际上这是不可能实现的，因此滤波器和滤波器之间会有带间干扰。(符合滤波器设计时的衰减参数即认为没问题了)

实时计算过程：
输入信号被馈送到由M个滤波器组成的滤波器组中，其中M取决于指定的带宽和采样率属性。滤波器组的每个分支乘以增益属性的相应元素的线性形式。对分支进行求和，并返回输出信号。
matlab中使用时，用eq = graphicEQ(‘Structure’,‘Parallel’)即可

级联滤波器的图形均衡器

图形均衡器的另一种实现方式使用级联均衡滤波器（峰值或陷波），实现为biquad滤波器。均衡器的传递函数可以写成单个biquad的传递函数的乘积。
parallelGraphicEQ = graphicEQ(‘Structure’,‘Parallel’);
visualize(parallelGraphicEQ);

会注意到，幅度响应不是平坦的。这是因为滤波器是独立设计的，每个滤波器都有一个幅度响应下降的过渡宽度。此外，由于非理想阻带，存在从一个滤波器的阻带到相邻滤波器的通带的泄漏。泄漏可能导致实际增益与预期增益不同。(比如设定增益为10dB，输出很难准确的保证是10dB增益，但是中心频率会接近这个增益)
级联结构可以在一定程度上减轻这种情况，因为增益是滤波器设计中固有的。
把 cascade biquad的增益设置成0就可以bypass掉，由于这种类型的结构中没有分支，这意味着在输入和输出之间没有增益路径(滤波器直接把增益包含进去了)。graphicEQ默认实现级联结构。
级联设计的缺点是，无论何时相应的增益发生变化，都需要重新设计biquad滤波器的系数。这对于并行实现来说是不需要的，因为增益只是每个并行分支的倍数。带通滤波器的并联也避免了在级联中发现的相位误差和量化噪声的累积。

滤波器组设计：级联实现使用multibandParametricEQ-来设计图形均衡器滤波器组。
——后补(multibandParametricEQ)
输入信号被馈送到M个biquad滤波器的级联中，其中M取决于指定的带宽和采样率属性。

————后补程序：

——参考：matlab graphicEQ