CN105656448A - 一种音频滤波器 - Google Patents
一种音频滤波器 Download PDFInfo
- Publication number
- CN105656448A CN105656448A CN201511018313.0A CN201511018313A CN105656448A CN 105656448 A CN105656448 A CN 105656448A CN 201511018313 A CN201511018313 A CN 201511018313A CN 105656448 A CN105656448 A CN 105656448A
- Authority
- CN
- China
- Prior art keywords
- omega
- wave filter
- phi
- theta
- psi
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03H—IMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
- H03H17/00—Networks using digital techniques
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03H—IMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
- H03H17/00—Networks using digital techniques
- H03H2017/0072—Theoretical filter design
Landscapes
- Soundproofing, Sound Blocking, And Sound Damping (AREA)
Abstract
本发明公开了一种音频滤波器,所述音频滤波器由N个分数阶Gammatone滤波器构成;所述音频滤波器的冲击响应为h(t),hi(t)为第i个Gammatone滤波器的冲击响应,其中,Ai,mi,bi,fi和φi分别表示第i个Gammatone滤波器冲击响应的幅度、阶数、衰减因子、频率和相位,u(t)代表阶跃函数;mi为分数。所述音频滤波器在时域和频域上更为集中,具有滤波效果好、性能稳定,对音频有更好的降噪效果。
Description
技术领域
本发明涉及一种音频滤波器。
背景技术
连续时间滤波器广泛应用于工程技术中,比如A/D、音频信号处理系统等,但是其稳定性需要通过计算冲击响应的绝对值积分是否有界来判别。连续时间有理数滤波器稳定性容易确定,在设计和应用中广受欢迎,但是AlanV.Oppenheim等人在《SignalsandSystems》中指出,由于连续时间有理数滤波器的频率响应由频率复数多项式构成,功率函数在频域上不够集中,滤波效果比连续时间无理数滤波器要差。2008年StefanStrahl等人在BrainResearch上发表题为“SparsegammatonesignalmodeloptimizedforEnglishspeechdoesnotmatchthehumanauditoryfilters”的文章,提出了连续时间整数阶Gammatone滤波器,其频率响应在频域上更为集中,因而广泛应用于音频信号处理,但是单独的连续时间整数阶Gammatone滤波器仍未能取得较好的频率选择性。
发明内容
本发明针对上述现有滤波器技术的不足,提出一种音频滤波器,在时域和频域上更为集中,具有滤波效果好、性能稳定,对音频有更好的降噪效果。
所述音频滤波器由N个分数阶Gammatone滤波器构成;所述音频滤波器的冲击响应为h(t),hi(t)为第i个Gammatone滤波器的冲击响应,其中,Ai,mi,bi,fi和φi分别表示第i个Gammatone滤波器冲击响应的幅度、阶数、衰减因子、频率和相位,u(t)代表阶跃函数;mi为分数。
进一步地,第i个Gammatone滤波器的频率响应为Hi(ω),令p为微分阶数,p为分数;
Dp为微分算子,
Γ(z)代表Gamma函数,
则第i个Gammatone滤波器的频率响应为
进一步地,Ai,mi,bi,fi和φi分别可以优化设计。
本发明的音频滤波器的冲击响应由多项式函数、指数函数、余弦函数构成,响应在时频域中更为集中,频率选择性好;此滤波器冲击响应是指数下降的,满足绝对可积的条件,具有稳定性;此滤波器冲击响应中的参数,Ai,mi,bi,fi和φi分别可以优化设计,具有灵活性。本发明的音频滤波器可以实现更好的音频降噪效果。
附图说明
图1是本发明的音频滤波器的频率响应图。
图2是原始飞机噪声信号的频谱图;
图3是采用Butterworth滤波器对图2的飞机噪声信号滤波后的频谱图;
图4是采用本发明的音频滤波器对图2的飞机噪声信号滤波后的频谱图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
所述音频滤波器由N个分数阶Gammatone滤波器构成;所述音频滤波器的冲击响应为h(t),hi(t)为第i个Gammatone滤波器的冲击响应,其中,Ai,mi,bi,fi和φi分别表示第i个Gammatone滤波器冲击响应的幅度、阶数、衰减因子、频率和相位,u(t)代表阶跃函数;mi为分数。
进一步地,第i个Gammatone滤波器的频率响应为Hi(ω),令p为微分阶数,p为分数;
Dp为微分算子,
Γ(z)代表Gamma函数,
则第i个Gammatone滤波器的频率响应为
假设:
vi=[xiyi]T
θi=[mibifi]T
可得:
假设:
Ri=xiψa(θi)+xiψb(θi)+yiψc(θi)-yiψd(θi)
Ii=-xiψc(θi)-xiψd(θi)+yiψa(θi)-yiψb(θi)
可得:
Hi(ω)=Ricosαi-Iisinαi+j(Risinαi+Iicosαi)
假设:
ψ1(θi)=[ψa(θi)+ψb(θi)ψc(θi)-ψd(θi)]T
ψ2(θi)=[-ψc(θi)-ψd(θi)ψa(θi)-ψb(θi)]T
可得:
假设:
从而求得:
音频滤波器的频率响应:
音频滤波器的幅频响应的平方:
音频滤波器的相频响应:
进一步地,Ai,mi,bi,fi和φi分别可以优化设计。
如图1所示,其是本发明的音频滤波器的频率响应图。
请对比图2、图3与图4,图2是原始飞机噪声信号的频谱图,图3是采用Butterworth滤波器对图2的飞机噪声信号滤波后的频谱图,图4是采用本发明的音频滤波器对图2的飞机噪声信号滤波后的频谱图。
可见,采用本发明的音频滤波器对噪声滤波后,滤波效果更好。
本发明的音频滤波器的冲击响应由多项式函数、指数函数、余弦函数构成,响应在时频域中更为集中,频率选择性好;此滤波器冲击响应是指数下降的,满足绝对可积的条件,具有稳定性;此滤波器冲击响应中的参数,Ai,mi,bi,fi和φi分别可以优化设计,具有灵活性。本发明的音频滤波器可以实现更好的音频降噪效果。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程,是可以通过计算机程序来指令相关的硬件来完成,所述的程序可存储于一计算机可读取存储介质中,该程序在执行时,可包括如上述各方法的实施例的流程。其中,所述的存储介质可为磁碟、光盘、只读存储记忆体(Read-OnlyMemory,ROM)或随机存储记忆体(RandomAccessMemory,RAM)等。
以上所述是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也视为本发明的保护范围。
Claims (3)
1.一种音频滤波器,其特征在于,所述音频滤波器由N个分数阶Gammatone滤波器构成;所述音频滤波器的冲击响应为h(t),hi(t)为第i个Gammatone滤波器的冲击响应,其中,Ai,mi,bi,fi和φi分别表示第i个Gammatone滤波器冲击响应的幅度、阶数、衰减因子、频率和相位,u(t)代表阶跃函数;mi为分数。
2.如权利要求1所述的音频滤波器,其特征在于,第i个Gammatone滤波器的频率响应为Hi(ω),令p为微分阶数,p为分数;
Dp为微分算子,
Γ(z)代表Gamma函数,
则第i个Gammatone滤波器的频率响应为
3.如权利要求1或2所述的音频滤波器,其特征在于,Ai,mi,bi,fi和φi分别可以优化设计。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201511018313.0A CN105656448B (zh) | 2015-12-28 | 2015-12-28 | 一种音频滤波器 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201511018313.0A CN105656448B (zh) | 2015-12-28 | 2015-12-28 | 一种音频滤波器 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105656448A true CN105656448A (zh) | 2016-06-08 |
CN105656448B CN105656448B (zh) | 2018-12-25 |
Family
ID=56478395
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201511018313.0A Expired - Fee Related CN105656448B (zh) | 2015-12-28 | 2015-12-28 | 一种音频滤波器 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105656448B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104485114A (zh) * | 2014-11-27 | 2015-04-01 | 湖南省计量检测研究院 | 一种基于听觉感知特性的语音质量客观评估的方法 |
CN104658544A (zh) * | 2013-11-20 | 2015-05-27 | 大连佑嘉软件科技有限公司 | 一种语音中瞬态噪声抑制的方法 |
CN104778948A (zh) * | 2015-04-29 | 2015-07-15 | 太原理工大学 | 一种基于弯折倒谱特征的抗噪语音识别方法 |
-
2015
- 2015-12-28 CN CN201511018313.0A patent/CN105656448B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104658544A (zh) * | 2013-11-20 | 2015-05-27 | 大连佑嘉软件科技有限公司 | 一种语音中瞬态噪声抑制的方法 |
CN104485114A (zh) * | 2014-11-27 | 2015-04-01 | 湖南省计量检测研究院 | 一种基于听觉感知特性的语音质量客观评估的方法 |
CN104778948A (zh) * | 2015-04-29 | 2015-07-15 | 太原理工大学 | 一种基于弯折倒谱特征的抗噪语音识别方法 |
Non-Patent Citations (1)
Title |
---|
胡峰松等: ""基于 Gammatone 滤波器组的听觉特征提取"", 《计算机工程》 * |
Also Published As
Publication number | Publication date |
---|---|
CN105656448B (zh) | 2018-12-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
WO2021148052A9 (zh) | Pdcch的监控方法、装置、用户设备及存储介质 | |
CN104706350A (zh) | 一种心电信号降噪方法 | |
CN105393553A (zh) | 虚拟采样率增加的音频滤波 | |
CN104251934A (zh) | 谐波分析方法和装置以及确定谐波间杂波的方法和装置 | |
CN109446649B (zh) | 致密油藏体积压裂水平井三维渗流模型的建立方法 | |
CN102881293A (zh) | 一种适用于语音压缩感知的过完备字典构造方法 | |
CN203416230U (zh) | 滤波器 | |
CN105656448A (zh) | 一种音频滤波器 | |
CN104091591B (zh) | 一种音频处理方法及装置 | |
CN104375184B (zh) | 一种高效的地震数据随机噪声衰减方法 | |
CN109088616A (zh) | 一种基于混沌振子的信号去噪方法 | |
CN101355367B (zh) | 干扰消除的方法和装置 | |
CN104091595A (zh) | 一种音频处理方法及装置 | |
CN105046654A (zh) | 一种基于粒子群优化的心电信号自适应非局部均值降噪方法 | |
GB201102075D0 (en) | Data indexing | |
CN206732265U (zh) | 一种刻环锯片基体 | |
CN105212963B (zh) | 实现超声彩色多普勒血流成像的壁滤波方法及系统 | |
Belitser et al. | On the empirical Bayes approach to adaptive filtering | |
Corwin et al. | American Modern: Documentary Photographs by Abbott, Evans, and Bourke-White | |
CN101719812B (zh) | 减少无线传感器网络节点信号冗余的方法 | |
CN105656452A (zh) | 一种音频滤波器的设计方法 | |
Yoong | Ritzenthaler Publishes ‘The Ribboncutters’ and ‘A Good Place Gone’ | |
Cserne et al. | Introduction: Law and Economics as Interdisciplinary Exchange | |
Širca et al. | Transformations of Functions and Signals | |
De Alencar et al. | 8 Operational Amplifier Model |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20181225 Termination date: 20211228 |