摘? 要：張量分解作為一種高維數(shù)據(jù)分析工具能夠結(jié)合多個模態(tài)的信息從而獲取具有判別信息的特征，但是在高維空間上進行張量分解存在計算復(fù)雜度高的問題。為了解決該問題，研究借助隨機奇異值分解速度快的特點，提出基于隨機奇異值分解的張量Tucker分解（張量R-Tucker分解），并將其用于BCICIV2b數(shù)據(jù)集的特征提取和分類中。實驗結(jié)果顯示：相比張量Tucker分解，張量R-Tucker分解特征提取速度提升22%，并且平均分類準(zhǔn)確率達(dá)到80.93%，與現(xiàn)有基于矩陣的方法相比提高10.12%。

關(guān)鍵詞：張量分解；隨機奇異值分解；運動想象；腦機接口

中圖分類號：TP301.6? ? 文獻標(biāo)識碼：A? 文章編號：2096-4706（2023）02-0001-07

Research on BCI Data Classification Based on Tensor R-Tucker Decomposition

ZHANG Shuai

（Taiyuan Normal University， Jinzhong? 030600， China）

Abstract： As a high-dimensional data analysis tool， Tensor Decomposition can combine the information of multiple modalities to obtain features with discriminative information， but there is a problem of high computational complexity problem for Tensor decomposition on high dimensional space. In order to solve the problem， this research takes advantage of the b fast speed of Randomized Singular Value Decomposition， and proposes the Tensor Tucker Decomposition （Tensor R-Tucker Decomposition） based on Randomized Singular Value Decomposition， and uses it for feature extraction and classification of the BCICIV2b dataset. The experimental results show that compared with the Tensor-Tucker Decomposition， the feature extraction speed of the Tensor R-Tucker Decomposition is increased by 22%， and the average classification accuracy rate reaches 80.93%， which is 10.12% higher than the existing matrix-based methods.

Keywords： Tensor Decomposition;Randomized Singular Value Decomposition;motor imagery; Brain Computer Interface

0? 引? 言

腦機接口（Brain-Computer Interface， BCI）系統(tǒng)是一種新型的高效人機交互方式[1]，該系統(tǒng)通過信號分析將大腦神經(jīng)活動信息直接轉(zhuǎn)換成控制外部設(shè)備的命令，從而在大腦與外部設(shè)備之間建立了直接的交流通道，實現(xiàn)人與外界的交流或?qū)ν獠凯h(huán)境的控制。常見的用于BCI的腦電（Electroencephalogram， EEG）信號類型有視覺誘發(fā)電位（Visual Evoked Potential， VEP）[2]、事件相關(guān)電位（Eevent-related potential， ERP）[3]和運動想象（Motor Imagery， MI）[4]等。其中，基于運動想象的腦電信號作為一種經(jīng)典范例已被研究和開發(fā)數(shù)十年。其生理學(xué)基礎(chǔ)是事件相關(guān)同步/去同步（Event-Related Synchronization/Event- Related Desynchronization， ERS/ERD）現(xiàn)象[5]，即當(dāng)人體在想象單側(cè)肢體運動時，同側(cè)大腦皮層的電位活動活躍，特定頻率腦電信號對應(yīng)的振幅明顯增強，而對側(cè)大腦皮層的電位活動被抑制，特定頻率對應(yīng)的腦電信號振幅減弱。MI-EEG信號在時間上是非平穩(wěn)連續(xù)的，具有信噪比（Signal-to-Noiseratio， SNR）低的特點，容易受到多種生物學(xué)因素（如眨眼，心跳和肌肉活動產(chǎn)生的眼電，心電和肌電）或環(huán)境偽影（如市電和工頻噪聲）的影響；同時，MI-EEG個體樣本數(shù)據(jù)量較小，不同受試者間個體差異較大，以上因素均導(dǎo)致MI-EEG分類模型性能較差。

EEG信號是非平穩(wěn)的時域信號且主要特征反映在頻域上，因此EEG信號的時頻分析成為一種常見的特征提取方式[6，7]。它把時域信號在時間-頻率平面中展開，將以時間為自變量的信號表示成以時間和頻率兩個參數(shù)為自變量的函數(shù)，以便表現(xiàn)出信號不同時間點的頻率成分。由此，多通道EEG信號成為一種張量形式的數(shù)據(jù)。然而在傳統(tǒng)EEG分類中，一般會將該張量展開為矩陣的形式進行處理，這種展開不可避免地丟掉了一些維度內(nèi)與維度間可能存在的信息和交互作用。使用張量分解對EEG進行分析可以從時間、空間、頻率等多個模態(tài)上提取和分類有關(guān)的投影模式，并通過投影計算得到特征，由于同時考慮多階信息，該特征具有更好的判別能力。因此很多研究人員將張量分解用于EEG數(shù)據(jù)的分析[8，9]，Zhang等使用動態(tài)腦功能網(wǎng)絡(luò)在時間模態(tài)建立張量模型，使用帶正交約束和部分對稱約束的張量分解方法對構(gòu)造的張量進行了分析，并使用分解得到的核心張量作為特征進行分類，分類結(jié)果表明該方法優(yōu)于基于矩陣的方法[10]。Roáková使用張量分解檢測和分析了腦卒中偏癱患者在運動想象的神經(jīng)反饋訓(xùn)練中運動感覺區(qū)域的腦電節(jié)律激活情況，結(jié)果驗證該方法對于檢測偏癱患者特定運動想象相關(guān)的腦電節(jié)律具有較好的性能[11]。Krishnan等使用短時傅里葉變換對基于運動想象的EEG信號建立張量結(jié)構(gòu)，并在此基礎(chǔ)使用張量分解進行分析，結(jié)果表明此方法可以得到EEG信號可區(qū)分的分量且具有更好的分類性能[12]。李潔提出判別張量分解算法并將其用于基于運動想象的EEG信號的特征提取中，結(jié)果表明該方法由于同時考慮多階信息，對于獲取分類特征，識別與任務(wù)相關(guān)的判別信息非常有效[13]。上述文章均對時域EEG信號建立了張量結(jié)構(gòu)，然后使用張量分解進行分析且取得不錯的效果，然而在高維空間進行張量分解，存在計算復(fù)雜度高的問題。

隨機奇異值分解（Randomized Singular Value Decomposition， RSVD）是一種高效精確的奇異值分解方法。該方法利用隨機投影的思想，通過使用隨機向量生成變換矩陣，將高維的數(shù)據(jù)映射到低維空間中實現(xiàn)降維，同時在低維空間中的數(shù)據(jù)能夠保持高維矩陣的結(jié)構(gòu)特性，通過對低維空間矩陣的分解，從而實現(xiàn)原始大規(guī)模矩陣的高效分解[14]。因此，本文提出一種基于RSVD的張量Tucker分解方法，在張量每個階矩陣化后使用RSVD構(gòu)造因子矩陣，以提高了張量Tucker分解的速度，最后將其應(yīng)用于MI-EEG信號的特征提取，并與張量Tucker分解進行比較，驗證了張量R-Tucker分解的有效性。

1? 基本原理

1.1? 隨機奇異值分解

為了解決大規(guī)模矩陣進行高效精確的奇異值分解（Singular Value Decomposition， SVD）的問題，Halko等人提出了RSVD算法[15]，其基本思想是通過生成一個高斯隨機矩陣，將原始的大規(guī)模矩陣變換為一個小矩陣，通過對該小矩陣的QR分解得到原始矩陣的一組標(biāo)準(zhǔn)正交基，再通過該正交基的變換實現(xiàn)原始矩陣的奇異值分解。Arvind等對該算法進行了詳細(xì)的誤差分析，結(jié)果證明該方法具有很高的精度[16]。RSVD的計算過程如下：

（1）構(gòu)造一個服從高斯分布的隨機投影矩陣P∈Rn×r（r＜n）對原始矩陣X∈Rm×n的列空間進行采樣：

Z=XP

然后對矩陣Z進行低秩QR分解：

Z=QR

其中Q是一個正交矩陣，R是一個上三角矩陣。

（2）利用低秩正交基Q，將X投影到更小的空間：

Y=QTX

然后計算矩陣Y的奇異值分解：

Y=UY∑VT

因為Q是標(biāo)準(zhǔn)正交的并且近似于X的列空間，矩陣∑和V對于Y和X是一樣的。

（3）利用UY和Q可以重構(gòu)出原始矩陣X的左奇異向量U：

U=QUY

得到最終的結(jié)果：

X=U∑V

1.2? 張量Tucker分解

張量Tucker分解是主成分分析方法在高階張量的推廣，它可以將一個N階高維張量? 分解為一個核心張量g和每階乘以一個因子矩陣（A，B，C）的形式[17]，如：一個三階張量X∈RI×J×K的Tucker分解可表示為：

其中 （R1， R2， R3為張量的Tucker秩[17]）是核心張量，表示每一階成分之間的聯(lián)系，通常小于原始張量，A，B，C被稱為因子矩陣，表示每一階成分之間的聯(lián)系，它們通常是正交的，其分解形式如圖1所示。

上述分解形式可通過基于高階奇異值分解（High Order Singular Value Decomposition， HOSVD）和高階正交迭代（High Order Orthogonal Iteration， HOOI）聯(lián)立實現(xiàn)[17，18]。其中，HOSVD的思想是利用矩陣分析中的奇異值分解算法對張量的每一階展開進行分解，并采用低階近似處理，過濾掉一些較小的奇異值；HOOI的基本思想是使用交替投影的方法多次迭代得到最優(yōu)結(jié)果。由于HOSVD算法只對張量的每階進行了一次奇異值分解，因此不能很好地擬合原始數(shù)據(jù)，因此一般用HOSVD算法求出的投影矩陣作為HOOI算法的初值條件，然后使用HOOI算法完成張量Tucker分解。

1.3? 張量R-Tucker分解

張量Tucker分解需要計算高階張量每個階矩陣化的SVD來構(gòu)造因子矩陣，然后使用因子矩陣來構(gòu)造核心張量。但是對于大型的、高階張量的SVD是非常耗時的。為了解決該問題，本文提出了張量R-Tucker分解，通過在張量每個階矩陣化[18]后使用RSVD構(gòu)造因子矩陣，并保持核心張量的計算方法保持不變。即在HOOI階段通過RSVD構(gòu)造左奇異向量對HOSVD階段得到的因子矩陣進行不斷迭代，提高張量Tucker分解的計算速度。具體方法如下：

在HOSVD階段，設(shè)? 為一個張量，其在第n（n=1，2…N）階進行張量矩陣化記作X（n），對X（n）進行RSVD得到因子矩陣 ，由因子矩陣即可得張量得核心張量g，計算公式為：

算法1是該算法的具體實現(xiàn)過程。一般將HOSVD算法得到的因子矩陣作為HOOI算法的初始化，在HOOI算法階段，張量R-Tucker分解可通過求解下式得到最優(yōu)解：

在最優(yōu)解時，核心張量滿足：

上述最小值問題可等價為以下最大值問題，文獻[19，20]對該過程進行了詳細(xì)的說明：

上述目標(biāo)函數(shù)寫作矩陣形式，即：

上式中使用RSVD將A（n）設(shè)為W的n個左奇異向量， 表示Kronecker積，通過求解以上優(yōu)化問題，即可得到R-Tucker分解的最終結(jié)果。

算法2 R-HOSVD

preocedure HOSVD（X， R1， R2，…RN）

for n=1， …N‘do

generate a gaussian random matrix

Z=X（n）P

QR decompositionZ=QR

Y=QTX

SVD decomposition

end for

return g， A（1）， A（2）， …A（N）

end procedure

算法2描述了該算法的具體實現(xiàn)過程。

算法2 R-HOOI

preocedure HOOI（X， R1， R2，…RN）

initialize? for n=1，…Nusing HOSVD

repeat

for n=1， …N‘do

generate a gaussian random matrix

Z=Y（N）P

QR decompositionZ=QR

Y=QTY（N）

SVD decomposition

end for

until fit ceases to improve or maximum iterations exhausted

return g， A（1）， A（2）， …A（N）

end procedure

2? R-Tucker分解在BCI信號分類中的應(yīng)用

MI-EEG數(shù)據(jù)分類研究的整體步驟如圖2所示。主要分為三個模塊：數(shù)據(jù)預(yù)處理模塊、特征提取模塊和特征分類模塊。內(nèi)容有：

（1）預(yù)處理模塊。主要包括讀取數(shù)據(jù)、濾波、重參考和分段；

（2）特征提取模塊。對預(yù)處理后的二維EEG信號（channel×time）進行時頻分析并建立起張量結(jié)構(gòu)（channel×frequency×time），在此基礎(chǔ)上使用基于隨機奇異值分解的Tucker分解進行特征提?。?/p>

（3）特征分類模塊。使用支持向量機對提取的特征進行分類。

Ⅰ是預(yù)處理模塊。Ⅱ是特征提取模塊。Ⅲ是特征分類模塊。輸入的運動想象原始信號有3個采樣通道，750個時間采樣點。

2.1? 數(shù)據(jù)集

本文的實驗數(shù)據(jù)來自2008年國際BCI競賽的數(shù)據(jù)集BCI Competition 2008-Grazdata set B（BCICIV2b）[21]，分為訓(xùn)練數(shù)據(jù)和測試數(shù)據(jù)，整個實驗共有9個受試者參與，所有受試者都擁有正常視力且為右利手。每個受試者做5組實驗，其中前2組實驗是無反饋的，后3組實驗是帶反饋的；前3組是訓(xùn)練數(shù)據(jù)，后2組是測試數(shù)據(jù)。每個實驗包括120個試次（60次左手運動想象，60次右手運動想象）。實驗采用AgCl電極，數(shù)據(jù)從國際標(biāo)準(zhǔn)10-20導(dǎo)聯(lián)系統(tǒng)的C3、Cz和C4三個通道獲得，采樣頻率為250 Hz，得到的腦電數(shù)據(jù)已經(jīng)進行帶工頻陷波處理和帶通濾波（0.5～100 Hz）。

2.2? 數(shù)據(jù)預(yù)處理

原始的EEG時域圖如圖3所示。預(yù)處理可以從包含大量噪聲的原始EEG信號中獲取較為純凈的MI-EEG成分，本文使用MNE-Python[22]對原始EEG數(shù)據(jù)進行預(yù)處理，首先利用0.5～45 Hz的前置帶通濾波去除包括線性漂移，心電，肌電在內(nèi)的大部分偽跡，濾波后的EEG時域圖如圖4所示。然后對連續(xù)的數(shù)據(jù)進行分段（Epoch），提取出執(zhí)行運動想象任務(wù)時（任務(wù)開始0～3 s）的MI-EEG信號，由于MI-EEG信號產(chǎn)生的腦電節(jié)律變化主要集中在mu節(jié)律（8～13 Hz）和beta節(jié)律（13～30 Hz），所以接下來對分段后的數(shù)據(jù)進行8～30 Hz的后置帶通濾波。

2.3? 構(gòu)建張量

EEG時域信號是以矩陣的形式存儲的時間、空間上的二維數(shù)據(jù)（channel×time=3×750），為了包含空間、時間和頻率的多階信息，須將EEG信號通過時頻變換建立起張量結(jié)構(gòu)。本文中使用連續(xù)小波變換（Continuous Wavelet Transform， CWT）[23]進行時頻變換。

對于二階EEG信號樣本X（c，t）（代表通道c和時間點t上的EEG信號值），通過計算和小波函數(shù)ω（ f，t）的卷積的幅值，可得三階（channel×frequency×time=3×23×750）的張量數(shù)據(jù)X c×f×t，X（c， f， t）的含義是在通道c、頻率f時間點t上的能量。計算公式為：

X（c， f， t）=‖ω（ f， t）×X（c，t）‖

使用小波函數(shù)為復(fù)Morlet小波，該函數(shù)為：

選擇中心頻率Ω=1，帶寬參數(shù)σ=2作為其母小波函數(shù)，因為它被廣泛應(yīng)用在EEG數(shù)據(jù)的處理上并且取得了很好的效果[24]。

2.4? 提取特征

對于每個MI-EEG數(shù)據(jù)X c×f×t進行基于RSVD的Tucker分解，分解形式如下所示：

其中X為原始運動想象EEG構(gòu)造的張量數(shù)據(jù)，A，B，C每個數(shù)據(jù)樣本在每一階的投影矩陣，矩陣之間相互正交，也可以被認(rèn)為是數(shù)據(jù)在每階的主分量。g為核心張量作為提取的特征，由原始EEG數(shù)據(jù)通過投影矩陣得到。

2.5? 特征分類

支持向量機（Support Vector Machines， SVM）是一種基于統(tǒng)計學(xué)習(xí)理論的模式識別方法，SVM建立在VC（Vapnik-Chervonenkis）維理論和結(jié)構(gòu)風(fēng)險最小化原理的基礎(chǔ)上，根據(jù)有限樣本信息在模型復(fù)雜度和學(xué)習(xí)能力之間尋求最佳折中，能較好地解決小樣本、非線性、高維度分類問題[25]。因此，本文采用支持SVM進行特征分類。

同一被試者數(shù)據(jù)的80%作為訓(xùn)練集，剩余20%作為測試集，在訓(xùn)練樣本上使用10折交叉驗證進行SVM分類器的參數(shù)選擇。針對MI-EEG信號非線性的特點，選擇高斯徑向基函數(shù)（Gaussian Radial Basis Function， RBF）核函數(shù)，并使用網(wǎng)格搜索法對懲罰因子C和核參數(shù)gamma進行參數(shù)尋優(yōu)。

2.6? 評估指標(biāo)

本文以基于RSVD的張量Tucker分解的時間和提取特征的分類準(zhǔn)確率為評價指標(biāo)。計算公式為：

Time=end time-start time

其中，Time是完成Tucker分解的用時，end time是Tucker分解結(jié)束時間，start time是Tucker分解開始時間。TP是真陽性數(shù)，TN是真陰性數(shù)，F(xiàn)P是假陽性數(shù)，F(xiàn)N是假陰性數(shù)。

3? 結(jié)果分析與討論

EEG信號是非平穩(wěn)的時域信號且主要特征反映在頻域上，因此對于多通道EEG信號，為了包含空間、時間和頻率的多維信息，本文首先通過時頻分析對EEG信號建立三階（channel×frequency×time）張量結(jié)構(gòu)，然后使用R-Tucker分解得到其核心張量將其作為特征，最后使用SVM進行特征分類。在對BCICIV2b數(shù)據(jù)集的實驗驗證了本文方法的有效性，取得了較為可信和穩(wěn)定的結(jié)果。

為了驗證張量R-Tucker分解對于分解速度的有效性，本文對張量Tucker分解和張量R-Tucker分解在不同Tucker秩情況下對分解時間進行比較。如表1所示，列出了BCICIV2b數(shù)據(jù)集在一組實驗中每個被試者張量Tucker分解和張量R-Tucker分解特征提取時間。由表1可知當(dāng)Tucker秩為[3，23，100]時，R-Tucker分解速度提升了63%；當(dāng)Tucker秩為[3，23，250]時，R-Tucker分解的平均時間提升了41%；當(dāng)Tucker秩為[3，23，500]時，R-Tucker分解的平均時間提升了22%。而且由表可知，對于每一個被試者，在相同Tucker秩的情況下，R-Tucker分解的效率均高于Tucker分解。

造成上述結(jié)果的原因是當(dāng)Tucker秩為[3，23，100]進行Tucker分解時，Tucker分解會對張量矩陣化后69×750大小的矩陣進行奇異值分解，然后選取分解結(jié)果的前100個左奇異向量構(gòu)成因子矩陣；而R-Tucker分解通過生成一個750×100大小服從高斯分布的隨機投影矩陣，由RSVD算法可知，通過投影矩陣與張量數(shù)據(jù)進行矩陣乘法會生成一個69×100大小的矩陣，然后按照RSVD算法的步驟得到的左奇異向量構(gòu)成的矩陣即為因子矩陣。由于Tucker分解需要對一個69×750大小的矩陣進行分解，而R-Tucker分解僅需對一個69×100大小的矩陣分解，因此R-Tucker分解的效率高于Tucker分解的效率，且分解時間提升了63%。同理，當(dāng)Tucker秩為[3，23，250]和[3，23，500]時，通過RSVD生成69×250和69×500大小的矩陣相比于69×750大小的矩陣，R-Tucker分解的的效率更高，同時這也導(dǎo)致Tucker秩越大，R-Tucker分解的提升效率越低。

通過上述分析可得如下結(jié)論，Tucker分解和R-Tucker分解相比，R-Tucker分解速度更快；對于同一個張量數(shù)據(jù)，它的Tucker秩越小，R-Tucker分解提升的效率越高。

為了驗證張量R-Tucker分解提取特征的分類性能，本文對張量R-Tucker分解提取的特征使用SVM進行了單被試分類實驗，由于Tucker秩為[3，23，100]和[3，23，250]的張量Tucker分解生成的核心張量過于小，丟失的與分類相關(guān)的信息較多，因此選取Tucker秩為[3，23，500]進行張量R-Tucker分解并進行分類。將其與連續(xù)小波變換得到的張量作為特征（表中為CWT）和文獻[26]基于傳統(tǒng)矩陣方法的分類結(jié)果進行比較，結(jié)果如表2所示。

表2列出不同方法每個被試者的準(zhǔn)確率以及BCICIV2b數(shù)據(jù)集的平均準(zhǔn)確率，使用CWT特征的單被試平均準(zhǔn)確率為62.07%，而張量R-Tucker分解后特征的單被試平均分類準(zhǔn)確率為80.93%，平均分類準(zhǔn)確率提高了18.86%；文獻[26]方法的單被試平均準(zhǔn)確率為70.81%，由表2可知，除了被試者S4和S5的分類準(zhǔn)確率略有下降，其余被試者的分類準(zhǔn)確率均顯著高于文獻[26]提出的方法，且平均準(zhǔn)確率提高了10.12%。

造成上述結(jié)果的原因如下：由R-Tucker分解得到的特征的分類準(zhǔn)確率高于分解前特征的準(zhǔn)確率可知，R-Tucker分解前張量數(shù)據(jù)的大小為3×23×750，R-Tucker分解后得到的核心張量的大小為3×23×500，說明經(jīng)過R-Tucker分解后得到的核心張量可以獲取原始數(shù)據(jù)潛藏的更緊湊的模式表示，有效對EEG信號進行降噪，從而提高了數(shù)據(jù)的信噪比。由R-Tucker分解得到的特征的分類準(zhǔn)確率高于文獻[26]可知，傳統(tǒng)方法通常將張量結(jié)構(gòu)的數(shù)據(jù)展開為矩陣的形式進行處理，這種展開不可避免地丟掉了一些維度內(nèi)與維度間可能存在的信息和交互作用，而張量R-Tucker分解能夠獲取EEG信號時間、空間和頻率多階的信息以及它們之間的關(guān)系，可以更有效的提取隱含在數(shù)據(jù)中的多階特征信息，從而獲取具有判別能力的特征。

通過上述分析可得如下結(jié)論，張量R-Tucker分解可以融合多個模態(tài)的信息使其提取的特征更具有判別能力，而且對于信噪比低的EEG數(shù)據(jù)可以提升其信噪比。

為了進一步對上述結(jié)果進行分析，本文對R-Tucker和CWT獲取的特征進行了t分布隨機領(lǐng)域嵌入（T-SNE）分析。T-SNE可以將高維數(shù)據(jù)投影到二維散點圖上，被廣泛用于評價特征向量的判別能力。評價標(biāo)準(zhǔn)如下：一個類的實例可以從其他類的實例中分離出的越多，相關(guān)特征的表現(xiàn)就越好。

如圖5和圖6所示，分別為受試者S1到S9在CWT和R-Tucker分解后的類別標(biāo)簽映射T-SNE圖中，其中藍(lán)色和綠色分別代表MI-EEG信號的2類別標(biāo)簽（左手和右手）。由圖可知，受試者S2、S3、S6和S9在T-SNE映射后兩類標(biāo)簽重疊部分較多，這導(dǎo)致CWT特征在這幾位受試者分類準(zhǔn)確率較低，而這幾位受試者在R-Tucker分解后經(jīng)過T-SNE映射后，兩類標(biāo)簽重疊部分有一定程度的降低，在經(jīng)過SVM核函數(shù)的映射之后，在高維空間的分類準(zhǔn)確率更高。而且通過觀察圖5和圖6，9位受試者在R-Tucker分解后兩類標(biāo)簽的重疊部分均明顯少于CWT，這也導(dǎo)致R-Tucker分解特征的平均分類準(zhǔn)確率高于CWT特征的分類準(zhǔn)確率。

如圖7所示，對BCICIV2b數(shù)據(jù)集受試者使用不同方法的分類結(jié)果使用T檢驗進行顯著性差異分析?；赥ucker分解提取特征的分類準(zhǔn)確率與文獻[26]存在顯著性差異（t=3.364 0，p＜0.01）。而基于Tucker分解提取特征相比于Tucker分解前的特征的分類準(zhǔn)確率也存在顯著性差異（t=13.986 2，p＜0.001）。而且由圖可知，基于R-Tucker分解方法的單被試分類準(zhǔn)確率方差小于其他兩種方法，說明其分類結(jié)果更為平穩(wěn)。

綜上所述，基于RSVD的張量R-Tucker分解不僅可以提高分解效率，而且能夠獲取EEG信號時間、空間和頻率多個模態(tài)的信息以及它們之間的關(guān)系，從而獲取更緊湊、具有判別能力的特征，而且通過比較本文方法和文獻[26]方法的平均準(zhǔn)確率以及方差，說明R-Tucker分解提取的特征具有更好的判別能力以及魯棒性。

4? 結(jié)? 論

本文對于兩類（左手和右手）MI-EEG數(shù)據(jù)的分類主要有三個模塊組成，分別是預(yù)處理模塊、特征提取模塊和特征分類模塊。預(yù)處理模塊主要完成了數(shù)據(jù)的濾波和分段；特征提取模塊對預(yù)處理后的數(shù)據(jù)建立張量結(jié)構(gòu)并使用張量Tucker分解提取特征；分類模塊使用SVM對特征提取模塊提取的特征進行分類。本文的主要創(chuàng)新是使用RSVD改進了基于HOSVD-HOOI的張量Tucker分解，并將其用于兩類（左手和右手）MI-EEG數(shù)據(jù)的特征提取。實驗結(jié)果表明，張量R-Tucker分解與張量Tucker分解相比，張量R-Tucker分解速度得到顯著的提升，并且使用張量R-Tucker分解提取的特征擁有更好的分類性能。

參考文獻：

[1] LOTTE F，BOUGRAIN L，CICHOCKI A，et al. A Review of Classification Algorithms for EEG-Based Brain–Computer Interfaces：A 10 Year Update [J].Journal of Neural Engineering，2018，15（3）：031005.

[2] NORCIA A M，APPELBAUM L G，ALES J M，et al. The Steady-State Visual Evoked Potential in Vision Research：A Review [J].Journal of Vision，2015，15（6）：4.

[3] LUCK S J.An Introduction to the Event-Related Potential Technique [M].Bradford：A Bradford Book，2014.

[4] PFURTSCHELLER G，NEUPER C. Motor Imagery and Direct Brain-Computer Communication [J].Proceedings of the IEEE，2001，89（7）：1123-1134.

[5] PFURTSCHELLER G.Event-Related EEG/MEG Synchronization and Desynchronization：Basic Principles [J].Neurophysiol Clin，1999，110（11）：1842-1857.

[6] MA Y，ZHENG L S，YI Z K，et al. Short-Time Fourier Transform Covariance and Selection，A Feature Extraction Method for Binary Motor Imagery Classification [C]//2021 IEEE International Conference on Real-Time Computing and Robotics （RCAR）.Xining：IEEE，2021：1316-1322.

[7] YOU Y，CHEN W Z，ZHANG T. Motor Imagery EEG Classification Based on Flexible Analytic Wavelet Transform [J].Biomedical Signal Processing and Control，2020，62（615）：102069.

[8] PHAN A H，CICHOCKI A. Tensor Decompositions for Feature Extraction and Classification of High Dimensional Datasets [J].Nonlinear Theory and Its Applications，IEICE，2010，1（1）：37-68.

[9] CONG F Y，LIN Q H，KUANG L D，et al. Tensor Decomposition of EEG Signals：A Brief Review [J].Journal of Neuroscience Methods，2015，248：59-69.

[10] ZHANG Q Z，GUO B，KONG W Z，et al. Tensor-Based Dynamic Brain Functional Network for Motor Imagery Classification [J].Biomed Signal Process Control，2021，69.

[11] RO?KOV? Z，ROSIPAL R，SEIFPOUR S. Tucker Tensor Decomposition of Multi-session EEG Data [C]//29th International Conference on Artificial Neural Networks.Bratislava：Springer，2020：115-126.

[12] KRISHNAN K K，SOMAN K P. Canonical Polyadic Decomposition of EEG Image Tensor for BCI Applications [M]//ICT Systems and Sustainability.Singapore：Springer，2022：819-826.

[13] 李潔.多模態(tài)腦電信號分析及腦機接口應(yīng)用 [D].上海：上海交通大學(xué)，2009.

[14] BRUNTON S L，KUTZ J N. Data-Driven Science and Engineering：Machine Learning，Dynamical Systems，and Control [M].Cambridge：Cambridge University Press，2019.

[15] HALKO N，MARTINSSON P-G，SHKOLNISKY Y，et al.An Algorithm for the Principal Component Analysis of Large Data Sets [J].SIAM Journal on Scientific Computing，2011，33（5）：2580-2594.

[16] ZHANG J N，SAIBABA A K，KILMER M E，et al. A Randomized Tensor Singular Value Decomposition Based on the T-Product [J/OL].arXiv：1609.07086 [math.NA].[2022-08-09].https：//arxiv.org/abs/1609.07086.

[17] KOLDA T G，BADER B W. Tensor Decompositions and Applications [J].SIAM Review，2009，51（3）：455-500.

[18] DE L L，VANDEWALLE J，DE M B. A Multilinear Singular Value Decomposition [J].SIAM Journal on Matrix Analysis and Applications，2000，21（4）：1253-1278.

[19] ANDERSSON C A，BRO R. Improving the Speed of Multi-Way Algorithms：Part I.Tucker3 [J].Chemometrics and Intelligent Laboratory Systems，1998，42（1-2）：93-103.

[20] LATHAUWER L D，MOOR B D，VANDEWALLE J. On the Best Rank-1 and Rank-（R1，R2，...，Rn） Approximation of Higher-Order Tensors [J].SIAM Journal on Matrix Analysis and Applications，2006，21（4）：1324-1342.

[21] LEEB R，BRUNNER C，MULLER-PUTZ G R，et al.BCI Competition 2008 Graz data set B[EBOL].[2022-08-12].https：//www.bbci.de/competition/iv/.

[22] ALEXANDRE G. MEG and EEG Data Analysis with MNE-Python [J].Front Neurosci，2013，7（7）：267.

[23] COHEN M X. A Better Way to Define and Describe Morlet Wavelets for Time-Frequency Analysis [J].Neuroimage，2019，199：81-86.

[24] LE S，GORDON E，GYSELS E. Phase Synchrony Rate for the Recognition of Motor Imagery in Brain-Computer Interface [C]//Advances in Neural Information Processing Systems.Vancouver：DBLP，2005：1265-1272.

[25] BURGES C J C. A Tutorial on Support Vector Machines for Pattern Recognition [J].Data Mining and Knowledge Discovery，1998（2）：121-167.

[26] RAZA H，CECOTTI H，PRASAD G. A Combination of Transductive and Inductive Learning for Handling Non-Stationarities in Motor Imagery Classification [C]//2006 International Joint Conference on Neural Networks（IJCNN）.IEEE：Vancouver，2016：763-770.

作者簡介：張帥（1995—），男，漢族，山西忻州人，碩士研究生在讀，研究方向：智能數(shù)據(jù)分析與應(yīng)用。

收稿日期：2022-09-01