一种基于毫米波感知的非接触式声纹生物认证方法技术

本发明专利技术公开了一种基于毫米波感知的非接触式声纹生物认证方法。包括以下步骤：毫米波探针感知喉咙的震动，对用户的喉咙区域发射调频连续波，然后接收用户喉咙区域反射的带有声音震动信息的调频连续波；毫米波探针将接收的调频连续波的周期性线性调频信号转化为中频信号，根据中频信号计算声音的震动位移信息；去除震动位移信息的背景杂波和动态杂波；从距离多普勒矩阵中提取声源特征和声轨特征；将提取的声源特征和声轨特征输入细粒度认证模型中得到用户的分类结果。本发明专利技术基于低成本、便携式、低能耗的硬件系统实现了安全非接触式的声纹认证系统，不但可以从本质上抵御重放攻击，而且可以方便的迁移到智能手机等移动设备的认证系统中。

【技术实现步骤摘要】
一种基于毫米波感知的非接触式声纹生物认证方法
本专利技术涉及生物认证和无线
，尤其涉及一种基于毫米波(millimeterwave，mmWave)感知的非接触式声纹生物认证方法。
技术介绍
声纹是一种混合生理特征和行为特征的生物特征，像指纹一样被认为是每个人都独一无二的生物特征。由于声控设备和服务的发展，基于声音的用户认证技术流行起来。亚马逊的Echo和GoogleHome等商用声控设备都配备了语音认证功能来保证用户的信息安全。然而基于声纹的生物认证技术有很大的隐患。声纹认证是基于无向感知的方式，由于这个本质的弱点，声纹认证极易遭受重放攻击的威胁。
技术实现思路
为了解决
技术介绍
中的问题，本专利技术提供了一种基于毫米波感知的非接触式声纹生物认证方法，通过对毫米波探针收集到的调频连续波(FMCW)进行特征提取和分类后实现用户身份的认证本专利技术采用的技术方案：一种基于毫米波感知的非接触式声纹生物认证方法包括以下步骤步骤一：声音震动的毫米波感知：使用一个高精度的毫米波探针感知喉咙的震动，这种感知是非接触式的，对用户的喉咙区域发射调频连续波(FMCW)，然后接收用户喉咙区域反射的带有声音震动信息的调频连续波；毫米波探针将接收的调频连续波的周期性线性调频信号转化为中频(IF)信号，通过计算得到多个线性调频周期的中频信号快速傅里叶变换后的结果，并计算声音的震动位移信息；步骤二：杂波抑制处理：将步骤一中多个线性调频周期的中频信号快速傅里叶变换后的结果进行第二次快速傅里叶变换得到距离多普勒矩阵，通过更新距离多普勒矩阵完成背景杂波的分离和动态杂波的去除；步骤三：特征提取：从步骤二去除杂波后的距离多普勒矩阵中提取声源特征和声轨特征；步骤四：用户识别：构建细粒度认证模型，将步骤三中提取的声源特征和声轨特征输入细粒度认证模型中得到用户的分类结果。所述步骤一中的根据中频信号计算声音震动位移的具体过程为：1.1)设一个线性调频周期为Tr，根据中频信号获得第m个线性调频周期d(mTr)的震动位移的计算公式为；其中，m为线性调频周期的编号，m＝0，1，2...，M-1；M为线性调频周期的数量；f0是载波频率；Δψm是对第m个线性调频周期的中频信号进行快速傅里叶变换(FFT)后的结果；最终得到M个线性调频周期的震动位移。Δψm包含有第m个线性调频周期的距离信息。1.2)通过平移变化对M个线性调频周期的震动位移进行身体运动补偿，即对线性调频周期的距离信息进行校准：身体的随机运动是影响获取准确的声音震动数据的主要障碍，当人体运动幅度大于距离轮廓分辨率的一半时，将很难准确的计算出实际的中频信号傅里叶变换结果Δψm，我们使用一种细粒度的距离分箱对准方法来解决这个问题，傅里叶变换后的信号包含了物体的距离信息，细粒度的距离分箱对准方法通过对原始的距离信息进行平移来消除身体运动补偿。定义Sm(l)为获取的第m个线性调频周期的距离信息，的值代表了第m个线性调频周期校准后的距离信息；其中，l＝0，1，2，..，L-1，L为第m个线性调频周期内采样点的数量，l为采样点的编号；定义第m个线性调频周期的相对距离信息为Qm(l)，具体为：其中，Qm-1(l)为第m-1个线性调频周期的相对距离信息；所述相对距离信息Qm(l)为校准后的距离信息与步骤1.1)获取的距离信息Sm(l)之间的相对距离。定义Xm为Sm(l)的平移距离，为最优平移距离，根据相对距离信息Qm(l)计算最优平移距离最优平移距离通过最大化下列公式得到：其中，Sm-1(l-Xm)表示第m-1个线性调频周期平移后的距离信息；第m个线性调频周期校准后的距离信息为：其中，j为单位复数，Δ是向量[0，1，...，L-1]；第m个线性调频周期的中频信号快速傅里叶变换(FFT)后的结果为第m个线性调频周期校准后的距离信息；从而得到M个线性调频周期的中频信号快速傅里叶变换后的结果，通过M个线性调频周期的中频信号快速傅里叶变换后的结果得到身体运动补偿后的震动位移信息。所述步骤二具体为：我们将背景杂波看作是由大量小物体影响产生的，这些背景杂波的振幅和相位都是随机的，但他们的频谱包络a遵循韦伯分布：因此他们可以用如下方法进行杂波分离。2.1)分离背景杂波：将步骤1得到的M个线性调频周期的中频信号快速傅里叶变换后的结果进行第二次快速傅里叶变换得到距离多普勒矩阵(RDM)，然后对距离多普勒矩阵进行对数正则化；将设定的杂波分离率代入杂波分离率的定义公式中计算阈值μ0，杂波的分离率通过人为设定，杂波的分离率定义公式为：其中，γ是欧拉常数，pc为杂波分离率；定义弹性矩阵弹性矩阵中每个元素定义如下：其中，代表距离多普勒矩阵，为弹性矩阵中第i行的第j个元素，Rij为距离多普勒矩阵中第i行第j个元素，为无偏估计，对距离多普勒矩阵进行更新，具体过程如下：J为01矩阵，ο为哈达玛积，sgn是符号函数；最终得到更新后的距离多普勒矩阵，从而完成杂波分离；2.2)去除动态杂波：由于行人、车辆等造成的杂波不遵循韦伯分布，上述步骤无法去除这些动态杂波。本专利技术使用跨越多个RDM的信息来消除行人和车辆移动引起的动态杂波。持续发射调频连续波(FMCW)，得到n个连续的距离多普勒矩阵，每个距离多普勒矩阵包含M个线性调频周期，计算多个距离多普勒矩阵中物体的平均移动速度，确定单位时间移动距离大于距离轮廓分辨率ΔRES的物体，并移除物体在距离多普勒矩阵中对应的值(该值为物体的速度)，进而实现对行人和车辆造成的动态杂波的去除。二次傅里叶变换得到的距离多普勒矩阵包含了速度信息；距离多普勒矩阵中的物体为调频连续波探测到的所有障碍物。所述步骤三中的声源特征为声门流倒谱系数(glottalflowcepstrumcoefficients，GFCC)、残余相位倒谱系数(residualphaseCepstrumcoefficients,RPCC)、Teager相位倒谱系数(Teagerphasecepstrumcoefficients,TPCC)；所述步骤三中的声轨特征为频谱特征，梅尔频率倒谱系数(Melfrequencycepstralcoefficients,MFCC)，线性预测系数(LPC),线性预测倒谱系数(LPCC)和线频谱频率(LSP)。3.1)声源特征提取：声源信号代表了声带的肌肉结构和张力，以及相关的声门脉冲参数，例如声门的闭合瞬时率，张开持续时间和张开程度。声带的振动模式不仅为语音生产提供了声音来源，而且还具有独特的非线性流动模式。本专利技术使用声门流倒谱系数(glottalflowcepstrumcoefficients，GFCC)表示说话人声门激励模式的频谱幅度特征。使用残余相位倒谱系数(residu本文档来自技高网...

【技术保护点】
1.一种基于毫米波感知的非接触式声纹生物认证方法，其特征在于，包括以下步骤：/n步骤一：声音震动的毫米波感知：使用毫米波探针感知喉咙的震动，对用户的喉咙区域发射调频连续波，然后接收用户喉咙区域反射的带有声音震动信息的调频连续波；毫米波探针将接收的调频连续波的周期性线性调频信号转化为中频信号，通过计算得到多个线性调频周期的中频信号快速傅里叶变换后的结果，并计算声音的震动位移信息；/n步骤二：杂波抑制处理：将步骤一中多个线性调频周期的中频信号快速傅里叶变换后的结果进行第二次快速傅里叶变换得到距离多普勒矩阵，通过更新距离多普勒矩阵完成背景杂波的分离和动态杂波的去除；/n步骤三：从步骤二去除杂波后的距离多普勒矩阵中提取声源特征和声轨特征；/n步骤四：用户识别：构建细粒度认证模型，将步骤三中提取的声源特征和声轨特征输入细粒度认证模型中得到用户的分类结果。/n

【技术特征摘要】
1.一种基于毫米波感知的非接触式声纹生物认证方法，其特征在于，包括以下步骤：
步骤一：声音震动的毫米波感知：使用毫米波探针感知喉咙的震动，对用户的喉咙区域发射调频连续波，然后接收用户喉咙区域反射的带有声音震动信息的调频连续波；毫米波探针将接收的调频连续波的周期性线性调频信号转化为中频信号，通过计算得到多个线性调频周期的中频信号快速傅里叶变换后的结果，并计算声音的震动位移信息；
步骤二：杂波抑制处理：将步骤一中多个线性调频周期的中频信号快速傅里叶变换后的结果进行第二次快速傅里叶变换得到距离多普勒矩阵，通过更新距离多普勒矩阵完成背景杂波的分离和动态杂波的去除；
步骤三：从步骤二去除杂波后的距离多普勒矩阵中提取声源特征和声轨特征；
步骤四：用户识别：构建细粒度认证模型，将步骤三中提取的声源特征和声轨特征输入细粒度认证模型中得到用户的分类结果。


2.根据权利要求1所述的一种基于毫米波感知的非接触式声纹生物认证方法，其特征在于，所述步骤一中的根据中频信号计算声音震动位移的具体过程为：
1.1)设一个线性调频周期为Tr，根据中频信号获得第m个线性调频周期的震动位移d(mTr)的计算公式为：



其中，m为线性调频周期的编号，m＝0,1,2…,M-1；M为线性调频周期的数量；f0是载波频率；Δψm是对第m个线性调频周期的中频信号进行快速傅里叶变换后的结果；
从而根据第m个线性调频周期的液晶屏幕位移d(mTr)得到M个线性调频周期的震动位移。
1.2)为了避免身体的随机运动对步骤1.1)结果的干扰，通过平移变化对步骤1.1)获得的M个线性调频周期的震动位移进行身体运动补偿：
定义Sm(l)为获取的第m个线性调频周期的距离信息，的值为第m个线性调频周期校准后的距离信息；其中，l＝0,1,2,..,L-1，L为第m个线性调频周期内采样点的数量，l为采样点的编号；
定义第m个线性调频周期的相对距离信息为Qm(l)，具体为：



其中，Qm-1(l)为第m-1个线性调频周期的相对距离信息；
定义Xm为Sm(l)的平移距离，为最优平移距离，根据相对距离信息Qm(l)计算最优平移距离最优平移距离通过最大化下列公式得到：



其中，Sm-1(l-Xm)表示第m-1个线性调频周期平移后的距离信息；
第m个线性调频周期校准后的距离信息为：



其中，j为单位复数，Δ是向量[0,1,…,L-1]；
第m个线性调频周期的中频信号快速傅里叶变换后的结果为第m个线性调频周期校准后的距离信息；
根据第m个线性调频周期的中频信号快速傅里叶变换后的结果得到M个线性调频周期的中频信号快速傅里叶变换后的结果，通过M个线性调频周期的中频信号快速傅...

【专利技术属性】
技术研发人员：许文曜，林峰，李勤，李正雄，陈百成，
申请(专利权)人：杭州环木信息科技有限责任公司，
类型：发明
国别省市：浙江;33