1.基于模糊聚类优化高斯混合模型的高频振荡节律检测方法，其特征在于：包括以下步骤：S101：获取癫痫脑电信号的四个特征：模糊熵、短时能量、功率比和频谱质心构成一个特征向量，在不同时刻得到多组所述四个特征，构成多个特征向量；

S102：采用期望最大化高斯混合模型聚类算法对所述多个特征向量进行分类，得到聚类结果；

分类过程中，采用模糊c‐均值聚类算法对所述期望最大化高斯混合模型聚类算法的参数进行初始化；初始化的具体步骤为：S201：随机产生初始聚类中心vc，c＝1,2,…,C，C是聚类数目；

S202：根据聚类中心vc，计算第i个特征向量属于第c类的隶属度值μic：其中，μic满足条件： C是聚类数目，p是加权模糊参数，Vi是所述D维特征向量，且由提取的所述模糊熵、所述短时能量、所述功率比和所述频谱质心构成；

S203：根据公式(1)更新聚类中心vc：

其中，Vi是所述D维特征向量，由所述模糊熵、所述短时能量、所述功率比和所述频谱质心构成，Nfl是所述D维特征向量Vi的长度，μic是第i个特征向量属于第c类的隶属度值；

S204：当满足max1≤c≤C‖vc,new-vc,old‖＜ε时，所述模糊c‐均值聚类算法结束；其中，vc,new是更新后的聚类中心，vc,old是更新前的聚类中心，ε是设定的精度误差；

根据所述模糊c-均值聚类算法，通过公式(2)、公式(3)和公式(4)分别得到所述基于期望最大化高斯混合模型的聚类算法的初始化参数：权重Wc、聚类中心uc和协方差矩阵Σc：uc＝vc     (3)

其中，Nc是第c类特征向量的数目之和， 是属于第c类的特征向量构成的矩阵；

S103：根据所述期望最大化高斯混合模型聚类算得到的聚类结果，采用箱形图的方法对不同类别特征进行统计学分析，检测出癫痫脑电信号的高频振荡节律。

2.如权利要求1所述的基于模糊聚类优化高斯混合模型的高频振荡节律检测方法，其特征在于：步骤S101中，所述模糊熵的具体计算步骤如下：S301：获取癫痫脑电信号时间序列X＝[x(1),x(2),...,x(N)]；

S302：将所述癫痫脑电信号时间序列按顺序构造m维向量其中，i＝1,2,...,N-m+1，m为正整数，且m≥1，I是元素全为1的m维行向量，S303：计算出向量 与 之间的距离 其

中，i,j∈1,2,...,N-m+1，j≠i，k＝0,1,…,m-1；

S304：根据所述距离 采用模糊隶属度函数 计算所述 与所述 之间的相似度 其中， 为 与 之间的相似度，r为相似容忍度；

S305：根据相似度 得到函数 其中， 为

与 之间的相似度，N为所述癫痫脑电信号时间序列的长度，m为所述向量I的维数；

S306：将所述向量I的维数增加到m+1，得到函数Φm+1(n,r)；

S307：根据所述函数Φm(n,r)和所述函数Φm+1(n,r)，得到所述模糊熵：FuzzyEn(n,r,N)m m+1＝lnΦ(n,r)-lnΦ (n,r)；

*

采用均方法计算得到所述短时能量的估计值： 其中，E (t)为所述短时能量的估计值，x(k)为脑电信号的幅值；

所述功率比R的计算公式为： 其中，P[80-200]是疑似高频振荡节律在ripples带宽内的功率，P[250-500]是疑似高频振荡节律在fast ripples带宽内的功率；

所述频谱质心fc的计算公式为： 其中，T是采样周期，L是窗长，M(k)是多窗口功率谱密度估计值，且 w(l)是长度为L的汉明窗。

3.如权利要求1所述的基于模糊聚类优化高斯混合模型的高频振荡节律检测方法，其特征在于：在步骤102中，基于期望最大化高斯混合模型的聚类算法进行分类的具体步骤为：S401：采用模糊c-均值聚类算法得到所述基于期望最大化高斯混合模型的聚类算法的初始化参数；

S402：根据所述初始化参数，通过公式(5)计算D维特征向量Vi属于第c类的概率其中，t是迭代次数，T≥1， 和 分别是第t次迭代的权重、聚类中心和协方差矩t阵，P(Vi) 是第t次迭代所述D维特征向量Vi对应的概率；

计算所述特征向量Vi对应的概率的公式是：

其中，C是聚类数目，Wc是第c个高斯成分的权重，并且满足 uc为聚D D×D

类中心，uc∈R，Σc为协方差矩阵，Σc∈R ；

S403：根据计算出的所述D维特征向量Vi属于第c类的概率，通过公式(7)、公式(8)和公式(9)，更新高斯混合模型的参数，分别得到第t+1次迭代的权重 聚类中心 和协方差矩阵其中， 是所述D维特征向量Vi属于第c类的概率， 是所述D维特征向量Vi属于第c类的总数，且S404：根据第t+1次迭代的权重 聚类中心 和协方差矩阵 通过公式(5)得到第t+1次迭代所述D维特征向量Vi对应的概率P(vi)t+1，再通过公式(10)计算得到第t+1次迭代的对数似然函数：其中，P(vi)t+1为第t+1次迭代时所述D维特征向量Vi对应的概率；

t+1 t

当|L -L|＜E且t

4.如权利要求1所述的基于模糊聚类优化高斯混合模型的高频振荡节律检测方法，其特征在于：在步骤S103中，使用箱形图得到中位数和四分位距两个指标，对不同类别的统计学特征进行分析，确定出每个类别随时间变化的波形状态，从而检测出高频振荡节律。