1.一种基于粒子群算法的神经元群模型参数自适应优化方法，其特征在于，其步骤包括：

步骤S1：初始化粒子群，设置粒子群算法基本参数，设置每个粒子的参数组合搜索范围，设置神经元群模型基本参数；

步骤S2：计算所述初始化粒子群中每个粒子的适应度值，根据所述初始化粒子群中粒子适应度值，初始化粒子个体极值和全局极值，具体步骤为：步骤S2-1：建立神经元群模型微分方程：

其中，He为兴奋性平均突触增益、Hi为抑制性平均突触增益、τe为膜平均时间常数、τi为树突平均时间常数；C1＝α1·C，C2＝α2·C，C3＝C4＝α3·C＝α4·C；p表示外部输入高斯白噪声，p＝mean+randn*sigma，mean＝220，sigma＝20，randn为随机数；y0为锥体神经元平均膜电位即锥体神经元子群输出；y1为兴奋性中间神经元平均膜电位即兴奋性神经元子群输出，y2为抑制性中间神经元平均膜电位即抑制性神经元子群输出；y3，y4，y5分别为y0，y1，y2的导数；Sigm()为非线性函数，具体为非线性函数模块 v为神经元群的平均膜电位，vmax为神经元群的最大激发值，v0为达到50％最大激发值时的突触后电位，r为该函数在 的斜率；

根据每个粒子的四个参数：兴奋性平均突触增益He、抑制性平均突触增益Hi、膜平均时间常数τe、树突平均时间常数τi，采用4阶Runge-Kutta方法对神经元群模型微分方程进行求解，针对每个粒子都获得y1和y2，根据y1和y2生成相对于每个粒子的仿真脑电信号波形为：eeg＝y1-y2，其中，eeg为锥体神经元突触后膜平均膜电位，代表模型的输出，用于仿真脑电信号，y1为兴奋性中间神经元平均膜电位即兴奋性神经元子群输出，y2为抑制性中间神经元平均膜电位即抑制性神经元子群输出；

步骤S2-2：对相对于所述粒子的仿真脑电信号波形进行1Hz高通滤波以去除0频直流信号，然后对滤波后的信号进行傅里叶变换，获取频谱图中的幅峰值及其对应的频率点fc；

步骤S2-3：每个粒子适应度值计算式为：fitness＝|fc-fg|，其中，fc为当前由频谱图获取的幅峰值对应频率，fg为设定的目标频率；

步骤S2-4：重复步骤S2-1至S2-3计算获得粒子群中每个粒子的粒子适应度值，初始化粒子个体极值和全局极值；

步骤S3：更新粒子群中粒子的速度和位置，形成新的粒子群；

粒子群中每个粒子更新粒子速度和位置所使用的计算公式为：

式中， 是粒子i在第k+1次迭代中第d维的速度； 是粒子i在第k次迭代中第d维的速度，c1，c2分别是第一学习因子和第二学习因子，第一学习因子调节全局最优飞行的最大步长，第二学习因子调节个体最优飞行的最大步长；rand1和rand2是[0，1]之间的随机数；

是粒子i在第k次迭代中第d维的位置， 是粒子i在第k+1次迭代中第d维的位置； 是粒子i在第k次迭代中第d维的个体极值点的位置； 是整个种群在第k次迭代中第d维的全局极值点的位置；wk是惯性权重，其计算式如下：式中，ws为初始惯性权重，we为迭代到最大次数时的惯性权重，k为当前迭代次数，Kmax为最大迭代次数；

当粒子群中所有粒子的位置更新完成后，形成了新的粒子群；

步骤S4：采用所述步骤S2中计算粒子适应度值的方法，对新的粒子群计算每个粒子适应度值；再根据每个粒子适应度值更新粒子个体极值和全局极值，具体为如果粒子的适应度值优于该粒子当前的个体极值，则 被设置为该粒子的位置，并更新个体极值；如所有粒子的个体极值中最优的优于当前的全局极值，则 被设置为该粒子的位置，且更新全局极值；

步骤S5：判断是否满足最大迭代次数，是则输出全局最优粒子，否则返回步骤S3；

步骤S6：根据步骤S5输出的全局最优粒子，获得脑电相应频段的最优参数组合，再根据脑电相应频段的最优参数组合生成这个频段的仿真脑电频段波形。

2.根据权利要求1所述的一种基于粒子群算法的神经元群模型参数自适应优化方法，其特征在于，步骤还包括：脑电包括5个频段δ、θ、α、lowβ、highβ和γ，对脑电5个频段分别设置目标频率fg的值，根据相应的目标频率fg执行步骤S1～S6，最后分别获得脑电5个频段的仿真脑电频段波形。

3.根据权利要求2所述的一种基于粒子群算法的神经元群模型参数自适应优化方法，其特征在于，步骤包括：脑电5个频段的目标频率fg分别设为：δ＝3Hz，θ＝5Hz，α＝12Hz，lowβ＝17Hz，highβ＝

25Hz，γ＝41Hz。

4.根据权利要求1所述的一种基于粒子群算法的神经元群模型参数自适应优化方法，其特征在于，所述步骤S1具体包括：步骤S1-1：设置粒子群算法基本参数，所述粒子群算法基本参数包括第一学习因子c1和第二学习因子c2、粒子维数d、最大迭代次数Kmax、种群规模sizepop、种群粒子搜索范围bound、初始权重we与达到最大迭代时的权重ws；

步骤S1-2：每个粒子所对应的参数组合包括：兴奋性平均突触增益He、抑制性平均突触增益Hi、膜平均时间常数τe和树突平均时间常数τi；设置每个粒子的参数组合搜索范围；

步骤S1-3：设置神经元群模型的基本参数，包括兴奋性反馈环上平均突触连接第一平均概率α1和第二平均概率α2；抑制性反馈环上平均突触连接第三平均概率α3和第四平均概率α4；设置非线性函数Sigm的参数vmax、v0和r，非线性函数 中vmax为神经元群的最大激发值，v0为达到50％最大激发值时的突触后电位，r为该函数在 的斜率；设置外部输入高斯白噪声p＝mean+randn*sigma，mean＝220，sigma＝20，randn为随机数；子群间的最大突触平均连接数目C。

5.根据权利要求4所述的一种基于粒子群算法的神经元群模型参数自适应优化方法，其特征在于，所述步骤S1-2中还包括：每个粒子所对应的参数组合搜索范围为：兴奋性平均突触增益He为2.6～9.75mV、抑制性平均突触增益Hi为17.6～110mV、膜平均时间常数τe为2～150ms、树突平均时间常数τi为2～150ms。