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

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

Info

Publication number
CN113630690A
CN113630690A CN202111037509.XA CN202111037509A CN113630690A CN 113630690 A CN113630690 A CN 113630690A CN 202111037509 A CN202111037509 A CN 202111037509A CN 113630690 A CN113630690 A CN 113630690A
Authority
CN
China
Prior art keywords
loudspeaker
loudspeaker system
determining
equivalent circuit
tube
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.)
Granted
Application number
CN202111037509.XA
Other languages
English (en)
Other versions
CN113630690B (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
Priority claimed from CN202111037509.XA external-priority 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

Images

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

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。
具体的，基于等效电路模型构建与扬声器系统的工作变量相关的电力学方程组按照如下公式表示：
Figure BDA0003247857250000061
本实施例中，扬声器系统包括如下至少一个单元：扬声器单元、闭箱、倒相箱及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表示如下：
Figure BDA0003247857250000071
其中，
Figure BDA0003247857250000072
Figure BDA0003247857250000081
Figure BDA0003247857250000082
Figure BDA0003247857250000083
Figure BDA0003247857250000084
Figure BDA0003247857250000085
Figure BDA0003247857250000086
Figure BDA0003247857250000087
Figure BDA0003247857250000088
Figure BDA0003247857250000089
Figure BDA00032478572500000810
Ts为当前时刻与下一时刻的时间间隔。
具体的，在获得给定频率信号下各时刻的工作变量矩阵后，截取各时刻的工作变量矩阵中处于中间稳态部分的数据，基于中间稳态部分的数据计算该给定频率信号下的声压级SPL、阻抗IMP、扬声器位移DIS_DRV、总谐波失真THD、二阶谐波失真2HD、3阶谐波失真3HD、无缘辐射器位移DISPR和倒相管中空气流速v_air。
其中，基于中间稳态部分的数据按照如下公式计算所述扬声器系统在各频率信号下的非线性特性值：
其中，t为时域求解中的时长；r为声压级测试中麦克风到扬声器的垂直距离；N为求解的最大谐波失真阶数；pnf为基频为f时的第n阶谐波声压，由时域声压经过傅里叶变换之后得到，Stube为导向管横截面积，ρ为空气密度，Sd为扬声器有效振动面积，Sp为无缘辐射器(PR)的有效振动面积。
Figure BDA0003247857250000091
以频率为循环参数，对每一个频率进行上述计算之后，便能得到扬声器系统声压级、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；
获取所述工作变量矩阵的时间变化规律；
基于所述时间变化规律和当前时刻的工作变量矩阵确定各时刻的工作变量矩阵；
基于各时刻的工作变量矩阵确定所述扬声器系统在各频率信号下的非线性特性值。
可选的，基于所述等效电路模型构建与所述扬声器系统的工作变量相关的电力学方程组按照如下公式表示：
Figure BDA0003247857250000121
可选的，所述扬声器系统包括如下至少一个单元：扬声器单元、闭箱、倒相箱及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 x1v1 x2 v2]′，F和G表示如下：
Figure BDA0003247857250000131
其中，
Figure BDA0003247857250000132
Figure BDA0003247857250000133
Figure BDA0003247857250000134
Figure BDA0003247857250000135
Figure BDA0003247857250000136
Figure BDA0003247857250000137
Figure BDA0003247857250000138
Figure BDA0003247857250000139
Figure BDA00032478572500001310
Figure BDA00032478572500001311
Figure BDA00032478572500001312
其中，Ts为当前时刻与下一时刻的时间间隔。
可选的，非线性特性包括：声压级SPL、阻抗IMP、扬声器位移DIS_DRV、总谐波失真THD、二阶谐波失真2HD、3阶谐波失真3HD、无缘辐射器位移DISPR和倒相管中空气流速v_air；非线性特性值确定模块240，还用于；
截取各时刻的工作变量矩阵中处于中间稳态部分的数据；
基于中间稳态部分的数据按照如下公式计算所述扬声器系统在各频率信号下的非线性特性值：
Figure BDA0003247857250000141
其中，t为时域求解中的时长；r为声压级测试中麦克风到扬声器的垂直距离；N为求解的最大谐波失真阶数；pnf为基频为f时的第n阶谐波声压，由时域声压经过傅里叶变换之后得到，Stube为导向管横截面积，Sd为扬声器有效振动面积，Sp为无缘辐射器(PR)的有效振动面积。
上述装置可执行本发明前述所有实施例所提供的方法，具备执行上述方法相应的功能模块和有益效果。未在本实施例中详尽描述的技术细节，可参见本发明前述所有实施例所提供的方法。
注意，上述仅为本发明的较佳实施例及所运用技术原理。本领域技术人员会理解，本发明不限于这里所述的特定实施例，对本领域技术人员来说能够进行各种明显的变化、重新调整和替代而不会脱离本发明的保护范围。因此，虽然通过以上实施例对本发明进行了较为详细的说明，但是本发明不仅仅限于以上实施例，在不脱离本发明构思的情况下，还可以包括更多其他等效实施例，而本发明的范围由所附的权利要求范围决定。

Claims (9)

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

Priority Applications (1)

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

Applications Claiming Priority (1)

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

Publications (2)

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

Family

ID=

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116055951A (zh) * 2022-07-20 2023-05-02 荣耀终端有限公司 信号处理方法和电子设备

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105246008A (zh) * 2015-10-26 2016-01-13 上海创米科技有限公司 共振型扬声器
CN105530586A (zh) * 2015-12-25 2016-04-27 矽力杰半导体技术(杭州)有限公司 扬声器振膜的保护方法和扬声器控制装置
US20160134982A1 (en) * 2014-11-12 2016-05-12 Harman International Industries, Inc. System and method for estimating the displacement of a speaker cone
CN106068007A (zh) * 2016-06-07 2016-11-02 瑞声科技(新加坡)有限公司 扬声器非线性系统辨识方法
US20170318390A1 (en) * 2016-04-29 2017-11-02 Cirrus Logic International Semiconductor Ltd. Audio signals
WO2017220816A1 (en) * 2016-06-24 2017-12-28 Harman Becker Automotive Systems Gmbh Method for simulating total harmonic distortion of a loudspeaker
CN109906617A (zh) * 2017-10-06 2019-06-18 株式会社特瑞君思半导体 扬声器驱动装置、扬声器装置和程序
CN110542793A (zh) * 2019-08-07 2019-12-06 华南理工大学 一种无源器件智能等效电路模型与参数测量装置及工作方法
US10542361B1 (en) * 2018-08-07 2020-01-21 Samsung Electronics Co., Ltd. Nonlinear control of loudspeaker systems with current source amplifier
CN112637734A (zh) * 2019-09-24 2021-04-09 陈文剑 一种扬声器系统控制方法及电子设备
CN113286233A (zh) * 2021-05-20 2021-08-20 维沃移动通信有限公司 扬声器仿真方法、仿真装置、扬声器以及电子设备

Patent Citations (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 上海创米科技有限公司 共振型扬声器
CN105530586A (zh) * 2015-12-25 2016-04-27 矽力杰半导体技术(杭州)有限公司 扬声器振膜的保护方法和扬声器控制装置
US20170318390A1 (en) * 2016-04-29 2017-11-02 Cirrus Logic International Semiconductor Ltd. Audio signals
CN106068007A (zh) * 2016-06-07 2016-11-02 瑞声科技(新加坡)有限公司 扬声器非线性系统辨识方法
WO2017220816A1 (en) * 2016-06-24 2017-12-28 Harman Becker Automotive Systems Gmbh Method for simulating total harmonic distortion of a loudspeaker
CN109906617A (zh) * 2017-10-06 2019-06-18 株式会社特瑞君思半导体 扬声器驱动装置、扬声器装置和程序
US10542361B1 (en) * 2018-08-07 2020-01-21 Samsung Electronics Co., Ltd. Nonlinear control of loudspeaker systems with current source amplifier
CN110542793A (zh) * 2019-08-07 2019-12-06 华南理工大学 一种无源器件智能等效电路模型与参数测量装置及工作方法
CN112637734A (zh) * 2019-09-24 2021-04-09 陈文剑 一种扬声器系统控制方法及电子设备
CN113286233A (zh) * 2021-05-20 2021-08-20 维沃移动通信有限公司 扬声器仿真方法、仿真装置、扬声器以及电子设备

Non-Patent Citations (7)

* Cited by examiner, † Cited by third party
Title
A. DUSTOR等: "Speaker Identification And Verification Based On Cepstral Features And Fuzzy Nonlinear Classif", 《PROCEEDINGS OF THE INTERNATIONAL CONFERENCE MIXED DESIGN OF INTEGRATED CIRCUITS AND SYSTEM》 *
刘成等: "微型扬声器蠕变效应与小信号参数研究", 《声学技术》, vol. 28, no. 6, pages 91 - 93 *
周东宪等: "四阶带通式无源辐射器扬声器系统", 《电子技术》, pages 114 - 115 *
孔晓鹏: "电动扬声器分数阶建模及非线性失真分析", 《中国博士学位论文电子期刊网》 *
孔晓鹏: "电动扬声器分数阶建模及非线性失真分析", 《中国博士学位论文电子期刊网信息科技辑》, no. 11, pages 70 - 102 *
薛政: "动圈扬声器非线性参数辨识", 《中国优秀硕士学位论文全文数据库》 *
黄坤朋等: "使用逆滤波器方法获得扬声器的目标传输函数", 《电声技术》, vol. 41, no. 1 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116055951A (zh) * 2022-07-20 2023-05-02 荣耀终端有限公司 信号处理方法和电子设备
CN116055951B (zh) * 2022-07-20 2023-10-20 荣耀终端有限公司 信号处理方法和电子设备

Similar Documents

Publication Publication Date Title
CN103634726B (zh) 一种扬声器自动均衡方法
CN104980877B (zh) 识别和补偿机电换能器中的非线性振动的装置和方法
CN106068007B (zh) 扬声器非线性系统辨识方法
CN1951148B (zh) 用于限制扬声器位移的系统
CN101231847B (zh) 一种结构振动噪声主动控制方法及其装置
Rivet et al. Broadband low-frequency electroacoustic absorbers through hybrid sensor-/shunt-based impedance control
US9210506B1 (en) FFT bin based signal limiting
CN102843633A (zh) 扬声器输出的控制
US8054983B2 (en) Method for parameter identification and parameter optimization of microspeakers
JP5364271B2 (ja) 変換器パラメータの最適推定装置および方法
CN103686530A (zh) 通过对电动扬声器的总响应建模来处理音频信号的方法
CN103916810B (zh) 一种时域声能量对比控制方法和系统
EP3225038B1 (en) Low frequency active acoustic absorber by acoustic velocity control through porous resistive layers
Zhou et al. A frequency band constrained filtered–x least mean square algorithm for feedback active control systems
CN113630690A (zh) 一种扬声器系统的非线性特性的确定方法及装置
CN113630690B (zh) 一种扬声器系统的非线性特性的确定方法及装置
CN112804626B (zh) 一种动态控制扬声器振幅的方法及系统、移动终端
Bai et al. Robust control of a sensorless bass-enhanced moving-coil loudspeaker system
TW202026635A (zh) 多孔材料的預測方法及其系統
Payal et al. Equalization of excursion and current-dependent nonlinearities in loudspeakers
CN113395639A (zh) 扬声器振膜振动位移的计算方法、扬声器保护装置及介质
CN106559722A (zh) 基于人耳听觉特性的音频重放系统均衡方法
Leishman et al. Evaluation of moving-coil loudspeaker and passive radiator parameters using normal-incidence sound transmission measurements: Theoretical developments
Schurer et al. Identification and compensation of the electrodynamic transducer nonlinearities
US20130243221A1 (en) Method and system of equalization pre-preocessing for sound receivng 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