馬 鵬,王澤宇,鐘衛(wèi)東,王緒安

(武警工程大學 網(wǎng)絡(luò)與信息安全武警部隊重點實驗室,西安 710086)

0 概述

密碼設(shè)備在數(shù)據(jù)加密過程中通常會有功耗[1]、時間[2]和電磁輻射[3]等側(cè)信道信息泄漏。側(cè)信道攻擊(Side Channel Attack,SCA)即利用泄漏的側(cè)信道信息攻擊密碼設(shè)備,采用數(shù)學統(tǒng)計分析方法計算泄漏信息與加密數(shù)據(jù)(或者解密數(shù)據(jù))之間的關(guān)系以破解密鑰。功耗攻擊是側(cè)信道攻擊的一種,其中,BRIER等人在2004年提出的相關(guān)功耗攻擊(Correlation Power Attack,CPA)[4]因具有較強攻擊性和密鑰破解高效性被廣泛應用。

在實施側(cè)信道攻擊前,通常要對密碼設(shè)備泄漏信息進行采樣,采樣數(shù)據(jù)的純凈度與攻擊效率和密鑰破解正確率密切相關(guān)。在實際功耗攻擊中,采集功耗信息會受到環(huán)境噪聲、熱噪聲和算法噪聲等干擾,而噪聲在一定程度上會降低功耗攻擊效率與密鑰破解正確率,甚至導致攻擊失敗。為消除噪聲影響,提高采樣信號質(zhì)量,研究人員采用多種方法對數(shù)據(jù)降噪。文獻[5]用估計的4階累積量代替原始信號來執(zhí)行傳統(tǒng)CPA和差分功耗攻擊(Differential Power Attack,DPA)以提高攻擊性能。文獻[6]利用小波變換對功耗曲線進行對準,增強功耗的信噪比信息,提高了DPA的攻擊效率。文獻[7]提出一種基于信息論的特殊降噪閾值,但其在功耗數(shù)據(jù)降噪方面不具備普遍性。文獻[8]通過仿真驗證了小波變換的降噪效率要優(yōu)于高階累積量。基于小波變換的去噪過程能較好地表征以低頻信息為主的信號,但其忽略高頻信息,不能很好地分解包含大量細節(jié)信息的信號,在表征數(shù)據(jù)細節(jié)信息方面存在一定缺陷[9]。針對小波變換無法處理高頻信號的問題,文獻[10]在進行數(shù)據(jù)預處理時使用小波包閾值法對數(shù)據(jù)降噪,從而提高功耗攻擊效率。在實際采用小波包閾值對不同功耗數(shù)據(jù)降噪時,閾值選取通常隨功耗數(shù)據(jù)特征而變化,并會影響到小波包降噪性能。此外,小波包閾值只對低頻信息降噪,忽略了高頻信息中仍存在噪聲含量,從而降低功耗數(shù)據(jù)質(zhì)量。

針對上述問題,本文提出一種改進的小波包分解降噪方法。采用奇異譜分析(Singular Spectrum Analysis,SSA)[11]通過奇異值分解將復雜信號分為不同的子序列。奇異譜分析與主成分分析[12]不同,可基于特殊的矩陣結(jié)構(gòu)處理單次功耗數(shù)據(jù),因此將奇異譜分析用于處理小波包分解的低頻部分和高頻部分,并依據(jù)奇異熵的波動趨勢[13]自適應地從各部分提取功耗信息,以提高功耗數(shù)據(jù)的質(zhì)量與純凈度。在此基礎(chǔ)上,使用原始功耗數(shù)據(jù)以及改進前后的小波包降噪功耗數(shù)據(jù)對硬件實現(xiàn)的SM4算法進行選擇明文攻擊,對所得相關(guān)功耗攻擊的攻擊效率和密鑰破解準確率進行對比分析。

1 基礎(chǔ)知識

1.1 相關(guān)功耗攻擊的原理與流程

側(cè)信道攻擊中相關(guān)功耗攻擊方法是利用真實功耗數(shù)據(jù)與模擬功耗數(shù)據(jù)之間線性關(guān)系來破解密鑰,該方法包括以下4個步驟:

1)選取N組不同明文(或者密文)通過密碼設(shè)備進行加密(或者解密)操作,采集密碼設(shè)備的功耗曲線,記為P。

2)選取中間值函數(shù)f(d,k),中間值函數(shù)的選取與部分密鑰和部分明文(或者密文)有關(guān),然后猜測密鑰,通過中間值函數(shù)計算對應中間值,根據(jù)中間值的漢明重量模型或者漢明距離模型,計算模擬功耗數(shù)據(jù),記為H。

3)根據(jù)式(1)將真實功耗P與模擬功耗H按列求出相關(guān)系數(shù)ρ,最后得到相關(guān)系數(shù)矩陣R。

(1)

4)觀察矩陣R中各值,其中最大值對應的猜測密鑰即為破解的正確密鑰[14]。

1.2 小波包分解

小波變換(Wavelet Transform,WT)方法是基于傅里葉變換、泛函數(shù)分析、數(shù)值分析等數(shù)學分析方法提出的一種信號分析與處理方法,廣泛應用于信號處理、圖像分析、語音處理等領(lǐng)域[15]。小波變換方法根據(jù)頻率不同,通過對小波基采取伸縮和平移等操作對信號進行多尺度精細化分析,從而實現(xiàn)對信號時域和頻域的局部變換,提取信號中有效信息。信號S的兩層小波分解流程如圖1所示。其中:A為信號低頻數(shù)據(jù)(即近似部分),為信號的主要信息;D為信號高頻數(shù)據(jù)(即細節(jié)部分),為信號的次要信息,通常被視為噪聲。由圖1可知,小波變換方法只是針對信號低頻數(shù)據(jù),而忽略信號高頻數(shù)據(jù)[16],這導致小波變換無法充分表征包含高頻數(shù)據(jù)的信號,因此小波變換方法去噪存在一定缺陷。

圖1 信號S的兩層小波分解流程

小波包變換(Wavelet Packet Transform,WPT)方法是對小波變換方法的提升與改進[17],其基本思想是讓信息能量集中,在細節(jié)中尋找有序性并篩選出其中規(guī)律,從而對信號進行精細分析。與小波變換方法僅對信號低頻數(shù)據(jù)分解不同,小波包變換方法對信號高頻數(shù)據(jù)也進行分解,并根據(jù)被分析信號的特征自適應地選擇相應頻帶,使之與信號頻譜相匹配,從而提高時頻分辨率。小波包變換在小波變換基礎(chǔ)上提供更多可使用的正交基,信號S的兩層小波包分解流程如圖2所示。

圖2 信號S的兩層小波包分解流程

小波包分解的函數(shù)方程如式(2)[18]所示:

(2)

其中,j為尺度指標(頻域參數(shù)),k為位置指標(時間參數(shù)),2j為分辨率,n=0,1,…,N為振蕩次數(shù)。

在使用小波包分析信號時,通常選擇n=0時函數(shù)φ(t)和ψ(t)作為正交尺度函數(shù)和小波函數(shù)進行分解和變換,表達式如下:

(3)

(4)

其中,h0,k和h1,k為濾波器系數(shù)。

當n=1,2,…,N時,對應的小波包函數(shù)為:

(5)

(6)

(7)

(8)

(9)

兩層小波包分解的重構(gòu)算法為:

(10)

其中,j=J-1,…,1,0;i=2j,…,2,1;J=lbN;h和g為小波重構(gòu)濾波器。

1.3 奇異譜分析

奇異譜分析方法是一種分析數(shù)據(jù)不同成分分布的方法,主要應用于非線性數(shù)據(jù)。該方法先將數(shù)據(jù)在軌跡矩陣重構(gòu)變換,經(jīng)奇異值分解(Singular Value Decomposition,SVD)[19]后將數(shù)據(jù)分組重構(gòu),最終使用不同成分數(shù)據(jù)取代原始數(shù)據(jù)。與主成分分析(Principal Component Analysis,PCA)[12]不同的是,由于奇異譜分析可以基于特殊的矩陣結(jié)構(gòu)處理單次功耗數(shù)據(jù),因此其用于處理小波包分解低頻部分和高頻部分。

奇異譜分析包括分解和重構(gòu)兩部分[20]:

1)分解。

假設(shè)待分析數(shù)據(jù)是長度為T的1維離散時間序列,選擇適當窗口寬度L(2≤L≤T),將待分析數(shù)據(jù)YT=(y1,y2,…,yT)轉(zhuǎn)化為多維軌跡矩陣X:

(11)

其中,K=T-L+1,窗口寬度L按照式(12)[21]選取最優(yōu)值。

L=?lg(T)c」,c∈[1.5,3]

(12)

計算XXT得到L個特征值λ1≥λ2≥…≥λL≥0,U1,U2,…,UL為相應的特征向量,另外選取d=max(i,λi>0)=R(A),將多維軌跡矩陣X分解如下:

X=X1+X2+…+Xd

(13)

2)重構(gòu)。

對式(13)中Xi進行變換生成相應時間序列,得到的每組數(shù)據(jù)均表征原始數(shù)據(jù)某方面的特征。處理后的數(shù)據(jù)通常存于幾個主要奇異值對應的時間序列中,例如M?{1,2,…,d}。選擇合適的主成分并根據(jù)式(14)重構(gòu)新的時間序列數(shù)據(jù)H,表達式為:

(14)

根據(jù)式(15)對XL矩陣求對角平均可計算得到H。

(15)

其中,x*為矩陣XL中的元素,且:

L*=min{L,K}

(16)

K*=max{L,K}

(17)

2 基于改進小波包分解的相關(guān)功耗攻擊

在使用小波包變換方法降噪的過程中,閾值選擇和閾值量化方法與數(shù)據(jù)降噪效果緊密相關(guān)[22]。在實際應用中,可選擇默認閾值或者不斷測試調(diào)整參數(shù)兩種方式進行降噪。其中,默認閾值的方式缺乏針對性,無法根據(jù)功耗數(shù)據(jù)特點進行準確降噪。而在不斷測試調(diào)整參數(shù)的方式中,參數(shù)只適用于當前功耗數(shù)據(jù)降噪,不具有普遍性。此外,小波包變換方法只針對每層高頻系數(shù)閾值進行降噪,忽略低頻系數(shù)中大量噪聲[23],從而降低功耗攻擊效率與密鑰破解準確率。

針對上述問題,本文將奇異譜分析添加到小波包分解降噪過程中,對功耗數(shù)據(jù)進行預處理以提高功耗攻擊效率和密鑰破解準確率。針對從實驗室采集到的功耗數(shù)據(jù),使用小波包對每條功耗高頻信息與低頻信息逐層分解形成小波包分解樹(本文選擇使用sym6小波,分解6層),選擇一個合適熵標準,使用Matlab自帶的小波包相關(guān)函數(shù)(besttree()函數(shù)計算小波包最佳樹,由leaves()函數(shù)獲得小波包樹所有節(jié)點,wpcoef()函數(shù)計算節(jié)點系數(shù)值)求解小波包分解的最佳小波包樹,如圖3所示。

圖3 最佳小波包分解樹

使用奇異譜分析計算小波包樹各節(jié)點最佳奇異譜值。某節(jié)點分解后低頻系數(shù)和高頻系數(shù)的奇異譜值分布如圖4所示。而根據(jù)奇異熵定義:

(18)

其中,λ為對應的奇異值,k為奇異熵階次,ΔEi為奇異熵在階次i處的增量,Ek為k階對應的奇異熵。由式(18)計算得到各奇異值對應的奇異熵,獲得節(jié)點奇異熵分布,如圖5所示。

圖4 節(jié)點低頻系數(shù)和高頻系數(shù)奇異譜值分布

圖5 節(jié)點低頻系數(shù)和高頻系數(shù)奇異熵分布

由圖5可以看出,節(jié)點低頻系數(shù)和高頻系數(shù)奇異熵均隨奇異值的增大而減小,最終趨于平緩,可認為處于平緩位置的奇異熵值所代表的節(jié)點信息為噪聲含量。根據(jù)奇異熵分布趨勢去除其中的噪聲,通過奇異譜分析重構(gòu)去除噪聲節(jié)點系數(shù),對處理后各節(jié)點進行小波包重構(gòu),可得到預處理后的功耗數(shù)據(jù)。

圖6為改進后的小波包降噪流程,具體步驟為:

1)利用小波包分解功耗數(shù)據(jù)。先挑選一個對稱性、緊支撐性和正交性較好的小波作為小波包分解的小波基,再確定小波包分解層數(shù)。

2)獲得最佳小波包樹。

3)使用奇異譜分析奇異熵分布趨勢,對小波包分解系數(shù)進行降噪處理。

4)小波包重構(gòu)。將上述步驟中處理的節(jié)點系數(shù)重新寫入小波包樹節(jié)點,通過小波包重構(gòu)得到降噪后的功耗數(shù)據(jù)。

圖6 改進后小波包降噪流程

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

為評估本文方法對相關(guān)功耗攻擊效率和密鑰破解準確率的提升效果,選擇基于硬件的SM4算法進行選擇明文攻擊。通過SM4對特定明文進行加密,選擇第1輪4個S盒作為泄漏點采集功耗數(shù)據(jù),其中采集的1條功耗數(shù)據(jù)如圖7所示,采樣速率為5 Gp/s。使用原始功耗數(shù)據(jù)、改進前小波包降噪功耗數(shù)據(jù)和改進后小波包降噪功耗數(shù)據(jù)進行相關(guān)功耗攻擊,分析不同功耗數(shù)據(jù)對應的攻擊效率和密鑰破解準確率。

圖7 功耗數(shù)據(jù)

3.1 攻擊性能分析

在相關(guān)功耗攻擊中,相關(guān)系數(shù)是衡量去噪方法性能的重要指標。使用原始功耗數(shù)據(jù)與經(jīng)過預處理的功耗數(shù)據(jù)進行相關(guān)功耗攻擊,通過比較正確密鑰攻擊結(jié)果的相關(guān)系數(shù)大小可確定去噪性能高低。圖8為使用原始功耗數(shù)據(jù)、改進前小波包降噪功耗數(shù)據(jù)和改進后小波包降噪功耗數(shù)據(jù)的相關(guān)功耗攻擊對比結(jié)果,可見正確密鑰(242)對應的相關(guān)系數(shù)均最大,表示3種功耗數(shù)據(jù)都能確保相關(guān)功耗攻擊成功破解密鑰。在原始功耗數(shù)據(jù)與改進前小波包降噪功耗數(shù)據(jù)下相關(guān)功耗攻擊中,出現(xiàn)與正確密鑰相關(guān)系數(shù)0.192 9、0.191 2對應峰形較接近的尖峰,通常被稱為“鬼峰”。鬼峰會對相關(guān)功耗攻擊結(jié)果造成干擾,導致功耗曲線條數(shù)增加。而使用改進后小波包降噪功耗數(shù)據(jù)在相關(guān)功耗攻擊下,其鬼峰與正確密鑰相關(guān)系數(shù)對應峰形距離較大,對相關(guān)功耗攻擊結(jié)果干擾較小,相關(guān)功耗攻擊結(jié)果可靠性得到提高。此外,原始功耗曲線、改進前小波包降噪功耗曲線和改進后小波包降噪功耗曲線對應的最大相關(guān)系數(shù)分別為0.194 5、0.242 1和0.412 7。由相關(guān)功耗攻擊原理可知,最大相關(guān)系數(shù)值越大則越易區(qū)分出正確密鑰對應的功耗,因此,使用改進后小波包降噪功耗數(shù)據(jù)的相關(guān)功耗攻擊性能較原始功耗數(shù)據(jù)和小波包降噪功耗數(shù)據(jù)分別提升了53%與41%。

圖8 3種數(shù)據(jù)相關(guān)功耗攻擊結(jié)果對比

3.2 攻擊效率分析

圖9～圖12分別為使用原始功耗數(shù)據(jù)、改進前小波包降噪功耗數(shù)據(jù)和改進后小波包降噪功耗數(shù)據(jù)對4個S盒進行相關(guān)功耗攻擊,得到相關(guān)系數(shù)與功耗曲線數(shù)量之間的關(guān)系?？梢钥闯?當原始功耗曲線數(shù)量分別為190條、120條、200條和130條時成功破解密鑰(此時正確密鑰與錯誤密鑰的`相關(guān)系數(shù)開始分離且后續(xù)不再重合);使用改進前小波包降噪對數(shù)據(jù)處理后,在破解密鑰上會導致功耗曲線數(shù)量增長,無法正確破解第3個S盒的密鑰;使用改進后小波包降噪對數(shù)據(jù)處理后,在破解密鑰上功耗曲線數(shù)量會出現(xiàn)不同程度地減少,破解4個S盒密鑰分別需要120條、100條、130條和120條功耗曲線,其攻擊效率較使用原始功耗數(shù)據(jù)分別提高37%、17%、35%和8%。

圖9 第1個S盒功耗曲線數(shù)量與相關(guān)系數(shù)的關(guān)系

圖10 第2個S盒功耗曲線數(shù)量與相關(guān)系數(shù)的關(guān)系

圖11 第3個S盒功耗曲線數(shù)量與相關(guān)系數(shù)的關(guān)系

圖12 第4個S盒功耗曲線數(shù)量與相關(guān)系數(shù)的關(guān)系

4 結(jié)束語

本文將奇異譜分析與小波包降噪相結(jié)合,提出一種基于改進小波包分解的相關(guān)功耗攻擊降噪方法。使用小波包分解求出最優(yōu)小波包樹,提取各節(jié)點中功耗數(shù)據(jù)的低頻系數(shù)及高頻系數(shù),運用奇異譜分析并根據(jù)各奇異熵分布趨勢去除節(jié)點系數(shù)的噪聲信息,并將處理后的節(jié)點系數(shù)寫入小波包樹,通過數(shù)據(jù)重構(gòu)獲得降噪后的功耗數(shù)據(jù)。對SM4算法進行選擇明文的相關(guān)功耗攻擊實驗表明,該方法能有效解決小波包分解中功耗數(shù)據(jù)預處理缺乏針對性的問題,保留高頻功耗數(shù)據(jù)信息并提高功耗質(zhì)量,較改進前小波包降噪方法的攻擊效率更高。雖然本文方法在攻擊性能和效率上有所提升,但其本質(zhì)是將功耗數(shù)據(jù)作為信號數(shù)據(jù)來進行預處理降噪,并未考慮功耗數(shù)據(jù)本身特征,導致降噪效果和攻擊效率不穩(wěn)定。后續(xù)將分析功耗數(shù)據(jù)特征并提出更具針對性的功耗數(shù)據(jù)預處理方法,以進一步提升相關(guān)功耗攻擊性能和效率。