一种数字信号处理方法、系统、终端及存储介质与流程

1.本技术涉及数字信号处理领域，特别是涉及一种数字信号处理方法、系统、终端及存储介质。


背景技术：

2.在实际数字信号采样中，由于采样的是离散数据，即等时间间隔获取对应的样本值，如采样频率为1000hz，即每隔1ms获取一次样本值，会有一部分有用样本值未被采集到，不利于对时域进行精密分析。在实际运用中有两种解决方法：一种处理方法是在采样时提高采样频率，这种方法得到的数据结果一般可以满足精密分析要求，但是数据量较大，会占用更大的储存空间，同时导致工作效率降低，不满足数字信号采集与分析的实际应用需求；另一种处理方法是用采集到的已知数字信号样本来估算未知所需数字信号样本的计算。而插值作为一种数值逼近的重要方法，可以通过有限个样本的取值情况估计出未知样本的近似值。目前常见的数字信号插值算法有拉格朗日插值、牛顿插值、厄米特插值等。这些插值方式均为线性插值，计算简单稳定，但是其插值误差较大，而且光滑性较差，会导致数字信号失真、峰值丢失、还引发虚假频率成分等问题。
3.因此，如何提高插值算法进行数字信号处理的插值精度，降低插值得到的数据和真实数据的误差，是本领域技术人员亟待解决的技术问题。


技术实现要素：

4.为解决上述技术问题，本技术提供一种数字信号处理方法，能够提高插值算法进行数字信号处理的插值精度，降低插值得到的数据和真实数据的误差。本技术还提供一种数字信号处理系统、终端及存储介质，具有相同的技术效果。
5.本技术的第一个目的为提供一种数字信号处理方法。
6.本技术的上述申请目的一是通过以下技术方案得以实现的：
7.一种数字信号处理方法，包括：
8.获取原始数字信号；
9.利用二阶导插值方法，运用函数二阶导数对所述原始数字信号进行插值处理，得到第一数据。
10.优选地，所述数字信号处理方法中，所述利用二阶导插值方法，运用函数二阶导数对所述原始数字信号进行插值处理，得到第一数据，包括：
11.在所述原始数字信号相邻两个数据节点组成的各区间，构建二阶导插值函数；
12.构建所述二阶导插值函数求解的矩阵方程组；
13.将所述原始数字信号的各数据节点代入所述矩阵方程组，求解得到二次微分值；
14.将所述二次微分值代入到所述二阶导插值函数，计算得出所述各区间的二阶导插值函数；
15.判断插值节点所在的区间，根据所述插值节点所在的区间的二阶导插值函数，计
算得到所述插值节点对应的函数值；
16.根据所述插值节点和所述插值节点对应的函数值，对所述原始数字信号进行插值处理，得到第一数据。
17.优选地，所述数字信号处理方法中，所述在所述原始数字信号相邻两个数据节点组成的各区间，构建二阶导插值函数，包括：
18.设所述原始数字信号包含n个数据节点，分别为(x1,y1),(x2,y2),
…
,(xn,yn)，n为大于1的正整数，得到所述原始数字信号相邻两个数据节点组成的各区间[xi,x
i+1
](i＝1,2,
…
,n-1)；
[0019]
构建所述各区间[xi,x
i+1
](i＝1,2,
…
,n-1)的二阶导插值函数si(x)＝ai+bi(x-xi)+ci(x-xi)2+di(x-xi)3(xi≤x≤x
i+1
)，式中，x表示插值函数自变量，ai,bi,ci,di为二阶导插值函数的各项系数，具体为：
[0020][0021]
式中，mi为二次微分值，步长hi＝x
i+1-xi(i＝1,2,
…
,n-1)。
[0022]
优选地，所述数字信号处理方法中，所述构建所述二阶导插值函数求解的矩阵方程组，包括：
[0023]
根据所述数据节点(x1,y1),(x2,y2),
…
,(xn,yn)构建矩阵方程式一：s(xi)＝yi(i＝1,2,
…
,n)，式中，y表示数据节点x对应的函数值，i表示数据节点的个数，n为大于1的正整数；
[0024]
根据所述s(x)、s
′
(x)、s
″
(x)在所述原始数字信号中具有二阶连续导数时，构建方矩阵方程式二：
[0025][0026]
式中，i表示数据节点的个数，n为大于1的正整数；
[0027]
根据所述s(x)在所述原始数字信号中为一个整周期数据时，构建矩阵方程式三：
[0028][0029]
式中，n为大于1的正整数。
[0030]
优选地，所述数字信号处理方法中，所述将所述数字信号的各数据节点代入所述矩阵方程组，求解得到二次微分值，包括：
[0031]
将所述数据节点代入所述矩阵方程式一、所述矩阵方程式二、所述矩阵方程式三，求解得到二次微分值mi＝si″
(xi)；
[0032]
相应地，所述将所述二次微分值代入到所述二阶导插值函数，计算得出所述各区间的二阶导插值函数，包括：
[0033]
根据所述二次微分值、所述数据节点、所述步长，求解得到ai,bi,ci,di，并根据所述ai,bi,ci,di，计算得出所述各区间的二阶导插值函数si(x)；
[0034]
相应地，所述判断插值节点所在的区间，根据所述插值节点所在的区间的二阶导插值函数，计算得到所述插值节点对应的函数值，包括：
[0035]
判断插值节点x
′k所在区间[xi,x
i+1
](i＝1,2,
…
,n-1)，根据所述插值节点x
′k所在的区间的二阶导插值函数si(x)计算得到所述插值节点x
′k对应的函数值y
′k，其中，k表示插值节点的数量；
[0036]
相应地，所述根据所述插值节点和所述插值节点对应的函数值，对所述原始数字信号行插值处理，得到第一数据，包括：
[0037]
根据所述插值节点x
′k和所述插值节点x
′k对应的函数值y
′k，对所述原始数字信号进行插值处理，得到第一数据。
[0038]
优选地，所述数字信号处理方法中，所述获取原始数字信号具体包括：
[0039]
根据采样频率fs，时间长度为t，时间序列为时间间隔δt＝1/fs，采样得到原始数字信号：其中，na表示采样点数，na＝tfs。
[0040]
优选地，所述数字信号处理方法中，所述利用二阶导插值方法，运用函数二阶导数对所述原始数字信号进行插值处理，得到第一数据具体包括：
[0041]
利用二阶导插值方法，对所述原始数字信号进行插值处理，得到第一数据：其中，nb表示插值处理后的采样点数，nb＝tf
′s，f
′s表示插值处理后的采样频率，所述第一数据对应时间序列为
[0042]
优选地，所述数字信号处理方法中，所述采样得到原始数字信号数据之后，还包括：
[0043]
对所述原始数字信号，按照分段时间间隔为t
dmin
进行分段处理，则插值后数据的采样频率f
′s根据所述分段时间间隔t
dmin
计算得到：
[0044][0045]
式中，m取大于等于2的正整数；
[0046]
选择截断时间点作为插值点，利用二阶导插值方法，将采样频率为fs的所述原始数字信号插值为采样频率为f
′s的数据，得到第一数据：(z
′1,z
′2,
…
,z
′n′
)，其中，下标n
′
值为：n
′
＝[tf
′s]，[]为向下取整符号，所述第一数据对应的时间序列为(t
′1,t
′2,
…
,t
′n′
)；
[0047]
对所述第一数据，按照所述分段时间间隔t
dmin
进行截断处理，得到第二数据。
[0048]
优选地，所述数字信号处理方法中，所述利用二阶导插值方法，运用函数二阶导数对所述原始数字信号进行插值处理，得到第一数据具体包括：
[0049]
对所述原始数字信号，利用二阶导插值方法，将采样频率为fs的所述原始数字信号插值为采样频率为f
′s的数据，得到第一数据：(z
′1,z
′2,
…
,z
′n′
)，其中，下标n
′
值为：n
′
＝[tf
′s]，[]为向下取整符号，f
′s大于fs，所述第一数据对应的时间序列为(t
′1,t
′2,
…
,t
′n′
)。
[0050]
优选地，所述数字信号处理方法中，所述得到第一数据之后，还包括：
[0051]
对所述第一数据进行抗混叠滤波处理，得到第三数据。
[0052]
本技术的第二个目的为提供一种数字信号处理系统。
[0053]
本技术的上述申请目的二是通过以下技术方案得以实现的：
[0054]
一种数字信号处理系统，其特征在于，包括：
[0055]
数据获取模块，用于获取原始数字信号；
[0056]
插值处理模块，用于利用二阶导插值方法，运用函数二阶导数对所述原始数字信号进行插值处理，得到第一数据。
[0057]
本技术的第三个目的为提供一种数字信号处理终端。
[0058]
本技术的上述申请目的三是通过以下技术方案得以实现的：
[0059]
一种数字信号处理终端，包括：存储介质和处理器；
[0060]
所述存储介质内存储计算机执行指令；
[0061]
所述处理器执行所述存储介质存储的计算机执行指令，以实现上述数字信号处理方法中任一所述方法。
[0062]
本技术的第四个目的为提供一种计算机可读存储介质。
[0063]
本技术的上述申请目的四是通过以下技术方案得以实现的：
[0064]
一种计算机可读存储介质，所述计算机可读存储介质中存储有计算机执行指令，所述计算机执行指令被处理器执行时用于实现上述数字信号处理方法中任一所述方法。
[0065]
上述技术方案，利用二阶导插值方法，对所述原始数字信号进行插值处理，实现了对原始数字信号的进一步处理，克服了线性插值存在信号失真、峰值丢失、引发假频率成分等问题，保证插值精度，使得插值得到的数据和真实数据的误差大幅度降低；同时，相比高阶次线性插值法还极大缩短了处理时间，提高了数字信号处理的效率，有利于在线实时监控中进行数据处理；综上可知，上述技术方案能够提高插值算法进行数字信号处理的插值精度，降低插值得到的数据和真实数据的误差。
附图说明
[0066]
为了更清楚地说明本技术实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本技术中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
[0067]
图1为本技术实施例中提供的一种数字信号处理方法的流程示意图；
[0068]
图2为本技术实施例中插值处理的流程示意图；
[0069]
图3为本技术实施例中原始数字信号的示意图；
[0070]
图4为本技术实施例中第一数据的示意图；
[0071]
图5为本技术实施例中提供的一种数字信号处理系统的结构示意图；
[0072]
图6为本技术实施例中提供的一种数字信号处理终端的结构示意图。
具体实施方式
[0073]
为了使本领域的技术人员更好地理解本技术中的技术方案，下面将对本技术实施
例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本技术的一部分实施例，而不是全部的实施例。基于本技术中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本技术保护的范围。
[0074]
在本技术所提供的实施例中，应该理解到，所揭露的方法和系统，可以通过其它的方式实现。以下所描述的系统实施例仅仅是示意性的，例如，模块的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，如：多个模块或组件可以结合，或可以集成到另一个系统，或一些特征可以忽略，或不执行。另外，所显示或讨论的各组成部分相互之间的耦合、或直接耦合、或通信连接可以是通过一些接口，设备或模块的间接耦合或通信连接，可以是电性的、机械的或其它形式的。
[0075]
本技术实施例采用递进的方式撰写。
[0076]
如图1所示，本技术实施例提供一种数字信号处理方法，包括：
[0077]
s101.获取原始数字信号；
[0078]
在s101中，具体地，可以采集不同信号源获取相应的数字信号作为所述原始数字信号，其中，数字信号可以是通过振动传感器获取振动数据，也可以是温度传感器获取温度数据等。其具体的获取方式不影响本实施例的实现。
[0079]
s102.利用二阶导插值方法，运用函数二阶导数对所述原始数字信号进行插值处理，得到第一数据。
[0080]
在s102中，利用所述二阶导插值法，可以运用函数二阶导数的思想构建插值函数，实现数据插值；具体地，可以利用二阶导插值函数，对所述原始数字信号进行插值处理。例如，可以将所述原始数据信号划分成多个区间，建立每个区间对应的二阶导插值函数并求解所述二阶导插值函数，最终得到每个区间对应的插值节点和插值节点对应的函数值，利用所述插值节点和插值节点对应的函数值，对所述原始数字信号进行插值处理，得到第一数据。
[0081]
插值作为一种数值逼近的重要方法，可以通过有限个样本的取值情况估计出未知样本的近似值。目前常见的数字信号插值算法有拉格朗日插值、牛顿插值、厄米特插值等。这些插值方式均为线性插值，计算简单稳定，但是其插值误差较大，而且光滑性较差，会导致数字信号失真、峰值丢失、还引发虚假频率成分等问题。
[0082]
上述实施例，利用二阶导插值方法，对所述原始数字信号进行插值处理，实现了对原始数字信号的进一步处理，克服了线性插值存在信号失真、峰值丢失、引发假频率成分等问题，保证插值精度，使得插值得到的数据和真实数据的误差大幅度降低；同时，相比高阶次线性插值法还极大缩短了处理时间，提高了数字信号处理的效率，有利于在线实时监控中进行数据处理；综上可知，上述实施例能够提高插值算法进行数字信号处理的插值精度，降低插值得到的数据和真实数据的误差。
[0083]
如图2所示，在基于上述实施例的基础上，所述利用二阶导插值方法，运用函数二阶导数对所述原始数字信号进行插值处理，得到第一数据的步骤的其中一种实现方式可以包括以下步骤：
[0084]
s201.在所述原始数字信号相邻两个数据节点组成的各区间，构建二阶导插值函数；
[0085]
在s201中，具体地，可以设所述原始数字信号包含n个数据节点，分别为(x1,y1),
(x2,y2),
…
,(xn,yn)，n为大于1的正整数，得到所述原始数字信号相邻两个数据节点组成的各区间[xi,x
i+1
](i＝1,2,
…
,n-1)；
[0086]
构建所述各区间[xi,x
i+1
](i＝1,2,
…
,n-1)的二阶导插值函数si(x)＝ai+bi(x-xi)+ci(x-xi)2+di(x-xi)3(xi≤x≤x
i+1
)，式中，x表示插值函数自变量，ai,bi,ci,di为二阶导插值函数的各项系数，具体为：
[0087][0088]
式中，mi为二次微分值，步长hi＝x
i+1-xi(i＝1,2,
…
,n-1)。
[0089]
s202.构建所述二阶导插值函数求解的矩阵方程组；
[0090]
在s202中，具体地，可以根据所述数据节点(x1,y1),(x2,y2),
…
,(xn,yn)构建矩阵方程式一：s(xi)＝yi(i＝1,2,
…
,n)，式中，y表示数据节点x对应的函数值，i表示数据节点的个数，n为大于1的正整数；
[0091]
根据所述s(x)、s
′
(x)、s
″
(x)在所述原始数字信号中具有二阶连续导数时，构建方矩阵方程式二：
[0092][0093]
式中，i表示数据节点的个数，n为大于1的正整数；
[0094]
根据所述s(x)在所述原始数字信号中为一个整周期数据时，构建矩阵方程式三：
[0095][0096]
式中，n为大于1的正整数。
[0097]
s203.将所述原始数字信号的各数据节点代入所述矩阵方程组，求解得到二次微分值；
[0098]
在s203中，具体地，可以将所述数据节点代入所述矩阵方程式一、所述矩阵方程式二、所述矩阵方程式三，求解得到二次微分值mi＝si″
(xi)。
[0099]
s204.将所述二次微分值代入到所述二阶导插值函数，计算得出所述各区间的二阶导插值函数；
[0100]
在s204中，具体地，可以根据所述二次微分值、所述数据节点、所述步长，求解得到ai,bi,ci,di，并根据所述ai,bi,ci,di，计算得出所述各区间的二阶导插值函数si(x)。
[0101]
s205.判断插值节点所在的区间，根据所述插值节点所在的区间的二阶导插值函数，计算得到所述插值节点对应的函数值；
[0102]
在s205中，具体地，可以判断插值节点x
′k所在区间[xi,x
i+1
](i＝1,2,
…
,n-1)，根据所述插值节点x
′k所在的区间的二阶导插值函数si(x)计算得到所述插值节点x
′k对应的
函数值y
′k，其中，k表示插值节点的数量。
[0103]
s206.根据所述插值节点和所述插值节点对应的函数值，对所述原始数字信号进行插值处理，得到第一数据。
[0104]
在s206中，具体地，可以根据所述插值节点x
′k和所述插值节点x
′k对应的函数值y
′k，对所述原始数字信号进行插值处理，得到第一数据。
[0105]
上述实施例提出的基于二阶导插值的数字信号处理方法，能够适用于不同场景的需要，例如，在基于上述实施例的另一实施例中，当需要对信号进行截断时，需要截断的时刻没有采样点，则可利用本方法插值后进行截取，如整周期截断时包括：
[0106]
s301.根据采样频率fs，时间长度为t，时间序列为时间间隔δt＝1/fs，采样得到原始数字信号：其中，na表示采样点数，na＝tfs。
[0107]
s302.对所述原始数字信号，按照分段时间间隔为t
dmin
进行分段处理，则插值后数据的采样频率f
′s根据所述分段时间间隔t
d min
计算得到：
[0108][0109]
式中，m取大于等于2的正整数，其可以根据实际使用情况选择。
[0110]
s303.选择截断时间点作为插值点，利用二阶导插值方法，将采样频率为fs的所述原始数字信号插值为采样频率为f
′s的数据，得到第一数据：(z
′1,z
′2,
…
,z
′n′
)，其中，下标n
′
值为：n
′
＝[tf
′s]，[]为向下取整符号，所述第一数据对应的时间序列为(t
′1,t
′2,
…
,t
′n′
)；
[0111]
其中，插值后时间间隔δt
′
＝1/f
′s，所述第一数据表示插值处理后得到的数据。
[0112]
s304.对所述第一数据，按照所述分段时间间隔t
dmin
进行截断处理，得到第二数据。
[0113]
在s304中，对插值处理后得到的第一数据进行截断处理，得到第二数据。
[0114]
在本实施例中，当对原始数字信号进行截断处理时，而截断点并没有采样点时，采用二阶导插值法对截断点进行插值之后，再进行数据截取，这样该段数据就是有起点有终点的完整数据。
[0115]
在基于上述实施例的另一实施例中，当需要对由于数字信号采样导致的峰值丢失进行恢复时，则可利用本方法插值后恢复丢失峰值，包括：
[0116]
s401.根据采样频率fs，时间长度为t，时间序列为时间间隔δt＝1/fs，采样得到原始数字信号：其中，na表示采样点数，na＝tfs。
[0117]
s402.对所述原始数字信号，利用二阶导插值方法，将采样频率为fs的所述原始数字信号插值为采样频率为f
′s的数据，得到第一数据：(z
′1,z
′2,
…
,z
′n′
)，其中，下标n
′
值为：n
′
＝[tf
′s]，[]为向下取整符号，f
′s大于fs，所述第一数据对应的时间序列为(t
′1,t
′2,
…
,t
′n′
)。
[0118]
其中，插值后时间间隔δt
′
＝1/f
′s，所述第一数据表示插值处理后得到的数据，其可以有效还原丢失的峰值。
[0119]
以如下具体实施例为例：
[0120]
原始数字信号为：幅值为1，频率为22hz的正弦信号；由采样频率为100hz得到的数字信号，如图3所示。
[0121]
信号中由于数字采样导致出现峰值丢失现象，最小峰值为0.77，丢失23％。
[0122]
通过本文所提二阶导插值方法，将采样频率为100hz数字信号插值为采样频率为400hz的第一数据，如图4所示。
[0123]
信号中由于数字采样导致出现峰值丢失的位置，通过插值将峰值还原为0.9852，峰值丢失从23％缩小为1.5％，极大解决了原始数字信号的峰值丢失问题。
[0124]
在基于上述实施例的另一实施例中，将二阶导插值方法应用于振动信号处理领域，提供另一种数字信号处理方法，包括：
[0125]
s501.根据采样频率fs，时间长度为t，时间序列为时间间隔δt＝1/fs，采样得到原始数字信号：其中，na表示采样点数，na＝tfs。
[0126]
在s501中，具体地，所述原始数字信号为振动信号，采集所述振动信号的振动源可以是旋转部件。
[0127]
s502.利用二阶导插值方法，对所述原始数字信号进行插值处理，得到第一数据：其中，nb表示插值处理后的采样点数，nb＝tf
′s，f
′s表示插值处理后的采样频率，所述第一数据对应时间序列为
[0128]
在s502中，具体地，根据所述原始数字信号，按照旋转周期f
n1
，每旋转周期内采样n1个点，利用二阶导插值方法，对所述原始数字信号进行插值处理，得到第一数据：所述第一数据表示插值处理后得到的数据，其中，f
′s＝f
n1
n1，所述旋转周期f
n1
可以是转频周期、发动机叶片旋转周期、保持架旋转周期等；插值处理后数据的时间序列为时间间隔δt
′
＝1/f
′s，对应的原始数字信号：其中，t1＝t
′1，z1＝z
′1，
[0129]
在本实施例中，当原始数字信号为振动信号时，按某个旋转的部件的旋转周期，如转频周期、发动机叶片旋转周期、保持架旋转周期等截取整周期时，通过上述二阶导插值方法实现对振动信号整周期截取。
[0130]
在基于上述实施例的另一实施例中，所述得到第一数据之后，还包括：
[0131]
s601.对所述第一数据进行抗混叠滤波处理，得到第三数据。
[0132]
在s601中，对所述第一数据，根据插值处理后的采样频率，进行抗混叠滤波，得到所述第三数据，具体地，所述抗混叠滤波处理可以采用低通滤波器完成，截止频率需小于插值处理后的采样频率的一半。抗混叠滤波处理目的是减少频谱混叠，最后得到所述第三数据。
[0133]
在本技术的其他实施例中，还可以在原始数字信号，设置采样频率为fs，通过上述二阶导插值方法，将采样频率变化为f
′s，实现对数字信号采样频率变化。
[0134]
在工程信号中，既含有频率与旋转机械的转速频率v成正比例的被检测信号，例如轴承、齿轮、车轮的故障信号以及轨道损伤与车轮作用的信号，它们最宜用转速跟踪采样的
方法检测，即等圆周角度间隔进行信号采集；也含有频率与旋转机械的转速频率无关的、基本固定频率的被检测信号，例如机械的广义共振信号、轨道波浪纹与车轮作用的信号，它们最宜用天文时钟采样的方法检测，即按等时间间隔进行信号采集。在本技术的其他实施例中，原始数字信号的采集可以采用天文时钟采样，根据采集得到的转速脉冲信号通过二阶导插值可以为转化为转速跟踪采样数据；在本技术的其他实施例中，也可以在原始数字信号采集时，采用转速跟踪采样，通过二阶导插值转化为为天文时钟采样数据，以此实现采样方式转换。
[0135]
本技术实施例提供的所述数字信号处理方法，不仅保证了数字信号处理的精度要求，还能够适用于不同场景的需要，不局限于上述实施例中功能的实现，可以通过二阶导插值法处理数字信号的所有功能。
[0136]
如图5所示，在本技术的另一实施例中，还提供一种数字信号处理系统，包括：
[0137]
数据获取模块10，用于获取原始数字信号；
[0138]
插值处理模块11，用于利用二阶导插值方法，运用函数二阶导数对所述原始数字信号进行插值处理，得到第一数据。
[0139]
如图6所示，在本技术的另一实施例中，还提供一种数字信号处理终端，包括：存储介质12和处理器13；
[0140]
所述存储介质12内存储计算机执行指令；
[0141]
所述处理器13执行所述存储介质12存储的计算机执行指令，以实现上述数字信号处理方法中任一所述方法。
[0142]
其中，处理器13可以包括一个或者多个处理核心。处理器13通过运行或执行存储在存储介质12内的指令、程序、代码集或指令集，调用存储在存储介质12内的数据，执行本技术的各种功能和处理数据。处理器13可以为特定用途集成电路、数字信号处理器、数字信号处理装置、可编程逻辑装置、现场可编程门阵列、中央处理器、控制器、微控制器和微处理器中的至少一种。可以理解地，对于不同的设备，用于实现上述处理器13功能的电子器件还可以为其它。
[0143]
其中，存储介质12可用于存储指令、程序、代码、代码集或指令集。存储介质12可以包括存储程序区和存储数据区，其中存储程序区可存储用于实现操作系统的指令、用于至少一个功能的指令以及用于实现上述任一数字信号处理方法的指令等；存储数据区可存储上述任一数字信号处理方法中涉及到的数据等。
[0144]
在本技术的另一实施例中，还提供一种计算机可读存储介质，所述计算机可读存储介质中存储有计算机执行指令，所述计算机执行指令被处理器执行时用于实现上述数字信号处理中任一所述方法。
[0145]
其中，所述计算机可读存储介质，可以为u盘、移动硬盘、只读存储器、随机存取存储器或者光盘等各种可以存储程序代码的介质。
[0146]
对所公开的实施例的上述说明，使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的，本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下，在其它实施例中实现。因此，本发明将不会被限制于本文所示的这些实施例，而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。