基于主成分分析的音乐事件相关脑电变化研究

目的:对在不同音乐背景下获取的脑电数据进行主成分分析.方法:在受测者闭眼欣赏喜欢音乐,闭眼欣赏不喜欢音乐,安静闭眼期间获得脑电数据,并对其进行主成分分析,提取第一主成分,主成分荷载,得出在不同条件和不同事件下主成分的变化规律以及影响主成分变化脑电极的变化规律.数据处理和统计分析采用独立设计的脑电图分析工具箱和主成分分析程序.结果:脑电图数据经过主成分分析后,第一主成分贡献率变化就有比较明显的规律,在闭眼欣赏喜欢音乐时,第一主成分贡献率明显比欣赏不喜欢音乐和安静闭眼时高.在前头部,欣赏喜欢音乐时和欣赏不喜欢音乐时,影响第一主成分的脑电极分别在右前头部和左前头部脑电信号比较强.结论:在不同音乐背景下,大脑左右前头部会受不同音乐所影响,同时也会对人的感觉,情感,情绪等会产生重要的影响.     

脑电信号的高阶奇异谱分析

奇异谱分析是脑电信号分析的一种新方法,脑电信号的奇异谱可以反映脑电的特征,它有助于研究大脑的动力学行为。奇异谱分析方法是基于二阶统计的方法,反映的是信号时间上和空间上的一种线性相关关系。而脑电信号属于非线性信号,其内在的非线性关系很难通过奇异谱得到真实的反映,从而会丢失某些有用的信息。提出一种新的基于高阶统计的脑电奇异谱分析方法,并将其运用于正常脑电和癫痫患者的脑电分析中。大量的实测信号样本仿真实验结果表明,正常脑电和癫痫脑电的奇异谱有明显的不同。此外,基于高阶统计的奇异谱和基于二阶统计的奇异谱相比更能反映出信号的细节。     

基于多特征融合的多分类运动想象脑电信号识别研究

针对目前多分类运动想象脑电识别存在特征提取单一、分类准确率低等问题,提出一种多特征融合的四分类运动想象脑电识别方法来提高识别率。对预处理后的脑电信号分别使用希尔伯特-黄变换、一对多共空间模式、近似熵、模糊熵、样本熵提取结合时频—空域—非线性动力学的初始特征向量,用主成分分析降维,最后使用粒子群优化支持向量机分类。该算法通过对国际标准数据集BCI2005 Data set IIIa中的k3b受试者数据经MATLAB仿真处理后获得93.30%的识别率,均高于单一特征和其它组合特征下的识别率。分别对四名实验者实验采集运动想象脑电数据,使用本研究提出的方法处理获得了72.96%的平均识别率。结果表明多特征融合的特征提取方法能更好的表征运动想象脑电信号,使用粒子群支持向量机可取得较高的识别准确率,为人脑的认知活动提供了一种新的识别方法。     

多通道脑电信号的盲分离

提出一种新的多通道脑电信号盲分离的方法,将小波变换和独立分量分析(independent component analysis,ICA)相结合,利用小波变换的滤噪作用,将混合在原始脑电的部分高频噪声滤除后,再重构原始脑电作为ICA的输入信号,有效地克服了现有ICA算法不能区分噪声的缺陷。实验结果表明,该方法对多通道脑电的盲分离是很有效的。     

基于小波变换的心电信号去噪算法

目的:去除在心电信号采集过程中混入的肌电干扰、工频干扰、基线漂移等噪声信号,避免噪声对心电信号特征点的识别和提取造成误判和漏判。方法:首先利用coif4小波对心电信号按Mallat算法进行分解,然后采用软、硬阈值折衷与小波重构的算法进行去噪。结果:采用MIT/BIH Arrhythmia Database中的心电信号进行仿真、验证,有效去除了三种常见的噪声信号。结论:本方法实时性好,为临床分析与诊断奠定了基础。     

基于独立元分析和非线性指数分析的脑电信号中伪迹分量的自动去除

脑电(electroencephalography,EEG)信号中不可避免地存在着眼动、心跳、肌电信号以及线性噪声等伪迹干扰,这些伪迹的存在极大地影响了脑电信号分析的准确性,因此在进行脑电信号分析前需要去除伪迹干扰。为了有效地去除伪迹,结合独立元分析和非线性指数分析,提出一种自动识别并去除脑电信号中伪迹分量的方法。该方法还可同时用于提取脑电信号中的基本节律如!波等。相应的模拟与实际脑电数据的实验结果表明所提议的方法具有很好的识别和去除脑电信号伪迹分量的性能。     

癫痫病人脑电信号的奇异谱

癫痫是一种常见的神经系统疾患,其唯一客观证据为脑电图的癫痫样发放。在癫痫发作间期,仅有偶发的很难辨别的癫痫样放电,为了正确诊断癫痫病,往往需要医生长时间监测病人的脑电信号,在对脑电信号进行相空间重构,进而对其进行奇异系统分析,发现癫痫病人无论在癫痫发作前、发作中、发作后,其脑电信号的奇异谱曲线不存在噪声平台,明显区别于正常人。是否可以认为脑电信号的奇异谱正代表着大脑的一种基本状态,癫痫患者在未发作时,大脑的基本状态已经处于异常。无论如休,奇异系统分析方法使得可以利用很短的一段脑电数据诊断癫痫。无疑为癫痫病人的临床诊断提供了一条简单、有效的途径。     

基于EMD的左右手运动想象脑电信号分析研究

基于经验模态分解（EMD）理论,提出一种左右手运动想象脑电信号分析方法。首先利用时间窗对脑电信号数据进行划分,对每段数据通过经验模态分解法将其分解为一组固有模态函数IMF,提取主要信号所在的IMF层去除信号中的噪声。对含有主要信号的几层IMF进行Hilbert变换,得到瞬时频率与对应的瞬时幅值。再提取左右手想象的特定频段mu节律和beta节律的能量信号作为特征,分别利用支持向量机（SVM）和Fisher进行了分类比较。对EMD和小波包在去噪和特征提取进行了比较。结果表明,EMD是一种很有效的去噪方法,经过EMD分解后提取的能量信号在区分左右手想象上更具有优势,识别率高。     

自由活动大鼠睡眠状态和呼吸节律的同步监测方法

目的:设计一种睡眠呼吸暂停综合征动物模型的监测方法.方法:通过多导生理仪同步监测自由活动大鼠6 h生理睡眠时间的皮层脑电、颈部肌电,以及呼吸波形.并用人机交互的方法进行数据处理.结果:监测信号稳定、清晰易辨,根据记录到的信号可进行不同清醒-睡眠状态分期,并判断大鼠的呼吸暂停情况.结论:该法操作简便,是小动物睡眠和呼吸监测的可靠方法.     

脑磁图神经活动源数目的估计

在脑磁图信号的分析中,正确估计出脑磁图神经活动源的数目是进一步分析脑磁图信号的前提。目前广泛采用的信息论方法和主成分分析方法都是根据特征值来确定源的数目,这两种方法在源数目较多、噪声较强的情况下,会导致误判。该文提出了一种噪声调节自动阈值的脑磁图源数目判断方法,利用基于噪声调节的主成分分析并结合聂曼- 皮尔逊准则对脑磁图源数目进行估计。同时,该方法采用了基于小波的噪声方差估计,实现了脑磁图信号中噪声方差的精确估计。通过对基于信息论方法、主成分分析方法以及该文所提议方法的实验结果的比较,表明该文所提议方法能更准确地估计脑磁图源数目,特别是在源数目较多、信噪比较小的情况下,仍能准确地估计脑磁图源数目,具有较大的实际意义。     

小波和主分量分析方法研究思维脑电

研究自发脑电和思维活动的关系．利用小波和主分量分析结合的WPCA算法对不同思维任务记录的六导脑电进行处理,并对思维特征的频谱能量和变化率等多指标进行综合分析和计算。结果表明WPCA算法不仅可以实现噪声的去除,而且能提高主分量的贡献率,降低输入矢量的维数。对脑电主分量的分析揭示了脑电与思维个体、思维种类、复杂度以及注意力的联系,思维任务的神经网络分类结果验证了WPCA方法研究脑电和思维的有效性,为进一步理解认知和思维过程,实现对思维的定位和分类提供了依据。     

Evidence for age-associated disinhibition of the wake drive provided by scoring principal components of the resting EEG spectrum in sleep-provoking conditions

Age-associated changes in different bandwidths of the human electroencephalographic (EEG) spectrum are well documented, but their functional significance is poorly understood. This spectrum seems to represent summation of simultaneous influences of several sleep–wake regulatory processes. Scoring of its orthogonal (uncorrelated) principal components can help in separation of the brain signatures of these processes. In particular, the opposite age-associated changes were documented for scores on the two largest (1st and 2nd) principal components of the sleep EEG spectrum. A decrease of the first score and an increase of the second score can reflect, respectively, the weakening of the sleep drive and disinhibition of the opposing wake drive with age. In order to support the suggestion of age-associated disinhibition of the wake drive from the antagonistic influence of the sleep drive, we analyzed principal component scores of the resting EEG spectra obtained in sleep deprivation experiments with 81 healthy young adults aged between 19 and 26 and 40 healthy older adults aged between 45 and 66 years. At the second day of the sleep deprivation experiments, frontal scores on the 1st principal component of the EEG spectrum demonstrated an age-associated reduction of response to eyes closed relaxation. Scores on the 2nd principal component were either initially increased during wakefulness or less responsive to such sleep-provoking conditions (frontal and occipital scores, respectively). These results are in line with the suggestion of disinhibition of the wake drive with age. They provide an explanation of why older adults are less vulnerable to sleep deprivation than young adults.     

Quantification of sleepiness through principal component analysis of the electroencephalographic spectrum

Although circadian and sleep research has made extraordinary progress in the recent years, one remaining challenge is the objective quantification of sleepiness in individuals suffering from sleep deprivation, sleep restriction, and excessive somnolence. The major goal of the present study was to apply principal component analysis to the wake electroencephalographic (EEG) spectrum in order to establish an objective measure of sleepiness. The present analysis was led by the hypothesis that in sleep-deprived individuals, the time course of self-rated sleepiness correlates with the time course score on the 2nd principal component of the EEG spectrum. The resting EEG of 15 young subjects was recorded at 2-h intervals for 32-50 h. Principal component analysis was performed on the sets of 16 single-Hz log-transformed EEG powers (1-16 Hz frequency range). The time course of self-perceived sleepiness correlated strongly with the time course of the 2nd principal component score, irrespective of derivation (frontal or occipital) and of analyzed section of the 7-min EEG record (2-min section with eyes open or any of the five 1-min sections with eyes closed). This result indicates the possibility of deriving an objective index of physiological sleepiness by applying principal component analysis to the wake EEG spectrum.     

Quantification of Sleepiness Through Principal Component Analysis of the Electroencephalographic Spectrum

Although circadian and sleep research has made extraordinary progress in the recent years, one remaining challenge is the objective quantification of sleepiness in individuals suffering from sleep deprivation, sleep restriction, and excessive somnolence. The major goal of the present study was to apply principal component analysis to the wake electroencephalographic (EEG) spectrum in order to establish an objective measure of sleepiness. The present analysis was led by the hypothesis that in sleep-deprived individuals, the time course of self-rated sleepiness correlates with the time course score on the 2nd principal component of the EEG spectrum. The resting EEG of 15 young subjects was recorded at 2-h intervals for 32–50?h. Principal component analysis was performed on the sets of 16 single-Hz log-transformed EEG powers (1–16?Hz frequency range). The time course of self-perceived sleepiness correlated strongly with the time course of the 2nd principal component score, irrespective of derivation (frontal or occipital) and of analyzed section of the 7-min EEG record (2-min section with eyes open or any of the five 1-min sections with eyes closed). This result indicates the possibility of deriving an objective index of physiological sleepiness by applying principal component analysis to the wake EEG spectrum. (Author correspondence: putilov@ngs.ru)     

Principal components of electroencephalographic spectrum as markers of opponent processes underlying ultradian sleep cycles

Sleep-wake regulation involves reciprocal interactions between sleep- and wake-promoting processes that inhibit one another. To uncover the signatures of the opponent processes underlying ultradian sleep cycles, principal component analysis was performed on the sets of 16 single-Hz log-transformed electroencephalographic (EEG) power densities (1-16?Hz frequency range). Data were collected during unrestricted night sleep followed by 9 20-min naps (14 women aged 17-55 yrs) and during 12 20-min naps after either restriction or deprivation of sleep (9 males and 9 males, respectively, aged 18-22 yrs). It was found that any subset of power spectra could be reduced to the invariant four-principal component structure. The time courses of scores on these four components might be interpreted as the spectral EEG markers of the kinetics of two pairs of opponent chronoregulatory processes. In a sequence of ultradian sleep cycles, the 1st and 2nd components represent the alternations between competing drives for sleep and wakefulness, respectively, whereas the 3rd and 4th components reflect the alternations between light and deep sleep, respectively. The results suggest that principal component structuring of EEG spectrum can be employed for derivation of the parameters of the quantitative models conceptualizing the three major aspects of sleep-wake regulation—homeostatic, circadian, and ultradian processes.     

High Resolution Topography of Age-Related Changes in Non-Rapid Eye Movement Sleep Electroencephalography

Sleeping brain activity reflects brain anatomy and physiology. The aim of this study was to use high density (256 channel) electroencephalography (EEG) during sleep to characterize topographic changes in sleep EEG power across normal aging, with high spatial resolution. Sleep was evaluated in 92 healthy adults aged 18–65 years old using full polysomnography and high density EEG. After artifact removal, spectral power density was calculated for standard frequency bands for all channels, averaged across the NREM periods of the first 3 sleep cycles. To quantify topographic changes with age, maps were generated of the Pearson’s coefficient of the correlation between power and age at each electrode. Significant correlations were determined by statistical non-parametric mapping. Absolute slow wave power declined significantly with increasing age across the entire scalp, whereas declines in theta and sigma power were significant only in frontal regions. Power in fast spindle frequencies declined significantly with increasing age frontally, whereas absolute power of slow spindle frequencies showed no significant change with age. When EEG power was normalized across the scalp, a left centro-parietal region showed significantly less age-related decline in power than the rest of the scalp. This partial preservation was particularly significant in the slow wave and sigma bands. The effect of age on sleep EEG varies substantially by region and frequency band. This non-uniformity should inform the design of future investigations of aging and sleep. This study provides normative data on the effect of age on sleep EEG topography, and provides a basis from which to explore the mechanisms of normal aging as well as neurodegenerative disorders for which age is a risk factor.     

Few-sweep estimation of evoked potential based on a generalized subspace approach

A generalized subspace approach is proposed for single channel brain evoked potential (EP) extraction from background electroencephalogram (EEG) signal. The method realizes the optimum estimate of EP signal from the observable noisy signal. The underlying principle is to project the signal and noise into signal and noise coefficient subspace respectively by applying projection matrix at first. Secondly, coefficient weighting matrix is achieved based on the autocorrelation matrices of the noise and the noisy signal. With the coefficient weighting matrix, we can remove the noise projection coefficients and estimate the signal ones. EP signal is then obtained by averaging the signals estimated with the reconstruction matrix. Given different signal-to-noise ratio (SNR) conditions, the algorithm can estimate the EP signal with only two sweeps observable noisy signals. Our approach is shown to have excellent capability of estimating EP signal even in poor SNR conditions. The interference of spontaneous EEG has been eliminated with significantly improved SNR. The simulation results have demonstrated the effectiveness and superior performance of the proposed method.     

Principal Components of Electroencephalographic Spectrum as Markers of Opponent Processes Underlying Ultradian Sleep Cycles

Sleep-wake regulation involves reciprocal interactions between sleep- and wake-promoting processes that inhibit one another. To uncover the signatures of the opponent processes underlying ultradian sleep cycles, principal component analysis was performed on the sets of 16 single-Hz log-transformed electroencephalographic (EEG) power densities (1–16?Hz frequency range). Data were collected during unrestricted night sleep followed by 9 20-min naps (14 women aged 17–55 yrs) and during 12 20-min naps after either restriction or deprivation of sleep (9 males and 9 males, respectively, aged 18–22 yrs). It was found that any subset of power spectra could be reduced to the invariant four–principal component structure. The time courses of scores on these four components might be interpreted as the spectral EEG markers of the kinetics of two pairs of opponent chronoregulatory processes. In a sequence of ultradian sleep cycles, the 1st and 2nd components represent the alternations between competing drives for sleep and wakefulness, respectively, whereas the 3rd and 4th components reflect the alternations between light and deep sleep, respectively. The results suggest that principal component structuring of EEG spectrum can be employed for derivation of the parameters of the quantitative models conceptualizing the three major aspects of sleep-wake regulation—homeostatic, circadian, and ultradian processes. (Author correspondence: putilov@ngs.ru)     

An automatic EEG-based sleep staging system with introducing NAoSP and NAoGP as new metrics for sleep staging systems

Different biological signals are recorded in sleep labs during sleep for the diagnosis and treatment of human sleep problems. Classification of sleep stages with electroencephalography (EEG) is preferred to other biological signals due to its advantages such as providing clinical information, cost-effectiveness, comfort, and ease of use. The evaluation of EEG signals taken during sleep by clinicians is a tiring, time-consuming, and error-prone method. Therefore, it is clinically mandatory to determine sleep stages by using software-supported systems. Like all classification problems, the accuracy rate is used to compare the performance of studies in this domain, but this metric can be accurate when the number of observations is equal in classes. However, since there is not an equal number of observations in sleep stages, this metric is insufficient in the evaluation of such systems. For this purpose, in recent years, Cohen’s kappa coefficient and even the sensitivity of NREM1 have been used for comparing the performance of these systems. Still, none of them examine the system from all dimensions. Therefore, in this study, two new metrics based on the polygon area metric, called the normalized area of sensitivity polygon and normalized area of the general polygon, are proposed for the performance evaluation of sleep staging systems. In addition, a new sleep staging system is introduced using the applications offered by the MATLAB program. The existing systems discussed in the literature were examined with the proposed metrics, and the best systems were compared with the proposed sleep staging system. According to the results, the proposed system excels in comparison with the most advanced machine learning methods. The single-channel method introduced based on the proposed metrics can be used for robust and reliable sleep stage classification from all dimensions required for real-time applications.Electronic supplementary materialThe online version of this article (10.1007/s11571-020-09641-2) contains supplementary material, which is available to authorized users.