非线性音频效果的数字建模方法、计算机设备及存储介质

1.本发明属于数字音频技术领域，特别涉及一种非线性音频效果的数字建模方法、计算机设备及存储介质。


背景技术：

2.非线性效果的数字建模一直是数字音频效果领域深度研究的内容。压缩(compressor)、激励(exciter)、过载(overdrive)和失真(distortion)等各种音频效果都属于非线性效果的范畴。但是，能够实现这些效果的模拟设备通常较为昂贵，且体积庞大、难以携带。因此，研究非线性效果的数字建模方法是有必要的。
3.非线性效果的数字建模主要有两种方式，即白盒建模与黑盒建模。白盒建模也被称为正向建模，其利用了目标系统的所有已知信息，包括目标系统的电路图以及电子元器件的具体参数以及非线性特性等。白盒建模的优势在于其能够提供非常高精度的建模结果，而其劣势在于必须要求获取目标系统的详细特征，并且，模拟电路复杂的拓扑结构与大量非线性元件可能造成实时处理的计算成本过高。黑盒建模也被称为反向建模，其利用抽象的非线性系统模型来表示目标系统，再通过数值计算方法计算建模结果。这类抽象非线性模型主要包括volterra模型、chebyshev模型与人工神经网络模型等。黑盒建模的优势在于其不需要获取目标系统的详细特征，只要目标系统符合抽象非线性模型的特征，就可以对其进行建模，而其劣势在于能够提供的建模精度相对白盒建模来说较低。
4.对于以非线性效果数字建模为核心业务的厂商来说，建模精度更高的白盒建模无疑是首选的建模方式，而对于普通用户来说，便携性和自由度则是更值得考虑的因素。因此，数字音频效果器厂商通常选择白盒建模作为产品的核心技术，而将黑盒建模作为产品的特色功能之一，供用户自行选择。现有技术中，黑盒建模的非线性模型通常较为复杂，模型参数的计算成本较高，所需的计算时间长，这导致配置有黑盒建模技术的产品价格普遍较高，且用户在对目标系统建模时等待时间较长，体验感较差。因此，亟需一种复杂度较低，且计算成本不高的建模方法。


技术实现要素：

5.本发明的主要目的在于提供一种非线性音频效果的数字建模方法，以解决上述现有技术中，黑盒建模的非线性模型通常较为复杂，模型参数的计算成本较高，所需的计算时间长，导致配置有黑盒建模技术的产品价格普遍较高，且用户在对目标音频效果建模时等待时间较长，体验感较差的问题。
6.为实现上述目的，本发明采用了如下技术方案：
7.第一方面，本发明提供一种非线性音频效果的数字建模方法，所述方法包括：
8.获取音频源输入信号的采样样本，所述音频源包括任何具有宽频带特征的音频信号；
9.将输入信号的采样样本发送至目标系统的输入，并获取目标系统输出信号的采样
样本；
10.基于w-h非线性模型，以输入信号的采样样本作为w-h非线性模型的输入，以目标系统输出信号的采样样本作为w-h非线性模型的期望输出，使用自适应矩估计计算w-h非线性模型的参数；
11.获取w-h非线性模型的期望输出信号的采样样本和真实输出信号的采样样本，对每个采样样本进行短时傅里叶变换，得到两个信号的幅度谱，计算两个幅度谱之间的差异，并采用级联的滤波器拟合该差异，得到最终的w-h非线性模型；
12.基于最终的w-h非线性模型输出处理后的信号。
13.进一步的，所述w-h非线性模型由滤波器单元与非线性单元级联构成。
14.进一步的，计算w-h非线性模型的参数具体包括：
15.首先按下式计算fir滤波器的输出：
16.y[n]＝x[n]*h[n]
[0017]
其中，x[n]表示输入信号；y[n]表示输出信号；h[n]表示滤波器传递函数；
[0018]
按下式计算非线性单元的输出：
[0019][0020]
其中，k
p
、kn、g
p
、gn参数用于调节非线性曲线的响应特征，该非线性单元用于模拟电子元器件的非线性特性；
[0021]
然后，计算目标系统的输出信号y
real
[n]与模型的输出信号y
model
[n]在时域上的均方误差j，如下式所示：
[0022][0023]
其中，y
model
[n]表示模型的输出信号；y
real
[n]表示目标系统的输出信号；n为信号样本采样点数；
[0024]
之后，计算在自适应矩估计中需要用到的偏导数，包括之后，计算在自适应矩估计中需要用到的偏导数，包括
[0025]
进一步的，还包括使用自适应矩估计对w-h模型的参数进行更新，计算步骤如下：
[0026]
(1)计算移动平均梯度v：
[0027][0028]
(2)计算误差修正后的移动平均梯度v
corrected
：
[0029][0030]
(3)计算移动平均平方梯度s：
[0031][0032]
(4)计算误差修正后的移动平均平方梯度s
corrected
：
[0033][0034]
(5)对参数θ进行更新：
[0035][0036]
其中，θ表示w-h模型中的任意参数；t表示模型当前的迭代次数；β1、β2用于确定移动平均数计算的窗口大小；ε用于防止分母项过小而导致梯度爆炸。
[0037]
进一步的，还包括在参数更新时使用自适应学习率衰减方法，具体包括：
[0038]
(1)计算当前迭代次数下均方误差j曲线的斜率大小；
[0039]
(2)设置阈值，判断当前迭代次数下斜率是否小于该阈值；若是，则认为模型需要衰减学习率，转至步骤(3)，若否，则学习率不变；
[0040]
(3)按下式衰减学习率：
[0041]
α＝α
×r[0042]
其中，α表示学习率，r表示衰减率。
[0043]
进一步的，基于最终的w-h非线性模型输出处理后的信号还包括：
[0044]
用户介入试听，向用户提供输入信号的增益调节，并基于用户指令对输入信号的增益量进行调节。
[0045]
第二方面，本发明提供一种计算机设备，包括存储器、处理器，所述存储器中存储有可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述第一方面所述的方法。
[0046]
第三方面，本发明提供一种计算机可读存储介质，所述计算机可读存储介质中存储有计算机执行指令，所述计算机执行指令被处理器执行时用于实现上述第一方面所述的方法。
[0047]
本发明的有益技术效果如下：
[0048]
本发明提供的非线性音频效果的数字建模方法、计算机设备及存储介质，通过获取输入信号，将其发送至目标系统，再获取目标系统的输出信号，基于w-h模型，结合自适应矩估计对目标系统进行数字建模，同时还引入了增益调节与自适应学习率衰减，加速了模型的收敛过程。解决了现有技术中，黑盒建模的非线性模型通常较为复杂，模型参数的计算成本较高，所需的计算时间长，导致配置有黑盒建模技术的产品价格普遍较高，且用户在对目标音频效果建模时等待时间较长，体验感较差的问题。
附图说明
[0049]
此处的附图被并入说明书中并构成本说明书的一部分，示出了符合本公开的实施例，并与说明书一起用于解释本公开的原理。
[0050]
图1是本发明实施例示出的非线性音频效果的数字建模方法的流程示意图。
[0051]
图2是以peavey 5150型号电子管音箱作为建模对象，用户将配置有本发明方法的装置与其他设备连接的方式示意图。
[0052]
图3是同一段电吉他信号分别由peavey 5150音箱输出的真实信号以及本发明方法建模得到的模型输出波形图。
[0053]
图4是本发明实施例提供的计算机设备的组成示意图。
[0054]
通过上述附图，已示出本公开明确的实施例，后文中将有更详细的描述。这些附图和文字描述并不是为了通过任何方式限制本公开构思的范围，而是通过参考特定实施例为本领域技术人员说明本公开的概念。
具体实施方式
[0055]
这里将详细地对示例性实施例进行说明，其示例表示在附图中。以下示例性实施例中所描述的实施方式并不代表与本发明相一致的所有实施方式。相反，它们仅是与权利要求书中所详述的、本技术的一些方面相一致的装置和方法的例子。
[0056]
现有技术中，黑盒建模的非线性模型通常较为复杂，模型参数的计算成本较高，所需的计算时间长，这导致配置有黑盒建模技术的产品价格普遍较高，且用户在对目标系统建模时等待时间较长，体验感较差。因此，亟需一种复杂度较低，且计算成本不高的建模方法。
[0057]
针对上述问题，本发明实施例提供一种非线性音频效果的数字建模方法、计算机设备及存储介质，以解决上述现有技术中的问题。
[0058]
下面以具体地实施例对本发明的技术方案以及本发明的技术方案如何解决上述技术问题进行详细说明。下面将结合附图，对本发明的实施例进行描述。
[0059]
图1为本发明实施例提供的一种非线性音频效果的数字建模方法的流程示意图，如图1所示，该方法包括：
[0060]
s101、获取音频源输入信号的采样样本。
[0061]
音频源可以是任何具有宽频带特征的音频信号，例如吉他的演奏信号、歌手的演唱信号以及正弦扫频信号等等。
[0062]
s102、将输入信号发送至目标系统的输入，并获取目标系统输出信号的采样样本；
[0063]
s103、基于w-h非线性模型，以输入信号的采样样本作为w-h非线性模型的输入，以目标系统输出信号的采样样本作为w-h非线性模型的期望输出，使用自适应矩估计计算w-h非线性模型的参数；
[0064]
本实施例中，w-h非线性模型由fir滤波器单元与非线性单元级联构成，在计算模型参数的过程中，首先按下式计算fir滤波器的输出：
[0065]
y[n]＝x[n]*h[n]
[0066]
其中，x[n]表示输入信号；y[n]表示输出信号；h[n]表示滤波器传递函数。
[0067]
按下式计算非线性单元的输出：
[0068][0069]
其中，k
p
、kn、g
p
、gn参数用于调节非线性曲线的响应特征，该非线性单元用于模拟
电子元器件的非线性特性。
[0070]
然后，计算目标系统的输出信号y
real
[n]与模型的输出信号y
model
[n]在时域上的均方误差j，如下式所示：
[0071][0072]
其中，y
model
[n]表示模型的输出信号；y
real
[n]表示目标系统的输出信号；n为信号样本采样点数。
[0073]
之后，计算在自适应矩估计中需要用到的偏导数，包括之后，计算在自适应矩估计中需要用到的偏导数，包括
[0074]
使用自适应矩估计对模型参数进行更新，计算步骤如下：
[0075]
(1)计算移动平均梯度v：
[0076][0077]
(2)计算误差修正后的移动平均梯度v
corrected
：
[0078][0079]
(3)计算移动平均平方梯度s：
[0080][0081]
(4)计算误差修正后的移动平均平方梯度s
corrected
：
[0082][0083]
(5)对参数θ进行更新：
[0084][0085]
其中，θ表示w-h非线性模型中的任意参数；t表示模型当前的迭代次数；β1、β2用于确定移动平均数计算的窗口大小；ε用于防止分母项过小而导致梯度爆炸。
[0086]
为了加速模型的收敛，减少计算时间，在参数更新时还需要衰减学习率。此处应用了一种自适应学习率衰减方法，其步骤如下：
[0087]
(1)计算当前迭代次数下均方误差j曲线的斜率大小；
[0088]
(2)设置阈值，判断当前迭代次数下斜率是否小于该阈值。若是，则认为模型需要衰减学习率，转至步骤(3)，若否，则学习率不变；
[0089]
(3)按下式衰减学习率：
[0090]
α＝α
×r[0091]
其中，α表示学习率，r表示衰减率。
[0092]
s104、获取w-h非线性模型的期望输出信号与真实输出信号的采样样本，对每个采样样本进行短时傅里叶变换，得到两个信号的幅度谱，计算两个幅度谱之间的差异，并采用级联的滤波器拟合该差异，得到最终的w-h非线性模型。
[0093]
优选的，本实施例中，滤波器采用iir滤波器。
[0094]
s105、基于最终的w-h非线性模型输出处理后的信号。
[0095]
此时，允许用户介入试听，并向用户提供输入信号的增益调节，由用户确定输入信号的增益量。
[0096]
示例性地，以peavey 5150型号电子管音箱作为建模对象，用户以图2所示的方式将配置有本方法的装置与其他设备连接，并输入一段若干秒时长的电吉他信号，该信号将被发送至peavey 5150电子管音箱中，同时音箱将返回其输出信号。若干秒后，建模完成，提示用户介入进行增益调节。最后，由用户确认并储存建模数据。图3显示了该示例中，同一段电吉他信号分别由peavey 5150音箱输出的真实信号以及本方法建模得到的模型输出。可以看到，由本方法建模得到的模型基本还原了真实的peavey 5150音箱的特征。
[0097]
图4为本发明实施例提供的一种计算机设备的结构示意图，如图4所示，计算机设备500包括：存储器501和处理器502；
[0098]
其中，存储器501用于存储处理器可执行的计算机指令；
[0099]
处理器502在执行计算机指令时实现上述实施例中方法中的各个步骤。具体可以参见前述方法实施例中的相关描述。
[0100]
可选地，上述存储器501既可以是独立的，也可以跟处理器502集成在一起。当存储器501独立设置时，该检测设备还包括总线，用于连接存储器501和处理器502。
[0101]
本技术实施例还提供了一种非临时性计算机可读存储介质，当该存储介质中的指令由电子设备的处理器执行时，使得电子设备能够执行上述实施例提供的方法。
[0102]
本技术实施例还提供了一种计算机程序产品，计算机程序产品包括：计算机程序，计算机程序存储在可读存储介质中，电子设备的至少一个处理器可以从可读存储介质读取计算机程序，至少一个处理器执行计算机程序使得电子设备执行上述任一实施例提供的方案。
[0103]
以上对本技术所述的一种非线性音频效果的数字建模方法、装置、设备及存储介质进行了详细地的介绍，以上的实例说明只是用于帮助理解本技术的方法及其核心思想而非对其进行限制，其他的任何未背离本技术的精神实质与原理下所作改变、修饰、替代、组合、简化，均应为等效的置换方式，都包含在本技术的保护范围之内。

技术特征：
1.一种非线性音频效果的数字建模方法，其特征在于，所述方法包括：获取音频源输入信号的采样样本，所述音频源包括任何具有宽频带特征的音频信号；将输入信号的采样样本发送至目标系统的输入，并获取目标系统输出信号的采样样本；基于w-h非线性模型，以输入信号的采样样本作为w-h非线性模型的输入，以目标系统输出信号的采样样本作为w-h非线性模型的期望输出，使用自适应矩估计计算w-h非线性模型的参数；获取w-h非线性模型的期望输出信号的采样样本和真实输出信号的采样样本，对每个采样样本进行短时傅里叶变换，得到两个信号的幅度谱，计算两个幅度谱之间的差异，并采用级联的滤波器拟合该差异，得到最终的w-h非线性模型；基于最终的w-h非线性模型输出处理后的信号。2.如权利要求1所述的非线性音频效果的数字建模方法，其特征在于，所述w-h非线性模型由滤波器单元与非线性单元级联构成。3.如权利要求2所述的非线性音频效果的数字建模方法，其特征在于，计算w-h非线性模型的参数具体包括：首先按下式计算fir滤波器的输出：y[n]＝x[n]*h[n]其中，x[n]表示输入信号；y[n]表示输出信号；h[n]表示滤波器传递函数；按下式计算非线性单元的输出：其中，k
p
、k
n
、g
p
、g
n
参数用于调节非线性曲线的响应特征，该非线性单元用于模拟电子元器件的非线性特性；然后，计算目标系统的输出信号y
real
[n]与模型的输出信号y
model
[n]在时域上的均方误差j，如下式所示：其中，y
model
[n]表示模型的输出信号；y
real
[n]表示目标系统的输出信号；n为信号样本采样点数；之后，计算在自适应矩估计中需要用到的偏导数，包括之后，计算在自适应矩估计中需要用到的偏导数，包括4.如权利要求3所述的非线性音频效果的数字建模方法，其特征在于，还包括使用自适应矩估计对w-h模型的参数进行更新，计算步骤如下：(1)计算移动平均梯度v：
(2)计算误差修正后的移动平均梯度v
corrected
：(3)计算移动平均平方梯度s：(4)计算误差修正后的移动平均平方梯度s
corrected
：(5)对参数θ进行更新：其中，θ表示w-h模型中的任意参数；t表示模型当前的迭代次数；β1、β2用于确定移动平均数计算的窗口大小；ε用于防止分母项过小而导致梯度爆炸。5.如权利要求4所述的非线性音频效果的数字建模方法，其特征在于，还包括在参数更新时使用自适应学习率衰减方法，具体包括：(1)计算当前迭代次数下均方误差j曲线的斜率大小；(2)设置阈值，判断当前迭代次数下斜率是否小于该阈值；若是，则认为模型需要衰减学习率，转至步骤(3)，若否，则学习率不变；(3)按下式衰减学习率：α＝α
×
r其中，α表示学习率，r表示衰减率。6.如权利要求1-5任一项所述的非线性音频效果的数字建模方法，其特征在于，基于最终的w-h非线性模型输出处理后的信号还包括：用户介入试听，向用户提供输入信号的增益调节，并基于用户指令对输入信号的增益量进行调节。7.一种计算机设备，其特征在于，包括存储器、处理器，所述存储器中存储有可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述权利要求1-6任一项中的方法。8.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质中存储有计算机执行指令，所述计算机执行指令被处理器执行时用于实现如权利要求1-6任一项中的方法。

技术总结
本申请提供了一种非线性音频效果的数字建模方法、计算机设备及存储介质，该方法包括：获取音频源输入信号的采样样本；将输入信号的采样样本发送至目标系统的输入，并获取目标系统输出信号的采样样本；基于W-H非线性模型，以输入信号的采样样本作为W-H非线性模型的输入，以目标系统输出信号的采样样本作为W-H非线性模型的期望输出，使用自适应矩估计计算W-H非线性模型的参数；获取W-H非线性模型的期望输出信号的采样样本和真实输出信号的采样样本，对每个采样样本进行短时傅里叶变换，得到两个信号的幅度谱，计算两个幅度谱之间的差异，并采用级联的滤波器拟合该差异，得到最终的W-H非线性模型；基于最终的W-H非线性模型输出处理后的信号。出处理后的信号。出处理后的信号。


技术研发人员：张宏鑫 毛怡帆 吴泓嘉
受保护的技术使用者：浙江大学
技术研发日：2023.05.30
技术公布日：2023/8/21