CN102447446A - 基于振动元件运动状态反馈的扬声器频响均衡方法和装置 - Google Patents

基于振动元件运动状态反馈的扬声器频响均衡方法和装置 Download PDF

Info

Publication number
CN102447446A
CN102447446A CN201110407787XA CN201110407787A CN102447446A CN 102447446 A CN102447446 A CN 102447446A CN 201110407787X A CN201110407787X A CN 201110407787XA CN 201110407787 A CN201110407787 A CN 201110407787A CN 102447446 A CN102447446 A CN 102447446A
Authority
CN
China
Prior art keywords
signal
sound
analog
frequency response
digital
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.)
Pending
Application number
CN201110407787XA
Other languages
English (en)
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.)
SHANGSHENG ELECTRONIC CO Ltd SUZHOU
Suzhou Sonavox Electronics Co Ltd
Original Assignee
SHANGSHENG ELECTRONIC CO Ltd SUZHOU
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 SHANGSHENG ELECTRONIC CO Ltd SUZHOU filed Critical SHANGSHENG ELECTRONIC CO Ltd SUZHOU
Priority to CN201110407787XA priority Critical patent/CN102447446A/zh
Publication of CN102447446A publication Critical patent/CN102447446A/zh
Pending legal-status Critical Current

Links

Images

Landscapes

  • Circuit For Audible Band Transducer (AREA)

Abstract

本发明公开了一种基于振动元件运动状态反馈的扬声器频响均衡方法和装置。均衡方法包括:1)利用传感器件拾取扬声器振动元件的运动状态信息并转换为模拟电信号;2)利用模数转换器件将反馈的模拟信号转换为数字信号;3)计算出多个级联的均衡器参数;4)利用级联均衡器对声源信号进行逐级均衡处理,消除系统频响起伏;5)将均衡后声源信号送至功率放大器,并由功放驱动扬声器发声。均衡装置包括声源、数字信号处理器、功率放大器、扬声器单元、传感器、模数转换器,各单元依次顺序连接。本发明能够实时有效的消除扬声器的频响起伏性,算法实现简单可靠,并且对系统频响变化具有很好的自适应均衡特性,提高了扬声器单元的智能化水平。

Description

基于振动元件运动状态反馈的扬声器频响均衡方法和装置
技术领域
本发明涉及一种扬声器的频响均衡方法和装置,具体地说,本发明特别涉及一种基于振动元件运动状态反馈的扬声器频响均衡方法和装置。
 
背景技术
在扬声器的生产制作过程中,因材料、结构及加工精度等各方面所存在的缺陷及障碍,会导致扬声器单元本身的频响曲线在幅度上存在着较多的峰谷同时相位上存在着一定的非线性,为了改善扬声器的频响特性,保证重放信号的真实还原,需要考虑对扬声器单元的幅度和相位响应进行均衡处理,提高幅度响应的平坦度和相位响应的线性度。
目前,中国专利ZL 200820047855.X和美国专利US 4130727、US 4739514、US 4130726中所涉及的传统的扬声器频响均衡方法, 都是基于预先测量出的扬声器频响曲线数据进行均衡器参数设计,但是在实际应用过程中扬声器自身的频响特性经常会发生改变,从而会导致先前测量的频响曲线与实际扬声器的工作频响曲线出现较大的偏差,从而造成基于先验测量频响数据所设计的均衡器在实用中出现均衡失效现象。这些基于固定参数的均衡处理方法,缺乏随系统频响变化自适应调整均衡器参数的能力,从而造成这些均衡方法在诸多实用环境中存在着较大的缺陷性。另外,这些传统的扬声器频响均衡方法,都是采用传声器接收待测扬声器单元的发射信号,进而获取扬声器单元的频响曲线数据,但是在测试过程中要扣除掉传声器本身的频响以及扬声器与传声器之间空间传播信道的响应特性对扬声器单元频响的影响因素,以便获得准确的扬声器单元频响曲线数据。这些测量实验需要使用昂贵的测量仪器和传声器单元,并且频响曲线的求取过程也较为麻烦。
为了克服传统扬声器频响曲线测量和计算过程的麻烦性以及传统方法在自适应均衡方面的缺陷性,美国专利US 5191619、US 4118600、US 5588065、US 5764781、US 6104817和US 20030194097A都提出了各种扬声器单元振动元件的运动状态拾取方法,并利用拾取的振动元件运动状态信息对系统频响进行补偿校正处理。这些方法都是将扬声器振动元件的运动状态信息作为反馈信号送至系统功放输入端,与功放输入端的源信号进行相减操作,从而达到改善扬声器系统整体频响特性的目的。这些基于运动状态反馈信息的系统频响补偿方法,并没有利用运动反馈信号进行更为深入的扬声器频响均衡处理研究,其频响补偿能力和频响补偿的自适应能力都非常有限。
针对现有的扬声器频响均衡方法所存在的缺陷和不足,需要寻找更为有效的扬声器频响均衡方法,以提高扬声器频响均衡方法的响应均衡效果同时提高均衡方法的自适应能力,并制作具有自适应频响均衡能力的扬声器系统装置。
发明内容
本发明的目的是克服现有的扬声器频响均衡方法所存在的缺陷和不足,提出一种基于振动元件运动状态反馈的扬声器频响均衡方法和装置。
为了达到上述目的,本发明一方面提供一种基于振动元件运动状态反馈的扬声器频响均衡方法,包括如下步骤:
    1)利用传感器件拾取扬声器振动元件的运动状态信息并转换为模拟电信号;
2)利用模数转换器将反馈的模拟信号转换为数字信号;
3)结合反馈的数字信号和迭代最小二乘准则,计算出多个级联的均衡器参数;
4)利用级联均衡器对声源信号进行逐级均衡处理,消除系统频响起伏;
5)将均衡后的声源信号送至功率放大器,并由功放驱动扬声器发声。
进一步地,步骤1)中所述传感器件的运动状态拾取方式主要分为四种:a)感应线圈拾取振动元件的运动状态信息;b)感应电容拾取振动元件的运动状态信息;c)激光拾取振动元件的运动状态信息;d)MEMS(Micro Electro Mechanical systems)加速度传感器拾取振动元件的运动状态信息。
感应线圈拾取振动元件的运动状态信息:感应线圈与音圈粘接在一起,感应线圈可以置于音圈上侧、中间或者下侧,其粘接位置可以根据扬声器的制作需要选择,感应线圈随音圈一起运动并在磁场中产生感应电流,从而拾取出音圈的运动状态信息。
感应电容拾取振动元件的运动状态信息:感应电容器的下极板固定在磁柱上,上极板粘接在扬声器的锥盆上,上极板随着锥盆一起运动,这样锥盆的运动状态将引起电容器上下极板间距的变化,从而引起电容器的电容值发生变化,从而可以基于电容容值的变化来获取扬声器的锥盆运动状态信息。
激光拾取振动元件的运动状态信息:激光拾振装置可以放置在扬声器锥盆的上方或者固定在扬声器磁柱上,其装置的粘接位置可以结合扬声器的制作工艺进行合理选择,该装置可以拾取出锥盆的运动状态信息。
MEMS加速度传感器拾取振动元件的运动状态信息:MEMS加速度传感器可以粘接在扬声器锥盆上表面、锥盆内部、锥盆的下表面、音圈上表面、音圈内部、音圈下表面等位置点处,MEMS加速度传感器的粘接位置可以根据扬声器的制作工艺进行合理选择。利用加速度传感器可以拾取出扬声器音圈或者锥盆的运动状态信息。
上述四种传感器件的运动状态拾取手段,不仅仅限于动圈扬声器的使用范围,可以推广应用于压电或者静电扬声器的振动元件运动状态拾取。
进一步地,步骤2)中所述的模数转换器件,其内部包含了对反馈电信号的调理电路,该调理电路可以对反馈信号进行适当的放大和滤波处理,在放大到足够幅度的基础上再进行模数转换处理。
进一步地,步骤3)中所述的级联的均衡器参数,其测算过程如下:
假设数字格式的输入声源信号的时域序列矢量为:
Figure 219614DEST_PATH_IMAGE001
 
其中,
Figure 381605DEST_PATH_IMAGE002
是声源信号时域离散序列的采样点数。假设待均衡的系统时域脉冲响应序列表达式:
Figure 708681DEST_PATH_IMAGE003
其中,
Figure 914404DEST_PATH_IMAGE004
为系统时域脉冲响应的序列长度。该系统脉冲响应由系统内部功率放大器和扬声器单元两部分脉冲响应合成后形成系统总体的脉冲响应,同时假定模数转换器为理想的无失真传输系统。假设声源信号经过系统功率放大器、扬声器单元、传感器和模数转换器之后获得的反馈信号时域序列矢量为:
在第1次迭代时,假定待求的均衡器时域脉冲响应序列的长度为,那么该响应序列矢量可以表示为:
反馈信号
Figure 443546DEST_PATH_IMAGE008
经由待求均衡器
Figure 15473DEST_PATH_IMAGE009
处理后,得到第1次均衡后的反馈信号表达式为:
Figure 581583DEST_PATH_IMAGE010
结合声源信号
Figure 70203DEST_PATH_IMAGE011
和第1次均衡后的反馈信号,根据最小二乘准则,计算出使声源信号
Figure 755579DEST_PATH_IMAGE011
和均衡后反馈信号
Figure 125380DEST_PATH_IMAGE012
之间的误差功率取最小值时的均衡器脉冲响应矢量,作为第1次迭代时设计的最优均衡器响应矢量的估计值,其表达式为: 
Figure 468506DEST_PATH_IMAGE013
其中,
Figure 203244DEST_PATH_IMAGE014
Figure 812079DEST_PATH_IMAGE015
维的数据矩阵,其表达式为:
Figure 906943DEST_PATH_IMAGE016
对第1次迭代后的均衡器响应矢量估计值
Figure 917625DEST_PATH_IMAGE017
进行归一化处理,得到:
Figure 823264DEST_PATH_IMAGE018
由归一化的均衡器响应矢量估计值
Figure 919396DEST_PATH_IMAGE017
,可以构造均衡器响应矩阵的估计值
Figure 820880DEST_PATH_IMAGE019
,其表达式为:
Figure 686068DEST_PATH_IMAGE020
其中
Figure 346036DEST_PATH_IMAGE021
维的数据矩阵。根据第1次迭代后均衡器的响应矩阵估计值,可以得出声源信号经过第1次均衡处理后的信号表达式为:
Figure 48282DEST_PATH_IMAGE022
作为本发明的一种优选,由于仅靠第1次迭代所获得均衡器参数值,对系统响应进行均衡处理,仍不能达到很好的均衡效果,为了提高均衡处理的响应修正精度,需要进行迭代均衡处理,迭代过程利用迭代最小二乘准则估计多个级联均衡器的响应,具体步骤包括:
假定迭代循环次数为
Figure 767976DEST_PATH_IMAGE023
,在第
Figure 15418DEST_PATH_IMAGE024
次迭代时,待求的预滤波器冲激响应的序列矢量定义为:
Figure 86142DEST_PATH_IMAGE025
利用第
Figure 326500DEST_PATH_IMAGE024
个均衡器,对经过
Figure 103963DEST_PATH_IMAGE026
次均衡处理后的反馈信号进行第
Figure 329594DEST_PATH_IMAGE024
次均衡处理,获得经
Figure 921112DEST_PATH_IMAGE024
次均衡处理后的反馈信号为:
Figure 553082DEST_PATH_IMAGE028
根据迭代最小二乘准则,联合声源信号和第次均衡处理后的反馈信号
Figure 835049DEST_PATH_IMAGE029
,计算出使两者之间误差功率最小时,对应的第
Figure 587104DEST_PATH_IMAGE024
个均衡器响应矢量,作为第
Figure 409567DEST_PATH_IMAGE024
次迭代过程中,所设计的第
Figure 863551DEST_PATH_IMAGE024
个均衡器响应矢量,其表达式如下:
Figure 62451DEST_PATH_IMAGE030
其中,
Figure 669013DEST_PATH_IMAGE031
Figure 662377DEST_PATH_IMAGE015
维的数据矩阵,其表达式为:
Figure 603657DEST_PATH_IMAGE032
对第
Figure 340668DEST_PATH_IMAGE024
次迭代后的均衡器响应矢量估计值进行归一化处理,得到:
Figure 231581DEST_PATH_IMAGE034
由归一化的均衡器响应矢量估计值
Figure 660157DEST_PATH_IMAGE035
,可以构造均衡器响应矩阵的估计值
Figure 935281DEST_PATH_IMAGE036
,其表达式为:
Figure 516435DEST_PATH_IMAGE037
其中
Figure 851601DEST_PATH_IMAGE036
Figure 770403DEST_PATH_IMAGE021
维的数据矩阵。根据第
Figure 849218DEST_PATH_IMAGE024
次迭代后均衡器响应矩阵估计值,可以得出声源信号经过第次均衡处理后的信号表达式为:
Figure 56525DEST_PATH_IMAGE038
                        
作为本发明的又一种优选,在整个迭代的均衡器参数估计过程中,所述的迭代的均衡器参数估计过程,对其迭代循环次数的控制方式为:
首先,在第
Figure 194114DEST_PATH_IMAGE024
次迭代过程结束之后,计算出经
Figure 76620DEST_PATH_IMAGE024
次均衡处理后归一化的反馈信号与归一化的声源信号之间的均方根误差值,其表达式如下:
Figure 309335DEST_PATH_IMAGE040
                      
然后,在第
Figure 934220DEST_PATH_IMAGE041
次迭代开始之前,通过比较第
Figure 354837DEST_PATH_IMAGE024
次迭代的均方根误差
Figure 765090DEST_PATH_IMAGE039
与用户设定的期望均方根误差
Figure 878539DEST_PATH_IMAGE042
之间的大小关系来控制迭代循环的运行;如果
Figure 990721DEST_PATH_IMAGE043
,那么继续执行次迭代;如果
Figure 214209DEST_PATH_IMAGE044
,则停止迭代。
进一步地,步骤5)中所述功率放大器,其具体实现分为模拟和数字两种情况。在模拟实现情况下,经迭代均衡后数字格式的声源信号,需要先经过数模转换器转换为模拟信号,再经过滤波和缓冲放大电路处理后,送至功率放大器的模拟输入端口进行信号放大。在数字实现情况下,经迭代均衡处理后数字格式的声源信号,直接送至功率放大器的数字输入端口进行信号放大。
进一步地,步骤5)中所述扬声器,其具体实现可以为动圈扬声器、压电扬声器和静电扬声器等各种电声转换装置。
本发明另一方面提供一种基于振动元件运动状态反馈的扬声器频响均衡装置,包括:
一声源,是系统待播放的信息;
一数字信号处理器,与声源的输出端相连接,用于计算多个级联均衡器的参数值,同时利用所计算的多个级联均衡器对声源信号进行均衡处理,提高系统幅度响应的平坦度和相位响应的线性度。
一功率放大器,与数字信号处理器的输出端相连接,用于对均衡处理后的声源信号进行功率放大,以达到驱动扬声器发声的功率要求;
一扬声器单元,与功率放大器的输出端相连接,用于对均衡处理和功率放大后的声源信号进行电声变换,在空气中还原出真实的声源信号;
一传感器,与扬声器单元的输出端相连接,用于拾取扬声器振动元件的运动状态信息;
一模数转换器,与传感器的输出端相连接,用于将传感器拾取的模拟电信号转换为数字信号,并送至数字信号处理器中。
进一步地,声源的输入信号可以为模拟格式信号和数字格式信号两种情况。在模拟格式输入情况下,声源需要通过模数转换器将模拟输入信号转换为指定字长和采样率要求的数字格式信号并送至数字信号处理器中;在数字格式输入情况下,声源需要将数字格式输入信号按照指定的字长和采样率要求转换为相应格式的数字信号并送至数字信号处理器中。
数字信号处理器内部的信号处理流程为:将声源信号和反馈信号输入到第1个LS(Least Squares)估计器,计算出第1次迭代过程的均衡器的参数,然后将声源信号和经第1次均衡处理后的反馈信号送入到第2个LS估计器,计算出第2次迭代过程的均衡器的参数,按照这一过程依次类推,可以一直计算出第K次迭代过程的均衡器K的参数。在每次迭代计算之前,都要检测声源信号和经迭代均衡处理后的反馈信号之间的误差信号大小,如果误差信号大于设定的误差值则继续进行下一次迭代处理,如果误差信号小于设定误差值,则停止迭代处理过程,将经过迭代均衡处理后的声源信号送至功率放大器的输入端。
进一步地,功率放大器可以分为模拟和数字两种实现方式,对于模拟实现方式,数字信号处理器输出的数字信号,先经过数模转换器转换为模拟信号,再经过模拟滤波器和缓冲放大器调理后,送至功率放大器进行信号的功率放大;对于数字实现方式,数字信号处理器输出的数字信号直接送至功率放大器进行信号的功率放大。
进一步地,扬声器单元的实现形式可以为动圈扬声器、压电扬声器、静电扬声器等各种电声转换装置。
进一步地,传感器可以为感应线圈、感应电容、激光拾振器、MEMS加速度传感器这四种实现形式。
进一步地,模数转换器内部包含模拟滤波器、缓冲放大器和模数转换器三部分组成。传感器拾取的模拟电信号先经过模拟滤波器进行滤波处理后,再送至缓冲放大器进行放大处理,最后再通过模数转换器转换为数字信号,将该数字信号送至数字信号处理器中。
本发明的所述一种基于振动元件运动状态反馈的扬声器频响均衡方法和装置,适用于实时对扬声器系统频响进行均衡处理,能够在系统频响变化情况下自适应均衡响应,通过多次迭代的均衡处理过程,能够取得较好的系统频响均衡效果。
与现有技术相比,本发明的优点在于:
A. 本发明通过直接反馈扬声器振动元件的运动状态信息进行均衡处理,不需要预先测量出扬声器系统的频响曲线数据,省去了较为麻烦的测量和计算过程,省去了昂贵的测量仪器和昂贵的传声器接收单元。
B. 本发明能够实时地对扬声器系统的频响进行均衡处理,能够跟随系统频响变化进行自适应的均衡处理。
C. 本发明通过多次迭代的均衡方法,能够有效的提高频响均衡的精度,而且这种最小二乘的迭代方法收敛速度较快,仅需要几次迭代过程就可以达到期望的均衡效果。
D. 本发明基于数字信号处理的方法,解决了扬声器系统的频响缺陷性,能够有效提高扬声器系统的频响性能和扬声器系统的智能化水平。
E. 与传统的均衡方法相比,本发明所提出的方法实现简单、运行可靠,不需要复杂的电路设计及制作,有效节省均衡装置的制作成本和时间。
F.   本发明所采用的振动拾取装置,避免了传统的基于传声器拾取声信号方法在线处理时所存在的声反馈问题,从而可以实时的进行在线的均衡处理。
附图说明
图1是根据本发明的基于振动元件运动状态反馈的扬声器频响均衡方法和装置的信号处理流程图;
图2是图1中的传感器所使用的感应线圈的装配结构示意图;
图3是图1中的传感器所使用的感应电容器的装配结构示意图;
图4是图1中的传感器所使用的激光拾振装置的装配结构示意图;
图5是图1中的传感器所使用的MEMS加速度传感器的装配结构示意图;
图6是根据本发明的均衡方法所采用的多次迭代均衡过程的信号处理流程图;
图7是根据本发明的均衡装置所采用的信号处理器所执行的多次迭代计算均衡器参数的信号处理流程图;
图8是根据本发明的均衡装置所采用的功率放大器的实现方式示意图;
图9是根据本发明的均衡装置所采用的模数转换器的实现方式示意图;
图10是本发明一实施例的输入声源信号的时域波形图;
图11是本发明一实施例的输入声源信号的频域幅度谱曲线图;
图12是本发明一实施例的扬声器系统时域脉冲响应的波形图;
图13是本发明一实施例的扬声器系统幅频响应的曲线图;
图14是本发明一实施例的系统未施加均衡处理和系统经过第1次迭代均衡处理之后两种情况下系统幅频响应的对比曲线图;
图15是本发明一实施例的迭代均衡过程中误差随迭代次数变化的曲线图;
图16是本发明一实施例的系统未施加均衡处理、系统经过第1次迭代均衡处理之后和系统经过第10次迭代均衡处理之后三种情况下系统幅频响应的对比曲线图。
 
具体实施方式
下面结合附图和具体实施方式对本发明作进一步详细描述:
本发明首先通过传感器拾取扬声器系统的运动状态信号并转换为模拟电信号,然后再通过模数转换器将反馈的模拟电信号转换为数字信号并送至数字信号处理器中,在数字信号处理器中,结合声源信号和反馈信号,利用迭代的最小二乘估计方法,计算出多个级联均衡器的参数值,最后将声源信号通过多个级联均衡器处理后再送至功率放大器进行放大后驱动扬声器单元发声。
如图1所示的依据本发明的基于振动元件运动状态反馈的扬声器频响均衡装置,其主体由声源1、数字信号处理器2、功率放大器3、扬声器单元4、传感器5、模数转换器6等组成。声源1与所述数字信号处理器2的输入端连接,用于将待播放信号转换为制定字长和采样率的数字信号;数字信号处理器2的输入端与声源1的输出端和模数转换器6的输出端连接,用于估计多次迭代的多个级联均衡器参数值,并利用所设计的多个级联均衡器对声源1进行均衡处理,并将均衡处理后的声源信号送至功率放大器3;功率放大器3对均衡后声源信号进行放大,并驱动扬声器单元4发声;扬声器单元4的输出端与传感器5的输入端连接,由传感器5拾取扬声器单元4 的运动状态信息并转换为模拟电信号;模数转换器6将传感器5送入的模拟电信号转换为数字信号,并送至数字信号处理器2中。
传感器件的运动状态拾取方式主要分为四种:a)感应线圈拾取振动元件的运动状态信息;b)感应电容拾取振动元件的运动状态信息;c)激光拾取振动元件的运动状态信息;d)MEMS(Micro Electro Mechanical systems)加速度传感器拾取振动元件的运动状态信息。
感应线圈拾取振动元件的运动状态信息如图2所示,感应线圈与音圈粘接在一起,感应线圈可以置于音圈上侧、中间或者下侧,其粘接位置可以根据扬声器的制作需要选择,感应线圈随音圈一起运动并在磁场中产生感应电流,从而拾取出音圈的运动状态信息。
感应电容拾取振动元件的运动状态信息如图3所示,感应电容器的下极板固定在磁柱上,上极板粘接在扬声器的锥盆上,上极板随着锥盆一起运动,这样锥盆的运动状态将引起电容器上下极板间距的变化,从而引起电容器的电容值发生变化,从而可以基于电容容值的变化来获取扬声器的锥盆运动状态信息。
激光拾取振动元件的运动状态信息如图4所示,激光拾振装置可以放置在扬声器锥盆的上方或者固定在扬声器磁柱上,其装置的粘接位置可以结合扬声器的制作工艺进行合理选择,该装置可以拾取出锥盆的运动状态信息。
MEMS加速度传感器拾取振动元件的运动状态信息如图5所示,MEMS加速度传感器可以粘接在扬声器锥盆上表面、锥盆内部、锥盆的下表面、音圈上表面、音圈内部、音圈下表面等位置点处,MEMS加速度传感器的粘接位置可以根据扬声器的制作工艺进行合理选择。利用加速度传感器可以拾取出扬声器音圈或者锥盆的运动状态信息。
级联均衡器的参数值,其测算过程如下:
假设数字格式的输入声源信号的时域序列矢量为:
Figure 498560DEST_PATH_IMAGE001
                         
其中,
Figure 89248DEST_PATH_IMAGE002
是声源信号时域离散序列的采样点数。假设待均衡的系统时域脉冲响应序列表达式:
Figure 117247DEST_PATH_IMAGE003
其中,
Figure 236513DEST_PATH_IMAGE004
为系统时域脉冲响应的序列长度。该系统脉冲响应由系统内部功率放大器和扬声器单元两部分脉冲响应合成后形成系统总体的脉冲响应,同时假定模数转换器为理想的无失真传输系统。假设声源信号经过系统功率放大器、扬声器单元、传感器和模数转换器之后获得的反馈信号时域序列矢量为:
Figure 426185DEST_PATH_IMAGE005
如图6所示,在第1次迭代时,假定待求的均衡器时域脉冲响应序列的长度为,那么该响应序列矢量可以表示为:
Figure 79070DEST_PATH_IMAGE007
                  
反馈信号
Figure 318421DEST_PATH_IMAGE008
经由待求均衡器
Figure 678995DEST_PATH_IMAGE009
处理后,得到第1次均衡后的反馈信号表达式为:
Figure 253065DEST_PATH_IMAGE010
结合声源信号
Figure 622866DEST_PATH_IMAGE011
和第1次均衡后的反馈信号
Figure 451145DEST_PATH_IMAGE012
,根据最小二乘准则,计算出使声源信号
Figure 248200DEST_PATH_IMAGE011
和均衡后反馈信号
Figure 43986DEST_PATH_IMAGE012
之间的误差功率取最小值时的均衡器脉冲响应矢量,作为第1次迭代时设计的最优均衡器响应矢量的估计值,其表达式为: 
Figure 155162DEST_PATH_IMAGE013
其中,
Figure 165843DEST_PATH_IMAGE014
Figure 323680DEST_PATH_IMAGE015
维的数据矩阵,其表达式为:
Figure 419812DEST_PATH_IMAGE016
对第1次迭代后的均衡器响应矢量估计值
Figure 69099DEST_PATH_IMAGE017
进行归一化处理,得到:
Figure 121237DEST_PATH_IMAGE018
由归一化的均衡器响应矢量估计值
Figure 260095DEST_PATH_IMAGE017
,可以构造均衡器响应矩阵的估计值
Figure 781206DEST_PATH_IMAGE019
,其表达式为:
Figure 30921DEST_PATH_IMAGE020
其中
Figure 203146DEST_PATH_IMAGE019
Figure 512904DEST_PATH_IMAGE021
维的数据矩阵。根据第1次迭代后均衡器的响应矩阵估计值,可以得出声源信号经过第1次均衡处理后的信号表达式为:
Figure 521312DEST_PATH_IMAGE022
作为本发明的一种优选,由于仅靠第1次迭代所获得均衡器参数值,对系统响应进行均衡处理,仍不能达到很好的均衡效果,为了提高均衡处理的响应修正精度,如图6所示,需要进行迭代均衡处理。
数字信号处理器内部的信号处理流程如图7所示,将声源信号和反馈信号输入到第1个LS(Least Squares)估计器,计算出第1次迭代过程的均衡器的参数,然后将声源信号和经第1次均衡处理后的反馈信号送入到第2个LS估计器,计算出第2次迭代过程的均衡器的参数,按照这一过程依次类推,可以一直计算出第K次迭代过程的均衡器K的参数。在每次迭代计算之前,都要检测声源信号和经迭代均衡处理后的反馈信号之间的误差信号大小,如果误差信号大于设定的误差值则继续进行下一次迭代处理,如果误差信号小于设定误差值,则停止迭代处理过程,将经过迭代均衡处理后的声源信号送至功率放大器的输入端。
功率放大器可以分为模拟和数字两种实现方式,模拟实现方式如图8所示,数字信号处理器输出的数字信号,先经过数模转换器转换为模拟信号,再经过模拟滤波器和缓冲放大器调理后,送至功率放大器进行信号的功率放大;对于数字实现方式,数字信号处理器输出的数字信号直接送至功率放大器进行信号的功率放大。
模数转换器的实施方式可以参照图9所示,模数转换器6内部包含模拟滤波器、缓冲放大器和模数转换器,传感器5拾取的模拟电信号先经过模拟滤波器进行滤波处理,再送至缓冲放大器进行放大处理,最后再通过模数转换器转换为数字信号,将该数字信号送至数字信号处理器2中。
下面结合附图和一实施例对本发明进行详细地说明。
在本实施例中,声源信号假定为均值为0、方差为1的白高斯噪声信号,其时域序列长度为2048点,采样间隔为19.531微秒,采样率为51.2 KHz,其时域波形和频域谱线如图10和11所示。假设扬声器系统的时域脉冲响应序列的长度为170点,其时域波形和频域谱线如图12和13所示。传声器的反馈信号为声源信号与扬声器系统脉冲响应的卷积结果。
假设均衡器的长度为400点,在第1次迭代均衡完成之后,对比系统均衡前后的频响曲线,如图14所示。从图中可以看出,经过均衡器处理之后,系统频响的幅度谱得到了很好的校正,系统整个频带的幅度谱保持平坦,具有较好的频响特性。
为了观察多次迭代均衡的效果,我们设定迭代次数为10次,观察多次迭代均衡处理后归一化的反馈信号与归一化的声源信号之间的均方根误差大小,图15给出了误差大小随迭代次数变化的曲线。观察该曲线可以看出,在第1次迭代之后,系统频响均衡已经取得了较为精确的结果,在随后的迭代处理过程中,系统频响的均衡效果又进一步得到提高,在第3次均衡之后,系统频响均衡结果趋于稳定,这说明该迭代算法具有较快的收敛速度。图16给出了未施加均衡、第1次迭代均衡和第10次迭代均衡三种情况下,系统频响曲线的对比图,从图中也可以看出第1次迭代均衡之后,系统频响曲线取得了较好的均衡效果,但是仍有小量的幅度起伏,在第10次迭代均衡之后,系统频响曲线的小量幅度起伏得到了消除,这说明了通过迭代均衡操作,系统频响的均衡效果得到了提高。
最后所应说明的是,以上实施例仅用以说明本发明的技术方案而非限制。尽管参照实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,对本发明的技术方案进行修改或者等同替换,都不脱离本发明技术方案的精神和范围,其均应涵盖在本发明的权利要求范围当中。

Claims (16)

1.一种基于振动元件运动状态反馈的扬声器频响均衡方法,包括如下步骤:
1)利用传感器件拾取扬声器振动元件的运动状态信息并转换为模拟电信号;
2)利用模数转换器将反馈的模拟信号转换为数字信号;
3)结合反馈的数字信号和迭代最小二乘准则,计算出多个级联的均衡器参数;
4)利用级联均衡器对声源信号进行逐级均衡处理,消除系统频响起伏;
5)将均衡后的声源信号送至功率放大器,并由功放驱动扬声器发声。
2.根据权利要求1所述的扬声器频响均衡方法,其特征在于:步骤1)中的传感器件拾取运动状态的方式包括以下几种:a)感应线圈拾取振动元件的运动状态信息;b)感应电容拾取振动元件的运动状态信息;c)激光拾取振动元件的运动状态信息;d)MEMS加速度传感器拾取振动元件的运动状态信息。
3.根据权利要求1所述的扬声器频响均衡方法,其特征在于:步骤2)中所述的模数转换器件,其内部包含了对反馈电信号的调理电路,该调理电路用于对反馈信号进行适当的放大和滤波处理,在放大到足够幅度的基础上再进行模数转换处理。
4.根据权利要求1所述的扬声器频响均衡方法,其特征在于:步骤3)中所述的级联的均衡器参数,其测算方法包括:数字格式的输入声源信号的时域序列矢量为                                                                          
其中,
Figure 406503DEST_PATH_IMAGE002
是声源信号时域离散序列的采样点数;
待均衡的系统时域脉冲响应序列表达式为,其中,
Figure 117287DEST_PATH_IMAGE004
为系统时域脉冲响应的序列长度;
该系统脉冲响应由系统内部功率放大器和扬声器单元两部分脉冲响应合成后形成系统总体的脉冲响应,同时假定模数转换器为理想的无失真传输系统;
假设声源信号经过系统功率放大器、扬声器单元、传感器和模数转换器之后获得的反馈信号时域序列矢量为:
Figure 161335DEST_PATH_IMAGE005
在第1次迭代时,待求的均衡器时域脉冲响应序列的长度为
Figure 855622DEST_PATH_IMAGE006
,该响应序列矢量为
Figure 444866DEST_PATH_IMAGE007
反馈信号
Figure 680063DEST_PATH_IMAGE008
经由待求均衡器
Figure 75272DEST_PATH_IMAGE009
处理后,得到第1次均衡后的反馈信号表达式为:
Figure 561748DEST_PATH_IMAGE010
结合声源信号
Figure 649790DEST_PATH_IMAGE011
和第1次均衡后的反馈信号
Figure 103774DEST_PATH_IMAGE012
,根据最小二乘准则,计算出使声源信号
Figure 302674DEST_PATH_IMAGE011
和均衡后反馈信号之间的误差功率取最小值时的均衡器脉冲响应矢量,作为第1次迭代时设计的最优均衡器响应矢量的估计值,其表达式为: 
Figure 902599DEST_PATH_IMAGE013
其中,
Figure 843880DEST_PATH_IMAGE014
Figure 580891DEST_PATH_IMAGE015
维的数据矩阵,其表达式为:
Figure 41960DEST_PATH_IMAGE016
对第1次迭代后的均衡器响应矢量估计值
Figure 471804DEST_PATH_IMAGE017
进行归一化处理,得到:
Figure 634801DEST_PATH_IMAGE018
由归一化的均衡器响应矢量估计值
Figure 113187DEST_PATH_IMAGE017
,可以构造均衡器响应矩阵的估计值
Figure 756658DEST_PATH_IMAGE019
,其表达式为:
Figure 219388DEST_PATH_IMAGE020
其中
Figure 948309DEST_PATH_IMAGE019
Figure 214074DEST_PATH_IMAGE021
维的数据矩阵;
根据第1次迭代后均衡器的响应矩阵估计值,得出声源信号经过第1次均衡处理后的信号表达式为:
Figure 712052DEST_PATH_IMAGE022
5.根据权利要求4所述的扬声器频响均衡方法,其特征在于:声源信号经过迭代均衡处理,迭代过程利用迭代最小二乘法则测算多个级联均衡器的均衡器参数值,包括:在第
Figure 155803DEST_PATH_IMAGE023
次迭代时,待求的预滤波器冲激响应的序列矢量定义为:
Figure 372020DEST_PATH_IMAGE024
,                   
利用第
Figure 441476DEST_PATH_IMAGE023
个均衡器,对经过
Figure 793960DEST_PATH_IMAGE025
次均衡处理后的反馈信号进行第
Figure 112126DEST_PATH_IMAGE023
次均衡处理,获得经
Figure 719694DEST_PATH_IMAGE023
次均衡处理后的反馈信号为:
根据迭代最小二乘准则,联合声源信号和第
Figure 903048DEST_PATH_IMAGE023
次均衡处理后的反馈信号
Figure 328955DEST_PATH_IMAGE028
,计算出使两者之间误差功率最小时对应的第
Figure 656031DEST_PATH_IMAGE023
个均衡器响应矢量,作为第
Figure 878065DEST_PATH_IMAGE023
次迭代过程中,所设计的第
Figure 290592DEST_PATH_IMAGE023
个均衡器响应矢量,其表达式如下:
Figure 239962DEST_PATH_IMAGE029
其中,维的数据矩阵,其表达式为:
Figure 714303DEST_PATH_IMAGE031
对第
Figure 467364DEST_PATH_IMAGE023
次迭代后的均衡器响应矢量估计值
Figure 503453DEST_PATH_IMAGE032
进行归一化处理,得到:
Figure 67289DEST_PATH_IMAGE033
由归一化的均衡器响应矢量估计值,可以构造均衡器响应矩阵的估计值
Figure 745581DEST_PATH_IMAGE035
,其表达式为:
Figure 901756DEST_PATH_IMAGE036
其中
Figure 636494DEST_PATH_IMAGE035
Figure 245330DEST_PATH_IMAGE021
维的数据矩阵;
根据第
Figure 608703DEST_PATH_IMAGE023
次迭代后均衡器响应矩阵估计值,可以得出声源信号经过第
Figure 619384DEST_PATH_IMAGE023
次均衡处理后的信号表达式为:
Figure 525023DEST_PATH_IMAGE037
6.根据权利要求5所述的扬声器频响均衡方法,其特征在于:迭代的均衡器参数测算过程,对其迭代循环次数的控制方式为:
首先,在第次迭代过程结束之后,计算出经
Figure 191814DEST_PATH_IMAGE023
次均衡处理后归一化的反馈信号与归一化的声源信号之间的均方根误差值,其表达式如下:
Figure 648389DEST_PATH_IMAGE039
;                     
然后,在第
Figure 169500DEST_PATH_IMAGE040
次迭代开始之前,通过比较第
Figure 419216DEST_PATH_IMAGE023
次迭代的均方根误差
Figure 591440DEST_PATH_IMAGE038
与用户设定的期望均方根误差
Figure 901199DEST_PATH_IMAGE041
之间的大小关系来控制迭代循环的运行;如果
Figure 909606DEST_PATH_IMAGE042
,那么继续执行
Figure 963012DEST_PATH_IMAGE040
次迭代;如果
Figure 992673DEST_PATH_IMAGE043
,则停止迭代。
7.根据权利要求1所述的扬声器频响均衡方法,其特征在于:步骤5)中所述的功率放大器为模拟实现,经迭代均衡后数字格式的声源信号,先经过数模转换器转换为模拟信号,再经过滤波和缓冲放大电路处理后,送至功率放大器的模拟输入端口进行信号放大。
8.根据权利要求1所述的扬声器频响均衡方法,其特征在于:功率放大器为数字实现,经迭代均衡处理后数字格式的声源信号,直接送至功率放大器的数字输入端口进行信号放大。
9.根据权利要求1所述的扬声器频响均衡方法,其特征在于,步骤5)中所述的扬声器为动圈扬声器、压电扬声器或静电扬声器。
10.一种基于振动元件运动状态反馈的扬声器频响均衡装置,其特征在于包括:
一声源(1),是系统待播放的信息;
一数字信号处理器(2),其与声源(1)的输出端相连接,用于计算多个级联均衡器的参数值,同时利用所计算的多个级联均衡器对声源信号进行均衡处理,提高系统幅度响应的平坦度和相位响应的线性度;
一功率放大器(3),其与数字信号处理器(2)的输出端相连接,用于对均衡处理后的声源信号进行功率放大,以达到驱动扬声器发声的功率要求;
一扬声器单元(4),其与功率放大器(3)的输出端相连接,用于对均衡处理和功率放大后的声源信号进行电声变换,在空气中还原出真实的声源信号;
一传感器(5),其与扬声器单元(4)的输出端相连接,用于拾取扬声器振动元件的运动状态信息;
一模数转换器(6),其与传感器(5)的输出端相连接,用于将传感器拾取的模拟电信号转换为数字信号,并送至数字信号处理器(2)中。
11.根据权利要求10所述的扬声器频响均衡装置,其特征在于:声源(1)的输入信号为模拟格式信号,声源(1)通过模数转换器将模拟输入信号转换为指定字长和采样率要求的数字格式信号并送至数字信号处理器(2)中。
12.根据权利要求10所述的扬声器频响均衡装置,其特征在于:声源(1)的输入信号为数字格式信号,声源(1)将数字格式输入信号按照指定的字长和采样率要求转换为相应格式的数字信号并送至数字信号处理器(2)中。
13.根据权利要求10所述的扬声器频响均衡装置,其特征在于:功率放大器(3)为模拟实现,数字信号处理器(2)输出的数字信号,先经过数模转换器转换为模拟信号,再经过模拟滤波器和缓冲放大器调理,然后送至功率放大器(3)进行信号的功率放大。
14.根据权利要求10所述的扬声器频响均衡装置,其特征在于:所述扬声器单元(4)的为动圈扬声器、压电扬声器或静电扬声器。
15.根据权利要求10所述的扬声器频响均衡装置,其特征在于:所述传感器(5)为感应线圈、感应电容、激光拾振器或MEMS加速度传感器。
16.根据权利要求10所述的扬声器频响均衡装置,其特征在于:模数转换器(6)内部包含模拟滤波器、缓冲放大器和模数转换器,传感器(5)拾取的模拟电信号先经过模拟滤波器进行滤波处理,再送至缓冲放大器进行放大处理,最后再通过模数转换器转换为数字信号,将该数字信号送至数字信号处理器(2)中。
CN201110407787XA 2011-12-09 2011-12-09 基于振动元件运动状态反馈的扬声器频响均衡方法和装置 Pending CN102447446A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201110407787XA CN102447446A (zh) 2011-12-09 2011-12-09 基于振动元件运动状态反馈的扬声器频响均衡方法和装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201110407787XA CN102447446A (zh) 2011-12-09 2011-12-09 基于振动元件运动状态反馈的扬声器频响均衡方法和装置

Publications (1)

Publication Number Publication Date
CN102447446A true CN102447446A (zh) 2012-05-09

Family

ID=46009597

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201110407787XA Pending CN102447446A (zh) 2011-12-09 2011-12-09 基于振动元件运动状态反馈的扬声器频响均衡方法和装置

Country Status (1)

Country Link
CN (1) CN102447446A (zh)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102883243A (zh) * 2012-10-15 2013-01-16 苏州上声电子有限公司 在线迭代的声重放系统频响均衡方法和装置
CN102903367A (zh) * 2012-10-15 2013-01-30 苏州上声电子有限公司 离线迭代的声重放系统频响均衡方法和装置
CN103634726A (zh) * 2013-08-30 2014-03-12 苏州上声电子有限公司 一种扬声器自动均衡方法
CN107453773A (zh) * 2016-05-25 2017-12-08 格罗方德半导体公司 将取样内插用于可调均衡器调适的系统、方法及软件程序
CN108366331A (zh) * 2017-01-24 2018-08-03 瑞昱半导体股份有限公司 音频处理装置与音频处理方法
CN108462921A (zh) * 2017-02-20 2018-08-28 北京东方广视科技股份有限公司 扬声器声场均衡设备
CN112261229A (zh) * 2020-09-11 2021-01-22 歌尔科技有限公司 骨传导通话设备测试方法、装置及系统
CN114584896A (zh) * 2022-02-25 2022-06-03 苏州清听声学科技有限公司 一种参量阵扬声器的频响曲线测试方法及系统
US11579165B2 (en) 2020-01-23 2023-02-14 Analog Devices, Inc. Method and apparatus for improving MEMs accelerometer frequency response

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070127738A1 (en) * 2003-12-15 2007-06-07 Sony Corporation Audio signal processing device and audio signal reproduction system
CN201243266Y (zh) * 2008-05-16 2009-05-20 华南理工大学 一种基于dsp的声频系统频响特性均衡处理装置
CN101645264A (zh) * 2009-07-17 2010-02-10 中国科学院声学研究所 一种基于迭代最小二乘-预滤波的最优主动时反聚焦方法
CN101702093A (zh) * 2009-10-09 2010-05-05 江苏大学 一种振动结构辐射噪声的有源控制系统及控制方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070127738A1 (en) * 2003-12-15 2007-06-07 Sony Corporation Audio signal processing device and audio signal reproduction system
CN201243266Y (zh) * 2008-05-16 2009-05-20 华南理工大学 一种基于dsp的声频系统频响特性均衡处理装置
CN101645264A (zh) * 2009-07-17 2010-02-10 中国科学院声学研究所 一种基于迭代最小二乘-预滤波的最优主动时反聚焦方法
CN101702093A (zh) * 2009-10-09 2010-05-05 江苏大学 一种振动结构辐射噪声的有源控制系统及控制方法

Cited By (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102883243A (zh) * 2012-10-15 2013-01-16 苏州上声电子有限公司 在线迭代的声重放系统频响均衡方法和装置
CN102903367A (zh) * 2012-10-15 2013-01-30 苏州上声电子有限公司 离线迭代的声重放系统频响均衡方法和装置
WO2014059890A1 (zh) * 2012-10-15 2014-04-24 苏州上声电子有限公司 在线迭代的声重放系统频响均衡方法和装置
WO2014059897A1 (zh) * 2012-10-15 2014-04-24 苏州上声电子有限公司 离线迭代的声重放系统频响均衡方法和装置
CN103634726A (zh) * 2013-08-30 2014-03-12 苏州上声电子有限公司 一种扬声器自动均衡方法
CN103634726B (zh) * 2013-08-30 2017-03-08 苏州上声电子有限公司 一种扬声器自动均衡方法
CN107453773A (zh) * 2016-05-25 2017-12-08 格罗方德半导体公司 将取样内插用于可调均衡器调适的系统、方法及软件程序
CN107453773B (zh) * 2016-05-25 2019-10-08 格罗方德半导体公司 将取样内插用于可调均衡器调适的系统、方法及软件程序
CN108366331A (zh) * 2017-01-24 2018-08-03 瑞昱半导体股份有限公司 音频处理装置与音频处理方法
CN108462921A (zh) * 2017-02-20 2018-08-28 北京东方广视科技股份有限公司 扬声器声场均衡设备
CN108462921B (zh) * 2017-02-20 2020-06-02 北京东方广视科技股份有限公司 扬声器声场均衡设备
US11579165B2 (en) 2020-01-23 2023-02-14 Analog Devices, Inc. Method and apparatus for improving MEMs accelerometer frequency response
CN112261229A (zh) * 2020-09-11 2021-01-22 歌尔科技有限公司 骨传导通话设备测试方法、装置及系统
CN114584896A (zh) * 2022-02-25 2022-06-03 苏州清听声学科技有限公司 一种参量阵扬声器的频响曲线测试方法及系统
CN114584896B (zh) * 2022-02-25 2023-10-31 苏州清听声学科技有限公司 一种参量阵扬声器的频响曲线测试方法及系统

Similar Documents

Publication Publication Date Title
CN102447446A (zh) 基于振动元件运动状态反馈的扬声器频响均衡方法和装置
US11736093B2 (en) Identifying mechanical impedance of an electromagnetic load using least-mean-squares filter
CN103634726B (zh) 一种扬声器自动均衡方法
US8831246B2 (en) MEMS microphone with programmable sensitivity
CN108174320B (zh) 一种自适应耳道主动降噪耳机及自适应耳道主动降噪方法
US20150189445A1 (en) Wide Dynamic Range Microphone
CN107742522A (zh) 基于麦克风阵列的目标语音获取方法及装置
US20130294616A1 (en) Method and system for speech enhancement in a room
CN108401204B (zh) 一种新型主动降噪耳机
JP2005094777A5 (zh)
WO2008091403A1 (en) System and method for calibrating phase and gain mismatches of an array microphone
CN102883243B (zh) 在线迭代的声重放系统频响均衡方法
CN110291718A (zh) 校准麦克风截止频率的系统和方法
CN102685637A (zh) 信号处理装置和信号处理方法
CN112289295B (zh) 主动降噪系统训练方法以及相关设备
US10560787B2 (en) Hearing device comprising an amplifier system for minimizing variation in an acoustical signal caused by variation in gain of an amplifier
EP2077549B1 (en) Coefficient measurement apparatus, effect impartment apparatus, and musical sound generating apparatus
CN109379652B (zh) 一种耳机有源噪声控制的次级通道离线辨识方法
CN105356861B (zh) 一种有源降噪的方法和系统
CN111354368B (zh) 补偿处理后的音频信号的方法
Skordilis et al. Multichannel speech enhancement using MEMS microphones
RU2716846C2 (ru) Способ коррекции аудиосигнала
US20140269198A1 (en) Beamforming Sensor Nodes And Associated Systems
KR102412405B1 (ko) 시스템 및 방법
TW201626814A (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
C12 Rejection of a patent application after its publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20120509