張睿琳 賈建 張瑞

摘要：該文提出一種基于可調(diào)Q因子小波變換的熵特征提取方法用于癲癇腦電的快速檢測(cè)。首先，利用TWQT算法對(duì)腦電信號(hào)進(jìn)行分解，計(jì)算分解后得到的小波子帶的中心頻率;隨后，選擇對(duì)應(yīng)癲癇發(fā)作時(shí)異常波段頻率范圍的子帶進(jìn)行小波重構(gòu)，并對(duì)重構(gòu)信號(hào)提取排序熵和規(guī)則性指數(shù)等熵特征;最后，選擇包括支持向量機(jī)、隨機(jī)森林、極端梯度提升在內(nèi)的多種分類器進(jìn)行分類實(shí)驗(yàn)。實(shí)驗(yàn)中采用德國(guó)波恩大學(xué)癲癇研究中心的公開數(shù)據(jù)集和UPenn and Mayo Clinic′s Seizure Detection Challenge數(shù)據(jù)集進(jìn)行驗(yàn)證，所提出的方法針對(duì)小樣本數(shù)據(jù)集可以達(dá)到99.3％的準(zhǔn)確率，100％的靈敏度以及98.6％的特異度，體現(xiàn)了該算法的有效性。

關(guān)鍵詞：腦電信號(hào);可調(diào)Q因子小波變換;排序熵;規(guī)則性指數(shù);支持向量機(jī);隨機(jī)森林; 極端梯度提升

中圖分類號(hào)：TP391.4

DOI：10.16152/j.cnki.xdxbzr.2020-04-008開放科學(xué)（資源服務(wù)）標(biāo)識(shí)碼（OSID）：

Rapid detection of epileptic EEG based on TWQT and entropy

ZHANG Ruilin1，2， JIA Jian1，2， ZHANG Rui1，2

（1.School of Mathematics， Northwest University， Xi′an 710127， China;

2.Medical Big Data Research Center， Northwest University， Xi′an 710127， China）

Abstract： In this paper， an entropy feature extraction method based on tunable Q factor wavelet transform is proposed. First， using TWQT algorithm to decompose EEG signals，calculate the center frequency of each wavelet subband， then choose corresponding sub-bands for reconstructing and extract reconstructed signals′ permutation entropy and order index as the characteris-tic. In the end， choose? a variety of classifieiers including support vector machine（SVM）， random forest and XGBoost to classify. The proposed classification method is applied to the detection of epileptic EEG signals， and verified by the open dataset of the Epilepsy Research Center of Bonn University in Germany and UPenn and Mayo Clinic′s Seizure Detection Challenge dataset. The accuracy is 99.3%， sensitivity is 100% and specificity is 98.6%， which shows the effectiveness of the algorithm.

Key words： EEG; TWQT; permutation entropy; order index; SVM; random forest; XGBoost

癲癇是由大腦內(nèi)神經(jīng)元異常放電引起的慢性腦部疾病，通常表現(xiàn)為感覺(jué)、知覺(jué)以及行為上不同程度的障礙，具有突發(fā)性、重復(fù)性、自限性、暫時(shí)性等特點(diǎn)［1］。癲癇不僅嚴(yán)重?fù)p害了患者的身心健康，對(duì)患者的日常生活造成了嚴(yán)重的影響，還給家庭和社會(huì)帶來(lái)了沉重的經(jīng)濟(jì)負(fù)擔(dān)。迄今為止，腦電圖（Electroencephalogram，EEG）仍然是研究癲癇信號(hào)的主要工具，被用于進(jìn)行癲癇相關(guān)的神經(jīng)紊亂疾病監(jiān)測(cè)、診斷和治療。腦電圖會(huì)監(jiān)測(cè)到大量癲癇發(fā)作時(shí)的異常腦電活動(dòng)數(shù)據(jù)，通常由醫(yī)生對(duì)患者腦電圖進(jìn)行肉眼觀察完成疾病診斷，但這種方式會(huì)受到腦電信號(hào)的非癲癇瞬變等干擾，因而效率較低且非常耗時(shí)。此外，在長(zhǎng)時(shí)程的主體分析過(guò)程中，主觀性也會(huì)對(duì)診斷結(jié)果的準(zhǔn)確性造成影響。因此，癲癇發(fā)作的自動(dòng)檢測(cè)技術(shù)研究具有重要的意義［2］。

傳統(tǒng)的信號(hào)處理方法包括快速傅里葉變換（FFT）、短時(shí)傅里葉變換（STFT）、幅頻分析等方法，可以在頻域準(zhǔn)確、高效地分析信號(hào)，但是在時(shí)域上缺乏分辨能力。而小波變換既包含STFT的局部化思想，又克服了其主要缺陷，如窗口大小不能隨頻率變化，成為一種對(duì)信號(hào)進(jìn)行時(shí)頻分析的有效新方法［3］。但進(jìn)行小波變換時(shí)需要選定合適的小波基函數(shù)，由單個(gè)基函數(shù)導(dǎo)出的小波函數(shù)無(wú)法準(zhǔn)確表現(xiàn)信號(hào)的局部特征，因此在重構(gòu)時(shí)將失去原有的時(shí)域特征［4］。

可調(diào)Q因子小波變換（tunable Q-factor wavelet transform，TWQT）是一種參數(shù)可以根據(jù)小波的振蕩特性與特征波形的振蕩特性進(jìn)行調(diào)節(jié)的小波變換方法［5］，可以可靠地表征信號(hào)在時(shí)域上的稀疏性和非平穩(wěn)性，適用于分析包括腦電信號(hào)在內(nèi)的振蕩信號(hào)［6］。幾十年來(lái)，通過(guò)提取腦電圖特征的癲癇發(fā)作檢測(cè)方法已經(jīng)得到了廣泛的應(yīng)用。Patidar等采用可調(diào)Q因子小波變換分解腦電信號(hào)，以Kraskov熵為特征，利用最小二乘支持向量機(jī)（LS-SVM）進(jìn)行分類［7］;Peker等采用雙樹復(fù)小波變換（DTCWT）對(duì)癲癇腦電信號(hào)進(jìn)行處理和分類［8］;Sharma等采用時(shí)頻靈活分析小波變換（ATFFWT）和分形維數(shù)（FD）進(jìn)行癲癇發(fā)作檢測(cè)［9］;Jia等基于完全集合經(jīng)驗(yàn)?zāi)Ｊ椒纸猓–EEMD）提取譜特征并利用隨機(jī)森林進(jìn)行分類檢測(cè)［10］;Mutlu采用希爾伯特震動(dòng)分解（HVD）提取最大能量分量的瞬時(shí)頻率值作為L(zhǎng)S-SVM的特征以檢測(cè)癲癇發(fā)作［11］。以上方法基于傳統(tǒng)的機(jī)器學(xué)習(xí)方法進(jìn)行實(shí)驗(yàn)，均取得了十分有效的成果。除此之外，許多研究通過(guò)深度學(xué)習(xí)實(shí)現(xiàn)癲癇發(fā)作自動(dòng)檢測(cè)，深度學(xué)習(xí)會(huì)自動(dòng)發(fā)現(xiàn)并學(xué)習(xí)分類輸入所需的特征，并基于不同的深度神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)實(shí)現(xiàn)，例如完全連接的神經(jīng)網(wǎng)絡(luò)（FCNN）［12］，卷積神經(jīng)網(wǎng)絡(luò)（CNN）［13］以及循環(huán)神經(jīng)網(wǎng)絡(luò)（RNN）［14］。這種深度神經(jīng)網(wǎng)絡(luò)模型的實(shí)現(xiàn)需要依賴大量的數(shù)據(jù)，并且需要對(duì)模型的參數(shù)進(jìn)行調(diào)整和選擇，對(duì)于小樣本數(shù)據(jù)集的分類檢測(cè)具有一定的局限性。因此，本文研究了針對(duì)小樣本數(shù)據(jù)集的癲癇腦電快速檢測(cè)方法。

基于TWQT對(duì)癲癇腦電信號(hào)的特征提取，本文主要研究以下3個(gè)方面的內(nèi)容：①利用TWQT對(duì)腦電信號(hào)進(jìn)行分解，得到相應(yīng)的子帶;②根據(jù)癲癇發(fā)作時(shí)異常波的頻率對(duì)應(yīng)范圍，選擇相適應(yīng)的子帶進(jìn)行小波重構(gòu)，并且在重構(gòu)后的腦電信號(hào)上提取熵特征，采用排序熵（permutation entropy，PE）和規(guī)則性指數(shù)（order index，OI）研究癲癇患者的腦電信號(hào)復(fù)雜度;③通過(guò)多種分類器對(duì)腦電信號(hào)進(jìn)行分類，包括支持向量機(jī)、隨機(jī)森林、極端梯度提升。在采用德國(guó)波恩大學(xué)的癲癇研究中心的公共數(shù)據(jù)集和UPenn and Mayo Clinic′s Seizure Detection Challenge數(shù)據(jù)集的驗(yàn)證結(jié)果表明本文所述分類方法針對(duì)小樣本數(shù)據(jù)集可以達(dá)到快速準(zhǔn)確的檢測(cè)結(jié)果，其分類準(zhǔn)確率分別為99.3％和99％，并引入靈敏度、特異度和AUC值作為衡量指標(biāo)與現(xiàn)有方法進(jìn)行對(duì)比。

1 可調(diào)Q因子小波變換和熵特征

1.1 可調(diào)Q因子小波變換

可調(diào)Q因子小波變換是一種完全離散的小波變換［15］，其主要參數(shù)是品質(zhì)因子（Q）、冗余度（r）和分解層數(shù)（J），如表1所示。

Q因子主要影響子波的振蕩特性，Q的值越大，小波振蕩頻次越高，Q值越小，小波振蕩頻次越低;參數(shù)r是TWQT在無(wú)限多層J計(jì)算時(shí)的冗余度，在實(shí)驗(yàn)中，為了能夠使得分解后的小波具有局域化時(shí)域特征的能力，通常選擇r≥3;小波變換的階數(shù)（或?qū)訑?shù)）用J表示，TQWT的分解原理與傳統(tǒng)的小波變換相似，采用迭代的雙通道濾波器對(duì)信號(hào)進(jìn)行分解，每個(gè)濾波器組的低通輸出作為輸入到后續(xù)的過(guò)濾器組。每個(gè)輸出信號(hào)構(gòu)成小波變換的一個(gè)子帶，最終得到J+1個(gè)子帶，其中包括高通濾波器輸出的每個(gè)濾器組的信號(hào)，以及最終濾波器組的低通濾波器輸出的信號(hào)，它表征了信號(hào)中所含低頻分量的細(xì)化程度。例如，一個(gè)3階TWQT如圖1所示［16］。

，其中，T為葉結(jié)點(diǎn)的數(shù)量。結(jié)構(gòu)部分q（X^）的輸出是葉結(jié)點(diǎn)編號(hào)d，它的作用是將輸入W^映射到編號(hào)為d的葉結(jié)點(diǎn)。葉結(jié)點(diǎn)權(quán)重部分就是每個(gè)葉結(jié)點(diǎn)的值，它的作用是輸出編號(hào)為d的葉結(jié)點(diǎn)的值wd。

XGBoost同樣體現(xiàn)了集成學(xué)習(xí)的思想，并且在處理過(guò)擬合問(wèn)題時(shí)有著良好的應(yīng)對(duì)方法：一種方法稱作Shrinkage，目的在于適當(dāng)?shù)厝趸靠脴涞挠绊?，?duì)于決策樹的每次迭代，添加一個(gè)縮減因子以降低葉子結(jié)點(diǎn)的權(quán)重比例，這樣使得之后決策樹的產(chǎn)生可以伴隨著更大的空間優(yōu)化可能。第二種方法稱作column subsampling，與隨機(jī)選擇的方法類似，在森林中選取部分特征值來(lái)構(gòu)建決策樹，這種選取方式可以是分層隨機(jī)抽樣，在每一層的分類結(jié)點(diǎn)之前先計(jì)算每一個(gè)隨機(jī)選取部分特征值的信息增益；也可以是在一棵決策樹的生成前先隨機(jī)選取部分特征值，然后計(jì)算每一個(gè)分類結(jié)點(diǎn)處的信息增益。

2 實(shí)驗(yàn)結(jié)果

在癲癇的發(fā)作檢測(cè)中，腦電圖檢查對(duì)其病因和發(fā)作類型的判斷具有很大價(jià)值，其中最常用的是頭皮腦電圖，它將電極夾放置在頭皮上從而記錄腦細(xì)胞在活動(dòng)時(shí)的電位變化。

本文所使用的腦電圖公共數(shù)據(jù)來(lái)源于波恩大學(xué)的癲癇研究中心。它由5組數(shù)據(jù)集組成，即Z，O，N，F(xiàn)，S，每組含100個(gè)屬于單個(gè)實(shí)驗(yàn)者的單道EEG片段，每一單道段持續(xù)23.6 s?；趪?guó)際標(biāo)準(zhǔn)的10-20腦電圖電極放置系統(tǒng)，設(shè)置集Z和O記錄清醒和放松的信號(hào)，集F和N包括電極置于致癲癇區(qū)時(shí)產(chǎn)生的信號(hào)。集S中的數(shù)據(jù)由顱內(nèi)的、以及植入顳部和基底部的大腦皮層區(qū)域放置電極的信號(hào)組成。集S中的腦電圖數(shù)據(jù)對(duì)應(yīng)于癲癇發(fā)作，而在集F和N則對(duì)應(yīng)于未發(fā)作期。每個(gè)腦電圖信號(hào)段采用具有12位分辨率和173.61 Hz數(shù)字化的128通道放大器系統(tǒng)采樣。最后，每個(gè)片段的采樣長(zhǎng)度為173.61×23.6≈4 097，對(duì)應(yīng)的帶寬為86.8 Hz。本文利用腦電圖信號(hào)的子集S，F(xiàn)評(píng)估所建議方法的表現(xiàn)，所提出方法的檢測(cè)流程如圖2所示。

首先采用TWQT對(duì)原始的EEG數(shù)據(jù)進(jìn)行分解，分解時(shí)的參數(shù)設(shè)置為Q=1.5，r=3，J=6。分

解后的結(jié)果如圖3所示，最上方為某一原始EEG信號(hào)，下方的信號(hào)為分解后得到的7條小波子帶，可以看到，子帶中包含了不同的特征波形。

對(duì)分解后的子帶選擇具有顯著特征的子帶進(jìn)行小波重構(gòu)。這里挑選的準(zhǔn)則為每層小波子帶的中心頻率。TWQT分解第j層對(duì)應(yīng)的中心頻率fc表示為

fc=αj2-β4αfs。（20）

計(jì)算得到每層小波子帶所對(duì)應(yīng)的中心頻率如表2所示。

因在癲癇發(fā)作期大腦異常波段對(duì)應(yīng)的頻率范圍是3～25 Hz，因此選擇包含在這一頻率范圍之內(nèi)的第3～6層小波子帶對(duì)EEG信號(hào)進(jìn)行重構(gòu)。重構(gòu)后的腦電信號(hào)長(zhǎng)度與原始信號(hào)相同，之后計(jì)算其排序熵和規(guī)則性指數(shù)作為特征。嵌入維度和延遲時(shí)間為m=5，r=1。

圖4為重構(gòu)信號(hào)的排序熵與規(guī)則性指數(shù)的箱型圖。可以看出，癲癇發(fā)作期的PE值低于癲癇發(fā)作間期，OI值高于癲癇發(fā)作間期。這表明癲癇發(fā)作期的腦電序列復(fù)雜度更低。

接下來(lái)利用所提取特征進(jìn)行分類，這里考慮3種不同的輸入方式，一是以PE作為特征進(jìn)行分類，二是以O(shè)I作為特征進(jìn)行分類，三是以PE和OI構(gòu)成的特征向量進(jìn)行分類，其分類準(zhǔn)確率如表3所示。

由表3可見(jiàn)，由PE和OI構(gòu)成的特征向量作為分類指標(biāo)時(shí)準(zhǔn)確率（Accuracy）更高，通過(guò)比較不同分類器的分類結(jié)果可以看到，XGBoost分類器在以不同特征作為輸入時(shí)均產(chǎn)生了更好的結(jié)果。為了進(jìn)一步評(píng)估所述方法的實(shí)驗(yàn)性能，引入靈敏度（Sensitivity）、特異度（Specificity）和AUC（area under the curve）作為評(píng)估指標(biāo)［25］，表達(dá)式如下：

1） 靈敏度

Sensitivity=TPTP+FN。（21）

2） 特異度

Specificity=TNTN+FP。（22）

其中：TP代表正確分類的正類樣本數(shù)；TN代表正確分類的負(fù)類樣本數(shù)；FP代表被錯(cuò)誤分類的負(fù)類樣本數(shù)；FN代表被錯(cuò)誤分類的正類樣本數(shù)。

以XGBoost作為分類器，PE和OI共同作為特征輸入時(shí)，達(dá)到99.3％的最高分類準(zhǔn)確率，98.6％的特異度以及100％的靈敏度，0.996的AUC值也表明該分類方法具有較好的分類效果，對(duì)比實(shí)驗(yàn)數(shù)據(jù)如表4所示。

將本文提出的方法與現(xiàn)有的研究方法進(jìn)行對(duì)比，針對(duì)波恩大學(xué)癲癇研究中心的公開數(shù)據(jù)庫(kù)展開的部分研究數(shù)據(jù)如表5所示。本文所提方法得到了最高的準(zhǔn)確率99.3％和最高靈敏度100％，盡管Jia［10］提出方法的特異度略高于本文方法，但其準(zhǔn)確度卻低于本文方法?？梢钥闯鲎鳛橐环N基于小樣本和僅采用熵特征進(jìn)行分類的方法，本文所采用的方法在各度量標(biāo)準(zhǔn)上均取得了較好的性能。

為進(jìn)一步驗(yàn)證本文方法的可行性，在UPenn and Mayo Clinic′s Seizure Detection Challenge數(shù)據(jù)集上再次進(jìn)行實(shí)驗(yàn)，選取其中一名實(shí)驗(yàn)者的腦電數(shù)據(jù)，其中，包括70段癲癇發(fā)作時(shí)的EEG信號(hào)以及70段癲癇發(fā)作間隙的EEG片段，這些片段均為1 s段，采樣頻率為500Hz，對(duì)1s段進(jìn)行TWQT分解與小波重構(gòu)，參數(shù)設(shè)置為Q=1.5，r=1，J=9，即進(jìn)行9階TWQT分解，采用與上述實(shí)驗(yàn)相同的方式計(jì)算各個(gè)子帶的中心頻率，結(jié)果見(jiàn)表6。

選取包含于癲癇異常波段頻率范圍內(nèi)的第7～9條子帶進(jìn)行小波重構(gòu)，并在重構(gòu)后的信號(hào)上提取排序熵和規(guī)則性指數(shù)作為特征，參數(shù)設(shè)置為m=5，r=1。采用準(zhǔn)確率、特異度、靈敏度和AUC值作為評(píng)價(jià)指標(biāo)，并且將結(jié)果與先前在同樣數(shù)據(jù)集上實(shí)驗(yàn)的方法進(jìn)行比較，最終得到的結(jié)果如表7所示。 可以看出， 本文所采用方法的準(zhǔn)確度低于應(yīng)用CNN進(jìn)行分類的方法，高于應(yīng)用RNN進(jìn)行分類的方法，并且在靈敏度、特異度以及AUC值上取得了更好的結(jié)果。

相比神經(jīng)網(wǎng)絡(luò)這種需要大量數(shù)據(jù)來(lái)確定最優(yōu)參數(shù)和復(fù)雜網(wǎng)絡(luò)結(jié)構(gòu)的深度學(xué)習(xí)方法，本文通過(guò)TWQT分解和提取熵特征的方法取得了較好的分類結(jié)果，在利用PE和OI共同作為特征輸入時(shí)達(dá)到99％的最高準(zhǔn)確率以及97.9％的靈敏度，而特異度達(dá)到了100％，0.992的AUC值也表明該分類方法具有較好的分類效果。

3 結(jié) 論

近年來(lái)，在癲癇發(fā)作的自動(dòng)診斷中，通過(guò)機(jī)器學(xué)習(xí)的方法對(duì)正常腦電信號(hào)和癲癇發(fā)作間期腦電信號(hào)進(jìn)行分類取得了顯著的成果，而基于人工神經(jīng)網(wǎng)絡(luò)的深度學(xué)習(xí)方法提供了一種癲癇發(fā)作檢測(cè)的新模式，在不用提取特征的條件下便能得到很高的分類準(zhǔn)確率，但在應(yīng)用時(shí)需要通過(guò)大量數(shù)據(jù)確定最優(yōu)參數(shù)和網(wǎng)絡(luò)結(jié)構(gòu)。本文基于小樣本數(shù)據(jù)集和僅采用熵特征作為分類依據(jù)，引入靈活的TWQT對(duì)腦電信號(hào)進(jìn)行分解，并且根據(jù)癲癇異常波段對(duì)應(yīng)的頻率范圍選擇相應(yīng)的小波子帶進(jìn)行重構(gòu)［16］，提取重構(gòu)信號(hào)的排序熵和規(guī)則性指數(shù)作為特征，采用多種分類器進(jìn)行分類。選擇波恩大學(xué)癲癇研究中心的公共數(shù)據(jù)集進(jìn)行驗(yàn)證，最終的分類準(zhǔn)確率達(dá)到了99.3％。作為進(jìn)一步闡釋本文方法可行性的依據(jù)，采用UPenn and Mayo Clinic′s Seizure Detection Challenge數(shù)據(jù)集再次實(shí)驗(yàn)，最終達(dá)到99％的準(zhǔn)確率。

作為一種實(shí)現(xiàn)癲癇發(fā)作自動(dòng)檢測(cè)的方法，本研究不僅取得了有效的分類結(jié)果，并且表明癲癇發(fā)作間期和癲癇發(fā)作期的腦電信號(hào)在復(fù)雜度上具有顯著差異，與發(fā)作間期相比，癲癇發(fā)作期的EEG序列的PE值降低，OI值升高，說(shuō)明其腦活動(dòng)的復(fù)雜度降低。這種差異可能表現(xiàn)為癲癇患者的認(rèn)知障礙。本文提出的熵特征算法可以作為檢測(cè)癲癇發(fā)作的指標(biāo)，對(duì)于癲癇患者的早期診斷和有效治療都具有一定的指導(dǎo)意義。

參考文獻(xiàn)：

［1］ 孫殿榮， 鄒曉毅. 癲癇發(fā)病機(jī)制的研究進(jìn)展［J］.華西醫(yī)學(xué)， 2004，19（3）：524-525.

SUN D R，ZOU X Y.Advances in pathogenesis of epilepsy［J］.West China Medical Journal，2004，19（3）：524-525.

［2］ 黃瑞梅， 杜守洪， 陳子怡， 等. 基于支持向量機(jī)的癲癇腦電信號(hào)模式識(shí)別研究［J］.生物醫(yī)學(xué)工程學(xué)雜志， 2013，30（5）：919-924.

HUANG R M，DU S H，CHEN Z Y，et al.Study on pattern recognition of epileptic EEG based on support vector machine［J］.Journal of Biomedical Engineering，2013，30（5）：919-924.

［3］ 趙建林， 周衛(wèi)東， 劉凱， 等. 基于 SVM和小波分析的腦電信號(hào)分類方法［J］.計(jì)算機(jī)應(yīng)用與軟件， 2011， 28（5）：114-116.

ZHAO J L，ZHOU W D，LIU K，et al.Classification of EEG signals based on SVM and wavelet analysis［J］.Computer Applications and Software，2011，28（5）：114-116.

［4］ 陳萬(wàn)忠， 王曉旭， 張濤. 基于可調(diào) Q因子小波變換的識(shí)別左右手運(yùn)動(dòng)想象腦電模式研究［J］.電子與信息學(xué)報(bào)， 2019， 41（3）：530-536.

CHEN W Z，WANG X X，ZHANG T.Study on EEG pattern of left and right hand movement recognition based on tunable Q factor wavelet transform［J］.Journal of Electronics and Information Technology， 2019， 41（3）：530-536.

［5］ HE W P，ZI Y Y，CHEN B Q，et al.Automatic fault feature extraction of mechanical anomaly on induction motor bearing using ensemble super-wavelet transform［J］.Mechanical Systems and Signal Processing，2015， 54-55（5）：457-480.

［6］ HE W P， ZI Y Y， CHEN B Q， et al. Tunable Q-factor wavelet transform denoising with neighboring coefficients and its application to rotating machinery fault diagnosis［J］.Science China Technological Sciences，2013，56（8）：1956-1965.

［7］ PATIDAR S，PANIGRAHI T.Detection of epileptic seizure using Kraskov entropy applied on tunable-Q wavelet transform of EEG signals［J］.Biomedical Signal Processing and Control， 2017，34（4）：74-80.

［8］ PEKER M， SEN B H， DELEN D. A novel method for automated diagnosis of epilepsy using complex-valued classifiers［J］.IEEE Journal of Biomedical and Health Informatics， 2016，20（1）：108-118.

［9］ SHARMA M， PACHORI R B， RAJENDRA ACHARYA U. A new approach to characterize epileptic seizures using analytic time-frequency flexible wavelet transform and fractal dimension［J］.Pattern Recognition Letters， 2017， 94（15）：172-179.

［10］JIA J， GOPARAJU B， SONG J L， et al. Automated identification of epileptic seizures in EEG signals based on phase space representation and statistical features in the CEEMD domain［J］.Biomedical Signal Processing and? Control， 2017， 38：148-157.

［11］MUTLU A Y. Detection of epileptic dysfunctions in EEG signals using Hilbert vibration decomposition［J］.Biomedical Signal Processing and Control， 2018， 40（2）：33-40.

［12］JANG H J， CHO K O. Dual deep neural network-based classifiers to detect experimental seizures［J］.The Korean Journal of Physiology and Pharmacology， 2019，23（2）： 131-139.

［13］ULLAH I，HUSSAIN M， QAZIE U， et al.An automated system for epilepsy detection using EEG brain signals based on deep learning approach［J］.Expert Systems with Applications，2018，107（10）：61-71.

［14］HUSSEIN R， PALANGI H， WARD R K， et al. Optimized deep neural network architecture for robust detection of epileptic seizures using EEG signals［J］.Clinical Neurophysiology， 2019，130（1）：25-37.

［15］SELESNICK I W.Wavelet transform with tunable Q-factor［J］.IEEE Transactions on Signal Processing， 2011，59（8）：3560-3575.

［16］賀王鵬， 楊琳， 王芳， 等. 基于 TQWT的癲癇腦電信號(hào)的識(shí)別［J］.生物醫(yī)學(xué)工程研究， 2017， 36（4）：346-350.

HE W P，YANG L，WANG F，et al.Recognition of epileptic EEG based on TQWT［J］.Journal of Biomedical Engineering Research， 2017， 36（4）：346-350.

［17］CAO Y H， TUNG W W， GAO J B， et al. Detecting dynamical changes in time series using the permutation entropy［J］.Physical Review E， 2004， 70（4）：046217.

［18］LI D， LI X L， LIANG Z H， et al. Multiscale permutation entropy analysis of EEG recordings during sevoflurane anesthesia［J］.Journal of Neural Engineering， 2010， 7（4）：046010.

［19］BANDT C， POMPE B. Permutation entropy： A natural complexity measure for time series［J］.Physical Review Letters， 2002，88（17）：174102.

［20］OUYANG G X， DANG C Y， RICHARDS D A， et al. Ordinal pattern based similarity analysis for EEG recordings［J］.Clinical Neurophysiology， 2010，121（5）：694-703.

［21］李航. 統(tǒng)計(jì)學(xué)習(xí)方法［M］.北京： 清華大學(xué)出版社， 2012： 111-112.

［22］BREIMAN L.Random forests［J］.Machine Learning，2001，45（1）：5-32.

［23］SVETNIK V，LIAW A，TONG C，et al.Random forest： A classification and regression tool for compound classification and QSAR modeling［J］.Journal of Chemical in? Formation and Computer Sciences， 2003，43（6）：1947-1958.

［24］CHEN T Q， GUESTRIN C. XGBoost： A scalable tree boosting system［EB/OL］.2016： arXiv：1603.02754［cs.LG］.https：∥arxiv.org/abs/1603.02754.

［25］HASSAN A R，SUBASI A.Automatic identification of epileptic seizures from EEG signals using linear programming boosting［J］.Computer Methods and Programs in Biomedicine，2016，136：65-77.

［26］YUAN Q， ZHOU W D， LI S F，et al. Epileptic EEG classification based on extreme learning machine and nonlinear features［J］.Epilepsy Research， 2011，96（1）：29-38.

［27］ZHU G H， LI Y， WEN P. Epileptic seizure detection in EEGs signals using a fast weighted horizontal visibility algo-rithm［J］.Computer Methods and Programs in Biomedicine， 2014，115（2）：64-75.

［28］SIUL Y，LI Y，WEN P P.Clustering technique-based least square support vector machine for EEG signal classification［J］.Computer Methods and? Programs in Biomedicine，2011，104（3）：358-372.

［29］KUMAR Y， DEWAL M L， ANAND R S. Epileptic seizure detection using DWT based fuzzy approximate entropy and support vector machine［J］.Neurocomputing， 2014，133（8）：271-279.

［30］SONG J L，ZHANG R.Application of extreme learning machine to epileptic seizure detection based on lagged Poincare plots［J］. Multidimensional Systems and Signal Processing， 2017， 28（3）：945-959.

［31］ORSHEA A， LIGHTBODY G， BOYLAN G， et al. Investigating the impact of CNN depth on neonatal seizure detection performance［C］∥2018 40th Annual International Conference of the IEEE Engineering in Medicine and Biology Society （EMBC）. Honolulu：IEEE， 2018： 5862-5865.

［32］ZHOU M N， TIAN C， CAO R， et al. Epileptic seizure detection based on EEG signals and CNN［J］.Frontiers in Neuroinformatics， 2018， 12：95.

［33］CAO Y Z， GUO Y X， YU H， et al. Epileptic seizure auto-detection using deep learning method［C］∥2017 4th International Conference on Systems and Informat-ics （ICSAI）.Hangzhou：IEEE， 2017：1076-1081.

［34］CHO K O，JANG H J.Comparison of different input modalities and network structures for deep learning-based seizure detection［J］.Scientific Reports， 2020， 10（1）：122.

［35］STROBLC，BOULESTEIXAL，ZEILEISA，et al.Bias in random forest variable importance measures： Illustrations， sources and a solution［J］.BMC Bioinformatics， 2007，8（25）：25.

（編 輯 李 靜）