CN103632009B - 一种模拟反馈有源降噪耳机的设计方法 - Google Patents
一种模拟反馈有源降噪耳机的设计方法 Download PDFInfo
- Publication number
- CN103632009B CN103632009B CN201310674036.3A CN201310674036A CN103632009B CN 103632009 B CN103632009 B CN 103632009B CN 201310674036 A CN201310674036 A CN 201310674036A CN 103632009 B CN103632009 B CN 103632009B
- Authority
- CN
- China
- Prior art keywords
- frequency characteristic
- order
- function
- controller
- secondary path
- 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.)
- Active
Links
Landscapes
- Soundproofing, Sound Blocking, And Sound Damping (AREA)
Abstract
本发明公布了一种模拟反馈有源降噪耳机的设计方法,其特征为:基于H∞技术综合优化的控制器频率特性;先数字滤波器实现优化的控制器频率特性,再拟合数字滤波器的频率特性为低阶的连续域传递函数;用广义二阶GIC电路实现所得的低阶连续域传递函数。其显著优点在于:(1)基于H∞技术综合优化的控制器频率特性,保证了降噪性能的优化和系统的鲁棒性;(2)先数字滤波器实现优化的控制器频率特性,再拟合数字滤波器的频率特性为低阶的连续域传递函数,保证了控制器的最小相位特性;(3)应用广义二阶GIC电路可实现任意零极点位置的二阶稳定系统传递函数,且具有低灵敏度特性,物理可实现性好。
Description
一、技术领域
本发明涉及反馈有源噪声控制系统的控制设计方法,主要是对有源降噪耳机提出了一种模拟反馈设计方法。
二、背景技术
有源耳机作为一种听力保护装置,其结构如附图1所示,由传声器、电路(包括前置放大电路、控制电路和功率放大电路)、次级源(扬声器)组成。其降噪设计多基于反馈技术实现,通常是给定设计目标,应用一定的优化设计方法得到优化的控制器频率特性,采用数字反馈或模拟反馈实现。
数字反馈具有灵活性高,易实现的优点。Rafaely等应用H2性能标准和H∞约束设计了数字反馈控制器(Rafaely B.and Elliott S.J..H2/H∞ Active Control of Sound in a Headrest:Design and Implementation[J].IEEE Transactions on Control Systems Technology,1999(7):79-84.),Zhang等应用水床展平方法设计了数字反馈控制器(Zhang L.,Wu L.and Qiu X.AnIntuitive Approach for Feedback Active Noise Controller Design,Applied Acoustics,2013(74):160-168.)。然而由于数字反馈的时延会导致降噪性能下降和稳定性问题,且成本较模拟反馈的高,多数有源耳机采用模拟反馈实现。模拟反馈设计的关键是如何用模拟电路实现优化的控制器频率特性。
专利200510033638.6采用二阶低通滤波器、放大器和全通移相器实现耳机反馈闭环消噪。专利200710001442.8采用二阶低通滤波和放大电路实现有源耳机降噪。专利200910238545.5采用一阶电阻电容滤波和反馈回路增益单元实现耳机的主动消噪。专利200910253182.2采用两级陷波和放大电路实现耳机的反馈消噪。上述专利中均未给出如何得到优化的控制器频率特性和所采用电路中的参数是如何设计的。
Yu等采用H2范数最小化平均声功率设计控制器并采用4阶模拟滤波电路实现(Yu S.H.and Hu J.S..Controller Design for Active Noise Cancellation Headphones Using ExperimentalRaw Data[J].IEEE/Asme Transcations on Mechatronics,2001(6):483-490.),但并未给出具体的4阶模拟反馈电路形式。Pawelczyk采用多参数优化与大量搜索方法设计了4阶模拟反馈控制器,并采用一双二阶开关电容状态变量滤波器MF10实现(Pawelczyk M..Analogue Active NoiseControl of Acoustic Noise at a Virtual Location[J].IEEE Transactions on Control SystemsTechnology,2009(17):456-472.)。Bai等基于H∞技术与线性变换设计了连续域控制器,将其转换为一阶、二阶级联形式并采用基于单运放的一阶单元,二阶单元串接实现了前面五阶电路(Bai M.and Lee D.,Implementation of an Active Headset by Using the H∞ Robust ControlTheory,Journal ofAcoustical Society of America,1997(102):2184-2190.)。以上文献中未详细给出如何把优化的控制器频率特性转换为低阶的连续域传递函数形式,而这一步是模拟电路实现优化控制器的关键。此外,给定的模拟电路单元不能实现任意零极点位置的二阶稳定系统传递函数。若二阶稳定系统传递函数的零极点位置不当,会导致指定的单元电路实际电路参数为负值,物理上不可实现。在实际应用中根据优化的控制器得到的模拟反馈控制器的传递函数零极点位置是任意的,这就需要一个可以实现任意零极点位置的模拟电路单元。
三、发明内容
1、发明目的:提出一种模拟反馈有源降噪耳机的设计方法,这种方法包含了优化控制器设计、低阶的连续域传递函数转换和模拟电路实现三部分内容,且给出的广义二阶GIC(通用阻抗变换,General Impedance Converter)电路单元可实现任意零极点位置的二阶稳定系统传递函数,并具有低灵敏度特性,因此该方法具有设计性能优化,鲁棒性强和物理可实现性好的优点。
2、技术方案:为实现上述发明目的,本发明提出基于H∞技术、频率特性拟合和广义二阶GIC电路实现有源降噪耳机模拟反馈的设计,包括以下步骤:
(1)测试耳机正常佩戴时次级路径的频率特性,称为标称次级路径。
次级路径为附图1所示的B点相对于A点的路径,包括功率放大电路、扬声器、扬声器至传声器的声学路径、传声器和前置放大电路。测试次级路径的方法是将测试设备的输出信号接功率放大电路的输入端,前置放大电路输出接测试设备的输入端。对测试设备的输出信号与输入信号做FFT分析得到次级路径的频率特性。
(2)应用控制论中的状态空间模型拟合标称次级路径,得到标称次级路径模型G。
(3)设定性能指标加权函数W1,|W1(jω)|>1对应噪声抑制,|W1(jω)|<1对应水床放大。
(4)设定输入加权函数W2,其功能是去除控制输入中的高频成分。
(5)设定模型摄动加权函数W3,加权函数W3增益越大,稳定裕量相应地越大。
(6)由标称次级路径模型和加权函数得到增广对象模型P,表示为
(7)设计鲁棒H∞控制器,得到优化的控制器频率特性F。
(8)根据优化的控制器频率特性F设计数字滤波器D。
(9)应用控制论中的状态空间模型拟合数字滤波器D的频率特性为多个二阶串接的连续域传递函数。
(10)用多级广义二阶GIC电路实现多个二阶串接的连续域传递函数,其中广义二阶GIC电路如附图2所示,其对应的传递函数为
式中s为拉普拉斯变换因子,ω0为极点频率,Q为品质因数,a,b,c为系数,取值在0至1之间,设计时先设定a为0至1之间的某一值,根据所得的二阶连续域传递函数与(2)式对应项相等求得ω0,Q,b,c,若求得的b,c取值不在0至1之间,可重新设定a值并重新计算各参数直至b,c取值在0至1之间为止。
(11)将计算得到的噪声衰减量和设计目标相比较,若计算值没有达到设计目标,则返回步骤(2)至(10),通过调整加权函数改进设计,如此不断反复,直到计算值达到设计目标。
3、有益效果:本发明显著优点在于:(1)基于H∞技术综合优化的控制器频率特性,保证了降噪性能的优化和系统的鲁棒性。(2)先数字滤波器实现优化的控制器频率特性,再拟合数字滤波器的频率特性为低阶的连续域传递函数,保证了控制器的最小相位特性。(3)应用广义二阶GIC电路可实现任意零极点位置的二阶稳定系统传递函数,且具有低灵敏度特性,物理可实现性好。
四、附图说明
图1是有源降噪耳罩系统结构图。
图2是广义二阶GIC电路图。
图3是多种次级路径频率特性(a)幅频特性(b)相频特性。
图4是控制电路实现图。
图5是标称次级路径下降噪性能仿真与实验测试。
图6是多种次级路径下降噪性能仿真。
图7是次级路径2对应的降噪性能仿真与实验测试。
五、具体实施方式
本发明提出的模拟反馈有源降噪耳机设计方法的特征为:基于H∞技术综合优化的控制器频率特性;先数字滤波器实现优化的控制器频率特性,再拟合数字滤波器的频率特性为低阶的连续域传递函数;用广义二阶GIC电路实现所得低阶的连续域传递函数。下面以实现一模拟反馈有源耳机低于1000Hz以下频段最大降噪量18dB作为实例对该设计方法做详细说明。
(1)将耳机佩戴在人工头上,测试所得次级路径频率特性如附图3所示,其中标称路径为耳机居中佩戴时的次级路径;次级路径1是耳机向上偏离中间位置1cm佩戴时的次级路径,相对标称路径高频端频率特性变化较大;次级路径2是在耳罩与人工头之间横向夹入一只笔时测试所得次级路径,相对标称路径整个频带内频率特性均变化较大。
(2)用Matlab中fitfrd和zpk函数对标称次级路径建模,得到8阶标称次级路径模型G。
(3)设定性能指标加权函数W1为4个品质因数为0.6,中心频率为300Hz,中心增益为2.58的二阶带通滤波器的串接。
(4)设定输入加权函数W2为1e-5,即相当于W2=[]。
(5)设定模型摄动加权函数W3为6.5。
(6)由标称次级路径模型和加权函数得到增广对象模型P。
(7)用Matlab中hinfsyn函数设计鲁棒H∞控制器,用zpk和frd函数可以得到需求控制器的频率特性F。
(8)根据F用Matlab中的fdesign.arbmagnphase函数设计对应的250阶数字滤波器D。
(9)用Matlab中fitfrd函数拟合数字滤波器D的频率特性为低阶状态空间模型,并用zpk函数得到4阶连续域传递函数如下:
(10)用广义二阶GIC电路实现所得的4阶连续域传递函数,控制电路实现如附图4所示,包括一级增益调节电路,两级广义二阶GIC电路,其中第一级广义二阶GIC电路实现的传递函数为
第二级广义二阶GIC电路实现的传递函数为
说明广义二阶GIC电路可实现的零极点位置可以是在轴上的,也可以是在复平面上,物理可实现性好。增益调节电路的增益为-15.37,由0.9881除以0.1154和0.5571得到。
对于标称次级路径,降噪性能如附图5所示,其中黑色实线为仿真结果,红色虚线为实测结果,仿真结果与实测结果曲线图形相似,吻合较好,最大降噪量约19.0dB,降噪带宽为62Hz-1420Hz,最大噪声放大量在4210Hz处约9.5dB,达到了预期目标。
采用标称次级路径设计所得的模拟控制器,仿真附图3中多种次级路径下的降噪性能,如附图6所示,最大噪声放大量不超过10dB。附图7是次级路径2对应的降噪性能仿真与实验测试结果,两者吻合较好,说明所提出的方法具有很好的鲁棒性。
Claims (1)
1.一种模拟反馈有源降噪耳机的设计方法,其特征在于它包括以下步骤:
(1)测试耳机正常佩戴时次级路径的频率特性,称为标称次级路径;
(2)应用控制论中的状态空间模型拟合标称次级路径,得到标称次级路径模型G;
(3)设定性能指标加权函数W1,|W1(jω)|>1对应噪声抑制,|W1(jω)|<1对应水床放大;
(4)设定输入加权函数W2,设定为le-5;
(5)设定模型摄动加权函数W3,设定为一常数;
(6)由标称次级路径模型和加权函数得到增广对象模型P,表示为
式中I为单位矩阵;
(7)设计鲁棒H∞控制器,得到优化的控制器频率特性F;
(8)根据优化的控制器频率特性F设计数字滤波器D;
(9)应用控制论中的状态空间模型拟合数字滤波器D的频率特性为多个二阶串接的连续域传递函数;
(10)用多级广义二阶GIC电路实现多个二阶串接的连续域传递函数,其中广义二阶GIC电路对应的传递函数为
式中s为拉普拉斯变换因子,ω0为极点频率,Q为品质因数,a,b,c为系数,系数a,b,c取值在0至1之间,设计时先设定a为0至1之间的某一值,根据所得的二阶连续域传递函数与(2)式对应项相等求得ω0,Q,b,c,若求得的b,c取值不在0至1之间,可重新设定a值并重新计算各参数直至b,c取值在0至1之间为止;
(11)将通过步骤(2)至(10)计算得到的噪声衰减量和设计目标相比较,若计算值没有达到设计目标,则返回步骤(2)至(10),通过调整性能指标加权函数W1、输入加权函数W2和模型摄动加权函数W3改进设计,如此不断反复,直到计算值达到设计目标。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310674036.3A CN103632009B (zh) | 2013-12-13 | 2013-12-13 | 一种模拟反馈有源降噪耳机的设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310674036.3A CN103632009B (zh) | 2013-12-13 | 2013-12-13 | 一种模拟反馈有源降噪耳机的设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103632009A CN103632009A (zh) | 2014-03-12 |
CN103632009B true CN103632009B (zh) | 2016-08-17 |
Family
ID=50213049
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310674036.3A Active CN103632009B (zh) | 2013-12-13 | 2013-12-13 | 一种模拟反馈有源降噪耳机的设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103632009B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106782487B (zh) * | 2016-12-20 | 2020-09-22 | 歌尔科技有限公司 | 反馈式主动降噪耳机的降噪量仿真方法和系统 |
CN107393545B (zh) * | 2017-07-17 | 2020-12-08 | 会听声学科技(北京)有限公司 | 一种具有柔性增益的反馈式主动降噪系统及方法 |
CN107808657B (zh) * | 2017-11-30 | 2020-11-20 | 会听声学科技(北京)有限公司 | 主动降噪耳机的设计系统及方法 |
CN108646598B (zh) * | 2018-05-04 | 2022-02-08 | 南京大学 | 一种适于小型化的低阶模拟反馈控制器及其设计方法 |
CN110430520B (zh) * | 2019-08-12 | 2021-07-13 | 会听声学科技(北京)有限公司 | 一种反馈滤波器的设计方法、设计装置及耳机 |
CN111081214B (zh) * | 2019-12-12 | 2022-08-16 | 西安讯飞超脑信息科技有限公司 | 主动降噪方法及主动降噪装置的反馈滤波器的优化方法 |
CN112312261B (zh) * | 2020-11-02 | 2022-04-22 | 南京大学 | 一种高时延耳罩的制作方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102300140A (zh) * | 2011-08-10 | 2011-12-28 | 歌尔声学股份有限公司 | 一种通信耳机的语音增强方法、装置及降噪通信耳机 |
-
2013
- 2013-12-13 CN CN201310674036.3A patent/CN103632009B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102300140A (zh) * | 2011-08-10 | 2011-12-28 | 歌尔声学股份有限公司 | 一种通信耳机的语音增强方法、装置及降噪通信耳机 |
Non-Patent Citations (3)
Title |
---|
反馈有源控制系统水床效应的;张丽敏等;《应用声学》;20101130;第29卷(第6期);第401-408页 * |
反馈有源降噪耳机的高频约束条件研究;严馨叶等;《电声技术》;20130530;第37卷(第5期);第19-23页 * |
有源头靠系统中次级声源和控制电路的参数优化研究;祖峰磊等;《应用声学》;20090530;第28卷(第3期);第177-183页 * |
Also Published As
Publication number | Publication date |
---|---|
CN103632009A (zh) | 2014-03-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103632009B (zh) | 一种模拟反馈有源降噪耳机的设计方法 | |
CN102016984B (zh) | 用于动态声音传送的系统和方法 | |
CN110417414B (zh) | 可切换次级回放路径 | |
CN104935293B (zh) | 大功率变压器自适应有源降噪控制方法及控制系统 | |
US9398374B2 (en) | Systems and methods for nonlinear echo cancellation | |
Milani et al. | On maximum achievable noise reduction in ANC systems | |
US20130136282A1 (en) | System and Method for Spectral Personalization of Sound | |
WO2021189309A1 (zh) | 主动降噪的方法、系统、电子设备和芯片 | |
CN102638755A (zh) | 基于频率压缩搬移的数字助听器响度补偿方法 | |
KR102004460B1 (ko) | 블루투스 회로와 디지털 신호 처리를 이용한 디지털 히어링 디바이스 | |
CN103428609A (zh) | 用于去除噪声的设备和方法 | |
CN109195050B (zh) | 一种前馈滤波器设计方法及降噪耳机 | |
US20150092951A1 (en) | Headphone with active noise cancelling and auto-calibration method thereof | |
CN103945076A (zh) | 用于移动设备的实时音频回声和背景噪声的减小 | |
CN112562624B (zh) | 主动降噪滤波器设计方法、降噪方法、系统及电子设备 | |
CN103067629A (zh) | 回声消除装置 | |
Huang et al. | Practically efficient nonlinear acoustic echo cancellers using cascaded block RLS and FLMS adaptive filters | |
CN105794228A (zh) | 自适应残余反馈抑制 | |
CN104270539A (zh) | 基于系数差的成比例仿射投影回声消除方法 | |
US20160277860A1 (en) | Apparatus, Method and Computer Program for Providing an Audio Signal | |
KR20190035460A (ko) | 비대칭 다중 채널 오디오 동적 범위 처리 방법 | |
RU2614345C1 (ru) | Способ расширения динамического диапазона в радиотехнических системах | |
CN110728970B (zh) | 一种数字辅助隔音处理的方法及装置 | |
CN102576560B (zh) | 电子音频设备 | |
CN112562625B (zh) | 主动降噪滤波器设计方法、降噪方法、系统及电子设备 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |