CN102903367A - 离线迭代的声重放系统频响均衡方法和装置 - Google Patents

离线迭代的声重放系统频响均衡方法和装置 Download PDF

Info

Publication number
CN102903367A
CN102903367A CN2012103882599A CN201210388259A CN102903367A CN 102903367 A CN102903367 A CN 102903367A CN 2012103882599 A CN2012103882599 A CN 2012103882599A CN 201210388259 A CN201210388259 A CN 201210388259A CN 102903367 A CN102903367 A CN 102903367A
Authority
CN
China
Prior art keywords
sound
signal
frequency response
feedback signal
iteration
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
CN2012103882599A
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 CN2012103882599A priority Critical patent/CN102903367A/zh
Publication of CN102903367A publication Critical patent/CN102903367A/zh
Priority to PCT/CN2013/085000 priority patent/WO2014059897A1/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • H—ELECTRICITY
    • H03—ELECTRONIC CIRCUITRY
    • H03G—CONTROL OF AMPLIFICATION
    • H03G5/00—Tone control or bandwidth control in amplifiers
    • H03G5/16—Automatic control
    • H03G5/165—Equalizers; Volume or gain control in limited frequency bands
    • H—ELECTRICITY
    • H03—ELECTRONIC CIRCUITRY
    • H03G—CONTROL OF AMPLIFICATION
    • H03G5/00—Tone control or bandwidth control in amplifiers
    • H03G5/005—Tone control or bandwidth control in amplifiers of digital signals
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04R—LOUDSPEAKERS, MICROPHONES, GRAMOPHONE PICK-UPS OR LIKE ACOUSTIC ELECTROMECHANICAL TRANSDUCERS; DEAF-AID SETS; PUBLIC ADDRESS SYSTEMS
    • H04R3/00—Circuits for transducers, loudspeakers or microphones
    • H04R3/02—Circuits for transducers, loudspeakers or microphones for preventing acoustic reaction, i.e. acoustic oscillatory feedback
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04R—LOUDSPEAKERS, MICROPHONES, GRAMOPHONE PICK-UPS OR LIKE ACOUSTIC ELECTROMECHANICAL TRANSDUCERS; DEAF-AID SETS; PUBLIC ADDRESS SYSTEMS
    • H04R3/00—Circuits for transducers, loudspeakers or microphones
    • H04R3/04—Circuits for transducers, loudspeakers or microphones for correcting frequency response

Landscapes

  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Otolaryngology (AREA)
  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Acoustics & Sound (AREA)
  • Signal Processing (AREA)
  • Circuit For Audible Band Transducer (AREA)

Abstract

本发明公开了一种离线迭代的声重放系统频响均衡方法和装置。该方法包括：（1）利用测量仪器获取待均衡系统的脉冲响应序列；（2）选取噪声信号，利用噪声信号与脉冲响应序列的卷积获取反馈信号；（3）结合噪声信号和反馈信号，利用逐次迭代的最小二乘准则依次估计出各级均衡器的参数；（4）级联各级均衡器，形成合成均衡器；（5）利用合成均衡器的参数更新FIR滤波器系数，完成系统频响均衡。该装置包括：声源、数字信号处理器、功率放大器、扬声器；各单元依次顺序连接。本发明通过增加级联均衡器数量,可以明显提升均衡器的响应均衡能力,同时通过离线估计方式,便于处理复杂计算的均衡任务，同时也节省了均衡器的硬件制作成本。

Description

离线迭代的声重放系统频响均衡方法和装置
技术领域
本发明涉及一种声重放系统的频响均衡方法和装置，特别涉及一种离线迭代的声重放系统频响均衡方法和装置。
背景技术
近年来，随着大规模集成电路和数字信号处理技术的迅速发展，基于数字信号处理技术的声重放系统响应均衡问题也逐渐受到国内外多家研究机构及企业的关注，并有几家公司推出了几款带有响应均衡功能的声学产品。乌普萨拉大学（Uppsala University）旗下的Dirac公司推出了Dirac HD Sound技术，用于解决自由场环境下扬声器单元的频响曲线峰谷点均衡，同时还推出了Dirac Live技术，用于解决房间内声重放系统的频响曲线峰谷点均衡。 香港正然传讯有限公司推出了CONEQ技术，用于解决房间内扬声器系统的频响起伏均衡，该技术利用单个传声器在扬声器单元前方按照蛇形走线路径逐点采集扬声器系统的脉冲响应数据，然后通过特定的响应数据处理算法，实行对扬声器系统的频响均衡。丹麦林道夫（LYNGDORF）公司推出了房间均衡技术——Room Perfect，该技术利用单个传声器在房间内的多个位置点采集扬声器系统的响应数据，并利用这些多点响应信息完成听者位置的响应曲线均衡。美国KRK公司也推出了房间内声重放系统响应均衡产品——Ergo（Enhanced Room Geometry Optimization），该产品也是利用单个传声器在听者位置四周采集响应数据，并处理这些响应数据获得均衡器的参数，完成对听者位置点的响应均衡。
针对自由场及混响场内扬声器系统均衡方法的研究文献较多，其中一些代表性的研究成果如下：
文献1——Stephen T. Neely and Jont B. Allen, “Invertibility of a room impulse response ,” J. Acoust. Soc. Am, Vol. 66, No. 1, pp. 165- 169, July 1979. ——提出了房间响应中最小相位响应部分的计算方法，并对最小相位响应的零极点倒置，获得了最小相位响应的逆滤波响应，基于最小相位部分的逆滤波器参数对房间响应进行均衡处理，并通过实验证实了算法的有效性。 
文献2——Yoichi Haneda , Shoji Makino, and Yutaka Kaneda, “Common acoustical pole and zero modeling of room transfer functions,” IEEE Transaction on Speech and Audio Processing, Vol. 2, No. 2, pp. 320-328, April 1994. ——针对房间内多个位置点采集的多组传递函数，提出了基于共声学极点模型的房间传递函数融合算法，减少了房间模型的估计参数，提高了计算速度。
文献3——Aki Harma, Matti Karjalainen, Lauri Savioja, Vesa Valimaki, Unto K. Laine, and Jyri Huopaniemi, “Frequency-warped signal processing for audio application,” J. Audio Eng. Soc., Vol. 48, No. 11, pp. 1011-1029, November 2000.——提出基于弯折频率域的响应均衡方法，利用线性预测编码（Linear Predictive Coding——LPC）方法计算弯折频率域内的均衡器参数，从而提高对低频区域的均衡能力。
传统的声重放系统均衡方法，都是基于对系统脉冲响应函数进行分析，拟合这些脉冲响应函数的零极点模型，再通过零极点倒置找出系统的逆滤波器响应，从而获得了声重放系统的均衡器参数。这些方法的参数估计过程都是依赖于最小均方误差（Least Mean Squares——LMS）算法或者线性预测编码（LPC）算法的单次估计来计算均衡器的参数，这种基于单次估计方法所获得的逆滤波器参数与理想的逆滤波器参数之间仍然存在着一定程度的偏差，这些偏差将造成均衡后的声重放系统频响曲线在一些频带内仍有较为明显的峰谷起伏性，仍未达到较为理想的频响平直特性。这些均衡后频响曲线的峰谷起伏特性是由均衡器的参数误差所造成的，为了削弱均衡后频响曲线的峰谷起伏特征，需要进一步提高均衡器的参数估计精度，因此需要寻找更为精确有效的均衡器参数估计方法。
针对现有基于单次LMS或者单次LPC参数估计算法，在系统均衡器参数估计方面所存在的一定误差缺陷，需要考虑采用多次迭代估计的方法，通过多次迭代操作使所估计的级联均衡器响应逐步逼近理想的系统逆滤波器响应，减少均衡器的参数估计误差，从而保证均衡后频响曲线具有更好的平直特性。
发明内容
本发明的目的在于克服基于单次LMS或者单次LPC参数估计算法，在系统均衡器参数估计方面所存在的一定误差缺陷，提供一种离线迭代的声重放系统频响均衡方法和装置。
为了达到上述目的，本发明采用的技术方案如下：一种离线迭代的声重放系统频响均衡方法，如图1所示，包括如下步骤：
（1）利用测量仪器获取待均衡系统的脉冲响应序列；
（2）选取噪声信号，利用噪声信号与脉冲响应序列的卷积获取反馈信号； 
（3）结合噪声信号和反馈信号，利用逐次迭代的最小二乘准则依次估计出各级（均衡器的参数；
（4）级联各级均衡器，生成合成均衡器；
（5）利用合成均衡器的参数更新有限脉冲响应（Finite Impulse Response——FIR）滤波器系数，完成系统频响均衡。
进一步地，步骤（1）中的利用测量仪器获取待均衡系统的脉冲响应序列，具体步骤是通过测量仪器发出噪声信号，用传声器采集和记录由噪声激励所产生的系统脉冲响应序列。按照声重放系统的均衡环节和目标的差异，如图2所示，可以将脉冲响应的测量分为以下三种情况：
情况1：如果声重放系统仅需要均衡信号处理和功放所组成电路部分的系统响应，则在脉冲响应测试时，仅测试和记录由噪声激励所产生的功放输出端的脉冲响应序列。
情况2：如果声重放系统需要均衡电路部分和扬声器部分所组成系统的响应，则在脉冲响应测量时，将传声器置于空间期望位置点，测量和记录由噪声激励所产生的该点处的脉冲响应序列。
情况3：如果声重放系统需要均衡电路部分、扬声器部分和外部环境三者所组成系统的响应，则在脉冲响应测量时，将传声器依次置于空间多个位置点，记录多个位置点上的脉冲响应序列，并按照一定的融合方式对多位置点的脉冲响应进行平均化处理，获得平均后的脉冲响应序列。
进一步地，步骤（2）中的选取噪声信号，这种噪声信号可以为指定带宽的白噪声序列或者最大长度序列（Maximum Length Sequence——MLS）所产生的噪声信号，这种信号呈现平坦的功率谱特性，以用于训练出自由场或混响场内扬声器到传声器位置点的传递函数；根据实际测量所获得脉冲响应序列的长度                                               
Figure 2012103882599100002DEST_PATH_IMAGE002
来选取噪声信号的序列长度
Figure 2012103882599100002DEST_PATH_IMAGE004
，。
进一步地，步骤（2）中的利用噪声信号与脉冲响应序列的卷积获取反馈信号，其具体实现如下：
a．假设噪声信号的时域序列矢量为：
Figure 2012103882599100002DEST_PATH_IMAGE008
，                         
其中，
Figure 2012103882599100002DEST_PATH_IMAGE010
是噪声信号时域离散序列的采样点数；假设经测量仪器测量获得的声重放系统的时域脉冲响应序列表达式：
Figure 2012103882599100002DEST_PATH_IMAGE012
，
其中，为系统时域脉冲响应的序列长度；联合噪声序列和声重放系统的时域脉冲响应序列，获得反馈信号的时域序列矢量表达式为：
Figure 2012103882599100002DEST_PATH_IMAGE016
，
其中“
Figure 2012103882599100002DEST_PATH_IMAGE018
”代表两个序列矢量之间进行卷积操作。
更进一步地，在进行噪声序列和脉冲响应序列的卷积操作时，已经将卷积后所获得反馈信号序列的头部切除掉长度为的序列部分，从而获得了长度为的反馈信号序列。
进一步地，步骤（3）中的结合噪声信号和反馈信号，利用逐次迭代的最小二乘准则依次估计出各级均衡器的参数，其具体实现流程如下：
如图3所示，在第1次迭代时，假定待求的第1级均衡器的时域脉冲响应序列的长度为
Figure 2012103882599100002DEST_PATH_IMAGE020
，则该响应序列矢量表示为：
Figure 2012103882599100002DEST_PATH_IMAGE022
，
反馈信号 经待求的第1级均衡器处理后的表达式为：
Figure 2012103882599100002DEST_PATH_IMAGE026
，
结合声源信号
Figure 2012103882599100002DEST_PATH_IMAGE028
和第1级均衡器处理后的反馈信号
Figure 2012103882599100002DEST_PATH_IMAGE030
，根据最小二乘准则，计算出使声源信号
Figure 939465DEST_PATH_IMAGE028
和均衡后反馈信号
Figure 982245DEST_PATH_IMAGE030
之间的误差功率取最小值时的均衡器脉冲响应矢量，作为第1次迭代所设计的第1级最优均衡器响应矢量的估计值，其表达式为：
Figure 2012103882599100002DEST_PATH_IMAGE032
，
其中，
Figure 2012103882599100002DEST_PATH_IMAGE034
是
Figure 2012103882599100002DEST_PATH_IMAGE036
维的数据矩阵，其表达式为：
Figure 2012103882599100002DEST_PATH_IMAGE038
；
此时第1级均衡器的参数估计值为：，对第1级均衡器的参数估计值
Figure 2012103882599100002DEST_PATH_IMAGE042
进行归一化处理，得到：
Figure 2012103882599100002DEST_PATH_IMAGE044
其中，是矢量
Figure 343651DEST_PATH_IMAGE042
的转置矢量；反馈信号 经由第1级均衡器
Figure 2012103882599100002DEST_PATH_IMAGE048
处理后，得到经第1次均衡后的反馈信号表达式为：
Figure 2012103882599100002DEST_PATH_IMAGE050
，
结合第1次均衡后的反馈信号
Figure 2012103882599100002DEST_PATH_IMAGE052
和声源信号
Figure 974274DEST_PATH_IMAGE028
，计算出经第1级均衡器处理后的反馈信号与声源信号之间的均方根误差值
Figure 2012103882599100002DEST_PATH_IMAGE054
，其表达式如下：
Figure 2012103882599100002DEST_PATH_IMAGE056
，
b．按照步骤a中均衡器参数的估计过程，继续完成第1、2、…、
Figure 2012103882599100002DEST_PATH_IMAGE058
个均衡器的参数估计，假设在完成第
Figure 189005DEST_PATH_IMAGE058
个均衡器参数估计的基础上，如图3所示，第个均衡器参数的估计过程如下：
经前级均衡器滤波处理后的反馈信号
Figure 2012103882599100002DEST_PATH_IMAGE062
表示为：
Figure 2012103882599100002DEST_PATH_IMAGE064
，
经前
Figure 611207DEST_PATH_IMAGE058
级均衡器滤波处理后的反馈信号再经由待求的第
Figure 144181DEST_PATH_IMAGE060
级均衡器
Figure 2012103882599100002DEST_PATH_IMAGE066
处理后，得到经前
Figure 544594DEST_PATH_IMAGE060
级均衡器处理后的反馈信号表达式为：
Figure 2012103882599100002DEST_PATH_IMAGE068
，
结合声源信号
Figure 585100DEST_PATH_IMAGE028
和经前
Figure 40746DEST_PATH_IMAGE060
级均衡器处理后的反馈信号
Figure 2012103882599100002DEST_PATH_IMAGE070
，根据最小二乘准则，计算出使声源信号
Figure 787991DEST_PATH_IMAGE028
和均衡后反馈信号
Figure 435004DEST_PATH_IMAGE070
之间误差功率取最小值时的均衡器参数矢量，作为第
Figure 762474DEST_PATH_IMAGE060
次迭代所设计的第
Figure 894247DEST_PATH_IMAGE060
级最优均衡器参数矢量的估计值，其表达式为： 
Figure 2012103882599100002DEST_PATH_IMAGE072
，
其中，
Figure 2012103882599100002DEST_PATH_IMAGE074
是
Figure 697597DEST_PATH_IMAGE036
维的数据矩阵，其表达式为：
Figure 2012103882599100002DEST_PATH_IMAGE076
；
此时第
Figure 954660DEST_PATH_IMAGE060
个均衡器的参数估计值为：
Figure 2012103882599100002DEST_PATH_IMAGE078
，对第
Figure 454911DEST_PATH_IMAGE060
个均衡器的参数估计值
Figure 2012103882599100002DEST_PATH_IMAGE080
进行归一化处理，得到：
Figure 2012103882599100002DEST_PATH_IMAGE082
，
其中，
Figure 2012103882599100002DEST_PATH_IMAGE084
是矢量
Figure 439310DEST_PATH_IMAGE080
的转置矢量；经前级均衡器滤波处理后的反馈信号再经由第
Figure 351027DEST_PATH_IMAGE060
级均衡器
Figure 2012103882599100002DEST_PATH_IMAGE090
处理后，得到经前
Figure 78857DEST_PATH_IMAGE060
级均衡器处理后的反馈信号，其表达式为：
Figure 2012103882599100002DEST_PATH_IMAGE092
；
结合经前
Figure 318602DEST_PATH_IMAGE060
级均衡器处理后的反馈信号
Figure 418276DEST_PATH_IMAGE070
和声源信号，计算出经前
Figure 327905DEST_PATH_IMAGE060
级均衡器处理后的反馈信号与声源信号之间的均方根误差值
Figure 2012103882599100002DEST_PATH_IMAGE094
，其表达式如下：
Figure 2012103882599100002DEST_PATH_IMAGE096
；
c．按照步骤b中均衡器参数的估计过程，继续完成第、
Figure 2012103882599100002DEST_PATH_IMAGE100
、…、
Figure 2012103882599100002DEST_PATH_IMAGE102
个均衡器的参数估计，在完成第
Figure 755081DEST_PATH_IMAGE102
个均衡器的参数估计后，均衡根误差值
Figure 2012103882599100002DEST_PATH_IMAGE104
小于用户设定的期望均方根误差，算法不再继续计算下一级均衡器的参数，停止迭代。
更进一步地，在估计均衡器参数的迭代计算过程中，算法一直监视均衡根误差值
Figure 970030DEST_PATH_IMAGE094
的大小，在每次迭代完成后都与用户设定的期望均方根误差
Figure 528444DEST_PATH_IMAGE106
进行比较，并控制迭代循环的运行。例如，在第次迭代开始之前， 如果
Figure 2012103882599100002DEST_PATH_IMAGE108
，那么继续执行
Figure 2012103882599100002DEST_PATH_IMAGE110
次迭代；如果
Figure 2012103882599100002DEST_PATH_IMAGE112
，则停止迭代。
进一步地，步骤（4）中的级联各级均衡器，生成合成均衡器，其实现过程如下：
由
Figure 558552DEST_PATH_IMAGE102
次迭代估计所获得的这
Figure 62346DEST_PATH_IMAGE102
个均衡器的参数估计值分别为：
Figure 737916DEST_PATH_IMAGE048
、
Figure 2012103882599100002DEST_PATH_IMAGE114
、…、
Figure 2012103882599100002DEST_PATH_IMAGE116
，由这
Figure 463602DEST_PATH_IMAGE102
个均衡器级联所形成的合成均衡器
Figure 2012103882599100002DEST_PATH_IMAGE118
的表达式为：
Figure 2012103882599100002DEST_PATH_IMAGE120
，
其中“
Figure 289869DEST_PATH_IMAGE018
”代表时域矢量序列之间的卷积操作。
进一步地，步骤（5）中的利用合成均衡器的参数更新FIR滤波器系数，完成系统频响均衡，其具体实现如下：
在获得合成均衡器的基础上，将合成均衡器的参数估计值用于更新有限脉冲响应滤波器的系数，从而实现系统频响的均衡操作，如图4所示。假设声重放系统输入的声源信号时域序列矢量为：
Figure 2012103882599100002DEST_PATH_IMAGE124
经合成均衡器处理后的声源信号可以表示为：
Figure DEST_PATH_IMAGE126
。
经合成均衡器处理后的声源信号，通过功率放大器放大后送至扬声器端，从而驱动扬声器辐射声波，通过这种均衡处理操作，从扬声器到听者位置点处的声场传递函数得到了均衡，其传递函数响应曲线的峰谷得到了抑制，从而提高了声重放信号的质量。
本发明提供的另一技术方案为：一种离线迭代的声重放系统频响均衡装置，如图5所示，它包括：
声源，是系统待重放的声信息；
数字信号处理器，与所述声源的输出端相连接，按照合成均衡器的参数估计值对用户选定的输入信号进行均衡处理，并将均衡后信号送至功放输入端；
功率放大器，与所述数字信号处理器的输出端相连接，用于对均衡处理后的信号进行功率放大，以驱动扬声器发声；
扬声器，与所述功率放大器的输出端相连接，用于电声转换以便将声源信号重放到空气中。
进一步地，声源为来自于各种模拟装置所产生的模拟声源信号，或者为各种数字装置所产生的数字编码信号，或者为无线网络传输信号，所述的无线网络传输信号为无线发射装置传送来的广播信号并通过无线接收器进行接收和解调获得用户指定的声源信号。当声源为模拟声源信号时，需要借助模数转换器，将模拟信号转换为系统指定的数字输入格式；当声源为数字编码信号时，需要在数字信号处理器内部将该数字编码信号转换为系统指定的数字输入格式；当声源为无线网络传输信号时，需要将无线接收器解调的信号转换为系统制定的数字输入格式。
进一步地，功率放大器，其输入接口可以分为两种类型：即数字输入接口和模拟输入接口。当所述功率放大器具有数字输入接口时，则可以直接对数字信号处理器送来的数字信号进行功率放大处理后再送至扬声器端，因此功率放大器直接与数字信号处理器相连接；如果功率放大器仅具有模拟输入接口，则需要依靠数模转换器，将数字信号处理器送来的数字信号转换为模拟信号后再进行功率放大处理，最后送至扬声器端，因此功率放大器与数字信号处理器之间连接有数模转换器。
进一步地，扬声器并不局限于单一的扬声器单元，扬声器的实现形式可以为单个扬声器单元；也可以为多个扬声器单元组成的扬声器阵列，并且该阵列的形状可以根据扬声器单元数量和实际应用需求进行排列，组成适合于实际应用需求的各种阵列形状。
通过采用上述技术方案，本发明与现有技术相比，其优点在于：
1. 与传统的基于单次LMS或者单次LPC参数估计算法相比，本发明所提出的基于多次迭代的声重放系统频响均衡方法，通过增加级联均衡器的数量和增加每级均衡器的阶数，可以明显提升合成均衡器的通道响应均衡能力，使均衡后的系统整体频响曲线更为平直。
2. 本发明所提出的离线迭代均衡方法，仅需要通过测量仪器获取待均衡系统的脉冲响应曲线，其基于多次迭代的多级均衡器的参数估计过程和合成均衡器的生成过程都是在个人计算机（PC机）上通过软件进行离线计算来获得的。在实际使用均衡器的时候，仅需要在系统硬件中利用合成均衡器的参数估计值去更新FIR滤波器的系数就可以完成系统整体的频响均衡操作，而且这种FIR滤波器可以在低成本的DSP和FPGA等数字信号处理器件中实现，这种离线均衡方式的均衡器硬件制作成本较低，实现较为简单。这种离线计算在线均衡的方法，节省了硬件资源的消耗，从而减少了均衡器的实现成本。
3. 本发明所提出的迭代均衡方法，能够通过增加迭代均衡的次数，即增加级联均衡器的数量，来实现对用户期望的整个宽频带内频响进行更为平直的均衡处理，其对于低频带内频响的均衡能力要明显优于传统的均衡方法。
4. 传统的均衡器参数估计方法，需要通过时频域的变换以获取最小相位响应分量，实现复杂繁锁。与传统的均衡器参数估计方法相比，本发明所提出的迭代均衡方法，直接在时域内分析噪声信号和反馈信号，并直接在时域内完成均衡器的参数估计，其信号处理流程及硬件实现较为简单。
5. 本发明通过对多个已估计的均衡器进行级联生成单一的合成均衡器，在实际应用中，使用单个合成均衡器来执行通道的均衡操作，实现简单可靠。
6. 本发明所采用的离线计算方式，依赖计算机强大的信息处理能力，在解决混响环境下的声重放均衡任务时，能够快速的完成空间多位置点的脉冲响应融合处理，便于处理复杂的计算任务，从而可以解决更为复杂环境下的声场均衡任务。
附图说明
图1表示本发明的一种离线迭代的声重放系统频响均衡方法和装置的信号处理流程图；
图2表示本发明的一种离线迭代的声重放系统频响均衡方法和装置的三种不同环节的均衡处理示意图；
图3表示本发明的一种离线迭代的声重放系统频响均衡方法和装置的多级均衡器的参数估计示意图，其中为白噪声声源，
Figure DEST_PATH_IMAGE130
为经均衡处理后的声源信号；
图4表示本发明的一种离线迭代的声重放系统频响均衡方法和装置的合成均衡器的实现过程示意图，其中为输入声源信号，
Figure DEST_PATH_IMAGE134
为经合成均衡器处理后的输入信号；
图5表示本发明的一种离线迭代的声重放系统频响均衡装置的各组成模块的示意图；
图6表示本发明实施例中待均衡声重放系统的脉冲响应曲线图；
图7表示本发明实施例中选定噪声源信号的时域波形图；
图8表示本发明实施例中利用噪声序列和脉冲响应序列卷积所生成反馈信号的时域波形图；
图9表示本发明实施例中系统未施加均衡、经1次迭代均衡和经10次迭代均衡后系统频响曲线对比图。
图中标号为：
1、声源；2、数字信号处理器；3、功率放大器；4、扬声器。
具体实施方式
下面结合附图对本发明的较佳实施例进行详细阐述，以使本发明的优点和特征能更易于被本领域技术人员理解，从而对本发明的保护范围做出更为清楚明确的界定。
目前，传统的声重放系统均衡方法，都是基于对系统脉冲响应函数进行分析，拟合这些脉冲响应函数的零极点模型，再通过零极点倒置找出系统的逆滤波器响应，从而获得了声重放系统的均衡器参数。这些方法的参数估计过程都是依赖于最小均方误差（LMS）算法或者线性预测编码（LPC）算法的单次估计来计算均衡器的参数，这种基于单次估计方法所获得的逆滤波器参数与理想的逆滤波器参数之间仍然存在着一定程度的偏差，这些偏差将造成均衡后的声重放系统频响曲线在一些频带内仍有较为明显的峰谷起伏性，仍未达到较为理想的频响平直特性。为了克服基于单次LMS或者单次LPC参数估计算法，在系统均衡器参数估计方面所存在的一定误差缺陷，本发明提出了一种离线迭代的声重放系统频响均衡方法和装置，通过采用多次迭代的最小二乘估计方法逐级计算出多个级联均衡器的参数值，这些各级均衡器级联所形成的合成均衡器能够更好的逼近于系统理想的逆滤波器响应，从而减少了均衡器的参数估计误差，保证了均衡后系统频响曲线具有更好的平直特性。本发明通过增加级联均衡器数量和增加各级均衡器的阶数,可以明显提升均衡器的响应均衡能力,同时通过离线计算在线均衡的方式,依赖计算机强大的信息处理能力，能够解决复杂环境下的声场均衡任务。
 如图5所示，制作一个依据本发明的离线迭代的声重放系统频响均衡装置，其主体由声源1、数字信号处理器2、功率放大器3、扬声器4等组成。
声源1，为用户指定的待重放信号。
数字信号处理器2，与所述声源1的输出端相连接，在硬件实现上可以由DSP或者FPGA作为核心处理器来实现，基于合成均衡器的参数估计值对有限脉冲响应滤波器的系数进行更新，从而实现对声重放系统的均衡处理。
功率放大器3，与所述数字信号处理器2的输出端相连接，对数字信号处理器2送入的数字信号进行数模变换和功率放大处理。
扬声器4，与所述功率放大器3的输出端相连接，实现电声变换，在空气中重放声信号。扬声器4为口径为3.5寸、额定功率为10瓦、直流电阻为4欧姆、置于封闭箱体内的扬声器。
在本实施例中，我们利用测量仪器测试和记录扬声器前方轴线上1 米远处的脉冲响应曲线，如图6所示，测量仪器记录的脉冲响应序列的长度为200点、采样频率为23.8KHz。选定噪声序列为白噪声序列，其采样频率也为23.8KHz，比特位数为16，序列长度为2048点，其时域信号波形如图7所示。利用噪声序列和脉冲响应序列的卷积获取反馈信号，其时域信号波形如图8所示。假设待估计的各级均衡器的阶数均为600，设置迭代均衡的次数为10。
图9给出了在未均衡、经1次迭代均衡和经10次迭代均衡三种情况下，系统频响曲线的对比图。对比这三组曲线可以看出，在未施加均衡器的情况下系统频响曲线在1.5KHz～4.5KHz的频带范围内存在着非常明显的峰值；在第1次迭代均衡处理后，系统在在1.5KHz～4.5KHz的频带范围内的峰值已经得到了消除，但是在1.5KHz频点附近的区域内系统频响曲线仍有小量的起伏，同时在100 Hz～200 Hz的频带内，系统频响曲线仍有较大程度的起伏；在经10次迭代均衡处理后，系统在1.5KHz频点附近区域内的少量起伏已经得到了消除，同时在100 Hz～200 Hz频带内的幅度峰值也得到了较大程度的抑制。对比第1次迭代均衡和第10次迭代均衡处理后的频响曲线，可以看出：通过增加迭代次数，本发明所提出的迭代均衡方法，能够明显改善系统均衡后频响曲线的平直程度，这说明本发明所提出的多次迭代均衡方法，较传统均衡方法相比，具有更好的均衡效果，经其均衡后的频响曲线将会更为平直。
上述实施例只为说明本发明的技术构思及特点，其目的在于让熟悉此项技术的人士能够了解本发明的内容并据以实施，并不能以此限制本发明的保护范围。凡根据本发明精神实质所作的等效变化或修饰，都应涵盖在本发明的保护范围之内。

Claims (13)

1. 一种离线迭代的声重放系统频响均衡方法，包括如下步骤：
（1）利用测量仪器获取待均衡系统的脉冲响应序列；
（2）选取噪声信号，利用噪声信号与脉冲响应序列的卷积获取反馈信号； 
（3）结合噪声信号和反馈信号，利用逐次迭代的最小二乘准则依次估计出各级均衡器的参数；
（4）级联各级均衡器，生成合成均衡器；
（5）利用合成均衡器的参数更新有限脉冲响应滤波器系数，完成系统频响均衡。
2. 根据权利要求1所述的离线迭代的声重放系统频响均衡方法，其特征在于：所述的步骤（1）中的利用测量仪器获取待均衡系统的脉冲响应序列，具体步骤为：通过测量仪器发出噪声信号，用传声器采集和记录由噪声激励所产生的系统脉冲响应序列；如果所述声重放系统仅需要均衡信号处理和功放所组成电路部分的系统响应，则在脉冲响应测试时，仅测试和记录由噪声激励所产生的功放输出端的脉冲响应序列；如果声重放系统需要均衡电路部分和扬声器部分所组成系统的响应，则在脉冲响应测量时，将传声器置于空间期望位置点，测量和记录由噪声激励所产生的该点处的脉冲响应序列；如果声重放系统需要均衡电路部分、扬声器部分和外部环境三者所组成系统的响应，则在脉冲响应测量时，将传声器依次置于空间多个位置点，记录多个位置点上的脉冲响应序列，并按照一定的融合方式对多位置点的脉冲响应进行平均化处理，获得平均后的脉冲响应序列。
3. 根据权利要求1所述的离线迭代的声重放系统频响均衡方法，其特征在于：所述步骤（2）中的噪声信号为指定带宽的白噪声序列或者最大长度序列所产生的噪声信号，根据实际测量所获得脉冲响应序列的长度                                               来选取上述噪声信号的序列长度
Figure 2012103882599100001DEST_PATH_IMAGE004
，
Figure 2012103882599100001DEST_PATH_IMAGE006
。
4. 根据权利要求1所述的离线迭代的声重放系统频响均衡方法，其特征在于：所述步骤（2）中的利用噪声信号与脉冲响应序列的卷积获取反馈信号，其具体实现如下：
假设噪声信号的时域序列矢量为：
Figure 2012103882599100001DEST_PATH_IMAGE008
，
其中，
Figure 2012103882599100001DEST_PATH_IMAGE010
是噪声信号时域离散序列的采样点数；假设经测量仪器测量获得的声重放系统的时域脉冲响应序列表达式：
Figure 2012103882599100001DEST_PATH_IMAGE012
，
其中，
Figure 2012103882599100001DEST_PATH_IMAGE014
为系统时域脉冲响应的序列长度；联合噪声序列和声重放系统的时域脉冲响应序列，获得反馈信号的时域序列矢量表达式为：
Figure 2012103882599100001DEST_PATH_IMAGE016
，
其中“
Figure 2012103882599100001DEST_PATH_IMAGE018
”代表两个序列矢量之间进行卷积操作。
5. 根据权利要求4所述的离线迭代的声重放系统频响均衡方法，其特征在于：在进行噪声序列和脉冲响应序列的卷积操作时，已经将卷积后所获得反馈信号序列的头部切除掉长度为
Figure 355162DEST_PATH_IMAGE014
的序列部分，获得了长度为
Figure 806872DEST_PATH_IMAGE010
的反馈信号序列。
6. 根据权利要求1所述的离线迭代的声重放系统频响均衡方法，其特征在于：所述步骤（3）中的结合噪声信号和反馈信号，利用逐次迭代的最小二乘准则依次估计出各级均衡器的参数，其具体实现流程如下：
a．在第1次迭代时，假定待求的第1级均衡器的时域脉冲响应序列的长度为，该响应序列矢量则表示为：
Figure 2012103882599100001DEST_PATH_IMAGE022
，
反馈信号
Figure 2012103882599100001DEST_PATH_IMAGE024
经待求的第1级均衡器处理后的表达式为：
Figure 2012103882599100001DEST_PATH_IMAGE026
，
结合声源信号
Figure 2012103882599100001DEST_PATH_IMAGE028
和第1级均衡器处理后的反馈信号
Figure 2012103882599100001DEST_PATH_IMAGE030
，根据最小二乘准则，计算出使声源信号和均衡后反馈信号
Figure 543895DEST_PATH_IMAGE030
之间的误差功率取最小值时的均衡器脉冲响应矢量，作为第1次迭代所设计的第1级最优均衡器响应矢量的估计值，其表达式为：
Figure DEST_PATH_IMAGE032
，
其中，
Figure DEST_PATH_IMAGE034
是
Figure DEST_PATH_IMAGE036
维的数据矩阵，其表达式为：
Figure DEST_PATH_IMAGE038
；
此时第1级均衡器的参数估计值为：
Figure DEST_PATH_IMAGE040
，对第1级均衡器的参数估计值
Figure DEST_PATH_IMAGE042
进行归一化处理，得到：
Figure DEST_PATH_IMAGE044
其中，是矢量
Figure 824091DEST_PATH_IMAGE042
的转置矢量；反馈信号
Figure 130308DEST_PATH_IMAGE024
经由第1级均衡器
Figure DEST_PATH_IMAGE048
处理后，得到经第1次均衡后的反馈信号表达式为：
Figure DEST_PATH_IMAGE050
，
结合第1次均衡后的反馈信号和声源信号
Figure 226307DEST_PATH_IMAGE028
，计算出经第1级均衡器处理后的反馈信号与声源信号之间的均方根误差值
Figure DEST_PATH_IMAGE054
，其表达式如下：
；
b．按照步骤a中均衡器参数的估计过程，继续完成第1、2、…、个均衡器的参数估计，假设在完成第
Figure 673600DEST_PATH_IMAGE058
个均衡器参数估计的基础上，第个均衡器参数的估计过程如下：
经前级均衡器滤波处理后的反馈信号表示为：
Figure DEST_PATH_IMAGE064
，
经前
Figure 407956DEST_PATH_IMAGE058
级均衡器滤波处理后的反馈信号再经由待求的第
Figure 404917DEST_PATH_IMAGE060
级均衡器
Figure DEST_PATH_IMAGE066
处理后，得到经前
Figure 964598DEST_PATH_IMAGE060
级均衡器处理后的反馈信号表达式为：
Figure DEST_PATH_IMAGE068
，
结合声源信号
Figure 559921DEST_PATH_IMAGE028
和经前
Figure 419293DEST_PATH_IMAGE060
级均衡器处理后的反馈信号
Figure DEST_PATH_IMAGE070
，根据最小二乘准则，计算出使声源信号
Figure 277396DEST_PATH_IMAGE028
和均衡后反馈信号
Figure 499824DEST_PATH_IMAGE070
之间误差功率取最小值时的均衡器参数矢量，作为第
Figure 635139DEST_PATH_IMAGE060
次迭代所设计的第
Figure 727728DEST_PATH_IMAGE060
级最优均衡器参数矢量的估计值，其表达式为： 
Figure DEST_PATH_IMAGE072
，
其中，
Figure DEST_PATH_IMAGE074
是维的数据矩阵，其表达式为：
Figure DEST_PATH_IMAGE076
；
此时第
Figure 300845DEST_PATH_IMAGE060
个均衡器的参数估计值为：
Figure DEST_PATH_IMAGE078
，对第
Figure 415301DEST_PATH_IMAGE060
个均衡器的参数估计值
Figure DEST_PATH_IMAGE080
进行归一化处理，得到：
，
其中，
Figure DEST_PATH_IMAGE084
是矢量
Figure 183186DEST_PATH_IMAGE080
的转置矢量；经前
Figure DEST_PATH_IMAGE086
级均衡器滤波处理后的反馈信号再经由第
Figure 78199DEST_PATH_IMAGE060
级均衡器
Figure DEST_PATH_IMAGE090
处理后，得到经前
Figure 20923DEST_PATH_IMAGE060
级均衡器处理后的反馈信号，其表达式为：
Figure DEST_PATH_IMAGE092
；
结合经前
Figure 930498DEST_PATH_IMAGE060
级均衡器处理后的反馈信号
Figure 551841DEST_PATH_IMAGE070
和声源信号，计算出经前级均衡器处理后的反馈信号与声源信号之间的均方根误差值
Figure DEST_PATH_IMAGE094
，其表达式如下：
Figure DEST_PATH_IMAGE096
；
c．按照步骤b中均衡器参数的估计过程，继续完成第
Figure DEST_PATH_IMAGE098
、
Figure DEST_PATH_IMAGE100
、…、
Figure DEST_PATH_IMAGE102
个均衡器的参数估计，在完成第
Figure 782597DEST_PATH_IMAGE102
个均衡器的参数估计后，均衡根误差值
Figure DEST_PATH_IMAGE104
小于用户设定的期望均方根误差，算法不再继续计算下一级均衡器的参数，停止迭代。
7. 根据权利要求6所述的离线迭代的声重放系统频响均衡方法，其特征在于：在估计均衡器参数的迭代计算过程中，算法一直监视均衡根误差值
Figure 637159DEST_PATH_IMAGE094
的大小，在每次迭代完成后都与用户设定的期望均方根误差
Figure 385059DEST_PATH_IMAGE106
进行比较，并控制迭代循环的运行，在下一次迭代开始之前， 如果
Figure DEST_PATH_IMAGE108
，那么继续执行下一次迭代；如果
Figure DEST_PATH_IMAGE110
，则停止迭代。
8. 根据权利要求1所述的离线迭代的声重放系统频响均衡方法，其特征在于：所述步骤（4）中的级联各级均衡器，生成合成均衡器，其实现过程如下：
由
Figure 664993DEST_PATH_IMAGE102
次迭代估计所获得的这
Figure 483914DEST_PATH_IMAGE102
个均衡器的参数估计值分别为：
Figure 528618DEST_PATH_IMAGE048
、、…、
Figure DEST_PATH_IMAGE114
，由这
Figure 682256DEST_PATH_IMAGE102
个均衡器级联所形成的合成均衡器
Figure DEST_PATH_IMAGE116
的表达式为：
，
其中“
Figure 380261DEST_PATH_IMAGE018
”代表时域矢量序列之间的卷积操作。
9. 根据权利要求1所述的离线迭代的声重放系统频响均衡装置，其特征在于：所述步骤（5）中的利用合成均衡器的参数更新有限脉冲响应滤波器系数，完成系统频响均衡，其具体实现如下：
在获得合成均衡器
Figure DEST_PATH_IMAGE120
的基础上，将合成均衡器的参数估计值用于更新有限脉冲响应滤波器的系数，实现系统频响的均衡操作；假设声重放系统输入的声源信号时域序列矢量为：
，
经合成均衡器
Figure 883048DEST_PATH_IMAGE120
处理后的声源信号则表示为：
Figure DEST_PATH_IMAGE124
；
经合成均衡器处理后的声源信号，通过功率放大器放大后送至扬声器端，驱动扬声器辐射声波。
10. 一种离线迭代的声重放系统频响均衡装置，其特征在于：它包括声源（1）、与所述声源（1）的输出端相连接并按照合成均衡器的参数估计值对用户选定的输入信号进行均衡处理并将均衡后信号送至功放输入端的数字信号处理器（2）、与所述数字信号处理器（2）的输出端相连接并用于对均衡处理后的信号进行功率放大以驱动扬声器（4）发声的功率放大器（3）、与所述功率放大器（3）的输出端相连接并用于电声转换以便将声源信号重放到空气中的扬声器（4）所述的声源（1）是系统待重放的声信息。
11. 根据权利要求10所述的离线迭代的声重放系统频响均衡装置，其特征在于：所述的声源（1）为来自于各种模拟装置所产生的模拟声源信号，或者为各种数字装置所产生的数字编码信号，或者为无线网络传输信号，所述的无线网络传输信号为无线发射装置传送来的广播信号并通过无线接收器进行接收和解调获得用户指定的声源信号。
12. 根据权利要求10所述的离线迭代的声重放系统频响均衡装置，其特征在于：当所述功率放大器（3）具有数字输入接口时，则所述的功率放大器（3）直接与数字信号处理器（2）连接；当所述功率放大器（3）仅具有模拟输入接口时，则所述的功率放大器（3）与数字信号处理器（2）之间连接有将数字信号处理器（2）送来的数字信号转换为模拟信号的数模转换器。
13. 根据权利要求10所述的离线迭代的声重放系统频响均衡装置，其特征在于：所述的扬声器（4）为单个扬声器单元、或者为多个扬声器单元组成的根据扬声器单元数量和实际应用需求进行排列的扬声器阵列。
CN2012103882599A 2012-10-15 2012-10-15 离线迭代的声重放系统频响均衡方法和装置 Pending CN102903367A (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN2012103882599A CN102903367A (zh) 2012-10-15 2012-10-15 离线迭代的声重放系统频响均衡方法和装置
PCT/CN2013/085000 WO2014059897A1 (zh) 2012-10-15 2013-10-10 离线迭代的声重放系统频响均衡方法和装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2012103882599A CN102903367A (zh) 2012-10-15 2012-10-15 离线迭代的声重放系统频响均衡方法和装置

Publications (1)

Publication Number Publication Date
CN102903367A true CN102903367A (zh) 2013-01-30

Family

ID=47575570

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2012103882599A Pending CN102903367A (zh) 2012-10-15 2012-10-15 离线迭代的声重放系统频响均衡方法和装置

Country Status (2)

Country Link
CN (1) CN102903367A (zh)
WO (1) WO2014059897A1 (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2014059890A1 (zh) * 2012-10-15 2014-04-24 苏州上声电子有限公司 在线迭代的声重放系统频响均衡方法和装置
WO2014059897A1 (zh) * 2012-10-15 2014-04-24 苏州上声电子有限公司 离线迭代的声重放系统频响均衡方法和装置
CN109873781A (zh) * 2017-12-01 2019-06-11 晨星半导体股份有限公司 符合同轴电缆多媒体联盟标准的信号接收装置及其信号处理方法
CN117676418A (zh) * 2023-12-06 2024-03-08 广州番禺职业技术学院 一种用于混合相位系统中的声场均衡方法及系统

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5025472A (en) * 1987-05-27 1991-06-18 Yamaha Corporation Reverberation imparting device
WO2003015082A1 (en) * 2001-08-07 2003-02-20 Dspfactory Ltd. Sound intelligibilty enchancement using a psychoacoustic model and an oversampled fiolterbank
WO2010014663A2 (en) * 2008-07-29 2010-02-04 Dolby Laboratories Licensing Corporation Method for adaptive control and equalization of electroacoustic channels
CN101645264A (zh) * 2009-07-17 2010-02-10 中国科学院声学研究所 一种基于迭代最小二乘-预滤波的最优主动时反聚焦方法
CN102447446A (zh) * 2011-12-09 2012-05-09 苏州上声电子有限公司 基于振动元件运动状态反馈的扬声器频响均衡方法和装置
CN102711010A (zh) * 2012-05-29 2012-10-03 苏州上声电子有限公司 利用二次剩余序列的扬声器阵列宽带声场控制方法和装置

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102903367A (zh) * 2012-10-15 2013-01-30 苏州上声电子有限公司 离线迭代的声重放系统频响均衡方法和装置

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5025472A (en) * 1987-05-27 1991-06-18 Yamaha Corporation Reverberation imparting device
WO2003015082A1 (en) * 2001-08-07 2003-02-20 Dspfactory Ltd. Sound intelligibilty enchancement using a psychoacoustic model and an oversampled fiolterbank
WO2010014663A2 (en) * 2008-07-29 2010-02-04 Dolby Laboratories Licensing Corporation Method for adaptive control and equalization of electroacoustic channels
CN101645264A (zh) * 2009-07-17 2010-02-10 中国科学院声学研究所 一种基于迭代最小二乘-预滤波的最优主动时反聚焦方法
CN102447446A (zh) * 2011-12-09 2012-05-09 苏州上声电子有限公司 基于振动元件运动状态反馈的扬声器频响均衡方法和装置
CN102711010A (zh) * 2012-05-29 2012-10-03 苏州上声电子有限公司 利用二次剩余序列的扬声器阵列宽带声场控制方法和装置

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2014059890A1 (zh) * 2012-10-15 2014-04-24 苏州上声电子有限公司 在线迭代的声重放系统频响均衡方法和装置
WO2014059897A1 (zh) * 2012-10-15 2014-04-24 苏州上声电子有限公司 离线迭代的声重放系统频响均衡方法和装置
CN109873781A (zh) * 2017-12-01 2019-06-11 晨星半导体股份有限公司 符合同轴电缆多媒体联盟标准的信号接收装置及其信号处理方法
CN117676418A (zh) * 2023-12-06 2024-03-08 广州番禺职业技术学院 一种用于混合相位系统中的声场均衡方法及系统
CN117676418B (zh) * 2023-12-06 2024-05-24 广州番禺职业技术学院 一种用于混合相位系统中的声场均衡方法及系统

Also Published As

Publication number Publication date
WO2014059897A1 (zh) 2014-04-24

Similar Documents

Publication Publication Date Title
CN109074816B (zh) 远场自动语音识别预处理
CN102883243B (zh) 在线迭代的声重放系统频响均衡方法
JP6572894B2 (ja) 情報処理装置、情報処理方法及びプログラム
CN103634726B (zh) 一种扬声器自动均衡方法
CN106535076B (zh) 一种立体声音响系统的空间校准方法及其移动终端设备
KR20150128753A (ko) 등화 필터 계수 결정기, 장치, 등화 필터 계수 프로세서, 시스템 및 방법
GB2519676A (en) Method for optimizing the performance of a loudspeaker to compensate for low frequency room modes
CN105792072B (zh) 一种音效处理方法、装置及终端
US9860641B2 (en) Audio output device specific audio processing
CN112562624B (zh) 主动降噪滤波器设计方法、降噪方法、系统及电子设备
CN112584277B (zh) 一种室内音频均衡的方法
JP2013516826A (ja) 音声処理動作の構成の決定
CN102903367A (zh) 离线迭代的声重放系统频响均衡方法和装置
WO2022256577A1 (en) A method of speech enhancement and a mobile computing device implementing the method
US10582299B1 (en) Modeling room acoustics using acoustic waves
CN101516055A (zh) 一种能够模拟不同音响环境的三维回声音效方法和设备
TW201118719A (en) System and method for modifying an audio signal
CN113921007B (zh) 提升远场语音交互性能的方法和远场语音交互系统
JP2017083566A (ja) 雑音抑圧装置、雑音抑圧方法、およびプログラム
Deng et al. Modeling and estimating acoustic transfer functions of external ears with or without headphones
WO2022209171A1 (ja) 信号処理装置、信号処理方法およびプログラム
WO2023287782A1 (en) Data augmentation for speech enhancement
US11830471B1 (en) Surface augmented ray-based acoustic modeling
JP6721010B2 (ja) 機械学習方法および機械学習装置
JP4229435B2 (ja) 音場シミュレーション装置、音場シミュレーション方法、コンピュータプログラム、プログラム記録媒体

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: 20130130