本发明涉及语音检测的技术领域，更具体地，涉及一种基于自回归模型系数的合成语音检测方法。
背景技术：
自动扬声器验证(asv)被部署在越来越多的不同应用和服务中，例如移动电话、智能扬声器和呼叫中心，以便为个人身份验证提供一种低成本和灵活的生物识别解决方案。尽管asv系统的性能在近年来逐渐提高，但asv系统易受欺骗攻击。其中，在人们特别关注的音频篡改方面，主要存在两种欺骗攻击：语音合成攻击和语音转换攻击，这两种欺骗攻击对asv系统构成重大威胁。语音合成(texttospeech),简称tts，是将文字转化为语音的一种技术，类似于人类的嘴巴，通过不同的音色说出想表达的内容，tts系统可以生成完全人工的语音信号。语音转换(voiceconversion)，简称vc，通过输入一条语音，在保持说话内容不变的情况下，让它听起来像是另一个人说的，vc系统在自然语音上运行。ss和vc技术都可以产生高质量的语音信号，模仿特定目标的语音。
目前，在合成语音检测方面的研究大多须先从语音里提取特征，常用的语音特征包括mfcc、cqcc、spec。其中，mfcc是语音处理中最常用的基于量的特征之一，在mel尺度上对对数谱进行倒谱分析，适用于区分篡改语音和人类语音；cqcc是一种基于振幅的特征，它使用常数q变换(cqt)与传统的倒频谱分析相结合；相对于mfcc、cqcc来说，spec特征是更原始的特征，因为它通过计算汉明窗口上的stft值，然后计算每个部分的大小来得到。
2018年11月16日，中国专利(cn108831506a)中公开了一种基于gmm-bic的数字音频篡改点检测方法及系统，也属于语音检测的技术领域，此专利中提出的方法就是在分割出语音信号中的静音帧之后提取静音帧的mfcc特征，然后使用gmm-bic方法代替了传统的sgm-bic用于数字音频篡改点检测，数字音频篡改定位的自动化，自适应性好，鲁棒性高，也能保证检测准确率，但实际上，mfcc语音特征在语音检测过程中的计算量大，因此，此技术方案无法同时兼顾语音检测的检测效率和检测准确率。
自回归模型(autoregressivemodel，简称ar)是最常见的平稳时间序列模型之一，在统计上是一种处理时间序列的方法，ar被广泛运用在经济学、信息学、自然现象的预测上，主要对一维数据做时间序列分析。语音也
属于一维数据，可以通过ar线性预测模型评估语音序列之间的关系，因此，研究如何基于ar系数进行语音篡改检测，具有重要意义。
技术实现要素：
为解决现有直接应用已有的语音特征提取算法在进行语音检测时，无法同时兼顾语音检测的检测效率和检测准确率的问题，本发明提出一种基于自回归模型系数的合成语音检测方法，减少检测过程中的计算量，提高检测准确率。
为了达到上述技术效果，本发明的技术方案如下：
一种基于自回归模型系数的合成语音检测方法，至少包括：
s1.将已知数据库的训练集、验证集和测试集的语音片段固定为统一长度a；
s2.将训练集、验证集和测试集的语音均进行分段，提取分段后每段语音的不同阶数的ar系数，再对提取出的ar系数进行排列，组成二维ar语音特征，并构成二维ar语音特征的训练特征集、验证特征集和测试特征集；
s3.利用二维ar语音特征的训练特征集及验证特征集训练一个卷积神经网络分类器，保存卷积神经网络分类器的最优参数；
s4.利用训练好的卷积神经网络分类器对测试特征集进行分类，确认语音是否经过语音合成或语音转换篡改操作；
s5.对不同阶数的二维ar语音特征进行融合，确认语音是否经过语音合成或语音转换篡改操作。
在本技术方案中，已知数据库的训练集、验证集和测试集的语音分段后的每一帧语音信号均提取了ar系数，排列后的ar系数组成二维ar语音特征，然后训练卷积神经网络分类器进行分类，即可确认语音是否经过语音合成或语音转换篡改操作，不直接应用已有的语音特征提取算法，减少检测过程中的计算量，最后通过融合不同阶数的二维ar特征进一步提高可检测精度，兼顾了语音检测的检测效率和检测准确率。
优选地，步骤s1所述将已知数据库的训练集、验证集和测试集的语音片段固定为统一长度a的过程为：
s101.选定已知数据库，分别将已知数据库的训练集、验证集和测试集的语音片段长度固定为a个采样点；
s102.判断训练集、验证集和测试集中的任意一个在固定前的语音长度是否大于等于a，若是，进行截断操作，将语音长度固定为a；否则，通过复制语音对固定前的语音长度进行扩充，再固定语音长度为a。
优选地，步骤s2所述的构成二维ar语音特征的训练特征集、验证特征集和测试特征集的具体过程为：
s201.将训练集、验证集和测试集固定统一长度为a的语音片段均划分为b段；
s202.提取分段后每段语音的阶数为h的ar系数，将提取出的b维ar系数进行排列，组成二维ar语音特征；
s203.分别构成维度为b×h的二维ar语音特征的训练特征集、验证特征集和测试特征集。
优选地，步骤s202中h满足：
8≤h≤150，其中，h表示正整数，保证ar系数在应用时的效果。
优选地，二维ar语音特征的训练特征集中包括原始语音特征和篡改语音特征，步骤s3中，所述的卷积神经网络分类器通过梯度下降法训练得到；将验证特征集准确度指标最好时对应的卷积神经网络分类器的参数作为最优参数，并保存。
优选地，步骤s4中所述利用训练好的卷积神经网络分类器对测试特征集进行分类，确认语音是否经过语音合成或语音转换篡改操作的具体过程为：
s401.将维度为b×h的二维ar语音特征的测试特征集输入训练好的卷积神经网络分类器；
s402.根据语音分数cm的计算公式：
cm(f)＝log(p(bonafide|f；θ))-log(p(spoof|f；θ))
计算每条语音的分数，并保存分数文件，其中，f表示当前送入的语音特征，θ表示保存的最优参数；p(bonafide|f；θ)表示送入特征f为原始语音bonafide特征的概率，p(spoof|f；θ)送入特征f为合成篡改语音spoof特征的概率；
s403.判断cm(f)的值是否大于等于判定阈值t，若是，送入的语音特征为原始语音，否则，送入的语音特征为合成语音或转换篡改语音；
s403.根据保存的分数文件及已知数据库给定的分数文件，进一步确定当前测试特征集的t-dcf指标和等错误概率eer指标。
在此，cm(f)分数的值代表待测语音信号为原始语音的可能性，同时在步骤s403中引入t-dcf指标和等错误概率eer指标进一步评判二维ar语音特征的有效性。
优选地，步骤s5所述对不同阶数的二维ar语音特征进行融合，确认语音是否经过语音合成或语音转换篡改操作的过程为：
s501.将不同阶数的二维ar语音特征的测试特征集分别输入至训练好的卷积神经网络分类器，得到不同的分数文件并保存；
s502.将不同分数文件进行平均，得到融合后的分数文件；
s503.判断融合后的分数值是否大于等于判定阈值t，若是，融合后的语音特征为原始语音，否则，融合后的语音特征为合成语音或转换篡改语音；
s503.根据融合后的分数文件及已知数据库给定的分数文件，进一步确定融合后的t-dcf指标和等错误概率eer指标。
优选地，cm(f)分数越大，当前输入的二维ar语音特征为原始语音的概率越大；cm(f)分数越小，当前输入的二维ar语音特征为原始语音的概率越小。
优选地，t-dcf指标和等错误概率eer指标越小，表示二维ar语音特征越有效。
优选地，已知数据库为asvspoof2019语音数据集。
与现有技术相比，本发明技术方案的有益效果是：
本发明提出一种基于自回归模型系数的合成语音检测方法，首先将已知数据库中训练集、验证集和测试集的语音片段固定为统一长度，分段后提取每段语音信号的ar系数，组成二维ar语音特征，从而构建训练特征集、验证特征集和测试特征集，利用训练特征集和验证特征集训练卷积神经网络分类器，再将测试特征集通过训练好的卷积神经网络分类器进行分类，确认待测语音信号是否经过语音合成或语音转换篡改操作，与现有技术相比，不直接应用已有的语音特征提取算法，减少检测过程中的计算量，最后通过融合不同阶数的二维ar特征进一步提高可检测精度，兼顾了语音检测的检测效率和检测准确率。
附图说明
图1表示本发明实施例中提出的基于自回归模型系数的合成语音检测方法的流程图。
具体实施方式
附图仅用于示例性说明，不能理解为对本专利的限制；
为了更好地说明本实施例，附图某些部位会有省略、放大或缩小，并不代表实际尺寸；
对于本领域技术人员来说，附图中某些公知内容说明可能省略是可以理解的。
附图中描述位置关系的用于仅用于示例性说明，不能理解为对本专利的限制；
下面结合附图和实施例对本发明的技术方案做进一步的说明。
实施例1
如图1所示的基于自回归模型系数的合成语音检测方法的流程图，参见图1，所述方法包括：
s1.将已知数据库的训练集、验证集和测试集的语音片段固定为统一长度a；在本实施例中，已知数据库选取asvspoof2019语音数据集，统一长度为64000，具体的固定过程为：
s101.选定asvspoof2019语音数据集作为已知数据库，分别将已知数据库的训练集、验证集和测试集的语音片段长度固定为64000个采样点；
s102.判断训练集、验证集和测试集中的任意一个在固定前的语音长度是否大于等于64000，若是，进行截断操作，将语音长度固定为64000；否则，通过复制语音对固定前的语音长度进行扩充，再固定语音长度为64000，即训练集、验证集和测试集中的任意一个在固定前的语音长度存在超过64000的情况时，直接通过现有技术进行截断操作，使得语音长度满足64000的统一长度要求，若训练集、验证集和测试集中的任意一个在固定前的语音长度不足64000，则将语音片段进行复制，对原有的语音长度进行扩充，使得语音长度满足64000的统一长度要求；
s2.将训练集、验证集和测试集的语音均进行分段，提取分段后每段语音的不同阶数的ar系数，再对提取出的ar系数进行排列，组成二维ar语音特征，并构成二维ar语音特征的训练特征集、验证特征集和测试特征集；
其中，ar系数的提取可以通过不同的成熟技术手段实现，其它具体过程为：
s201.将训练集、验证集和测试集固定统一长度为64000的语音片段均划分为400段；
s202.提取分段后每段语音的阶数为h的ar系数，将提取出的400维ar系数进行排列，组成二维ar语音特征；
在本实施例中，为保证ar系数在应用时的效果，阶数h满足：8≤h≤150，其中，h表示正整数，即h可以取端点8或150，也可以取8～150之间的任意一个正整数，划分的400段语音片段中每段含有160个采样点；二维ar语音特征的训练特征集中包括原始语音特征和篡改语音特征；
s203.按照每段语音的提取的阶数为h的ar系数与划分的400个片段的排列，分别构成维度为400×h的二维ar语音特征的训练特征集、验证特征集和测试特征集。
s3.利用二维ar语音特征的训练特征集及验证特征集训练一个卷积神经网络分类器，保存卷积神经网络分类器的最优参数；在本实施例中，卷积神经网络分类器通过梯度下降法训练得到，不局限于某一种具体的卷积神经网络分类器，然后将验证特征集准确度指标最好时对应的卷积神经网络分类器的参数作为最优参数，并保存。
s4.利用训练好的卷积神经网络分类器对测试特征集进行分类，确认语音是否经过语音合成或语音转换篡改操作；
具体过程为：
s401.将维度为400×h的二维ar语音特征的测试特征集输入训练好的卷积神经网络分类器；
s402.根据语音分数cm的计算公式：
cm(f)＝log(p(bonafide|f；θ))-log(p(spoof|f；θ))
计算每条语音的分数，并保存分数文件，其中，f表示当前送入的语音特征，θ表示保存的最优参数；p(bonafide|f；θ)表示送入特征f为原始语音bonafide特征的概率，p(spoof|f；θ)送入特征f为合成篡改语音spoof特征的概率；
s403.判断cm(f)的值是否大于等于判定阈值t，若是，送入的语音特征为原始语音，否则，送入的语音特征为合成语音或转换篡改语音；
s403.根据保存的分数文件及已知数据库给定的分数文件，进一步确定当前测试特征集的t-dcf指标和等错误概率eer指标，其中cm(f)分数的值代表待测语音信号为原始语音的可能性，同时在步骤s403中引入t-dcf指标和等错误概率eer指标进一步评判二维ar语音特征的有效性，其中，t-dcf(atandemdetectioncostfunction)是串联代价函数指标，由两个系统(asv，cm)决定，计算公式如下：
其中，
等错误概率eer是指误拒绝率(falserejectionrate，frr)等于误接受率(falseacceptancerate，far)的点，此时的far与frr的值称为等错误率。
s5.对不同阶数的二维ar语音特征进行融合，确认语音是否经过语音合成或语音转换篡改操作，过程为：
s501.将不同阶数的二维ar语音特征的测试特征集分别输入至训练好的卷积神经网络分类器，得到不同的分数文件并保存；
s502.将不同分数文件进行平均，得到融合后的分数文件；
s503.判断融合后的分数值是否大于等于判定阈值t，若是，融合后的语音特征为原始语音，否则，融合后的语音特征为合成语音或转换篡改语音；
s503.根据融合后的分数文件及已知数据库给定的分数文件，进一步确定融合后的t-dcf指标和等错误概率eer指标。
在本实施例中，cm(f)分数越大，当前输入的二维ar语音特征为原始语音的概率越大；cm(f)分数越小，当前输入的二维ar语音特征为原始语音的概率越小，t-dcf指标和等错误概率eer指标越小，表示二维ar语音特征越有效。
具体的综合实施的过程如下：
首先固定已知数据库中训练集、验证集和测试集的语音长度为64000，对语音进行分段处理，每段采样点为160，对每段语音进行ar系数的提取，
其中，ar系数在8～150阶之间，选取
8，10，20，30，40，50，60，70，80，90，100，110，120，130，140，150阶的ar系数进行实验；构成的训练特征集和验证特征集训练卷积神经网络分类器，保存最优参数，并用测试特征集进行测试，除此之外，利用t-dcf和eer两个指标来评价特征效果，整体来看，在这些阶数上的效果会有所波动，但是都能体现良好的语音篡改检测的效果，实验结果如表1所示。
表1
其中，development表示验证集上的结果项，evaluation表示测试集上的结果项，表1中实验结果表明在10阶[ar(10)]时特征的效果最好，在表1的最后一行，将10阶特征和50阶特征进行融合，更是进一步提高了特征效果。
附图中描述位置关系的用于仅用于示例性说明，不能理解为对本专利的限制；
显然，本发明的上述实施例仅是为清楚地说明本发明所作的举例，而并非是对本发明的实施方式的限定。对于所属领域的普通技术人员来说，在上述说明的基础上还可以做出其它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等均应包含在本发明权利要求的保护范围之内。
本文用于企业家、创业者技术爱好者查询，结果仅供参考。
