全通滤波的级联实现EQ filter(一)

2023-10-22 21:30

本文主要是介绍全通滤波的级联实现EQ filter(一),希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!



上一篇文章引入了dynamic EQ这个概念,就是动态EQ的意思,就是EQfilter的系数是按照sample点实时更新的。目前EQ filter都是使用biquardfilter type包括lowpass/highpass/bandpass/lowshelf/highshelf/peaking filter,每种滤波器有cutofffrequency/Q/gain这三个参数可调。通常设置好各段EQ的参数后,初始化的时候统一计算各个滤波器的系数,之后EQ参数都是固定的滤波,这叫静态的EQ。动态EQ则会随着信号的大小实时的更新滤波器系数来达到不同信号幅度不同的滤波效果,如果按照业内的EQ filter计算公式,运算量会比较大,当EQ段数增加时更不可接受,而且我们发现动态EQ只改变了EQ filtergain参数。

最开始在dsp-related论坛上看到有人说可以采用全通滤波器级联的方式来实现EQ filter,后面发现自己经常翻的一本书DAFx[1]第二章里面有详细介绍,并查阅到文献[2]详细分析了全通滤波器级联可实现1阶的lowshelf/highshelf filter2阶的peakingfilter,大大简化了EQ filter的滤波器系数的计算。下面分析这两种滤波器怎样用全通级联来实现。

  1. 一阶lowshelf/highshelf

首先以lowshelf为例,文献[2]的结论得知lowshelffilter的系统函数可以用如下全通滤波器通过简单的级联即可实现:

         上式中K表示滤波器的增益,K2表示6dB增益,而全通滤波器的系数计算只需要计算a

         上式中Ω表示归一化截止频率,对应lowshelffilter的截止频率,a固定的情况下,全通滤波器A(z)固定,滤波的过程无需改变,滤波器增益的改变只需要改变K即可。如下图为a固定的情况下,改变K的滤波器频响图:

         上图发现确实K值对应了lowshelffilter的最大增益,但发现当增益为衰减的时候,频响有些不正常,如果K设置为20.5频响应该关于横轴对称。

         文献[1]中,一阶lowshelffilter用全通滤波实现的公式如下:

其中H0表示如下:

将上式带入(2.40)可得:

上式和文献[2]F(z)完全一样,其中V0对应于K,但是全通滤波的系数如下,B表示boostC表示cut

         上式中boost情况的系数和a一样,而cut的情况就不能和boost使用一样的公式计算全通滤波器的系数了,所以文献[2]boostcut的频响并不是横轴对称的,而下图文献[2]中是完全对称的:

因此,对于boostcutlowshelffilter计算还不一样,boost中全通滤波的计算不用考虑增益这个参数,所以gain更新全通无需更新;cut情况的滤波过程gain更新了还要更新全通滤波器的系数,当然不更新全通滤波器也可以,只是滤波器频响更陡一些。

 

参考文献:

[1] U. Z¨olzer.Digital Audio Signal Processing. John Wiley & Sons, Ltd, 2nd edition, 2005.

[2]PHILLIP A. REGALIA and SANJIT K. MITRA. Tunable Digital Frequency ResponseEqualization Filters.

这篇关于全通滤波的级联实现EQ filter(一)的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/263703

相关文章

SpringBoot集成redisson实现延时队列教程

《SpringBoot集成redisson实现延时队列教程》文章介绍了使用Redisson实现延迟队列的完整步骤,包括依赖导入、Redis配置、工具类封装、业务枚举定义、执行器实现、Bean创建、消费... 目录1、先给项目导入Redisson依赖2、配置redis3、创建 RedissonConfig 配

Python的Darts库实现时间序列预测

《Python的Darts库实现时间序列预测》Darts一个集统计、机器学习与深度学习模型于一体的Python时间序列预测库,本文主要介绍了Python的Darts库实现时间序列预测,感兴趣的可以了解... 目录目录一、什么是 Darts?二、安装与基本配置安装 Darts导入基础模块三、时间序列数据结构与

Python使用FastAPI实现大文件分片上传与断点续传功能

《Python使用FastAPI实现大文件分片上传与断点续传功能》大文件直传常遇到超时、网络抖动失败、失败后只能重传的问题,分片上传+断点续传可以把大文件拆成若干小块逐个上传,并在中断后从已完成分片继... 目录一、接口设计二、服务端实现(FastAPI)2.1 运行环境2.2 目录结构建议2.3 serv

C#实现千万数据秒级导入的代码

《C#实现千万数据秒级导入的代码》在实际开发中excel导入很常见,现代社会中很容易遇到大数据处理业务,所以本文我就给大家分享一下千万数据秒级导入怎么实现,文中有详细的代码示例供大家参考,需要的朋友可... 目录前言一、数据存储二、处理逻辑优化前代码处理逻辑优化后的代码总结前言在实际开发中excel导入很

SpringBoot+RustFS 实现文件切片极速上传的实例代码

《SpringBoot+RustFS实现文件切片极速上传的实例代码》本文介绍利用SpringBoot和RustFS构建高性能文件切片上传系统,实现大文件秒传、断点续传和分片上传等功能,具有一定的参考... 目录一、为什么选择 RustFS + SpringBoot?二、环境准备与部署2.1 安装 RustF

Nginx部署HTTP/3的实现步骤

《Nginx部署HTTP/3的实现步骤》本文介绍了在Nginx中部署HTTP/3的详细步骤,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友们下面随着小编来一起学... 目录前提条件第一步:安装必要的依赖库第二步:获取并构建 BoringSSL第三步:获取 Nginx

MyBatis Plus实现时间字段自动填充的完整方案

《MyBatisPlus实现时间字段自动填充的完整方案》在日常开发中,我们经常需要记录数据的创建时间和更新时间,传统的做法是在每次插入或更新操作时手动设置这些时间字段,这种方式不仅繁琐,还容易遗漏,... 目录前言解决目标技术栈实现步骤1. 实体类注解配置2. 创建元数据处理器3. 服务层代码优化填充机制详

Python实现Excel批量样式修改器(附完整代码)

《Python实现Excel批量样式修改器(附完整代码)》这篇文章主要为大家详细介绍了如何使用Python实现一个Excel批量样式修改器,文中的示例代码讲解详细,感兴趣的小伙伴可以跟随小编一起学习一... 目录前言功能特性核心功能界面特性系统要求安装说明使用指南基本操作流程高级功能技术实现核心技术栈关键函

Java实现字节字符转bcd编码

《Java实现字节字符转bcd编码》BCD是一种将十进制数字编码为二进制的表示方式,常用于数字显示和存储,本文将介绍如何在Java中实现字节字符转BCD码的过程,需要的小伙伴可以了解下... 目录前言BCD码是什么Java实现字节转bcd编码方法补充总结前言BCD码(Binary-Coded Decima

SpringBoot全局域名替换的实现

《SpringBoot全局域名替换的实现》本文主要介绍了SpringBoot全局域名替换的实现,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友们下面随着小编来一... 目录 项目结构⚙️ 配置文件application.yml️ 配置类AppProperties.Ja