基于隨機映射的聲紋模板保護方法*

朱華虹 賀前華 李艷雄 張雪源

(華南理工大學(xué) 電子與信息學(xué)院，廣東 廣州 510640)

近年來，生物識別技術(shù)的廣泛應(yīng)用對生物特征模板的存儲和傳輸安全提出了較高要求.生物特征模板保護旨在解決特征模板的可撤銷性、多樣性、安全性、性能保持等問題［1］，但由于生物特征數(shù)據(jù)的模糊性和生物識別技術(shù)的差異性［2］，尚未有任何一種算法能夠滿足所有生物特征模板保護的要求.

目前，指紋、臉型、虹膜等圖像特征的模板保護算法研究較多［3-5］，基于矢量量化(VQ)聲紋認證系統(tǒng)的模板保護算法也相繼提出.文獻［6］提出了傳統(tǒng)密碼體制下的聲紋模板保護方法，通過改進數(shù)字簽名方案并與認證協(xié)議相融合實現(xiàn)聲紋特征的保護，其缺點是模板庫中的聲紋特征仍然是原始特征，無法保證存儲的安全性.文獻［7］將經(jīng)典的模糊保險箱思想應(yīng)用于聲紋模板保護，使用雜湊點達到隱藏真實點的目的，此方法的缺陷是安全性依賴于雜湊點的數(shù)量，導(dǎo)致存儲效率不高、交叉匹配等問題.基于多子空間映射的可撤銷聲紋模板保護方法［8］將隨機映射應(yīng)用于基于VQ 的聲紋認證系統(tǒng)，通過分析變換前后特征點與特征點、特征點與碼本的距離證明方法的有效性，但缺乏安全性與身份認證性能保持的深入分析和實驗證明.

鑒于VQ 適用于小語料庫的說話人識別任務(wù)，而目前主流的文本無關(guān)說話人識別模型為高斯混合模型(GMM)，因此，文中在傳統(tǒng)聲紋認證系統(tǒng)基礎(chǔ)上采用隨機映射更為通用的形式化表示，提出一種可撤銷的聲紋模板保護方法，并就方法的安全性和認證性能保持之間的關(guān)系進行研究.最后通過仿真實驗驗證了方法的有效性.

1 可撤銷的聲紋模板保護方案

1.1 隨機映射的形式化表示

文獻［8］采用式(1)將原始特征x 隨機映射為變換后的特征y，并通過分析變換前后特征點間的歐氏距離，以其不變來證明身份認證性能得到保持.

式中:R 為正交矩陣，R∈d×n(d =n)，d、n 為隨機矩陣的行列數(shù).

模板保護技術(shù)要求在提高安全性的同時不會惡化身份認證性能［9］，但實際應(yīng)用中安全性和身份認證性能保持往往很難同時滿足.對于隨機映射來說，降維是提高安全性的有效途徑［10］.為了便于后文對算法的安全性和認證性能進行綜合評價，定義通用形式化表示如下:

相關(guān)研究［11］表明，元素服從高斯分布的隨機矩陣R 具有正交性，尤其在高維空間中，R 中元素服從高斯分布與正交是近似統(tǒng)一的.在高維空間或滿足Johnson-Lindenstrauss (JL)定理降維下限的條件下，可近似認為單位矩陣［10］.特別地，當(dāng)d=n 時，式(2)轉(zhuǎn)化為式(1).

1.2 基于隨機映射的聲紋認證系統(tǒng)實現(xiàn)

GMM 雖然具有識別效果好、魯棒性高等特點，但其識別精度對訓(xùn)練數(shù)據(jù)依賴性大，如果將變換后的聲紋特征直接存儲為模板，則認證時進行模型訓(xùn)練的時間較長，無法滿足應(yīng)用要求.因此，所提認證系統(tǒng)在模板庫中存儲的不是變換后的聲紋特征而是訓(xùn)練獲得的GMM 模型參數(shù)，以提高認證的實時性.圖1為文中提出的基于GMM 聲紋認證系統(tǒng)的模板保護算法框圖.在注冊階段，語音經(jīng)過預(yù)處理提取原始美爾頻率倒譜系數(shù)(MFCC)特征，并將原始特征映射到隨機空間后進行GMM 模型學(xué)習(xí)，存儲模型參數(shù)作為模板;在認證階段，提取的MFCC 特征和相同的隨機矩陣相乘后再與模型庫進行匹配.具體步驟如下.

圖1 聲紋認證系統(tǒng)框圖Fig.1 Block diagram of voiceprint verification system

注冊階段:

(1)提取注冊語音的MFCC 特征序列Xk，Xk=為d 維語音特征矢量.

(3)使用Yk訓(xùn)練GMM 模型，存儲模型參數(shù)作為模板.

認證階段:

(1)提取待認證語音的MFCC 特征序列Xq，為d 維語音特征矢量;

(2)使用注冊階段的矩陣R，用式(2)對聲紋矢量進行變換得到

(3)使用Yq進行模型匹配，并決策輸出.

訓(xùn)練GMM 需要較充分的訓(xùn)練語音，如果在實際應(yīng)用中訓(xùn)練數(shù)據(jù)較少，通用的做法是采用GMMUBM 框架［12］，即利用所有說話人的語音訓(xùn)練通用背景模型(UBM)，然后根據(jù)目標(biāo)說話人的數(shù)據(jù)從UBM 自適應(yīng)獲得目標(biāo)說話人模型.優(yōu)化GMM 性能的算法比較多［13］，鑒于文中重點研究隨機映射對于傳統(tǒng)GMM 聲紋模板的保護有效性，故所提算法中模型訓(xùn)練主要由目標(biāo)說話人的語音單獨訓(xùn)練得到.根據(jù)隨機映射的原理［10］，對于GMM-UBM 框架，文中也可將所有說話人的數(shù)據(jù)映射至與目標(biāo)說話人相同的隨機空間中進行訓(xùn)練，但應(yīng)用中目標(biāo)說話人各自所映射的空間不同，故需要在不同變換域中進行UBM 訓(xùn)練，并在目標(biāo)說話人對應(yīng)的變換域中進行識別.

研究［9］表明，要滿足所有生物特征模板保護的要求是很困難的.在實際應(yīng)用中，這些要求往往相互制約，因此需要根據(jù)相應(yīng)的應(yīng)用場合作出合理的選擇.下面就所提算法在身份認證性能保持、安全性、計算復(fù)雜度等方面進行分析.

2 算法分析

2.1 身份認證性能分析

基于GMM 的聲紋認證系統(tǒng)中認證話音與模型的匹配程度主要通過計算聲紋矢量和GMM 模型的對數(shù)似然度獲得［14］，如果在變換域聲紋特征的統(tǒng)計特性沒有改變，那么識別準(zhǔn)確性不變，即系統(tǒng)的身份認證性能可以保持.首先分析d=n 的條件下變換前后身份認證性能是否相等.

設(shè) 代表說話人的一組GMM 模型參數(shù)，={M，wi，Ui，Σi}，其中M 為混合度，wi、Ui、Σi分別為第i 個高斯的權(quán)重、均值和協(xié)方差矩陣;'代表變換后同一說話人的GMM 模型參數(shù)，' ={M，w'i，U'i，Σ'i}，其中w'i，U'i，Σ'i分別為變換后第i 個高斯的權(quán)值、均值和協(xié)方差矩陣.變換前后聲紋矢量序列Xq、Yq與GMM 模型的對數(shù)似然度分別為和，其計算方法如下:

其中，

故比較式(4)、(5)是否相等需研究變換后的模型參數(shù).聲紋特征經(jīng)映射后首先采用LBG 算法獲得初始的聚類中心，計算初始權(quán)值、均值和方差，再使用最大期望估計(EM)算法估計模型參數(shù).正交變換前后特征矢量之間的距離保持不變且特征矢量與碼字之間距離保持不變.從幾何空間上分析，即映射前后M 個類中的點集合是一致的，因而初始的權(quán)值(各個類中所包含的特征矢量的個數(shù)占總的特征矢量的百分比)也不變，而初始的均值U'i和協(xié)方差Σ'i分別變?yōu)?

其中Xl、Yl分別為變換前后第i 個高斯聚類中包含的特征矢量，m 為特征矢量個數(shù)，使用EM 算法進行迭代運算可得:

第i 個高斯的權(quán)重w'i為

第i 個高斯的均值U'i為

第i 個高斯的協(xié)方差矩陣Σ'i為

綜上所述，將式(14)-(16)代入式(4)、(5)，可推導(dǎo)出L(Xq)和L(Yq')相等，即變換前后身份認證性能不變.但此時變換前后的原始特征或模型參數(shù)存在一一對應(yīng)關(guān)系，R 一旦泄漏就可破解所有原始數(shù)據(jù).結(jié)合后續(xù)安全性的分析，需要研究采用降維進一步提高安全性，同時保持身份認證性能的可能性.根據(jù)JL 定理［15］，當(dāng)滿足降維下限時歐氏空間中兩點映射至低維空間，其距離變化保持在一個小范圍內(nèi)是可能的，同理GMM 均值間的距離也可以近似保持.另一方面，文獻［16］指出，d空間中的任何高斯分布映射至n 維子空間時，其超橢球體的離心率有可能變小.為便于分析，文中定義某個高斯分量的原始協(xié)方差參數(shù)間滿足β1≤β2≤…≤βd，可得變換前的離心率為變換后的協(xié)方差參數(shù)為可得:

其中rki∈R.可見變換后超橢球體的離心率小于等于從而更有利于聚類.故GMM 的認證性能相對于只考慮矢量間距離差異的算法(如VQ)可以保持得更好.

2.2 安全性分析

生物特征模板保護對安全性的要求體現(xiàn)在模板的變換或加密滿足不可逆性.不可逆性具有兩點重要特征［1］，一是識別在變換域進行;二是即使模板被盜(泄漏)，敵手也無法獲得原始數(shù)據(jù).文中算法將聲紋特征映射至變換域進行訓(xùn)練，獲得的模型參數(shù)和識別也均在變換域中，滿足不可逆性的第1 個特征.另一方面，所提方法的安全參數(shù)依賴于隨機矩陣，如果模板被盜，敵手在沒有獲得隨機矩陣的條件下很難獲取原始數(shù)據(jù)，滿足不可逆性的第2 個特征.因此，該條件下，方法滿足不可逆性.進一步放寬限制條件，假設(shè)敵手獲得了隨機矩陣R，根據(jù)QR 分解和不定方程最小二乘解的原理，只有降低R 的維度，不定方程y =RTx 才存在多個解而無法精確獲得所有原始數(shù)據(jù).下面分析這種情況下敵手破解原始數(shù)據(jù)的可能性.

定理1 如果矩陣Q 滿足“l(fā)-安全”［17］，則從Q 中選擇l+1 個非零行向量組成矩陣Λ，如果這些行向量線性無關(guān)，則線性方程組y = Λx 至少包含2l+1個變量.

證明 由于Λ 中l(wèi)+1 個向量線性無關(guān)，采用消元法求解方程組的解，可獲得階梯形矩陣，Λ 中l(wèi)+1階子式Dl+1=I.為了方便敘述，不妨設(shè)Dl+1在Λ 的左上角，表示為

根據(jù)矩陣Q 滿足“l(fā)-安全”的定義可知，Q 的每行中至少有l(wèi)+1 個非零元素，即Q 中除了Dl+1外的子矩陣中每行至少還有l(wèi) 個非零元素.因此，每行中的總元素為l+1 +l=2l+1.證畢.

從定理1 可知，當(dāng)d ≥2l +1 時，線性方程組y=Λx 無法求解向量中任一元素.令n =l +1，即當(dāng)時，x 中的任一元素都無法精確求出.即使R 泄漏，也無法通過逆變換求解原始模板，滿足無條件不可逆性［8］，即獲得保護性能的上界.

2.3 計算復(fù)雜度分析

文中提出的方法需要將原始數(shù)據(jù)先映射至隨機空間再進行訓(xùn)練和識別.相對于傳統(tǒng)的聲紋認證系統(tǒng)，計算復(fù)雜度的增加主要在于注冊和認證階段需將原始特征進行隨機變換，而此后在變換域中的模型訓(xùn)練和識別與傳統(tǒng)方法是相同的.因此，文中只需分析特征變換產(chǎn)生的新增計算復(fù)雜度.對于長度為p 的d 維語音特征矢量，進行隨機變換(R∈d×n，d≥n)的計算復(fù)雜度為O(dnp).如果采用GMMUBM 框架，還需增加對所有說話人的聲紋特征(長度為p'的d 維語音特征矢量)進行變換的計算復(fù)雜度O(dnp').如果采用降維(d ＞n)提高安全性，則還能進一步降低計算復(fù)雜度.

3 實驗結(jié)果與分析

實驗使用863 漢語普通話連續(xù)語音識別訓(xùn)練庫進行Matlab 仿真，該語料庫共有80 個說話人(男女各40 人)，其中每人各有1560 條不同文本內(nèi)容的語音［8］.所有語音是經(jīng)過挑選的標(biāo)準(zhǔn)普通話發(fā)音人在專業(yè)錄音環(huán)境下錄制，采用高質(zhì)量16 kHz 采樣、16位量化精度、單聲道WAV 格式存儲.在實驗中，隨機選取每人包含不同文本內(nèi)容的500 條訓(xùn)練語音作為訓(xùn)練集，每條訓(xùn)練語音用Cooledit Pro 2.0 去靜音后的長度為4～10 s，總長度約為1 h.從剩余的語音中選取長度為1～2 min 的語音作為測試集，平均每人85 條.特征采用典型的24 階MFCC 特征，對語音進行分幀處理，幀長32 ms，幀移16 ms.

首先研究文中提出的隨機映射方法對身份認證性能的影響.文中采用誤識率(FAR)和識真率(GAR)作為評價認證性能的指標(biāo)［10］.由于模板保護技術(shù)要求的特殊性，一般要求FAR 極低才能夠達到有效保護模板的目的［9］.表1 示出了不同映射條件下，變換前后各GMM 混合度的認證性能(FAR 均為0.62%下的GAR)，其中當(dāng)混合度為2 048 時，采用GMM-UBM 框架.由于文中重點研究隨機映射對傳統(tǒng)GMM 聲紋模板保護的有效性，故在此只初步給出GMM-UBM 在混合度為2048 時的認證性能.加有利于聚類學(xué)習(xí)，因此認證性能可以得到一定的保持.實驗進一步驗證了不同混合度下降維前后各高斯分量對應(yīng)的離心率變化情況，如圖2 所示;總體上降維后的離心率包絡(luò)均小于原始維度的離心率包絡(luò)，證明了理論分析的正確性.

表1 隨機映射前后各GMM 混合度下的系統(tǒng)認證性能Table 1 Verification performance of system before and after random projection with different GMM mixtures

由表1 可以看出，對于傳統(tǒng)GMM 模型，在不同高斯混合度下，經(jīng)正交變換后認證性能與映射前保持一致;當(dāng)維度降為20 時，變換后各混合度下GAR的下降幅度均小于0.2%，可視為近似保持;維度降為12 時，變換后各混合度下GAR 的下降幅度均較大，在8%左右.相同映射條件下，GMM-UBM 認證性能均高于GMM 系統(tǒng)，其主要原因是，GMM-UBM 用大量說話人的語音訓(xùn)練一個穩(wěn)定的高階GMM，可以更好地擬合數(shù)據(jù)集分布，因此具有較高的性能和魯棒性［12］.

下面進一步分析文中所提方法對GMM 模板保護的有效性.表2 示出了隨機映射前后各典型碼本大小下VQ 算法的認證性能(FAR 均為0.62%下的GAR).當(dāng)維度為20 時，變換后VQ 認證性能的下降幅度大于8%;當(dāng)維度降為12 時，變換后VQ 認證性能的下降幅度超過25%，均遠大于GMM 認證性能的下降幅度.可見在相同的降維條件下，GMM 的認證性能保持得更好，與第2 節(jié)中的理論分析一致:GMM 降維后離心率變小，即超橢球變得更圓從而更

表2 隨機映射前后不同碼本大小時VQ 算法的認證性能Table 2 Verification performance of VQ algorithm before and after random projection with different codebook sizes

圖2 不同混合度降維前后各高斯分量下的橢球體離心率Fig.2 Ellipsoid’s eccentricity corresponding to different Gaussian components before and after dimensionality reduction with different mixtures

綜上可知，對于GMM，當(dāng)維度不變或適當(dāng)降維時，認證性能可基本保持;當(dāng)維度下降滿足無條件不可逆時認證性能有一定下降.因此，在實際應(yīng)用中需要根據(jù)不同的場景要求，在安全性和認證性能之間進行折中選擇.

4 結(jié)語

文中基于隨機映射的優(yōu)點，結(jié)合目前主流的文本無關(guān)說話人識別模型GMM，提出了一種可撤銷的聲紋模板保護方法.理論分析和實驗結(jié)果表明，當(dāng)隨機矩陣為正交矩陣或適當(dāng)降維時，變換前后身份認證性能基本保持;當(dāng)滿足無條件不可逆時認證性能有一定的下降.所以說，沒有完美的生物特征模板保護技術(shù)可以滿足所有要求，需要根據(jù)相應(yīng)的應(yīng)用場合作出合理選擇.而且GMM 的優(yōu)化算法較多，隨機映射對優(yōu)化算法的有效性及針對優(yōu)化算法的模板保護是今后研究的主要目標(biāo).

［1］Breebaart J，Yang B，Buhan-Dulman I，et al.Biometric template protection—the need for open standards［J］.Datenschutz und Datensicherheit，2009，33(5):299-304.

［2］Hao F，Anderson R，Daugman J.Combining crypto with biometric effectively［J］.IEEE Transactions on Computers，2006，55(9):1081-1088.

［3］Jin Z，Beng Jin Teoh A，Ong T S，et al.Fingerprint template protection with minutiae-based bit-string for security and privacy preserving［J］.Expert Systems with Applications，2012，39(6):6157-6167.

［4］Lee H G，Beng Jin Teoh A，Jung H G，et al.A secure biometric discretization scheme for face template protection［J］.Future Generation Computer Systems，2012，28(1):218-231.

［5］Karthik N，Anil K J，Sharath P.Fingerprint-based fuzzy vault:implementation and performance ［J］.IEEE Transactions on Information Forensics and Security，2007，2(4):744-757.

［6］Xu W H，He Q H，Li Y X.Cancelable voiceprint templates based on knowledge signatures［C］∥Proceedings of 2008 International Symposium on Electronic Commerce and Security.Guangzhou:IEEE，2008:412-415.

［7］Xu W H，Chen M Y.Cancelable voiceprint template based on Chaff-Points-Mixture method［C］∥Proceedings of International Conference on Computational Intelligence and Security.Suzhou:IEEE，2008:263-266.

［8］徐文華，賀前華，李韜，等.基于MRP 的可撤銷模板設(shè)計及其分析［J］.電子學(xué)報，2009，37(12):2792-2796.Xu Wen-hua，He Qian-hua，Li Tao，et al.Design and analysis of MRP-based cancelable template［J］.Acta Electronica Sinica，2009，37(12):2792-2796.

［9］李鵬，田捷，楊鑫，等.生物特征模板保護［J］.軟件學(xué)報，2009，20(6):1553-1573.Li Peng，Tian Jie，Yang Xin，et al.Biometric template protection ［J］.Journal of Software，2009，20(6):1553-1573.

［10］Wand Y，Plataniotis K N.An analysis of random projection for changeable and privacy preserving biometric verification［J］.IEEE Transactions on Systems，Man，and Cybernetics，Part B:Cybernetics，2010，40(5):1280-1293.

［11］Hecht-Nielsen R.Context vectors:general purpose approximate meaning representations self-organized from raw data ［M］∥Zurada Z M.Computational Intelligence:Imitating Life.Piscataway:IEEE，1994:43-56.

［12］郭武.復(fù)雜信道下的說話人識別［D］.合肥:中國科學(xué)技術(shù)大學(xué)電子工程與信息科學(xué)系，2007:27-42.

［13］何俊，李艷雄，賀前華，等.變異特征加權(quán)的異常語音說話人識別算法［J］.華南理工大學(xué)學(xué)報:自然科學(xué)版，2012，40(3):106-111.He Jun，Li Yan-xiong，He Qian-hua，et al.Speaker recognition algorithm for abnormal speech based on abnormal feature weighting［J］.Journal of South China University of Technology:Natural Science Edition，2012，40(3):106-111.

［14］Stolcke A，Kajarekar S S，F(xiàn)errer L，et al.Speaker recognition with session variability normalization based on MLLR adaptation transforms［J］.IEEE Transactions on Audio，Speech，and Language Processing，2007，15(7):1987-1998.

［15］Johnson W B，Lindenstrauss J.Extension of Lipschitz mapping into Hilbert space ［J］.Contemporary Mathematics，1984，26(1):189-206.

［16］Dasgupta S.Learning mixtures of Gaussians［C］∥Proceedings of IEEE Symposium on Foundations of Computer Science.New York:IEEE，1999:634-644.

［17］Du W L，Chen S G，Han Y S.Privacy-preserving multivariate statistical analysis:linear regression and classification［C］∥Proceedings of the 4th SIAM International Conference on Data Mining.Lake Buena Vista:SIAM，2004:222-233.