[go: up one dir, main page]

CN113630690B - 一种扬声器系统的非线性特性的确定方法及装置 - Google Patents

一种扬声器系统的非线性特性的确定方法及装置 Download PDF

Info

Publication number
CN113630690B
CN113630690B CN202111037509.XA CN202111037509A CN113630690B CN 113630690 B CN113630690 B CN 113630690B CN 202111037509 A CN202111037509 A CN 202111037509A CN 113630690 B CN113630690 B CN 113630690B
Authority
CN
China
Prior art keywords
loudspeaker
equivalent circuit
variable matrix
determining
equivalent
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
Application number
CN202111037509.XA
Other languages
English (en)
Other versions
CN113630690A (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.)
Guoguang Electric Co Ltd
Original Assignee
Guoguang Electric Co Ltd
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 Guoguang Electric Co Ltd filed Critical Guoguang Electric Co Ltd
Priority to CN202111037509.XA priority Critical patent/CN113630690B/zh
Publication of CN113630690A publication Critical patent/CN113630690A/zh
Application granted granted Critical
Publication of CN113630690B publication Critical patent/CN113630690B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • 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
    • 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
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Acoustics & Sound (AREA)
  • Signal Processing (AREA)
  • Circuit For Audible Band Transducer (AREA)

Abstract

本发明实施例公开了一种扬声器系统的非线性特性的确定方法及装置。包括：对初始信号在频域进行预处理，并将预处理后的多个频率信号分别加载于扬声器系统的输入端；对于加载了各频率信号的扬声器系统，获取所述扬声器系统中的声学元件的结构属性参数；其中，所述结构属性参数包括如下至少一项：扬声器数量、扬声器的串并联属性、蠕变效应参数、热效应参数、倒相管数量及无源辐射器PR数量；根据所述结构属性参数构建所述扬声器系统的综合等效电路模型；基于所述等效电路模型确定所述扬声器系统在各频率信号下的非线性特性值。本发明实施例公开的扬声器系统的非线性特性的确定方法，可以提高确定非线性特性的准确性。

Description

一种扬声器系统的非线性特性的确定方法及装置
技术领域
本发明实施例涉及扬声器技术领域，尤其涉及一种扬声器系统的非线性特性的确定方法、装置及设备。
背景技术
扬声器在较大功率工作时，会引起振动系统的非线性振动。振动系统的非线性振动引起扬声器的非线性失真，表现出声音不纯净、不好听，因此研究和抑制扬声器振动中的非线性失真是中高档扬声器设计中的要点。
现有技术中，在确定非线性特性时，由于未考虑扬声器系统的结构属性参数(如：扬声器的串并联、蠕变效应参数、热效应参数、倒相管数量及无源辐射器数量)，使得非线性特性不够准确。
发明内容
本发明实施例提供一种扬声器系统的非线性特性的确定方法及装置，可以提高确定非线性特性的准确性。
第一方面，本发明实施例提供了一种扬声器系统的非线性特性的确定方法，包括：
对初始信号在频域进行预处理，并将预处理后的多个频率信号分别加载于扬声器系统的输入端；
对于加载了各频率信号的扬声器系统，获取所述扬声器系统中的声学元件的结构属性参数；其中，所述结构属性参数包括如下至少一项：扬声器数量、扬声器的串并联属性、蠕变效应参数、热效应参数、倒相管数量及无源辐射器PR数量；
根据所述结构属性参数构建所述扬声器系统的综合等效电路模型；
基于所述等效电路模型确定所述扬声器系统在各频率信号下的非线性特性值。
第二方面，本发明实施例还提供了一种扬声器系统的非线性特性的确定装置，包括：
信号预处理模块，用于对初始信号在频域进行预处理，并将预处理后的多个频率信号分别加载于扬声器系统的输入端；
结构属性参数获取模块，用于对于加载了各频率信号的扬声器系统，获取所述扬声器系统中的声学元件的结构属性参数；其中，所述结构属性参数包括如下至少一项：扬声器数量、扬声器的串并联属性、蠕变效应参数、热效应参数、倒相管数量及无源辐射器PR数量；
综合等效电路模型构建模块，用于根据所述结构属性参数构建所述扬声器系统的综合等效电路模型；
非线性特性值确定模块，用于基于所述等效电路模型确定所述扬声器系统在各频率信号下的非线性特性值。
本发明实施例公开了一种扬声器系统的非线性特性的确定方法及装置。对初始信号在频域进行预处理，并将预处理后的多个频率信号分别加载于扬声器系统的输入端；对于加载了各频率信号的扬声器系统，获取扬声器系统中的声学元件的结构属性参数；其中，结构属性参数包括如下至少一项：扬声器数量、扬声器的串并联属性、蠕变效应参数、热效应参数、倒相管数量及无源辐射器PR数量；根据结构属性参数构建扬声器系统的综合等效电路模型；基于等效电路模型确定扬声器系统在各频率信号下的非线性特性值。本发明实施例公开的扬声器系统的非线性特性的确定方法，根据结构属性参数构建扬声器系统的综合等效电路模型，并基于等效电路模型确定扬声器系统在各频率信号下的非线性特性值，可以提高确定非线性特性的准确性。
附图说明
图1是本发明实施例一中的一种扬声器系统的非线性特性的确定方法的流程图；
图2是本发明实施例一中的信号处理模块的结构示意图；
图3是本发明实施例一中的综合等效电路模型的结构示意图；
图4是本发明实施例二中的一种扬声器系统的非线性特性的确定装置的结构示意图。
具体实施方式
下面结合附图和实施例对本发明作进一步的详细说明。可以理解的是，此处所描述的具体实施例仅仅用于解释本发明，而非对本发明的限定。另外还需要说明的是，为了便于描述，附图中仅示出了与本发明相关的部分而非全部结构。
实施例一
图1为本发明实施例一提供的一种扬声器系统的非线性特性的确定方法的流程图，本实施例可适用于确定扬声器系统的非线性特性的情况，该方法可以由扬声器系统的非线性特性的确定装置来执行，如图1所示，该方法具体包括如下步骤：
步骤110，对初始信号在频域进行预处理，并将预处理后的多个频率信号分别加载于扬声器系统的输入端。
其中，扬声器系统可以是音箱系统。初始信号在整个频段内为恒定幅值的电压信号，其中频段可以是20Hz-2wHz。对初始信号在频域进行预处理可以理解为对初始信号在频域进行滤波和/或均衡处理。
本实施例中，可以采用信号处理模块对初始信号进行处理，图2为本实施例中的信号处理模块的结构示意图，如图2所示，该模块包括滤波器及均衡器等处理单元，在给定初始信号时，在频域中，将信号依次输入滤波器和均衡器(Equalizer，EQ)。具体的，对初始信号进行滤波、均衡及其他处理等一系列处理后，获得加载于扬声器系统输入的信号。每次抽取一个频率的信号加载于扬声器系统的输入端，从而获取扬声器系统各频率对应的非线性特性值。
步骤120，对于加载了各频率信号的扬声器系统，获取扬声器系统中的声学元件的结构属性参数。
其中，结构属性参数包括如下至少一项：扬声器数量、扬声器的串并联属性、蠕变效应参数、热效应参数、倒相管数量及无源辐射器(Passive Radiators，PR)数量。扬声器系统可以包括扬声器单元、闭箱、导向箱及PR箱。本实施例中，需要获取扬声器单元、闭箱、导向箱及PR箱中的声学元件以及各声学元件的结构属性参数。
步骤130，根据结构属性参数构建扬声器系统的综合等效电路模型。
本实施例中，可以基于电学-力学-声学类比方法构建扬声器系统的综合等效电路模型。具体的，根据结构属性参数构建扬声器系统的综合等效电路模型的过程可以是：根据结构属性参数确定扬声器系统对应的等效电路元件以及各等效电路元件的连接关系；基于连接关系连接各等效电路元件，获得综合等效电路模型。
其中，综合等效电路模型包括：加载于扬声器系统的频率信号对应的输入电压ue、扬声器单元直阻Re、音圈电感Le(x)、力电耦合因子Bl(x)、扬声器有效振动质量对应的电感Mms、扬声器悬置系统等效力阻对应的电阻Rms(v)、扬声器悬置系统等效劲度系数对应的电容Kms(x)、背腔等效力学劲度系数对应的电容Kbox、倒相管中空气力学等效质量对应的电感Mtube、倒相管中空气等效力阻对应的电阻Rtube、PR有效振动质量经过次级到初级的阻抗换算后对应的电感Mpr、PR等效力阻经过次级到初级的阻抗换算后对应的电阻Rpr及PR等效劲度系数经过次级到初级的阻抗换算后对应的电容Kpr。
图3为本实施例中的综合等效电路模型的结构示意图，如图3所示，ue表示扬声器系统两端的输入电压(单位：V)信号，Le(x)为音圈电感(单位：H)，Bl(x)为力电耦合因子(单位：T*m)；Mms为扬声器有效振动质量(单位：kg)；Rms(v)为扬声器悬置系统等效力阻(单位：kg/s)；Kms(x)为扬声器悬置系统等效劲度系数(单位：N/m)；Kbox为背腔等效力学劲度系数(单位：N/m)；Mtube为倒相管中空气力学等效质量(单位：kg)；Rtube为倒相管中空气等效力阻(单位：kg/s)；Mpr为PR有效振动质量经过次级到初级的阻抗换算后的值(单位：kg)；Rpr为PR等效力阻经过次级到初级的阻抗换算后的值(单位：kg/s)；Kpr为PR等效劲度系数经过次级到初级的阻抗换算后的值(单位：N/m)。
步骤140，基于等效电路模型确定扬声器系统在各频率信号下的非线性特性值。
本实施例中，基于等效电路模型确定扬声器系统在各频率信号下的非线性特性值的过程可以是：基于等效电路模型构建与扬声器系统的工作变量相关的电学力学方程组；对电力学方程组进行求解，获得当前时刻的工作变量矩阵；获取工作变量矩阵的时间变化规律；基于时间变化规律和当前时刻的工作变量矩阵确定各时刻的工作变量矩阵；基于各时刻的工作变量矩阵确定扬声器系统在各频率信号下的非线性特性值。
其中，工作变量包括音圈电流i、扬声器振动系位移x、扬声器振动系振动速度v、倒相管中空气柱位移经过次级到初级的阻抗换算后的值x1、倒相管中空气柱振动速度经过次级到初级的阻抗换算后的值v1、PR位移经过次级到初级的阻抗换算后的值x2、PR振动速度经过次级到初级的阻抗换算后的值v2。
具体的，基于等效电路模型构建与扬声器系统的工作变量相关的电力学方程组按照如下公式表示：
本实施例中，扬声器系统包括如下至少一个单元：扬声器单元、闭箱、倒相箱及PR箱。需要分别确定扬声器单元、闭箱、倒相箱及PR箱的非线性特性。
若确定扬声器单元的声学非线性特性，则需要计算扬声器单元对应的工作变量矩阵，此时，在求解上述方程组时，令Kbox＝0，1/Mtube＝0，Rtube＝0，1/Mpr＝0，Rpr＝0，Kpr＝0。
若确定闭箱的声学非线性特性，则需要计算闭箱对应的工作变量矩阵，此时，在求解上述方程组时，令1/Mtube＝0，Rtube＝0，1/Mpr＝0，Rpr＝0，Kpr＝0。
若确定倒相箱的声学非线性特性，则需要计算倒相箱对应的工作变量矩阵，此时，在求解上述方程组时，令1/Mpr＝0，Rpr＝0，Kpr＝0，且计算Kbox时的背腔体积包含倒相管内空气体积。
若确定PR箱的声学非线性特性，则需要计算PR箱对应的工作变量矩阵，此时，在求解上述方程组时，令1/Mtube＝0，Rtube＝0。
具体的，获取所述工作变量矩阵的时间变化规律按照如下公式表示：Y(n+1)＝F*Y(n)+G*ue，其中，F和G分别为系数矩阵；Y(n+1)下一时刻的工作变量矩阵，Y(n)为当前时刻的工作变量矩阵，矩阵Y为1列7行的矩阵，且Y＝[i x v x1 v1 x2 v2]′，F和G表示如下：
其中，
Ts为当前时刻与下一时刻的时间间隔。
具体的，在获得给定频率信号下各时刻的工作变量矩阵后，截取各时刻的工作变量矩阵中处于中间稳态部分的数据，基于中间稳态部分的数据计算该给定频率信号下的声压级SPL、阻抗IMP、扬声器位移DIS_DRV、总谐波失真THD、二阶谐波失真2HD、3阶谐波失真3HD、无缘辐射器位移DISPR和倒相管中空气流速vair。
其中，基于中间稳态部分的数据按照如下公式计算所述扬声器系统在各频率信号下的非线性特性值：
其中，t为时域求解中的时长；r为声压级测试中麦克风到扬声器的垂直距离；N为求解的最大谐波失真阶数；pnf为基频为f时的第n阶谐波声压，由时域声压经过傅里叶变换之后得到，Stube为导向管横截面积，ρ为空气密度，Sd为扬声器有效振动面积，Sp为无缘辐射器(PR)的有效振动面积。
以频率为循环参数，对每一个频率进行上述计算之后，便能得到扬声器系统声压级、THD、位移、阻抗、PR位移、导向管中空气流速等声学参数随频率的变化规律，即得到扬声器系统的声学非线性特性。本申请中，同时输出的声学非线性特性包括了频响、阻抗、扬声器位移、PR位移、导向管中空气流速、失真等，输出声学非线性特性更加全面。
本实施例的技术方案，对初始信号在频域进行预处理，并将预处理后的多个频率信号分别加载于扬声器系统的输入端；对于加载了各频率信号的扬声器系统，获取扬声器系统中的声学元件的结构属性参数；其中，结构属性参数包括如下至少一项：扬声器数量、扬声器的串并联属性、蠕变效应参数、热效应参数、倒相管数量及无源辐射器PR数量；根据结构属性参数构建扬声器系统的综合等效电路模型；基于等效电路模型确定扬声器系统在各频率信号下的非线性特性值。本发明实施例公开的扬声器系统的非线性特性的确定方法，根据结构属性参数构建扬声器系统的综合等效电路模型，并基于等效电路模型确定扬声器系统在各频率信号下的非线性特性值，可以提高确定非线性特性的准确性。
实施例二
图4是本发明实施例二提供的一种扬声器系统的非线性特性的确定装置的结构示意图。如图4所示，该装置包括：
信号预处理模块210，用于对初始信号在频域进行预处理，并将预处理后的多个频率信号分别加载于扬声器系统的输入端；
结构属性参数获取模块220，用于对于加载了各频率信号的扬声器系统，获取所述扬声器系统中的声学元件的结构属性参数；其中，所述结构属性参数包括如下至少一项：扬声器数量、扬声器的串并联属性、蠕变效应参数、热效应参数、倒相管数量及无源辐射器PR数量；
综合等效电路模型构建模块230，用于根据所述结构属性参数构建所述扬声器系统的综合等效电路模型；
非线性特性值确定模块240，用于基于所述等效电路模型确定所述扬声器系统在各频率信号下的非线性特性值。
可选的，信号预处理模块210，还用于：
对所述初始信号在频域进行滤波和/或均衡处理。
可选的，综合等效电路模型构建模块230，还用于：
根据所述结构属性参数确定所述扬声器系统对应的等效电路元件以及各等效电路元件的连接关系；
基于所述连接关系连接各所述等效电路元件，获得综合等效电路模型；
其中，所述综合等效电路模型包括：加载于扬声器系统的频率信号对应的输入电压ue、扬声器单元直阻Re、音圈电感Le(x)、力电耦合因子Bl(x)、扬声器有效振动质量对应的电感Mms、扬声器悬置系统等效力阻对应的电阻Rms(v)、扬声器悬置系统等效劲度系数对应的电容Kms(x)、背腔等效力学劲度系数对应的电容Kbox、倒相管中空气力学等效质量对应的电感Mtube、倒相管中空气等效力阻对应的电阻Rtube、PR有效振动质量经过次级到初级的阻抗换算后对应的电感Mpr、PR等效力阻经过次级到初级的阻抗换算后对应的电阻Rpr及PR等效劲度系数经过次级到初级的阻抗换算后对应的电容Kpr。
可选的，非线性特性值确定模块240，还用于：
基于所述等效电路模型构建与所述扬声器系统的工作变量相关的电学力学方程组；
对所述电力学方程组进行求解，获得当前时刻的工作变量矩阵；
其中，所述工作变量包括音圈电流i、扬声器振动系位移x、扬声器振动系振动速度v、倒相管中空气柱位移经过次级到初级的阻抗换算后的值x1、倒相管中空气柱振动速度经过次级到初级的阻抗换算后的值v1、PR位移经过次级到初级的阻抗换算后的值x2、PR振动速度经过次级到初级的阻抗换算后的值v2；
获取所述工作变量矩阵的时间变化规律；
基于所述时间变化规律和当前时刻的工作变量矩阵确定各时刻的工作变量矩阵；
基于各时刻的工作变量矩阵确定所述扬声器系统在各频率信号下的非线性特性值。
可选的，基于所述等效电路模型构建与所述扬声器系统的工作变量相关的电力学方程组按照如下公式表示：
可选的，所述扬声器系统包括如下至少一个单元：扬声器单元、闭箱、倒相箱及PR箱；非线性特性值确定模块240，还用于：
当计算扬声器单元对应的工作变量矩阵时，令Kbox＝0，1/Mtube＝0，Rtube＝0，1/Mpr＝0，Rpr＝0，Kpr＝0；
当计算闭箱对应的工作变量矩阵时，令1/Mtube＝0，Rtube＝0，1/Mpr＝0，Rpr＝0，Kpr＝0；
当计算倒相箱对应的工作变量矩阵时，令1/Mpr＝0，Rpr＝0，Kpr＝0，且计算Kbox时的背腔体积包含倒相管内空气体积；
当计算PR箱对应的工作变量矩阵时，令1/Mtube＝0，Rtube＝0；
相应的，基于各时刻的工作变量矩阵确定所述扬声器系统在各频率信号下的非线性特性值，包括：
基于各时刻的工作变量矩阵分别确定扬声器单元、闭箱、倒相箱及PR箱在各频率信号下的非线性特性值。
可选的，获取所述工作变量矩阵的时间变化规律按照如下公式表示：Y(n+1)＝F*Y(n)+G*ue，其中，F和G分别为系数矩阵；Y(n+1)下一时刻的工作变量矩阵，Y(n)为当前时刻的工作变量矩阵，矩阵Y为1列7行的矩阵，且Y＝[i x v x1 v1 x2 v2]′，F和G表示如下：
其中，
其中，Ts为当前时刻与下一时刻的时间间隔。
可选的，非线性特性包括：声压级SPL、阻抗IMP、扬声器位移DIS_DRV、总谐波失真THD、二阶谐波失真2HD、3阶谐波失真3HD、无缘辐射器位移DISPR和倒相管中空气流速vair；非线性特性值确定模块240，还用于；
截取各时刻的工作变量矩阵中处于中间稳态部分的数据；
基于中间稳态部分的数据按照如下公式计算所述扬声器系统在各频率信号下的非线性特性值：
其中，t为时域求解中的时长；r为声压级测试中麦克风到扬声器的垂直距离；N为求解的最大谐波失真阶数；pnf为基频为f时的第n阶谐波声压，由时域声压经过傅里叶变换之后得到，Stube为导向管横截面积，Sd为扬声器有效振动面积，Sp为无缘辐射器(PR)的有效振动面积。
上述装置可执行本发明前述所有实施例所提供的方法，具备执行上述方法相应的功能模块和有益效果。未在本实施例中详尽描述的技术细节，可参见本发明前述所有实施例所提供的方法。
注意，上述仅为本发明的较佳实施例及所运用技术原理。本领域技术人员会理解，本发明不限于这里所述的特定实施例，对本领域技术人员来说能够进行各种明显的变化、重新调整和替代而不会脱离本发明的保护范围。因此，虽然通过以上实施例对本发明进行了较为详细的说明，但是本发明不仅仅限于以上实施例，在不脱离本发明构思的情况下，还可以包括更多其他等效实施例，而本发明的范围由所附的权利要求范围决定。

Claims (7)

1.一种扬声器系统的非线性特性的确定方法，其特征在于，包括：
对初始信号在频域进行预处理，并将预处理后的多个频率信号分别加载于扬声器系统的输入端；
对于加载了各频率信号的扬声器系统，获取所述扬声器系统中的声学元件的结构属性参数；其中，所述结构属性参数包括如下至少一项：扬声器数量、扬声器的串并联属性、蠕变效应参数、热效应参数、倒相管数量及无源辐射器PR数量；
根据所述结构属性参数构建所述扬声器系统的综合等效电路模型；
基于所述等效电路模型确定所述扬声器系统在各频率信号下的非线性特性值；
其中，根据所述结构属性参数构建所述扬声器系统的综合等效电路模型，包括：
根据所述结构属性参数确定所述扬声器系统对应的等效电路元件以及各等效电路元件的连接关系；
基于所述连接关系连接各所述等效电路元件，获得综合等效电路模型；
其中，所述综合等效电路模型包括：加载于扬声器系统的频率信号对应的输入电压ue、扬声器单元直阻Re、音圈电感Le(x)、力电耦合因子Bl(x)、扬声器有效振动质量对应的电感Mms、扬声器悬置系统等效力阻对应的电阻Rms(v)、扬声器悬置系统等效劲度系数对应的电容Kms(x)、背腔等效力学劲度系数对应的电容Kbox、倒相管中空气力学等效质量对应的电感Mtube、倒相管中空气等效力阻对应的电阻Rtube、PR有效振动质量经过次级到初级的阻抗换算后对应的电感Mpr、PR等效力阻经过次级到初级的阻抗换算后对应的电阻Rpr及PR等效劲度系数经过次级到初级的阻抗换算后对应的电容Kpr；
其中，基于所述等效电路模型确定所述扬声器系统在各频率信号下的非线性特性值，包括：
基于所述等效电路模型构建与所述扬声器系统的工作变量相关的电学力学方程组；
对所述电学力学方程组进行求解，获得当前时刻的工作变量矩阵；
其中，所述工作变量包括音圈电流i、扬声器振动系位移x、扬声器振动系振动速度v、倒相管中空气柱位移经过次级到初级的阻抗换算后的值x1、倒相管中空气柱振动速度经过次级到初级的阻抗换算后的值v1、PR位移经过次级到初级的阻抗换算后的值x2、PR振动速度经过次级到初级的阻抗换算后的值v2；
获取所述工作变量矩阵的时间变化规律；
基于所述时间变化规律和当前时刻的工作变量矩阵确定各时刻的工作变量矩阵；
基于各时刻的工作变量矩阵确定所述扬声器系统在各频率信号下的非线性特性值。
2.根据权利要求1所述的方法，其特征在于，对初始信号在频域进行预处理，包括：
对所述初始信号在频域进行滤波和/或均衡处理。
3.根据权利要求1所述的方法，其特征在于，基于所述等效电路模型构建与所述扬声器系统的工作变量相关的电力学方程组按照如下公式表示：
4.根据权利要求3所述的方法，其特征在于，所述扬声器系统包括如下至少一个单元：扬声器单元、闭箱、倒相箱及PR箱；对所述电力学方程组进行求解，获得当前时刻的工作变量矩阵，包括：
当计算扬声器单元对应的工作变量矩阵时，令Kbox＝0，1/Mtube＝0，Rtube＝0，1/Mpr＝0，Rpr＝0，Kpr＝0；
当计算闭箱对应的工作变量矩阵时，令1/Mtube＝0，Rtube＝0，1/Mpr＝0，Rpr＝0，Kpr＝0；
当计算倒相箱对应的工作变量矩阵时，令1/Mpr＝0，Rpr＝0，Kpr＝0，且计算Kbox时的背腔体积包含倒相管内空气体积；
当计算PR箱对应的工作变量矩阵时，令1/Mtube＝0，Rtube＝0；
相应的，基于各时刻的工作变量矩阵确定所述扬声器系统在各频率信号下的非线性特性值，包括：
基于各时刻的工作变量矩阵分别确定扬声器单元、闭箱、倒相箱及PR箱在各频率信号下的非线性特性值。
5.根据权利要求1所述的方法，其特征在于，获取所述工作变量矩阵的时间变化规律按照如下公式表示：Y(n+1)＝F*Y(n)+G*ue，其中，F和G分别为系数矩阵；Y(n+1)为下一时刻的工作变量矩阵，Y(n)为当前时刻的工作变量矩阵，矩阵Y为1列7行的矩阵，且Y＝[i x v x1 v1x2 v2]′，F和G表示如下：
其中，
其中，Ts为当前时刻与下一时刻的时间间隔。
6.根据权利要求5所述的方法，其特征在于，非线性特性包括：声压级SPL、阻抗IMP、扬声器位移DIS_DRV、总谐波失真THD、二阶谐波失真2HD、3阶谐波失真3HD、无缘辐射器位移DISPR和倒相管中空气流速vair；基于各时刻的工作变量矩阵确定所述扬声器系统在各频率信号下的非线性特性值；
截取各时刻的工作变量矩阵中处于中间稳态部分的数据；
基于中间稳态部分的数据按照如下公式计算所述扬声器系统在各频率信号下的非线性特性值：
其中，t为时域求解中的时长；r为声压级测试中麦克风到扬声器的垂直距离；N为求解的最大谐波失真阶数；pnf为基频为f时的第n阶谐波声压，由时域声压经过傅里叶变换之后得到，Stube为导向管横截面积，Sd为扬声器有效振动面积，Sp为无缘辐射器(PR)的有效振动面积，ρ为空气密度。
7.一种扬声器系统的非线性特性的确定装置，其特征在于，包括：
信号预处理模块，用于对初始信号在频域进行预处理，并将预处理后的多个频率信号分别加载于扬声器系统的输入端；
结构属性参数获取模块，用于对于加载了各频率信号的扬声器系统，获取所述扬声器系统中的声学元件的结构属性参数；其中，所述结构属性参数包括如下至少一项：扬声器数量、扬声器的串并联属性、蠕变效应参数、热效应参数、倒相管数量及无源辐射器PR数量；
综合等效电路模型构建模块，用于根据所述结构属性参数构建所述扬声器系统的综合等效电路模型；
非线性特性值确定模块，用于基于所述等效电路模型确定所述扬声器系统在各频率信号下的非线性特性值；
其中，所述综合等效电路模型构建模块，还用于：
根据所述结构属性参数确定所述扬声器系统对应的等效电路元件以及各等效电路元件的连接关系；
基于所述连接关系连接各所述等效电路元件，获得综合等效电路模型；
其中，所述综合等效电路模型包括：加载于扬声器系统的频率信号对应的输入电压ue、扬声器单元直阻Re、音圈电感Le(x)、力电耦合因子Bl(x)、扬声器有效振动质量对应的电感Mms、扬声器悬置系统等效力阻对应的电阻Rms(v)、扬声器悬置系统等效劲度系数对应的电容Kms(x)、背腔等效力学劲度系数对应的电容Kbox、倒相管中空气力学等效质量对应的电感Mtube、倒相管中空气等效力阻对应的电阻Rtube、PR有效振动质量经过次级到初级的阻抗换算后对应的电感Mpr、PR等效力阻经过次级到初级的阻抗换算后对应的电阻Rpr及PR等效劲度系数经过次级到初级的阻抗换算后对应的电容Kpr；
其中，所述非线性特性值确定模块，还用于：
基于所述等效电路模型构建与所述扬声器系统的工作变量相关的电学力学方程组；
对所述电学力学方程组进行求解，获得当前时刻的工作变量矩阵；
其中，所述工作变量包括音圈电流i、扬声器振动系位移x、扬声器振动系振动速度v、倒相管中空气柱位移经过次级到初级的阻抗换算后的值x1、倒相管中空气柱振动速度经过次级到初级的阻抗换算后的值v1、PR位移经过次级到初级的阻抗换算后的值x2、PR振动速度经过次级到初级的阻抗换算后的值v2；
获取所述工作变量矩阵的时间变化规律；
基于所述时间变化规律和当前时刻的工作变量矩阵确定各时刻的工作变量矩阵；
基于各时刻的工作变量矩阵确定所述扬声器系统在各频率信号下的非线性特性值。
CN202111037509.XA 2021-09-06 2021-09-06 一种扬声器系统的非线性特性的确定方法及装置 Active CN113630690B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111037509.XA CN113630690B (zh) 2021-09-06 2021-09-06 一种扬声器系统的非线性特性的确定方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111037509.XA CN113630690B (zh) 2021-09-06 2021-09-06 一种扬声器系统的非线性特性的确定方法及装置

Publications (2)

Publication Number Publication Date
CN113630690A CN113630690A (zh) 2021-11-09
CN113630690B true CN113630690B (zh) 2024-06-28

Family

ID=78389162

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111037509.XA Active CN113630690B (zh) 2021-09-06 2021-09-06 一种扬声器系统的非线性特性的确定方法及装置

Country Status (1)

Country Link
CN (1) CN113630690B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116055951B (zh) * 2022-07-20 2023-10-20 荣耀终端有限公司 信号处理方法和电子设备
CN116709123A (zh) * 2023-05-31 2023-09-05 歌尔股份有限公司 音频信号处理方法、装置、设备及存储介质

Family Cites Families (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20160134982A1 (en) * 2014-11-12 2016-05-12 Harman International Industries, Inc. System and method for estimating the displacement of a speaker cone
CN105246008A (zh) * 2015-10-26 2016-01-13 上海创米科技有限公司 共振型扬声器
CN105530586B (zh) * 2015-12-25 2018-11-02 矽力杰半导体技术(杭州)有限公司 扬声器振膜的保护方法和扬声器控制装置
GB2556015B (en) * 2016-04-29 2018-10-17 Cirrus Logic Int Semiconductor Ltd Audio Signals
CN106068007B (zh) * 2016-06-07 2019-09-13 瑞声科技(新加坡)有限公司 扬声器非线性系统辨识方法
WO2017220816A1 (en) * 2016-06-24 2017-12-28 Harman Becker Automotive Systems Gmbh Method for simulating total harmonic distortion of a loudspeaker
JPWO2019069466A1 (ja) * 2017-10-06 2020-09-10 株式会社 Ｔｒｉｇｅｎｃｅ Ｓｅｍｉｃｏｎｄｕｃｔｏｒ スピーカ駆動装置、スピーカ装置およびプログラム
US10542361B1 (en) * 2018-08-07 2020-01-21 Samsung Electronics Co., Ltd. Nonlinear control of loudspeaker systems with current source amplifier
CN110542793B (zh) * 2019-08-07 2021-06-08 华南理工大学 一种无源器件智能等效电路模型与参数测量装置及工作方法
CN112637734A (zh) * 2019-09-24 2021-04-09 陈文剑 一种扬声器系统控制方法及电子设备
CN113286233B (zh) * 2021-05-20 2022-12-27 维沃移动通信有限公司 扬声器仿真方法、仿真装置、扬声器以及电子设备

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
四阶带通式无源辐射器扬声器系统;周东宪等;《电子技术》;第114-115页 *
微型扬声器蠕变效应与小信号参数研究;刘成等;《声学技术》;第28卷(第6期);第91-93页 *
电动扬声器分数阶建模及非线性失真分析;孔晓鹏;《中国博士学位论文电子期刊网信息科技辑》(第11期);第70-102页 *

Also Published As

Publication number Publication date
CN113630690A (zh) 2021-11-09

Similar Documents

Publication Publication Date Title
CN104980877B (zh) 识别和补偿机电换能器中的非线性振动的装置和方法
CN106068007B (zh) 扬声器非线性系统辨识方法
CN113630690B (zh) 一种扬声器系统的非线性特性的确定方法及装置
Olson Electronic control of noise, vibration, and reverberation
Rivet et al. Broadband low-frequency electroacoustic absorbers through hybrid sensor-/shunt-based impedance control
CN102638736B (zh) 扬声器输出的控制
JP4243021B2 (ja) 電気音響スピーカ用コンデンサの無いクロスオーバ・ネットワーク
Small Direct-radiator loudspeaker system analysis
CN101231847B (zh) 一种结构振动噪声主动控制方法及其装置
CN102742300A (zh) 扬声器输出的控制
CN103686530A (zh) 通过对电动扬声器的总响应建模来处理音频信号的方法
JPS61294992A (ja) 電気信号を音響信号に変換する装置
CN110188428B (zh) 一种基于有限元法的扬声器联合仿真方法
CN102037738A (zh) 声音输入装置及其制造方法、以及信息处理系统
CN106973348A (zh) 扬声器模组及电子设备
CN101247671A (zh) 变换器参数的最优估计
CN110248300B (zh) 一种基于自主学习的啸叫抑制方法及扩声系统
CN103916810B (zh) 一种时域声能量对比控制方法和系统
Huang et al. Analysis of total harmonic distortion of miniature loudspeakers used in mobile phones considering nonlinear acoustic damping
CN207744107U (zh) 耳机分频电路
Cao et al. Broadband sound absorption of a hybrid absorber with shunt loudspeaker and perforated plates
Bai et al. Optimization of microspeaker diaphragm pattern using combined finite element–lumped parameter models
TW202026635A (zh) 多孔材料的預測方法及其系統
CN2744105Y (zh) 一种低成本且高效率的抗噪耳机及配套装置
Bai et al. Robust control of a sensorless bass-enhanced moving-coil loudspeaker system

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant