本发明涉及信号处理， 特别地， 涉及谐波分析方法和装置以及确定谐波 间杂波的方法和装置。 背景技术

谐波分析是一种信号处理与分析的基本手段， 在信号识别、分类、编码、 增强等领域都有着广泛的应用。 谐波分析可以涵盖多种不同的分析目的， 例 如计算谐波的个数、 计算谐波的基频、 计算谐波的强弱、 计算谐波的帧内 / 帧间相关性、 计算谐噪比等等。 谐波分析可以基于时域信号进行分析， 也可 以基于频域信号进行分析。 一般来说， 谐波分析在频域信号上进行最常见， 效果也更好。

在现有的音频编码算法中， 一种常用的方法是在频域上釆用峰均比的方 法来确定谐波。 但是， 如果音频帧的频谱是 "模糊" 的， 即谐波间有较强的 杂波存在， 于是谐波的能量被削弱了， 此时含有谐波的子带的峰均比与不含 谐波的子带的峰均比的区别就会不那么明显， 用简单的将峰均比与门限做比 较的方法判断谐波是否存在时就会引入较大的 误差。

另一种常用的方法是根据波峰、 波谷的相对能量（或幅度） 关系得到波 峰的强度， 最后通过将这个强度与一个门限做比较来确定 是否存在谐波。 但 是， 如果音频帧的频谱是 "模糊" 的， 即谐波间有较强的杂波存在， 由于谐 波间杂波的存在，而可能把谐波附近的杂波的 波谷错误的当作是谐波的波谷， 从而导致计算出的谐波强度低于真实值而被漏 掉。

由此可见，在进行频域的谐波分析时，无论谐 波分析的具体目的是什么， 一个前提条件是要能得到一个尽量准确的谐波 结构 , 或者说一个尽量清晰的 谐波频语。 然而， 实际的时频变换所得到的谐波频谱通常都会包 含各种干扰 因素， 使得得到的谐波频谱被 "模糊" 了。 这种 "模糊" 可能使原本较强的 谐波被削弱， 也可能在谐波间引入杂波干扰。 在谐波分析时， 这些 "模糊" 可能会导致错误的或性能下降的分析结果。 因此， 如何得到一个尽可能 "清 晰" 的谐波频谱， 是决定谐波分析成败的关键因素。 发明内容

本发明提出了谐波分析方法和装置以及确定谐 波间杂波的方法和装置， 旨在解决如何获得清晰的和准确的谐波频谱的 问题。

第一方面， 提出了一种谐波分析的方法， 包括： 对输入的时域信号进行 时频转换， 以获得所述时域信号在频域上的频谱；对所述 频谱进行滤波处理， 以确定目标频谱； 对所述目标频谱的谐波间杂波做消除处理， 以便确定所述 目标频谱的谐波的结构； 利用所述目标频谱的谐波的结构对所述谐波进 行谐 波分析。

结合第一方面， 在第一方面的第一实施方式中， 所述对所述频谱进行滤 波处理， 以确定目标频谱， 包括： 对所述频谱进行 k次滤波处理， 并在所述 频谱经过所述 k次滤波处理的每次滤波处理之后记录波峰个� �以获得 k个波 峰个数， 其中 k为大于或等于 2的正整数； 确定所述 k个波峰个数中的最小 波峰个数； 确定得到所述最小波峰个数的滤波处理的频谱 为目标频谱。

结合第一方面， 在第一方面的第二实施方式中， 所述对所述频谱进行滤 波处理， 以确定目标频谱， 包括： 对所述频谱进行第一滤波处理， 并记录经 过所述第一滤波处理的频谱的第一波峰个数;� �所述频谱进行第二滤波处理， 并记录经过所述第二滤波处理的频谱的第二波 峰个数； 对所述频谱进行第三 滤波处理， 并记录经过所述第三滤波处理的频谱的第三波 峰个数； 当所述第 二波峰个数小于所述第一波峰个数， 且所述第二波峰个数小于所述第三波峰 个数， 则确定经过所述第二滤波处理得到的频谱为目 标频谱。

结合第一方面的第二实施方式， 在第一方面的第三实施方式中， 所述第 一滤波处理为第一阶滤波、 所述第二滤波处理为第二阶滤波， 以及所述第三 滤波处理为第三阶滤波， 所述对所述频谱进行滤波处理， 以确定目标频谱， 包括： 对所述频谱进行第一阶滤波， 并记录经过所述第一阶滤波的频谱的第 一波峰个数； 对所述频谱进行第二阶滤波， 并记录经过所述第二阶滤波的频 谱的第二波峰个数； 对所述频谱进行第三阶滤波， 并记录经过所述第三阶滤 波的频谱的第三波峰个数； 当所述第二波峰个数小于所述第一波峰个数， 且 所述第二波峰个数小于所述第三波峰个数， 则确定经过所述第二阶滤波的频 谱为目标频谱。

结合第一方面， 在第一方面的第四实施方式中， 所述对所述频谱进行滤 波处理， 以确定目标频谱， 包括： 釆用低通滤波器对所述频谱进行低通滤波 处理， 将所述经过所述低通滤波处理后的频语作为目 标频谱。

结合第一方面， 在第一方面的第五实施方式中， 所述对所述频谱进行滤 波处理， 以确定目标频谱， 包括： 对所述频谱进行第四阶滤波； 确定经过所 述第四阶滤波得到的频谱为目标频谱， 其中所述第四阶为固定值。

结合第一方面或第一至第五实施方式， 在第一方面的第六实施方式中， 所述对所述目标频谱的谐波间杂波做消除处理 ， 包括釆用以下方式中的至少 一种： 依据波瓣宽度对所述目标频谱的谐波间杂波做 消除处理， 其中所述波 瓣宽度是波峰两侧的两个最邻近的波谷之间的 宽度； 依据波峰的倾斜度对所 述目标频谱的谐波间杂波做消除处理， 其中所述波峰的倾斜度是所述波峰两 侧的两个最邻近的波谷之间的连线与水平方向 的夹角， 或者所述波峰两侧的 两个最邻近的波谷之间的连线的斜率， 或者所述波峰两侧的两个最邻近的波 谷之间的连线的夹角的三角函数。

结合第一方面的第六实施方式， 在第一方面的第七实施方式中， 所述依 据波瓣宽度对所述目标频谱的谐波间杂波做消 除处理， 包括： 确定所述目标 频谱的全部波峰和波谷； 根据确定的所述波峰和所述波谷， 确定所述目标频 谱中的波瓣宽度； 当所述波瓣宽度小于第一波瓣宽度， 则确定所述波瓣宽度 对应的波峰为杂波的波峰； 通过消去所述杂波的波瓣来消除所述目标频谱 的 杂波。

结合第一方面的第七实施方式， 在第一方面的第八实施方式中， 所述第 一波瓣宽度为预设的固定值， 或者所述第一波瓣宽度为平均波瓣宽度， 其中 所述平均波瓣宽度是所述目标频谱的宽度与所 述第二波峰个数的比值。

结合第一方面的第六实施方式， 在第一方面的第九实施方式中， 所述依 据波峰的倾斜度消除所述目标频谱的杂波， 包括： 确定所述目标频谱的全部 波峰和波谷； 根据确定的所述波峰和所述波谷， 确定所述波峰的倾斜度； 当 所述波峰的倾斜度大于第一阔值， 则确定所述波峰为杂波波峰； 通过消去所 述杂波的波瓣来消除所述目标频谱的杂波。 结合第一方面的第九实施方式， 在第一方面的第十实施方式中， 所述第 一阔值为预设的固定值。

结合第一方面的第七实施方式或第九实施方式 ， 在第一方面的第十一实 将所述杂波的波瓣位于两个波谷之间的频点幅 值以在两个波谷幅值之间的内 插来代替。

第二方面， 提出了一种确定谐波间杂波的方法， 包括： 确定所述频谱的 全部波峰和波谷； 根据确定的所述波峰和所述波谷， 确定所述频谱中的波瓣 宽度， 其中所述波瓣宽度是波峰两侧的两个最邻近的 波谷之间的宽度； 当所 述波瓣宽度小于第一波瓣宽度， 则确定所述波瓣宽度对应的波峰为杂波的波 峰。

结合第二方面， 在第二方面的第一实施方式中， 所述第一波瓣宽度为预 设的固定值。

结合第二方面， 在第二方面的第二实施方式中， 所述第一波瓣宽度为平 均波瓣宽度， 其中所述平均波瓣宽度是所述频谱的宽度与所 述全部波峰的个 数的比值。

第三方面， 提出了一种确定谐波间杂波的方法， 包括： 确定所述目标频 谱的全部波峰和波谷； 根据确定的所述波峰和所述波谷， 确定所述波峰的倾 斜度； 当所述波峰的倾斜度大于第一阔值， 则确定所述波峰为杂波的波峰。

结合第三方面， 在第三方面的第一实施方式中， 所述波峰的倾斜度是所 述波峰两侧的两个最邻近的波谷之间的连线与 水平方向的夹角； 或者所述波 峰的倾斜度是所述波峰两侧的两个最邻近的波 谷之间的连线的斜率： 或者所 述波峰的倾斜度是所述波峰两侧的两个最邻近 的波谷之间的连线的夹角的三 角函数。

结合第三方面或第三方面的第一实施方式， 在第三方面的第二实施方式 中， 所述第一阔值为预设的固定值。

第四方面， 提出了一种谐波分析的装置， 包括： 时频转换单元， 用于对 输入的时域信号进行时频转换， 以获得所述时域信号在频域上的频谱； 滤波 单元， 用于对所述时频转换单元转换得到的所述频谱 进行滤波处理， 以确定 目标频谱； 杂波消除单元， 用于对所述滤波单元确定的所述目标频谱的谐 波 间杂波做消除处理， 以便确定所述目标频谱的谐波的结构； 谐波分析单元， 用于利用所述目标频谱的谐波的结构对由所述 杂波消除单元消除了杂波的所 述谐波进行谐波分析。

结合第四方面，在第四方面的第一实施方式中 ，所述滤波单元具体用于： 对所述频谱进行 k次滤波处理， 并在所述频谱经过所述 k次滤波处理的每次 滤波处理之后记录波峰个数以获得 k个波峰个数， 其中 k为大于或等于 2的 正整数； 确定所述 k个波峰个数中的最小波峰个数； 确定得到所述最小波峰 个数的滤波处理的频谱为目标频谱。

结合第四方面，在第四方面的第二实施方式中 ，所述滤波单元具体用于： 对所述频谱进行第一滤波处理， 并记录经过所述第一滤波处理的频谱的第一 波峰个数； 对所述频谱进行第二滤波处理， 并记录经过所述第二滤波处理的 频谱的第二波峰个数； 对所述频谱进行第三滤波处理， 并记录经过所述第三 滤波处理的频谱的第三波峰个数； 当所述第二波峰个数小于所述第一波峰个 数， 且所述第二波峰个数小于所述第三波峰个数， 则确定经过所述第二滤波 处理得到的频谱为目标频谱。

结合第四方面的第二实施方式， 在第四方面的第三实施方式中， 所述第 一滤波处理为第一阶滤波、 所述第二滤波处理为第二阶滤波， 以及所述第三 滤波处理为第三阶滤波， 所述滤波单元具体用于： 对所述频谱进行第一阶滤 波， 并记录经过所述第一阶滤波的频谱的第一波峰 个数； 对所述频谱进行第 二阶滤波， 并记录经过所述第二阶滤波的频谱的第二波峰 个数； 对所述频谱 进行第三阶滤波， 并记录经过所述第三阶滤波的频谱的第三波峰 个数； 当所 述第二波峰个数小于所述第一波峰个数， 且所述第二波峰个数小于所述第三 波峰个数， 则确定经过所述第二阶滤波的频谱为目标频谱 。

结合第四方面，在第四方面的第四实施方式中 ，所述滤波单元具体用于： 釆用低通滤波器对所述频谱进行低通滤波处理 ， 将所述经过所述低通滤波处 理后的频语作为目标频谱。

结合第四方面，在第四方面的第五实施方式中 ，所述滤波单元具体用于： 对所述频谱进行第四阶滤波；确定经过所述第 四阶滤波到的频谱为目标频谱； 其中所述第四阶为固定值。

结合第四方面或第一至第五实施方式， 在第四方面的第六实施方式中， 所述杂波消除单元釆用以下方式中的至少一种 ： 依据波瓣宽度对所述目标频 谱的谐波间杂波做消除处理， 其中所述波瓣宽度是波峰两侧的两个最邻近的 波谷之间的宽度； 依据波峰的倾斜度对所述目标频谱的谐波间杂 波做消除处 理， 其中所述波峰的倾斜度是所述波峰两侧的两个 最邻近的波谷之间的连线 与水平方向的夹角， 或者所述波峰两侧的两个最邻近的波谷之间的 连线的斜 率， 或者所述波峰两侧的两个最邻近的波谷之间的 连线的夹角的三角函数。

结合第四方面的第六实施方式， 在第四方面的第七实施方式中， 所述杂 波消除单元具体用于： 确定所述目标频谱的全部波峰和波谷； 根据确定的所 述波峰和所述波谷， 确定所述目标频谱中的波瓣宽度； 当所述波瓣宽度小于 第一波瓣宽度， 则确定所述波瓣宽度对应的波峰为杂波的波峰 ； 通过消去所 述杂波的波瓣来消除所述目标频谱的杂波。

结合第四方面的第七实施方式， 在第四方面的第八实施方式中， 所述第 一波瓣宽度为预设的固定值， 或者所述第一波瓣宽度为平均波瓣宽度， 其中 所述平均波瓣宽度是所述目标频谱的宽度与所 述第二波峰个数的比值。

结合第四方面的第六实施方式， 在第四方面的第九实施方式中， 所述杂 波消除单元具体用于： 确定所述目标频谱的全部波峰和波谷； 根据确定的所 述波峰和所述波谷， 确定所述波峰的倾斜度； 当所述波峰的倾斜度大于第一 阔值， 则确定所述波峰为杂波波峰； 通过消去所述杂波的波瓣来消除所述目 标频谱的杂波。

结合第四方面的第九实施方式， 在第四方面的第十实施方式中， 所述第 一阔值为预设的固定值。

结合第四方面的第七或第九实施方式，在第四 方面的第十一实施方式中， 所述杂波消除单元具体用于： 将所述杂波的波瓣位于两个波谷之间的频点幅 值以在两个波谷幅值之间的内插来代替。

第五方面， 提出了一种确定谐波间杂波的装置， 包括： 第一确定单元， 用于确定所述频谱的全部波峰和波谷； 第二确定单元， 用于根据所述第一确 定单元确定的所述波峰和所述波谷， 确定所述频谱中的波瓣宽度， 其中所述 波瓣宽度是波峰两侧的两个最邻近的波谷之间 的宽度； 第三确定单元， 用于 当所述波瓣宽度小于第一波瓣宽度， 则确定由所述第二确定单元确定的所述 波瓣宽度对应的波峰为杂波的波峰。

结合第五方面， 在第五方面的第一实施方式中， 所述第一波瓣宽度为预 设的固定值。

结合第五方面， 在第五方面的第一实施方式中， 所述第一波瓣宽度为平 均波瓣宽度， 其中所述平均波瓣宽度是所述频谱的宽度与所 述全部波峰的个 数的比值。

第六方面， 提出了一种确定谐波间杂波的装置， 包括： 第四确定单元， 用于确定所述目标频谱的全部波峰和波谷； 第五确定单元， 用于根据所述第 四确定单元确定的所述波峰和所述波谷， 确定所述波峰的倾斜度； 第六确定 单元， 用于当所述波峰的倾斜度大于第一阔值， 则确定由所述第五确定单元 确定的所述波峰为杂波的波峰。

结合第六方面， 在第六方面的第一实施方式中， 所述波峰的倾斜度是所 述波峰两侧的两个最邻近的波谷之间的连线与 水平方向的夹角； 或者所述波 峰的倾斜度是所述波峰两侧的两个最邻近的波 谷之间的连线的斜率： 或者所 述波峰的倾斜度是所述波峰两侧的两个最邻近 的波谷之间的连线的夹角的三 角函数。

结合第六方面或第六方面的第一实施方式， 在第四方面的第二实施方式 中， 所述第一阔值为预设的固定值。

本发明实施例通过滤波确定谐波分析的目标频 谱， 并且能够有效地消除 杂波干扰， 从而获得清晰的和准确的谐波频谱， 以利于进行谐波分析。 附图说明

为了更清楚地说明本发明实施例的技术方案， 下面将对本发明实施例中 所需要使用的附图作简单地介绍， 显而易见地， 下面所描述的附图仅仅是本 发明的一些实施例， 对于本领域普通技术人员来讲， 在不付出创造性劳动的 前提下， 还可以根据这些附图获得其他的附图。

图 1是根据本发明实施例的谐波分析的方法的流� �图。 图 2是消除杂波的示意图。

图 3示出了通过波峰的倾斜度确定杂波的方法。

图 4是才艮据本发明实施例的确定谐波间杂波的� �法的流程图。

图 5是根据本发明另一实施例的确定谐波间杂波� �方法的流程图。

图 6是根据本发明实施例的谐波分析的装置的结� �示意图。

图 7是才艮据本发明实施例的确定谐波间杂波的� �置的结构示意图。

图 8是才艮据本发明另一实施例的确定谐波间杂� �的装置的结构示意图。 图 9是根据本发明实施例的编码器的结构示意图� � 具体实施方式

下面将结合本发明实施例中的附图， 对本发明实施例中的技术方案进行 清楚、 完整地描述， 显然， 所描述的实施例是本发明的一部分实施例， 而不 是全部实施例。 基于本发明中的实施例， 本领域普通技术人员在没有做出创 造性劳动的前提下所获得的所有其他实施例， 都应属于本发明保护的范围。 以下将结合图 1 , 描述根据本发明实施例的谐波分析的方法， 包括如下步骤。

11 , 对输入的时域信号进行时频转换， 以获得所述时域信号在频域上的 频谱。

例如， 快速傅里叶变换（FFT, Fast Fourier Transform )和离散余弦变换 ( DCT, Discrete Cosine Transform )等是常用的时频转换方法， 通过上述方 法对输入的时域信号进行时频转换， 以获得所述输入的时域信号在频域上的 频谱。

12, 对所述频谱进行滤波处理， 以确定目标频谱。

在一种实施方式中， 对所述频语进行 k次滤波处理， 并在所述频语经过 所述 k次滤波处理的每次滤波处理之后记录波峰个� �以获得 k个波峰个数， 其中 k为大于或等于 2的正整数；确定所述 k个波峰个数中的最小波峰个数； 确定得到所述最小波峰个数的滤波处理的频谱 为目标频谱。

在另一种实现方式中， 对所述频谱进行第一滤波处理， 并记录经过所述 第一滤波处理的频谱的第一波峰个数； 对所述频谱进行第二滤波处理， 并记 录经过所述第二滤波处理的频谱的第二波峰个 数； 对所述频谱进行第三滤波 处理， 并记录经过所述第三滤波处理的频谱的第三波 峰个数； 当所述第二波 峰个数小于所述第一波峰个数，且所述第二波 峰个数小于所述第三波峰个数， 则确定经过所述第二滤波处理的频谱为目标频 谱。

具体而言， 第一滤波处理可以为第一阶滤波、 第二滤波处理可以为第二 阶滤波， 以及第三滤波处理可以为第三阶滤波。 应理解， 第一阶滤波、 第二 阶滤波和第三阶滤波分别对应的滤波阶数是不 同的。 这样， 对所述频谱进行 第一阶滤波， 并记录经过所述第一阶滤波的频谱的第一波峰 个数； 对所述频 谱进行第二阶滤波， 并记录经过所述第二阶滤波的频谱的第二波峰 个数； 对 所述频谱进行第三阶滤波， 并记录经过所述第三阶滤波的频谱的第三波峰 个 数； 当所述第二波峰个数小于所述第一波峰个数， 且所述第二波峰个数小于 所述第三波峰个数， 则确定经过所述第二阶滤波的频谱为目标频谱 。

通常， 当第一阶滤波、 第二阶滤波和第三阶滤波均为低通滤波时， 更优 的滤波应该导致目标频谱具有更少的频谱波峰 个数。 可以理解， 上述滤波过 程是可以循环的。 例如， 在经过第一次滤波组合中的第一阶滤波、 第二阶滤 波和第三阶滤波之后， 若第二波峰个数并不同时小于第一波峰个数和 第三波 峰个数， 那么就再次经过第二次滤波组合中的第一阶滤 波、 第二阶滤波和第 三阶滤波， 直到再某次滤波组合中寻找到同时小于第一波 峰个数和第三波峰 个数的第二波峰个数。 也就是说， 这样的循环滤波寻找最小频谱波峰个数的 过程， 就是一个寻找最优滤波、 确定目标频谱的过程。

这样， 随着滤波的阶数由低向高增加， 每轮滤波后的频谱波峰个数会呈 现出一种先单调递减， 再单调递增的现象。 通常， 最优滤波就出现在波峰个 数折线由递减到递增的拐点处。 为了降低计算复杂度， 本实施方式中对确定 目标频谱的循环过程增加了一个提前终止判断 ， 即一旦检测到拐点， 则认为 已找到最优滤波， 不再执行剩下的循环。 本实施方式中每循环一次， 比较当 次循环后的目标频谱的波峰个数是否小于等于 前次循环后的目标频谱的波峰 个数。 若是， 则继续循环， 若否， 则终止循环， 以前次循环后的目标频语作 为后续步骤的目标频谱， 即将前次循环的滤波作为最优滤波。

在另一种实施方式中， 为了进一步降低计算的复杂度， 滤波阶数不再动 态的确定， 而是对每帧都以一个固定的阶数进行滤波。 经过对大的样本进行 统计可以发现， 对一给定的釆样率和帧长， 最优滤波阶数在某一固定阶数的 概率会显著高于其它阶数。 因此以该最大出现概率阶数作为固定的滤波阶 数 可以在节省计算复杂度的同时最大限度的逼近 最优滤波结果。 例如， 对所述 频谱进行第四滤波处理；确定经过所述第四滤 波处理得到的频谱为目标频谱。 具体地， 第四滤波处理为第四阶滤波， 则对所述频谱进行第四阶滤波， 其中 所述第四阶为固定值； 确定经过所述第四阶滤波的频谱为目标频谱。

在另一种实施方式中， 可以釆用任意形式的低通滤波器对所述频谱进 行 低通滤波处理， 以确定目标频谱为经过所述低通滤波处理后的 频谱。

13 , 对所述目标频谱的谐波间杂波做消除处理， 以便确定所述目标频谱 的谐波的结构。

例如， 可以釆用以下方式中的至少一种来对所述目标 频谱的谐波间杂波 做消除处理： 依据波瓣宽度对所述目标频谱的谐波间杂波做 消除处理； 依据 波峰的倾斜度对所述目标频谱的谐波间杂波做 消除处理。

依据波瓣宽度对所述目标频谱的谐波间杂波做 消除处理的方法包括如下 步骤： 首先， 确定所述目标频谱的全部波峰和波谷， 其中频谱中某一频点的 能量高于它的两个最邻近频点的能量时，则认 为该频点处是频谱的一个波峰， 而波谷的定义为两个相邻波峰间的能量最小的 频点； 然后， 根据确定的所述 波峰和所述波谷， 确定所述目标频谱中的波瓣宽度， 其中所述波瓣是一个波 峰两侧的两个最邻近波谷之间的频点所组成的 集合， 所述波瓣宽度是波峰两 侧的两个最邻近的波谷之间的宽度； 接着， 当所述波瓣宽度小于第一波瓣宽 度， 则确定所述波瓣宽度对应的波峰为杂波的波峰 ； 最后， 通过消去所述杂 波的波瓣来消除所述目标频谱的杂波， 例如将所述杂波的波瓣位于两个波谷 之间的频点能量以在两个波谷之间的内插来代 替。

这里， 第一波瓣宽度为预设的固定值， 或者所述第一波瓣宽度为平均波 瓣宽度， 其中所述平均波瓣宽度是所述目标频谱的宽度 与所述第二波峰个数 的比值。 确定所述目标频谱的全部波峰和波谷； 然后， 根据确定的所述波峰和所述波 谷， 确定所述波峰的倾斜度； 接着， 当所述波峰的倾斜度大于第一阔值， 则 确定所述波峰为杂波波峰； 最后， 通过消去所述杂波的波瓣来消除所述目标 频谱的杂波， 例如将所述杂波的波瓣位于两个波谷之间的频 点能量以在两个 波谷之间的内插来代替。

这里所述的波峰的倾斜度可以是所述波峰的倾 斜度是所述波峰两侧的两 个最邻近的波谷之间的连线与水平方向的夹角 ； 也可以是所述波峰两侧的两 个最邻近的波谷之间的连线的斜率； 也可以是所述波峰两侧的两个最邻近的 波谷之间的连线的夹角的三角函数。 第一阔值可以为预设的固定值， 也可以 与所述目标频谱的滤波阶数相关。

14, 针利用所述目标频谱的谐波的结构对所述谐波 进行谐波分析。

例如， 计算谐波的个数， 计算谐波的基频， 计算谐波的强弱， 计算谐波 的帧内 /帧间相关性， 计算谐噪比等等。

由此可见， 本发明实施例通过滤波确定谐波分析的目标频 谱， 并且能够 有效地消除杂波干扰， 从而获得清晰的和准确的谐波频谱， 以利于进行谐波 分析。

下面将结合具体实施例， 详细说明根据本发明实施例的谐波分析的方法 。

第一具体实施例， 以音频编码器为例。 一般地， 音频领域常见的时频域 变换包括 FFT, DCT, 改进型的离散余弦变换（MDCT, Modified Discrete Cosine Transform )等。

由于 MDCT是广泛应用于音频编码器的时频域变换工� �， 复用 MDCT 系数可以节省谐波分析的计算复杂度， 因此本具体实施例以 MDCT 变换为 例， 说明一个计算 32kHz釆样音频信号在 8k-16kHz频段上的谐波个数的谐 波分析方案， 输入音频帧为 20ms帧长。

第一步， 首先对输入编码器的时域音频帧 s(x)进行 MDCT变换， 得到输 入的时域音频帧的 MDCT系数 M(x), 其中 x=0,l, ...639。 通常， 时域音频帧 在进行 MDCT变换前都要先进行加窗操作，由于是公知� �术，这里不再详述， 默认 s(x)是加窗后的时域音频帧。

容易理解， 时域音频帧 s(x)的釆样个数取决于该时域音频帧的帧长和� � 样率。 在本具体实施例中， 根据 20ms帧长和 32Hz釆样率， 可以确定时域音 频帧 s(x)的釆样个数为 640。 同样容易理解， 时频变换后的频域音频帧的样 点个数取决于时域帧釆样个数和时频变换方法 。 在本具体实施例中， 根据

MDCT的变换方法， 可以确定频域音频帧 M(x)的样点个数也为 640。

然后， 根据输入的时域音频帧的 MDCT系数 M(x)得到输入的时域音频 帧在 8k-16kHz频段上的对数能量谱 S(i), 其中 i=0,l,...319。 为了方便后续描 述， 简称为频谱 S(i), 频谱 S(i)可以根据以下公式 1进行确定：

S(i) = log ₂ (M(i + 320) -M(i + 320)) 公式（1) 其中， i=0, 1, ...319。

第二步，对频谱 S (i)做 3阶滤波，得到滤波后的频谱 S ^[k] (j),这里 j=0, 1, ... 319, 其中 S ^[k )表示 k阶滤波后的频谱。 具体的滤波方法为以下公式（2)至 公式（4) 中之一或其组合：

∑ S(i) , j<(k-l)/2 公式（ 2 ) j + (k + l)/2

M j+{klc--\[))I/22

¹ X S(i) , (k-l)/2< j<319-(k-l)/2 公式（ 3 ) i-3-(k-\)/2

319

zrjz ∑ )， j>319-(k-l)/2 公式（ 4 ) 其中 k为滤波器的阶数， 例如 k=3。 实际得到的谐波频谱由于多种因素 的干扰可能是 "模糊" 的， 也就是， 除谐波外， 频谱上还会有很多干扰杂波。 一般而言， 相对于谐波， 干扰杂波可以被看作是一种更高频的成份。 通过滤 波， 可以有效地将部分杂波滤掉， 而保留住谐波成份。

第三步， 求频谱 S ^[3] (j)中波峰的个数 p ^[3] 。 当频谱中某一频点的能量高 于它的两个最邻近频点的能量时， 则认为该频点处是频谱的一个波峰。 此步 执行完后， 记录下求得的波峰个数1 ^[3] , 这里说的记录可以由软件实现， 例如 读取后存储在内容中， 也可以由硬件实现， 比如存储在寄存器中。 返回第二 步， 对频谱 S(i)再做 5阶滤波， 对 S ^[5] (j)再记录下其波峰个数 p ^[5] , 以此循 环往复再记录下其余的 P ^m 、 P ^[9] 等。

第四步， 在所有记录下的 p ^[k] 中找到最小值 ρ ^[τ] , 获取其对应的滤波后的 频谱， 即 S ^[T] (j), 下面称作目标频谱。 目标频谱 S ^[T] (j)的获取可以是从之前 在求 P ^[kI 的循环中记录下的频谱中直接获取，也可 以是根据 T值重新 由于每帧谐波频谱的特征不同， 对每帧使用相同的滤波方法不会在每帧都达 到最优的效果。 一般而言， 更优的滤波应该导致目标频谱具有更少的频谱 波 峰个数。 因此， 上述循环滤波寻找最小频谱波峰个数的过程， 就是一个寻找 最优滤波、 确定目标频谱的过程。

第五步， 搜索目标频谱 S ^[T] (j)的所有波峰和波谷， 记录下每个波峰和波 谷的位置和能量。 其中波谷的定义为两个相邻波峰间的能量最小 的频点。 特 峰与它所在端点之间的能量最小的频点。 通常， 波峰和波谷总是交替出现的。 波峰的能量和位置分别记作 E _p (m)和 Idx _p (m) , 其中 m=0, 1, ... L_l , L为目标频 谱的波峰个数 p ^m 。波谷的能量和位置分别记作 (1 )和 Idx _v (n) ,其中 n=0, 1 , ... L。 波谷个数比波峰个数多一个。

第六步， 依据波瓣宽度消除目标频谱的杂波。

通过滤波获取目标频谱的过程已经一定程度上 消除了部分杂波， 但一般 而言仍会有杂波未被消除。 此时， 由于杂波的存在， 目标频谱的波峰个数会 大于实际谐波的个数， 从而使目标频谱的平均波瓣宽度会小于谐波波 瓣宽度。 这时， 如果目标频谱的某个波峰的波瓣宽度小于平均 波瓣宽度， 则该波峰被 认为是杂波。 这里， 平均波瓣宽度的定义为： 目标频谱的带宽与波峰个数的 比值。

搜索出目标频谱 S ^[T] (j)的所有波峰和波谷后， 再计算每个波瓣宽度 Wp(m), 其中 m=0, 1, ... L-l。 其中， 波瓣的定义为由一个波峰两侧的两个最邻 近波谷之间的频点所组成的集合。 波瓣宽度的定义为该波瓣的两个波谷之间 的宽度， 即波瓣覆盖的频率范围大小。 波瓣宽度可以通过以下公式进行确定：

W _p {m) = Idx _v (m + \) - Idx _v (m) 公式（ 5 ) 被确认的杂波需要被从目标频谱中消去， 具体的消去方法是将杂波的波 瓣消去， 即将杂波波瓣位于两个波谷之间的频点能量以 在两个波谷之间的内 插来代替， 参见图 2。 图 1中位于中间的较小波峰即为杂波波峰， 经杂波消去 后， 该波峰即被虚线的频谱所代替。

其数 达为： - S ^[T] (Idx _v (m)) 公式 ( 6 ) 其中， 1=0, 1, ...Wp (m) , S ^[T] 'G)表示杂波波瓣被消除后的所在位置的频 谱。 为 了方便， 消去杂波后的目标频谱仍称作目标频谱。

第七步， 依据波峰的倾斜度消除目标频谱的杂波。

经第六步后的目标频谱可能仍是包含有杂波的 ， 还需利用其它方法进一 步消除。 图 3示出了一种典型的杂波频谱。 波峰 A是谐波波峰， 波峰 B和波 峰 C是杂波波峰。 虚线 M是穿过波峰 B的两个波谷的直线， 虚线 N是穿过波 峰 C的两个波谷的直线， 角 α和角 β分别是 M和 N与水平轴的夹角。 可见， 波峰 Β和波峰 C明显倾斜， 具体而言， 角 α和角 β的角度明显较大， 即角 α 和角 β的正切较大， 或直线 Μ和直线 Ν的斜率较大。

通过判断穿过某波峰的两个波谷的直线的斜率 ， 可以一定程度上确定该 波峰是否为杂波。 具体的， 将求得的斜率与一个门限 thrl做比较， 若大于该 门限， 则认为该波峰为杂波。 门限 thrl可以是一个预先设定的固定值， 也可 以是变量。 当门限是变量时， 可以根据第四步中最优滤波的阶数 T确定， 即 thrl=f (T) , 表示 thrl是滤波器的最优阶数 T的函数。 由于第六步的消除杂 波操作已经改变了目标频谱， 在依据波峰的倾斜度消除目标频谱的杂波之前 ， 需要重新搜索目标频谱的所有波峰和波谷。 确定出杂波后， 以与第六步中公 式（6 )表示的相同的方法消除， 参照图 2所示的方法将杂波消除。

第八步， 搜索谐波个数。 经过之前的步骤， 目标频谱的杂波已经在很大程度 上消除了， 谐波结构已变得清晰， 可以进行最后的谐波个数确认了。 重新搜 索目标频谱的所有波峰和波谷，计算每个波峰 的强度 I _P (y),其中 y=0,l,2...K, K为波峰个数。 波峰强度通过以下公式进行确定： ， ¹ ) 公式（7 ) 将每个波峰的强度与一个门限 thr2做比较。 若波峰的强度大于该门限， 则认为检测到一个谐波。 与门限 thrl类似， 门限 thr2可以是一个预先设定的 固定值， 也可以是变量。 当门限是变量时， 可以根据第四步中最优滤波的阶 数 T确定， 即 thr2=f(T), 表示 thr2是滤波器的最优阶数 T的函数。

第二具体实施例中， 其与第一具体实施例中的大部分步骤相似， 不同在 于对第一实施例中确定目标频谱的第二至四步 的简化处理。 具体而言， 第一具体实施例的第二至四步中的最优滤波器 的确定是通过 一个循环过程并取其中的最优滤波器来实现的 。 一般而言， 随着滤波的阶数 由低向高增加， 每轮滤波后的频谱波峰个数会呈现出一种先单 调递减， 再单 调递增的现象。 通常， 最优滤波就出现在波峰个数折线由递减到递增 的拐点 处。 为了降低计算复杂度， 第二具体实施例中对第一具体实施例中的循环 过 程增加了一个提前终止判断， 即一旦检测到拐点， 则认为已找到最优滤波， 不再执行剩下的循环。 每循环一次， 比较当次循环后的目标频谱的波峰个数 是否小于等于前次循环后的目标频谱的波峰个 数。 若是， 则继续循环， 若否， 则终止循环， 以前次循环后的目标频谱作为后续步骤的目标 频谱， 即将前次 循环的滤波作为最优滤波。

第三具体实施例中， 其与第一具体实施例中的大部分步骤相似， 不同在 于对第一实施例中确定目标频谱的第二至四步 的简化处理。

具体而言， 滤波阶数不再动态的确定， 而是对每帧都以一个固定的阶数 进行滤波。 经过对大的样本进行统计可以发现， 对一给定的釆样率和帧长， 最优滤波阶数在某一固定阶数的概率会显著高 于其它阶数。 因此以该最大出 现概率阶数作为固定的滤波阶数可以在节省计 算复杂度的同时最大限度的逼 近最优滤波结果。 以第一具体实施例中 32kHz釆样， 20ms帧长的输入信号 为例， 釆用的固定滤波阶数选为 5。

第四具体实施例中， 其与第一具体实施例中的大部分步骤相似， 不同在 于对第一实施例中确定目标频谱的第二至四步 的简化处理。

本发明实施例的输入信号可以不限制为音频信 号， 可以是任意一种可以 在频域分析的信号， 如图像信号等。

另外， 对原始频谱滤波获得目标频谱的过程所使用的 滤波器可以是任意 形式的低通滤波器。 例如， 以下公式 ( 8 )给出了一阶低通滤波器的一个实施 例。

^ ^] ( ) = ^( -1) + ^( 公式（8 ) 其中， i=0,l, ... ,319, /是常系数。

第五具体实施例中， 其与第一具体实施例中的大部分步骤相似， 不同在 于对第一实施例中确定目标频谱的第六步的依 据波瓣宽度消除杂波过程中波 瓣宽度的确定方法。

具体而言， 第一具体实施例在第六步依据波瓣宽度消除杂 波的步骤中， 杂波的确定是将某个波峰的波瓣宽度与平均波 瓣宽度做比较。 然而， 根据特 定的应用， 确定杂波时也可以将波瓣宽度与一个预先设定 的固定值做比较。 如要搜索 8000Hz〜 16000Hz频带内的非致密谐波（即在该频带内谐� �数少于 10的谐波）， 若当谐波数等于 10时，平均波瓣宽度应为 800Hz或 32个频点， 故非致密谐波的波瓣宽度应大于 32。所以可将门限设定为一个固定值 32,若 某波峰的波瓣宽度小于 32, 则认为是杂波。

下面将结合图 4和图 5分别说明根据本发明实施例的确定频谱的谐� �间 杂波的方法。

如图 4所示的确定频谱的谐波间杂波的方法包括如� �步骤：

41 , 确定频语的全部波峰和波谷；

42, 根据确定的所述波峰和所述波谷， 确定所述频谱中的波瓣宽度， 其 中所述波瓣宽度是波峰两侧的两个最邻近的波 谷之间的宽度；

43 , 当所述波瓣宽度小于第一波瓣宽度， 则确定所述波瓣宽度对应的波 峰为杂波的波峰。

其中， 第一波瓣宽度可以为预设的固定值。 或者， 第一波瓣宽度可以为 平均波瓣宽度， 其中所述平均波瓣宽度是所述频谱的宽度与所 述全部波峰的 个数的比值。

由此可见， 本发明实施例的确定频谱的谐波间杂波的方法 可以有效去除 频谱的谐波间杂波。

如图 5所示的确定频谱的谐波间杂波的方法包括如� �步骤：

51 , 确定所述目标频谱的全部波峰和波谷；

52, 根据确定的所述波峰和所述波谷， 确定所述波峰的倾斜度；

53 , 当所述波峰的倾斜度大于第一阔值，则确定所 述波峰为杂波的波峰。 其中， 波峰的倾斜度是所述波峰两侧的两个最邻近的 波谷之间的连线与 水平方向的夹角； 或者所述波峰的倾斜度是；

其中， 第一阔值可以为预设的固定值。

由此可见， 本发明实施例的确定频谱的谐波间杂波的方法 可以有效去除 频谱的谐波间杂波。

此外， 图 4与图 5所示的两种确定频谱的谐波间杂波的方法也� �以结合 实施， 这样可以更有效地去除频谱的谐波间杂波。

图 6示出了根据本发明实施例的谐波分析的装置� �结构示意图。 其中， 谐波分析的装置 60 包括时频转换单元 61、 滤波单元 62、 杂波消除单元 63 和谐波分析单元 64。 所述时频转换单元 61、 滤波单元 62、 杂波消除单元 63 和谐波分析单元 64 所对应的功能模块集成在数字信号处理（DSP, Digital Signal Processing ) 器件中。

其中， 时频转换单元 61用于对输入的时域信号进行时频转换， 以获得所 述时域信号在频域上的频谱。 滤波单元 62用于对所述时频转换单元 61转换 得到的所述频谱进行滤波处理， 以确定目标频谱。杂波消除单元 63用于对所 述滤波单元 62确定的所述目标频谱的谐波间杂波做消除处� ��，以便确定所述 目标频谱的谐波的结构。谐波分析单元 64用于利用所述目标频谱的谐波的结 构对由所述杂波消除单元 63消除了杂波的所述谐波进行谐波分析。

进一步， 滤波单元 62具体用于： 对所述频谱进行 k次滤波处理， 并在所 述频谱经过所述 k次滤波处理的每次滤波处理之后记录波峰个� �以获得 k个 波峰个数， 其中 k为大于或等于 2的正整数； 确定所述 k个波峰个数中的最 小波峰个数； 确定得到所述最小波峰个数的滤波处理的频谱 为目标频谱。

或者， 滤波单元 62用于： 对所述频谱进行第一滤波处理， 并记录经过所 述第一滤波处理的频谱的第一波峰个数； 对所述频谱进行第二滤波处理， 并 记录经过所述第二滤波处理的频谱的第二波峰 个数； 对所述频谱进行第三滤 波处理， 并记录经过所述第三滤波处理的频谱的第三波 峰个数； 当所述第二 波峰个数小于所述第一波峰个数， 且所述第二波峰个数小于所述第三波峰个 数， 则确定经过所述第二滤波处理得到的频谱为目 标频谱。

具体而言， 若所述第一滤波处理为第一阶滤波、 所述第二滤波处理为第 二阶滤波， 以及所述第三滤波处理为第三阶滤波， 滤波单元 62用于： 对所述 频谱进行第一阶滤波， 并记录经过所述第一阶滤波的频谱的第一波峰 个数； 对所述频谱进行第二阶滤波， 并记录经过所述第二阶滤波的频谱的第二波峰 个数； 对所述频谱进行第三阶滤波， 并记录经过所述第三阶滤波的频谱的第 三波峰个数； 当所述第二波峰个数小于所述第一波峰个数， 且所述第二波峰 个数小于所述第三波峰个数，则确定经过所述 第二阶滤波的频谱为目标频谱。

或者，滤波单元 62用于： 釆用任意形式的低通滤波器对所述频谱进行低 通滤波处理， 将所述经过所述低通滤波处理后的频语作为目 标频谱。

或者， 滤波单元 62用于： 对所述频谱进行第四阶滤波， 其中所述第四阶 为固定值； 确定经过所述第四阶滤波的频谱为目标频谱。

进一步地，杂波消除单元 63釆用以下方式中的至少一种：依据波瓣宽度 对所述目标频谱的谐波间杂波做消除处理， 其中所述波瓣宽度是波峰两侧的 两个最邻近的波谷之间的宽度； 依据波峰的倾斜度对所述目标频谱的谐波间 杂波做消除处理， 其中所述波峰的倾斜度是所述波峰两侧的两个 最邻近的波 谷之间的连线与水平方向的夹角， 或者所述波峰两侧的两个最邻近的波谷之 间的连线的斜率， 或者所述波峰两侧的两个最邻近的波谷之间的 连线的夹角 的三角函数。

具体而言，杂波消除单元 63用于：确定所述目标频谱的全部波峰和波谷� �� 根据确定的所述波峰和所述波谷， 确定所述目标频谱中的波瓣宽度； 当所述 波瓣宽度小于第一波瓣宽度，则确定所述波瓣 宽度对应的波峰为杂波的波峰； 为预设的固定值， 或者所述第一波瓣宽度为平均波瓣宽度， 其中所述平均波 瓣宽度是所述目标频谱的宽度与所述第二波峰 个数的比值。

或者， 具体地， 杂波消除单元 63用于： 确定所述目标频谱的全部波峰和 波谷； 根据确定的所述波峰和所述波谷， 确定所述波峰的倾斜度； 当所述波 峰的倾斜度大于第一阔值， 则确定所述波峰为杂波波峰； 通过消去所述杂波 的波瓣来消除所述目标频谱的杂波， 例如， 将所述杂波的波瓣位于两个波谷 之间的频点幅值以在两个波谷幅值之间的内插 来代替。 其中， 第一阔值可以 为预设的固定值。

由此可见， 本发明实施例通过滤波确定谐波分析的目标频 谱， 并且能够 有效地消除杂波干扰， 从而获得清晰的和准确的谐波频谱， 以利于进行谐波 分析。

图 7示出了根据本发明实施例的确定频谱的谐波� �杂波的装置。 如图所 示， 确定频谱的谐波间杂波的装置 70包括第一确定单元 71、 第二确定单元 72和第三确定单元 73。

其中，第一确定单元 71用于确定所述频谱的全部波峰和波谷。第二� ��定 单元 72用于根据所述第一确定单元 71确定的所述波峰和所述波谷， 确定所 述频谱中的波瓣宽度， 其中所述波瓣宽度是波峰两侧的两个最邻近的 波谷之 间的宽度。第三确定单元 73用于当所述波瓣宽度小于第一波瓣宽度，则� ��定 由所述第二确定单元 72确定的所述波瓣宽度对应的波峰为杂波的波� ��。

这里， 第一波瓣宽度可以为预设的固定值。 或者， 第一波瓣宽度可以为 平均波瓣宽度， 其中所述平均波瓣宽度是所述频谱的宽度与所 述全部波峰的 个数的比值。

由此可见， 本发明实施例可以通过波瓣宽度与预设值的比 较结果确定谐 波间的杂波， 以便再通过消除杂波来获得清楚的谐波结构。

图 8示出了根据本发明实施例的确定频谱的谐波� �杂波的装置。 如图所 示， 确定频谱的谐波间杂波的装置 80包括第四确定单元 81、 第五确定单元 82和第六确定单元 83。

其中，第四确定单元 81用于确定所述目标频谱的全部波峰和波谷。� ��五 确定单元 82用于根据所述第四确定单元 81确定的所述波峰和所述波谷， 确 定所述波峰的倾斜度。第六确定单元 83用于当所述波峰的倾斜度大于第一阔 值， 则确定由所述第五确定单元 82确定的所述波峰为杂波的波峰。

这里， 所述波峰的倾斜度是所述波峰两侧的两个最邻 近的波谷之间的连 线与水平方向的夹角； 或者所述波峰的倾斜度是所述波峰两侧的两个 最邻近 的波谷之间的连线的斜率： 或者所述波峰的倾斜度是所述波峰两侧的两个 最 邻近的波谷之间的连线的夹角的三角函数。 第一阔值可以为预设的固定值。 定谐波间的杂波， 以便再通过消除杂波来获得清楚的谐波结构。

容易理解， 可以将上两者相结合， 即通过确定频谱的谐波间杂波的装置 70确定杂波并消除杂波， 然后再通过确定频谱的谐波间杂波的装置 80进一 步确定杂波并消除杂波， 从而可以更有效地消除杂波。 也可以将确定频谱的 谐波间杂波的装置 70与确定频谱的谐波间杂波的装置 80集成在一起形成一 个确定频谱的谐波间杂波的装置。 这样， 第一确定单元 71 与第四确定单元 81的功能相同， 可以去除一个。 于是， 集成的确定频谱的谐波间杂波的装置 可以包括第一确定单元 71、 第二确定单元 72、 第三确定单元 73、 第五确定 单元 82和第六确定单元 83。

图 9示出了才艮据本发明实施例的编码器 90的结构示意图。其中，该编码 器 90包括处理器 91和存储器 92。

其中， 处理器 91实现根据本发明实施例的谐波分析的方法。 即， 处理器 91用于对输入的时域信号进行时频转换， 以获得所述时域信号在频域上的频 谱； 对所述频谱进行滤波处理， 以确定目标频谱； 对所述目标频谱的谐波间 杂波做消除处理， 以便确定所述目标频谱的谐波的结构； 针利用所述目标频 谱的谐波的结构对所述谐波进行谐波分析。 存储器 92用于存储处理器 91执 行的指令。

应理解， 本发明的每个权利要求所叙述的方案也应看做 是一个实施例， 并且是权利要求中的特征是可以结合的， 如本发明中的判断步骤后的执行的 不同分支的步骤可以作为不同的实施例。

本领域普通技术人员可以意识到， 结合本文中所公开的实施例描述的各 示例的单元及算法步骤， 能够以电子硬件、 或者计算机软件和电子硬件的结 合来实现。 这些功能究竟以硬件还是软件方式来执行， 取决于技术方案的特 定应用和设计约束条件。 专业技术人员可以对每个特定的应用来使用不 同方 法来实现所描述的功能， 但是这种实现不应认为超出本发明的范围。

所属领域的技术人员可以清楚地了解到， 为描述的方便和简洁， 上述描 述的系统、 装置和单元的具体工作过程， 可以参考前述方法实施例中的对应 过程， 在此不再赘述。

在本申请所提供的几个实施例中， 应该理解到， 所揭露的系统、 装置和 方法， 可以通过其它的方式实现。 例如， 以上所描述的装置实施例仅仅是示 意性的， 例如， 所述单元的划分， 仅仅为一种逻辑功能划分， 实际实现时可 以有另外的划分方式， 例如多个单元或组件可以结合或者可以集成到 另一个 系统， 或一些特征可以忽略， 或不执行。 另一点， 所显示或讨论的相互之间 的耦合或直接耦合或通信连接可以是通过一些 接口， 装置或单元的间接耦合 或通信连接， 可以是电性， 机械或其它的形式。 为单元显示的部件可以是或者也可以不是物理 单元， 即可以位于一个地方， 或者也可以分布到多个网络单元上。 可以根据实际的需要选择其中的部分或 者全部单元来实现本实施例方案的目的。

另外，在本发明各个实施例中的各功能单元可 以集成在一个处理单元中， 也可以是各个单元单独物理存在， 也可以两个或两个以上单元集成在一个单 元中。

所述功能如果以软件功能单元的形式实现并作 为独立的产品销售或使用 时， 可以存储在一个计算机可读取存储介质中。 基于这样的理解， 本发明的 技术方案本质上或者说对现有技术做出贡献的 部分或者该技术方案的部分可 以以软件产品的形式体现出来， 该计算机软件产品存储在一个存储介质中， 包括若干指令用以使得一台计算机设备（可以 是个人计算机， 服务器， 或者 网络设备等）执行本发明各个实施例所述方法 的全部或部分步骤。 而前述的 存储介质包括： U盘、 移动硬盘、 只读存储器（ROM, Read-Only Memory ), 随机存取存储器 （RAM, Random Access Memory )、 磁碟或者光盘等各种可 以存储程序代码的介质。

以上所述， 仅为本发明的具体实施方式， 但本发明的保护范围并不局限 于此， 任何熟悉本技术领域的技术人员在本发明揭露 的技术范围内， 可轻易 想到变化或替换， 都应涵盖在本发明的保护范围之内。 因此， 本发明的保护 范围应所述以权利要求的保护范围为准。