CN105656448B - 一种音频滤波器 - Google Patents

一种音频滤波器 Download PDF

Info

Publication number
CN105656448B
CN105656448B CN201511018313.0A CN201511018313A CN105656448B CN 105656448 B CN105656448 B CN 105656448B CN 201511018313 A CN201511018313 A CN 201511018313A CN 105656448 B CN105656448 B CN 105656448B
Authority
CN
China
Prior art keywords
filter
response
gammatone
tone
frequency
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.)
Expired - Fee Related
Application number
CN201511018313.0A
Other languages
English (en)
Other versions
CN105656448A (zh
Inventor
凌永权
张小志
李志光
方百立
梁卓銘
黄永恆
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Guangdong University of Technology
Original Assignee
Guangdong University of Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Guangdong University of Technology filed Critical Guangdong University of Technology
Priority to CN201511018313.0A priority Critical patent/CN105656448B/zh
Publication of CN105656448A publication Critical patent/CN105656448A/zh
Application granted granted Critical
Publication of CN105656448B publication Critical patent/CN105656448B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

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、音频信号处理系统等，但是其稳定性需要通过计算冲击响应的绝对值积分是否有界来判别。连续时间有理数滤波器稳定性容易确定，在设计和应用中广受欢迎，但是Alan V.Oppenheim等人在《Signals andSystems》中指出，由于连续时间有理数滤波器的频率响应由频率复数多项式构成，功率函数在频域上不够集中，滤波效果比连续时间无理数滤波器要差。2008年Stefan Strahl等人在Brain Research上发表题为“Sparse gammatone signal model optimized forEnglish speech does not match the human auditory filters”的文章，提出了连续时间整数阶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＝[xi yi]T
θi＝[mi bi fi]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(ω)＝Ri cosαi-Iisinαi+j(Risinαi+Ii cosα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-Only Memory，ROM)或随机存储记忆体(RandomAccessMemory，RAM)等。
以上所述是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也视为本发明的保护范围。

Claims (1)

1.一种音频滤波器，其特征在于，所述音频滤波器由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分别可以优化设计。
CN201511018313.0A 2015-12-28 2015-12-28 一种音频滤波器 Expired - Fee Related CN105656448B (zh)

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 CN105656448A (zh) 2016-06-08
CN105656448B true 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)

* Cited by examiner, † Cited by third party
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 太原理工大学 一种基于弯折倒谱特征的抗噪语音识别方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
Title
"基于 Gammatone 滤波器组的听觉特征提取";胡峰松等;《计算机工程》;20121130;第38卷(第21期);第2.1节 *

Also Published As

Publication number Publication date
CN105656448A (zh) 2016-06-08

Similar Documents

Publication Publication Date Title
JP2015519087A5 (zh)
CN105393554A (zh) 利用正弦函数的音频滤波器
WO2014158387A3 (en) System and method for attenuating noise in seismic data
CN104950335B (zh) Enpemf信号归一化stft‑wvd时频分析方法
CN103684621A (zh) 用于减少脉冲噪声干扰的方法和系统
CN105126187A (zh) 一种通过脑电波辅助睡眠的方法和系统
CN105656448B (zh) 一种音频滤波器
CN104133248A (zh) 一种高保真声波干扰压制方法
EP3739903A3 (en) Active room compensation in loudspeaker system
MX2013010518A (es) Procedimiento de señal acustica utilizando filtracion adaptable basada en modelos.
CN104242877B (zh) 基于模拟cr‑rc电路的核脉冲信号数字高斯成形方法
CN104836547B (zh) 一种短群延时数字滤波方法
CN109088616B (zh) 一种基于混沌振子的信号去噪方法
WO2015138485A3 (en) Hydrophone response compensation filter derivation, design and application
Devkate et al. Protein aggregation a review
CN105022917B (zh) 一种信号精确提取与处理方法
CN103441973B (zh) 一种差分信号抑制补偿方法、装置及系统
CN105656452B (zh) 一种音频滤波器的设计方法
CN104050363A (zh) 热工过程数据滤波方法及其装置
CN104067632A (zh) 指向性控制方法及装置
CN107112979A (zh) 具有在高解析音频的前响应频率下的群延迟的非线性滤波器
CN105811921B (zh) 一种抑制工频谐波干扰的方法及滤波器
CN106162435A (zh) 一种抑制谐振的麦克风单体
CN211529611U (zh) 一种用于语音识别系统的采样电路
WO2017054507A1 (zh) 音效模拟方法和装置及系统

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

Granted publication date: 20181225

Termination date: 20211228

CF01 Termination of patent right due to non-payment of annual fee