第一获得原始信号预处理中删除不需要的工件提出的信号。接下来,使用各种分解预处理信号分解和分析方法。在下一步中价值和信息特征提取信号的分解组件。这些提取的功能最终美联储的机器学习步骤分类模型开发分类特征实例两种情况下“关注”或“nonattention。”gydF4y2Ba
3.1。数据预处理gydF4y2Ba
获得原始ECG信号包含了一些固有的不必要的工件前需要处理任何分析可以执行。这些工件的原因,通常是频率噪音或基线的趋势,可能是由于多种原因,如受试者的运动导致运动工件,呼吸模式工件,松散皮肤接触的电极,和电气干扰(通常发现约55 Hz)。因此一个预处理步骤设计,确保信号是干净和工件自由之前分析。gydF4y2Ba
3.1.1。心电图预处理gydF4y2Ba
心电信号的预处理步骤如图所示gydF4y2Ba
3gydF4y2Ba。因为每个信号必须过滤基于固有噪声的类型不同,原始ECG信号首先过滤使用“SGolay”过滤方法。“SGolay”过滤器由Savitzky-Golay开发。这个过滤器是一个数字多项式滤波器基于最小二乘平滑机制。SGolay过滤器通常用于消除噪声信号的频率跨度大。他们执行比标准平均FIR滤波器,因为这些过滤器倾向于保留很大一部分信号高频内容而只去除噪声(gydF4y2Ba
14gydF4y2Ba]。gydF4y2Ba
心电图预处理。gydF4y2Ba
接下来,过滤后的心电图数据发送通过去除基线漂移的一步。通常基线漂移是观察心电图记录由于呼吸、肌肉收缩、电极阻抗变化由于主体的运动(gydF4y2Ba
15gydF4y2Ba]。首先去除基线漂移的回归线最适合的样品在一个窗口的大小等于确定了采样率。gydF4y2Ba
鉴于gydF4y2Ba
ngydF4y2Ba
点的心电图信号gydF4y2Ba
(gydF4y2Ba
xgydF4y2Ba
1gydF4y2Ba
,gydF4y2Ba
ygydF4y2Ba
1gydF4y2Ba
)gydF4y2Ba
,gydF4y2Ba
(gydF4y2Ba
xgydF4y2Ba
2gydF4y2Ba
,gydF4y2Ba
ygydF4y2Ba
2gydF4y2Ba
)gydF4y2Ba
,gydF4y2Ba
…gydF4y2Ba
,gydF4y2Ba
(gydF4y2Ba
xgydF4y2Ba
ngydF4y2Ba
,gydF4y2Ba
ygydF4y2Ba
ngydF4y2Ba
)gydF4y2Ba
,最适合线与这些点可以计算如下:gydF4y2Ba
(1)gydF4y2Ba
米gydF4y2Ba
=gydF4y2Ba
ngydF4y2Ba
(gydF4y2Ba
∑gydF4y2Ba
1gydF4y2Ba
ngydF4y2Ba
xgydF4y2Ba
ygydF4y2Ba
)gydF4y2Ba
- - - - - -gydF4y2Ba
(gydF4y2Ba
∑gydF4y2Ba
1gydF4y2Ba
ngydF4y2Ba
xgydF4y2Ba
)gydF4y2Ba
(gydF4y2Ba
∑gydF4y2Ba
1gydF4y2Ba
ngydF4y2Ba
ygydF4y2Ba
)gydF4y2Ba
ngydF4y2Ba
(gydF4y2Ba
∑gydF4y2Ba
1gydF4y2Ba
ngydF4y2Ba
xgydF4y2Ba
2gydF4y2Ba
)gydF4y2Ba
- - - - - -gydF4y2Ba
(gydF4y2Ba
∑gydF4y2Ba
1gydF4y2Ba
ngydF4y2Ba
xgydF4y2Ba
)gydF4y2Ba
2gydF4y2Ba
,gydF4y2Ba
bgydF4y2Ba
=gydF4y2Ba
∑gydF4y2Ba
1gydF4y2Ba
ngydF4y2Ba
ygydF4y2Ba
- - - - - -gydF4y2Ba
米gydF4y2Ba
(gydF4y2Ba
∑gydF4y2Ba
1gydF4y2Ba
ngydF4y2Ba
xgydF4y2Ba
)gydF4y2Ba
ngydF4y2Ba
,gydF4y2Ba
ygydF4y2Ba
=gydF4y2Ba
米gydF4y2Ba
xgydF4y2Ba
+gydF4y2Ba
bgydF4y2Ba
,gydF4y2Ba
在哪里gydF4y2Ba
ygydF4y2Ba是一个点,gydF4y2Ba
米gydF4y2Ba直线的斜率,gydF4y2Ba
bgydF4y2Ba是拦截。每个窗口的计算最适合线然后减去从原始信号窗口获得基线drift-free信号。gydF4y2Ba
原始ECG信号过滤后的噪声和基线漂移,然后信号被分成两个部分的基础上,采集和实验框架。信号的两个部分,即“有趣”和“noninteresting”从原始信号中提取使用时间戳,在信号采集记录和索引。分离和分析数据的两个部分分别促进监督学习机制在训练阶段的机器学习步骤。gydF4y2Ba
3.1.2。脑电图预处理gydF4y2Ba
EEG信号是由一个复杂和非线性组合的几个不同的波形也称为带组件。乐队的每个组件分类中存在的频率范围。个人意识的状态可能会使一个频率范围比其他人更明显(gydF4y2Ba
16gydF4y2Ba]。如图gydF4y2Ba
4gydF4y2Ba,不同的乐队组件从原始EEG信号中提取使用巴特沃斯带通滤波器。五个主要乐队EEG信号的提取,即δ(0.2 4赫兹),θ(4 - 8 Hz),α8—13赫兹(),β(13-30 Hz),和伽马(30-55赫兹)。gydF4y2Ba
脑电图预处理步骤。gydF4y2Ba
3.2。心电图分解:使用斯托克变换gydF4y2Ba
的gydF4y2Ba
年代gydF4y2Ba变换被斯托克和他的同事在1996年提出的。的区别gydF4y2Ba
年代gydF4y2Ba变换是它产生分解frequency-dependant分辨率的时频域而完全保留当地的相位信息。换句话说,gydF4y2Ba
年代gydF4y2Ba变换不仅当地的功率谱估计,而且当地的相位谱,这是非常可取的研究复杂的生理信号,如心电图。gydF4y2Ba
分析时动态频谱或当地的光谱性质不稳定等观察心电图的一些常见的方法包括短时傅里叶变换(STFT) [gydF4y2Ba
17gydF4y2Ba],伽柏变换[gydF4y2Ba
18gydF4y2Ba),复杂的解调(gydF4y2Ba
19gydF4y2Ba)产生的一系列带通滤过的声音,也与小波滤波器组理论等等。一些方法代表了变换在时间和频率域的组合如科恩类(gydF4y2Ba
20.gydF4y2Ba的广义时频分布(GTFD) Cone-Kernel分布(gydF4y2Ba
21gydF4y2Ba],Choi-Williams分布[gydF4y2Ba
22gydF4y2Ba)以及平滑伪魏格纳分布(PWD) [gydF4y2Ba
23gydF4y2Ba]。一个比较流行的方法分解和小波变换在时频域分析。离散小波变换和DWT执行分解的信号提供了出色的时间分辨率,同时保持关键光谱信息或频率分辨率gydF4y2Ba
24gydF4y2Ba,gydF4y2Ba
25gydF4y2Ba]。gydF4y2Ba
虽然gydF4y2Ba
年代gydF4y2Ba变换是类似于渐进分辨率,小波变换与小波变换不同,gydF4y2Ba
年代gydF4y2Ba变换保留绝对引用的相位信息。绝对引用阶段意味着计算的相位信息gydF4y2Ba
年代gydF4y2Ba变换是参考时间gydF4y2Ba
tgydF4y2Ba
=gydF4y2Ba
0gydF4y2Ba
,这也是对傅里叶变换的相位。唯一的区别是gydF4y2Ba
年代gydF4y2Ba变换提供了每个样本的绝对引用的相位信息的时频空间。gydF4y2Ba
3.2.1之上。数学公式的<斜体> < /斜体>变换gydF4y2Ba
有两个品种gydF4y2Ba
年代gydF4y2Ba变换、连续和离散。连续gydF4y2Ba
年代gydF4y2Ba变换(gydF4y2Ba
26gydF4y2Ba)本质上是STFT的延伸。它也可以被视为一种phase-corrected格式的连续小波变换(CWT)。gydF4y2Ba
STFT的信号gydF4y2Ba
h (t)gydF4y2Ba被定义为gydF4y2Ba
(2)gydF4y2Ba
STFTgydF4y2Ba
(gydF4y2Ba
τgydF4y2Ba
,gydF4y2Ba
fgydF4y2Ba
)gydF4y2Ba
=gydF4y2Ba
∫gydF4y2Ba
- - - - - -gydF4y2Ba
∞gydF4y2Ba
∞gydF4y2Ba
hgydF4y2Ba
(gydF4y2Ba
tgydF4y2Ba
)gydF4y2Ba
ggydF4y2Ba
(gydF4y2Ba
τgydF4y2Ba
- - - - - -gydF4y2Ba
tgydF4y2Ba
)gydF4y2Ba
egydF4y2Ba
- - - - - -gydF4y2Ba
jgydF4y2Ba
2gydF4y2Ba
πgydF4y2Ba
fgydF4y2Ba
tgydF4y2Ba
dgydF4y2Ba
tgydF4y2Ba
,gydF4y2Ba
在哪里gydF4y2Ba
τgydF4y2Ba
是光谱定位的时候,gydF4y2Ba
fgydF4y2Ba
傅里叶频率,gydF4y2Ba
ggydF4y2Ba
(gydF4y2Ba
tgydF4y2Ba
)gydF4y2Ba
表示一个窗口函数。gydF4y2Ba
的gydF4y2Ba
年代gydF4y2Ba变换可以从上面的STFT方程只需用窗函数gydF4y2Ba
ggydF4y2Ba
(gydF4y2Ba
tgydF4y2Ba
)gydF4y2Ba
高斯函数:gydF4y2Ba
(3)gydF4y2Ba
ggydF4y2Ba
(gydF4y2Ba
tgydF4y2Ba
)gydF4y2Ba
=gydF4y2Ba
|gydF4y2Ba
fgydF4y2Ba
|gydF4y2Ba
2gydF4y2Ba
πgydF4y2Ba
egydF4y2Ba
- - - - - -gydF4y2Ba
(gydF4y2Ba
tgydF4y2Ba
2gydF4y2Ba
fgydF4y2Ba
2gydF4y2Ba
)gydF4y2Ba
/gydF4y2Ba
2gydF4y2Ba
。gydF4y2Ba
因此,gydF4y2Ba
年代gydF4y2Ba变换是数学上定义如下:gydF4y2Ba
(4)gydF4y2Ba
年代gydF4y2Ba
(gydF4y2Ba
τgydF4y2Ba
,gydF4y2Ba
fgydF4y2Ba
)gydF4y2Ba
=gydF4y2Ba
∫gydF4y2Ba
- - - - - -gydF4y2Ba
∞gydF4y2Ba
∞gydF4y2Ba
hgydF4y2Ba
(gydF4y2Ba
tgydF4y2Ba
)gydF4y2Ba
|gydF4y2Ba
fgydF4y2Ba
|gydF4y2Ba
2gydF4y2Ba
πgydF4y2Ba
egydF4y2Ba
- - - - - -gydF4y2Ba
(gydF4y2Ba
(gydF4y2Ba
τgydF4y2Ba
- - - - - -gydF4y2Ba
tgydF4y2Ba
)gydF4y2Ba
2gydF4y2Ba
fgydF4y2Ba
2gydF4y2Ba
)gydF4y2Ba
/gydF4y2Ba
2gydF4y2Ba
egydF4y2Ba
- - - - - -gydF4y2Ba
jgydF4y2Ba
2gydF4y2Ba
πgydF4y2Ba
fgydF4y2Ba
tgydF4y2Ba
dgydF4y2Ba
tgydF4y2Ba
。gydF4y2Ba
自gydF4y2Ba
年代gydF4y2Ba变换本质上与高斯窗函数分解,它可以推断,一个更广泛的窗口在时间域的变换可以提供更高的分辨率较低的频率,和一个狭窄的高斯窗更高频率的分辨率是更好的突出。gydF4y2Ba
的应用程序gydF4y2Ba
年代gydF4y2Ba变换在这项研究中,连续gydF4y2Ba
年代gydF4y2Ba变换不被证明是一个实用的选择。仅仅因为收购ECG信号的本身进行离散采样和连续分解这个信号的频率可能会非常耗时,因此不计算务实。因此一个离散的版本的gydF4y2Ba
年代gydF4y2Ba变换采用分解的心电信号。gydF4y2Ba
离散gydF4y2Ba
年代gydF4y2Ba变换可以介绍如下。gydF4y2Ba
让gydF4y2Ba
hgydF4y2Ba
(gydF4y2Ba
kgydF4y2Ba
TgydF4y2Ba
]gydF4y2Ba
是离散时间序列信号的研究,gydF4y2Ba
kgydF4y2Ba
=gydF4y2Ba
0gydF4y2Ba
,gydF4y2Ba
1gydF4y2Ba
,gydF4y2Ba
…gydF4y2Ba
,gydF4y2Ba
NgydF4y2Ba
- - - - - -gydF4y2Ba
1gydF4y2Ba
,gydF4y2Ba
TgydF4y2Ba
是采样间隔的时间。傅里叶变换的离散格式可以如下所示:gydF4y2Ba
(5)gydF4y2Ba
HgydF4y2Ba
(gydF4y2Ba
ngydF4y2Ba
NgydF4y2Ba
TgydF4y2Ba
]gydF4y2Ba
=gydF4y2Ba
1gydF4y2Ba
NgydF4y2Ba
∑gydF4y2Ba
kgydF4y2Ba
=gydF4y2Ba
0gydF4y2Ba
NgydF4y2Ba
- - - - - -gydF4y2Ba
1gydF4y2Ba
hgydF4y2Ba
(gydF4y2Ba
kgydF4y2Ba
TgydF4y2Ba
]gydF4y2Ba
egydF4y2Ba
2gydF4y2Ba
jgydF4y2Ba
πgydF4y2Ba
ngydF4y2Ba
kgydF4y2Ba
/gydF4y2Ba
NgydF4y2Ba
。gydF4y2Ba
使用连续gydF4y2Ba
年代gydF4y2Ba变换方程和上面的方程,时间序列,gydF4y2Ba
hgydF4y2Ba
(gydF4y2Ba
kgydF4y2Ba
TgydF4y2Ba
]gydF4y2Ba
的gydF4y2Ba
年代gydF4y2Ba变换可以表示如下:(gydF4y2Ba
fgydF4y2Ba
→gydF4y2Ba
ngydF4y2Ba
/gydF4y2Ba
NgydF4y2Ba
TgydF4y2Ba
和gydF4y2Ba
τgydF4y2Ba
→gydF4y2Ba
jgydF4y2Ba
TgydF4y2Ba
)gydF4y2Ba
(6)gydF4y2Ba
年代gydF4y2Ba
(gydF4y2Ba
jgydF4y2Ba
TgydF4y2Ba
,gydF4y2Ba
ngydF4y2Ba
NgydF4y2Ba
TgydF4y2Ba
]gydF4y2Ba
=gydF4y2Ba
∑gydF4y2Ba
米gydF4y2Ba
=gydF4y2Ba
0gydF4y2Ba
NgydF4y2Ba
- - - - - -gydF4y2Ba
1gydF4y2Ba
HgydF4y2Ba
(gydF4y2Ba
米gydF4y2Ba
+gydF4y2Ba
ngydF4y2Ba
NgydF4y2Ba
TgydF4y2Ba
]gydF4y2Ba
egydF4y2Ba
2gydF4y2Ba
πgydF4y2Ba
2gydF4y2Ba
米gydF4y2Ba
2gydF4y2Ba
/gydF4y2Ba
ngydF4y2Ba
2gydF4y2Ba
egydF4y2Ba
2gydF4y2Ba
jgydF4y2Ba
πgydF4y2Ba
米gydF4y2Ba
jgydF4y2Ba
/gydF4y2Ba
NgydF4y2Ba
,gydF4y2Ba
ngydF4y2Ba
≠gydF4y2Ba
0gydF4y2Ba
,gydF4y2Ba
在哪里gydF4y2Ba
jgydF4y2Ba
,gydF4y2Ba
米gydF4y2Ba
,gydF4y2Ba
ngydF4y2Ba
=gydF4y2Ba
0gydF4y2Ba
,gydF4y2Ba
1gydF4y2Ba
,gydF4y2Ba
…gydF4y2Ba
,gydF4y2Ba
NgydF4y2Ba
- - - - - -gydF4y2Ba
1gydF4y2Ba
。gydF4y2Ba
3.2.2。应用<斜体> < /斜体>变换gydF4y2Ba
图gydF4y2Ba
5gydF4y2Ba显示了不同的步骤分解的心电信号使用gydF4y2Ba
年代gydF4y2Ba变换。首先,预处理通过窗口发送心电图信号机制。在这种机制中,预处理ECG信号分割成小窗户。这些窗户是不重叠的,包含心电图数据10秒的间隔(128 Hz * 10秒= 1280数据点/窗口)。gydF4y2Ba
年代gydF4y2Ba变换应用心电图信号。gydF4y2Ba
窗口后一步,每一个10秒的windows分解使用gydF4y2Ba
年代gydF4y2Ba变换。的输出gydF4y2Ba
年代gydF4y2Ba变换是一个复杂的二维矩阵的行代表频率和列代表时间值。的gydF4y2Ba
年代gydF4y2Ba变换算法应用于本研究优化生产逐步与步长是1赫兹频率范围和样本结果之间的时间间隔是1单位。gydF4y2Ba
的一个示例输出5秒钟后心电图数据的窗口gydF4y2Ba
年代gydF4y2Ba转换图gydF4y2Ba
6gydF4y2Ba。gydF4y2Ba
(一个)显示频谱的contour-based可视化,基于gydF4y2Ba
年代gydF4y2Ba变换信号的窗口。(b)显示了原始信号窗口。gydF4y2Ba
图gydF4y2Ba
6gydF4y2Ba显示准确的原始(图的点对点表示gydF4y2Ba
6 (b)gydF4y2Ba)信号gydF4y2Ba
年代gydF4y2Ba变换时频域。的gydF4y2Ba
年代gydF4y2Ba变换输出矩阵一直等值线图显示(图所示gydF4y2Ba
6(一)gydF4y2Ba)。gydF4y2Ba
3.2.3。特征提取gydF4y2Ba
每个窗口的输出是一个频率时间矩阵表示。矩阵的每个实例频率点和时间点(由行和列的位置,职责)。所以整个输出矩阵可以提出如下:gydF4y2Ba
圣gydF4y2Ba
(gydF4y2Ba
xgydF4y2Ba
,gydF4y2Ba
ygydF4y2Ba
)gydF4y2Ba
,gydF4y2Ba在哪里gydF4y2Ba
xgydF4y2Ba
频率(行)位置和吗gydF4y2Ba
ygydF4y2Ba
是时候(列)的位置。gydF4y2Ba
的提取特性推导出输出矩阵的圣执行两个步骤。在第一步的输出矩阵从两个维度减少到一个维度。这是通过计算某些统计措施以及频率维度gydF4y2Ba
xgydF4y2Ba
,同时保留时间维度的不连续性gydF4y2Ba
ygydF4y2Ba
为是。计算统计措施以及频率gydF4y2Ba
(gydF4y2Ba
fgydF4y2Ba
)gydF4y2Ba
如下:gydF4y2Ba
意思的频率gydF4y2Ba
(gydF4y2Ba
fgydF4y2Ba
)gydF4y2Ba
,gydF4y2Ba
的频率gydF4y2Ba
(gydF4y2Ba
fgydF4y2Ba
)gydF4y2Ba
,gydF4y2Ba
产品的频率gydF4y2Ba
(gydF4y2Ba
fgydF4y2Ba
)gydF4y2Ba
,gydF4y2Ba
标准偏差的频率gydF4y2Ba
(gydF4y2Ba
fgydF4y2Ba
)gydF4y2Ba
,gydF4y2Ba
范围gydF4y2Ba
(gydF4y2Ba
fgydF4y2Ba
)gydF4y2Ba
。gydF4y2Ba
在第一步得到的数组从频域特性如下:gydF4y2Ba
(7)gydF4y2Ba
FregydF4y2Ba
问gydF4y2Ba
场效应晶体管gydF4y2Ba
=gydF4y2Ba
(gydF4y2Ba
的意思是gydF4y2Ba
(gydF4y2Ba
fgydF4y2Ba
)gydF4y2Ba
,gydF4y2Ba
总和gydF4y2Ba
(gydF4y2Ba
fgydF4y2Ba
)gydF4y2Ba
,gydF4y2Ba
产品gydF4y2Ba
(gydF4y2Ba
fgydF4y2Ba
)gydF4y2Ba
,gydF4y2Ba
性病gydF4y2Ba
(gydF4y2Ba
fgydF4y2Ba
)gydF4y2Ba
,gydF4y2Ba
范围gydF4y2Ba
(gydF4y2Ba
fgydF4y2Ba
)gydF4y2Ba
]gydF4y2Ba
。gydF4y2Ba
下一步是计算沿时域统计特性。gydF4y2Ba
意思是:gydF4y2Ba
(8)gydF4y2Ba
的意思是gydF4y2Ba
(gydF4y2Ba
圣gydF4y2Ba
)gydF4y2Ba
=gydF4y2Ba
的意思是gydF4y2Ba
(gydF4y2Ba
fgydF4y2Ba
我gydF4y2Ba
)gydF4y2Ba
,gydF4y2Ba
在哪里gydF4y2Ba
fgydF4y2Ba
我gydF4y2Ba
∈gydF4y2Ba
FregydF4y2Ba
问gydF4y2Ba
场效应晶体管gydF4y2Ba
。gydF4y2Ba
总结:gydF4y2Ba
(9)gydF4y2Ba
总和gydF4y2Ba
(gydF4y2Ba
圣gydF4y2Ba
)gydF4y2Ba
=gydF4y2Ba
总和gydF4y2Ba
(gydF4y2Ba
fgydF4y2Ba
我gydF4y2Ba
)gydF4y2Ba
,gydF4y2Ba
在哪里gydF4y2Ba
fgydF4y2Ba
我gydF4y2Ba
∈gydF4y2Ba
FregydF4y2Ba
问gydF4y2Ba
场效应晶体管gydF4y2Ba
。gydF4y2Ba
自协方差的意思是:gydF4y2Ba
(10)gydF4y2Ba
的意思是gydF4y2Ba
(gydF4y2Ba
自协方差gydF4y2Ba
(gydF4y2Ba
圣gydF4y2Ba
)gydF4y2Ba
)gydF4y2Ba
=gydF4y2Ba
的意思是gydF4y2Ba
(gydF4y2Ba
自协方差gydF4y2Ba
(gydF4y2Ba
fgydF4y2Ba
我gydF4y2Ba
)gydF4y2Ba
)gydF4y2Ba
,gydF4y2Ba
在哪里gydF4y2Ba
fgydF4y2Ba
我gydF4y2Ba
∈gydF4y2Ba
FregydF4y2Ba
问gydF4y2Ba
场效应晶体管gydF4y2Ba
。gydF4y2Ba
互相关的总和:gydF4y2Ba
(11)gydF4y2Ba
总和gydF4y2Ba
(gydF4y2Ba
autocorrilationgydF4y2Ba
(gydF4y2Ba
圣gydF4y2Ba
)gydF4y2Ba
)gydF4y2Ba
=gydF4y2Ba
总和gydF4y2Ba
(gydF4y2Ba
自相关gydF4y2Ba
(gydF4y2Ba
fgydF4y2Ba
我gydF4y2Ba
)gydF4y2Ba
)gydF4y2Ba
,gydF4y2Ba
在哪里gydF4y2Ba
fgydF4y2Ba
我gydF4y2Ba
∈gydF4y2Ba
FregydF4y2Ba
问gydF4y2Ba
场效应晶体管gydF4y2Ba
。gydF4y2Ba
日志gydF4y2Ba2gydF4y2Ba方差:gydF4y2Ba
(12)gydF4y2Ba
日志gydF4y2Ba
2gydF4y2Ba
(gydF4y2Ba
方差gydF4y2Ba
(gydF4y2Ba
圣gydF4y2Ba
)gydF4y2Ba
)gydF4y2Ba
=gydF4y2Ba
日志gydF4y2Ba
2gydF4y2Ba
(gydF4y2Ba
方差gydF4y2Ba
(gydF4y2Ba
fgydF4y2Ba
我gydF4y2Ba
)gydF4y2Ba
)gydF4y2Ba
,gydF4y2Ba
在哪里gydF4y2Ba
fgydF4y2Ba
我gydF4y2Ba
∈gydF4y2Ba
FregydF4y2Ba
问gydF4y2Ba
场效应晶体管gydF4y2Ba
。gydF4y2Ba
两个额外的功能是计算从最初获得圣矩阵。gydF4y2Ba
的最大频率:gydF4y2Ba
(13)gydF4y2Ba
的意思是gydF4y2Ba
gydF4y2Ba
(gydF4y2Ba
马克斯gydF4y2Ba
gydF4y2Ba
(gydF4y2Ba
圣gydF4y2Ba
)gydF4y2Ba
)gydF4y2Ba
=gydF4y2Ba
的意思是gydF4y2Ba
gydF4y2Ba
(gydF4y2Ba
马克斯gydF4y2Ba
gydF4y2Ba
(gydF4y2Ba
圣gydF4y2Ba
1gydF4y2Ba
,gydF4y2Ba
ygydF4y2Ba
,gydF4y2Ba
圣gydF4y2Ba
2gydF4y2Ba
,gydF4y2Ba
ygydF4y2Ba
,gydF4y2Ba
…gydF4y2Ba
,gydF4y2Ba
圣gydF4y2Ba
xgydF4y2Ba
,gydF4y2Ba
ygydF4y2Ba
)gydF4y2Ba
)gydF4y2Ba
。gydF4y2Ba
平均绝对偏差的频率:gydF4y2Ba
(14)gydF4y2Ba
的意思是gydF4y2Ba
gydF4y2Ba
(gydF4y2Ba
腹肌gydF4y2Ba
(gydF4y2Ba
圣gydF4y2Ba
)gydF4y2Ba
)gydF4y2Ba
=gydF4y2Ba
的意思是gydF4y2Ba
gydF4y2Ba
(gydF4y2Ba
腹肌gydF4y2Ba
(gydF4y2Ba
圣gydF4y2Ba
- - - - - -gydF4y2Ba
的意思是gydF4y2Ba
(gydF4y2Ba
圣gydF4y2Ba
)gydF4y2Ba
)gydF4y2Ba
)gydF4y2Ba
。gydF4y2Ba
特征提取后,总特性集gydF4y2Ba
年代gydF4y2Ba转换步骤将包含(5(在步骤1中功能)* 5(步骤2)特性)+ 2(额外的noniterative特性)= 27(每个窗口特征列)。gydF4y2Ba
3.3。脑电图分解和分析:使用小波变换gydF4y2Ba
EEG信号表现出复杂的行为和非线性动力学。过去做了广泛的工作理解相关的复杂性与大脑通过多个窗口的数学、物理、工程、化学、生理学等等(gydF4y2Ba
27gydF4y2Ba,gydF4y2Ba
28gydF4y2Ba]。获取和分析脑电图的意图在这个研究是开发一种基准注意识别。本研究的关键是看心电图信号,可以收集到一个便携式臂章可以相对有效识别个人的注意力和焦点。gydF4y2Ba
小而复杂的不同频率结构在参与这项发现脑电图波形包含详细neuroelectric毫秒时间框架的底层信息处理系统,很多研究表明,波形结构在不同的尺度上保存重要的基础和临床信息(gydF4y2Ba
29日gydF4y2Ba,gydF4y2Ba
30.gydF4y2Ba]。小规模神经节奏,特别是与事件相关的振荡厄洛斯,被视为基本的知觉和认知gydF4y2Ba
29日gydF4y2Ba]。小波分析提供了一种有力的方法分离这样的节奏进行研究。有几个脑电图分析小波变换的应用。它已经被用于去除噪声小波系数以来从原始脑电图波形方便精确的噪声过滤机制通过归零或衰减系数主要与噪声在重构信号小波合成(gydF4y2Ba
31日gydF4y2Ba- - - - - -gydF4y2Ba
33gydF4y2Ba]。小波分析脑电图也被广泛用于信号处理应用程序在智能检测系统用于临床gydF4y2Ba
34gydF4y2Ba,gydF4y2Ba
35gydF4y2Ba]。小波变换也被用于压缩脑电图信号。小波压缩技术已经被证明改善neuroelectric数据压缩比与小信号信息的损失(gydF4y2Ba
36gydF4y2Ba,gydF4y2Ba
37gydF4y2Ba]。也可以看到为组件和事件检测以及峰值和瞬态检测脑电图波形。小波分析在许多研究已经证明非常有效(gydF4y2Ba
33gydF4y2Ba- - - - - -gydF4y2Ba
38gydF4y2Ba]。gydF4y2Ba
3.3.1。小波变换的数学公式gydF4y2Ba
小波变换本质上存在两种不同的类型:连续小波变换(CWT)和离散小波变换(DWT)。在这项研究中分析EEG信号的DWT方法已被采用。使用DWT的优点是它允许信号通过应用分析只有离散值的改变和扩展形成离散小波。同样,如果原始信号与一组合适的取样比例和变化值,可以重建整个连续信号的DWT(使用Inverse-DWT)。一个自然的方式设置参数gydF4y2Ba
一个gydF4y2Ba
(扩展)和gydF4y2Ba
bgydF4y2Ba
(转移)是使用对数的离散化gydF4y2Ba
“一个”gydF4y2Ba规模和链接,分别”之间的步长gydF4y2Ba
b”gydF4y2Ba位置或转变。链接”gydF4y2Ba
b”gydF4y2Ba“gydF4y2Ba
一个“gydF4y2Ba离散步骤采取每个位置”gydF4y2Ba
b。”gydF4y2Ba是成正比的”gydF4y2Ba
一个“gydF4y2Ba规模。这种类型的母小波可以以下表格所示。gydF4y2Ba
母亲离散小波表示:gydF4y2Ba
(15)gydF4y2Ba
ΨgydF4y2Ba
米gydF4y2Ba
,gydF4y2Ba
ngydF4y2Ba
(gydF4y2Ba
tgydF4y2Ba
)gydF4y2Ba
=gydF4y2Ba
1gydF4y2Ba
一个gydF4y2Ba
0gydF4y2Ba
米gydF4y2Ba
(gydF4y2Ba
tgydF4y2Ba
- - - - - -gydF4y2Ba
ngydF4y2Ba
bgydF4y2Ba
0gydF4y2Ba
一个gydF4y2Ba
0gydF4y2Ba
米gydF4y2Ba
一个gydF4y2Ba
0gydF4y2Ba
米gydF4y2Ba
)gydF4y2Ba
,gydF4y2Ba
在哪里gydF4y2Ba
整数的gydF4y2Ba
米gydF4y2Ba
和gydF4y2Ba
ngydF4y2Ba
分别控制小波变化和扩展,gydF4y2Ba
一个gydF4y2Ba
0gydF4y2Ba
指定固定膨胀一步参数集值大于1,gydF4y2Ba
bgydF4y2Ba
0gydF4y2Ba
是位置参数必须大于零。gydF4y2Ba
分析方程(DWT):gydF4y2Ba
(16)gydF4y2Ba
WgydF4y2Ba
米gydF4y2Ba
ngydF4y2Ba
=gydF4y2Ba
∫gydF4y2Ba
- - - - - -gydF4y2Ba
∞gydF4y2Ba
+gydF4y2Ba
∞gydF4y2Ba
xgydF4y2Ba
(gydF4y2Ba
tgydF4y2Ba
)gydF4y2Ba
ΨgydF4y2Ba
米gydF4y2Ba
ngydF4y2Ba
*gydF4y2Ba
(gydF4y2Ba
tgydF4y2Ba
)gydF4y2Ba
dgydF4y2Ba
tgydF4y2Ba
。gydF4y2Ba
合成方程(逆DWT):gydF4y2Ba
(17)gydF4y2Ba
xgydF4y2Ba
(gydF4y2Ba
tgydF4y2Ba
)gydF4y2Ba
=gydF4y2Ba
cgydF4y2Ba
∑gydF4y2Ba
米gydF4y2Ba
∑gydF4y2Ba
ngydF4y2Ba
WgydF4y2Ba
米gydF4y2Ba
ngydF4y2Ba
ΨgydF4y2Ba
米gydF4y2Ba
ngydF4y2Ba
(gydF4y2Ba
tgydF4y2Ba
)gydF4y2Ba
,gydF4y2Ba
在哪里gydF4y2Ba
cgydF4y2Ba
是一个常数与母亲小波。gydF4y2Ba
3.3.2。应用程序DWT的脑电图gydF4y2Ba
在这项研究中,应用离散小波变换和DWT的脑电图乐队组件中提取的预处理步骤。gydF4y2Ba
如图gydF4y2Ba
7gydF4y2Ba,每个提取的乐队组件发送通过“窗口”的一步。在这一步的有趣和无聊的部分带组件基于时间戳的原始脑电图提取并通过窗口发送机制。在这种机制中,每个乐队分量信号分割成小窗户。窗户是10秒长,不重叠的。脑电图信号获得的采样率为1000 Hz,所以每个窗口都有1000 Hz * 10秒= 10000数据点。gydF4y2Ba
脑电图使用小波变换分解和分析步骤。gydF4y2Ba
然后分解利用DWT每个窗口。小波变换的性能取决于所选择的母小波分解的信号。常见的启发式是选择一个感兴趣的信号的形状相似。对于组乐队组件从原始EEG信号中提取不同的母亲在分解小波适合不同波段的应用。gydF4y2Ba
如图gydF4y2Ba
8gydF4y2Baγ波分量的分析,母亲小波选择“bior3.9”双正交小波的家庭。δθ,α波组件分解使用“db4”作为他们的母亲Daubechies家族的小波的小波。最后β波分解使用“coif3”Coiflets母小波的小波家族。这些小波选择不仅基于形状和复杂性,也因为他们似乎是常用的等在相关研究中的应用。gydF4y2Ba
(一)“COIF3”小波,(b)“DB4小波,(c)“BOIR3.9”小波。gydF4y2Ba
小波变换的分解进程可以执行迭代分为几个层次。的数量水平为分解是特定于应用程序的选择,也取决于信号的复杂性。脑电图信号频带的窗口组件,5个层次分解似乎提供所有需要的有用的信息;进一步分解并没有产生一个更好的结果。的所有阶段的详细系数1到5和5级的近似系数保留特征提取步骤。gydF4y2Ba
3.3.3。特征提取步骤gydF4y2Ba
这些系数的特性计算如下。(在这里,gydF4y2Ba
(gydF4y2Ba
xgydF4y2Ba
1gydF4y2Ba
,gydF4y2Ba
xgydF4y2Ba
2gydF4y2Ba
,gydF4y2Ba
…gydF4y2Ba
,gydF4y2Ba
xgydF4y2Ba
ngydF4y2Ba
)gydF4y2Ba
代表每个系数的值从10秒窗口。)gydF4y2Ba
标准偏差:gydF4y2Ba
(18)gydF4y2Ba
性病gydF4y2Ba
=gydF4y2Ba
1gydF4y2Ba
ngydF4y2Ba
∑gydF4y2Ba
我gydF4y2Ba
=gydF4y2Ba
1gydF4y2Ba
ngydF4y2Ba
xgydF4y2Ba
我gydF4y2Ba
2gydF4y2Ba
。gydF4y2Ba
熵,熵是一个随机性的统计测量。它是非常有用的在评估中的信息信号:gydF4y2Ba
(19)gydF4y2Ba
熵gydF4y2Ba
=gydF4y2Ba
- - - - - -gydF4y2Ba
总和gydF4y2Ba
(gydF4y2Ba
pgydF4y2Ba
*gydF4y2Ba
日志gydF4y2Ba
gydF4y2Ba
2gydF4y2Ba
(gydF4y2Ba
pgydF4y2Ba
)gydF4y2Ba
)gydF4y2Ba
,gydF4y2Ba
在哪里gydF4y2Ba
pgydF4y2Ba
是信号的直方图。gydF4y2Ba
方差的日志:让每个元素的概率质量函数如下gydF4y2Ba
xgydF4y2Ba
1gydF4y2Ba
↦gydF4y2Ba
pgydF4y2Ba
1gydF4y2Ba
,gydF4y2Ba
…gydF4y2Ba
,gydF4y2Ba
xgydF4y2Ba
ngydF4y2Ba
↦gydF4y2Ba
pgydF4y2Ba
ngydF4y2Ba
,然后gydF4y2Ba
(20)gydF4y2Ba
方差gydF4y2Ba
=gydF4y2Ba
∑gydF4y2Ba
我gydF4y2Ba
=gydF4y2Ba
1gydF4y2Ba
ngydF4y2Ba
pgydF4y2Ba
我gydF4y2Ba
*gydF4y2Ba
(gydF4y2Ba
xgydF4y2Ba
我gydF4y2Ba
- - - - - -gydF4y2Ba
μgydF4y2Ba
)gydF4y2Ba
2gydF4y2Ba
,gydF4y2Ba
在哪里gydF4y2Ba
μgydF4y2Ba
是期望值,gydF4y2Ba
(21)gydF4y2Ba
μgydF4y2Ba
=gydF4y2Ba
∑gydF4y2Ba
我gydF4y2Ba
=gydF4y2Ba
1gydF4y2Ba
ngydF4y2Ba
pgydF4y2Ba
我gydF4y2Ba
*gydF4y2Ba
xgydF4y2Ba
我gydF4y2Ba
。gydF4y2Ba
因此gydF4y2Ba
,gydF4y2Ba
日志gydF4y2Ba
的gydF4y2Ba
方差gydF4y2Ba
=gydF4y2Ba
日志gydF4y2Ba
gydF4y2Ba
2gydF4y2Ba
gydF4y2Ba
(gydF4y2Ba
方差gydF4y2Ba
(gydF4y2Ba
xgydF4y2Ba
)gydF4y2Ba
)gydF4y2Ba
。gydF4y2Ba
意思的频率(离散傅里叶域中):gydF4y2Ba
(22)gydF4y2Ba
dftgydF4y2Ba
(gydF4y2Ba
xgydF4y2Ba
kgydF4y2Ba
)gydF4y2Ba
=gydF4y2Ba
∑gydF4y2Ba
kgydF4y2Ba
=gydF4y2Ba
1gydF4y2Ba
NgydF4y2Ba
- - - - - -gydF4y2Ba
1gydF4y2Ba
XgydF4y2Ba
(gydF4y2Ba
jgydF4y2Ba
)gydF4y2Ba
egydF4y2Ba
jgydF4y2Ba
(gydF4y2Ba
2gydF4y2Ba
πgydF4y2Ba
/gydF4y2Ba
NgydF4y2Ba
)gydF4y2Ba
kgydF4y2Ba
ngydF4y2Ba
,gydF4y2Ba
净的gydF4y2Ba
NgydF4y2Ba时间样本、dft (gydF4y2Ba
xgydF4y2Ba
kgydF4y2Ba
),表示正弦和余弦分量的大小的样本gydF4y2Ba
egydF4y2Ba
jgydF4y2Ba
(gydF4y2Ba
2gydF4y2Ba
πgydF4y2Ba
/gydF4y2Ba
NgydF4y2Ba
)gydF4y2Ba
kgydF4y2Ba
ngydF4y2Ba
:gydF4y2Ba
(23)gydF4y2Ba
的意思是gydF4y2Ba
的gydF4y2Ba
傅里叶gydF4y2Ba
域gydF4y2Ba
=gydF4y2Ba
的意思是gydF4y2Ba
(gydF4y2Ba
dftgydF4y2Ba
(gydF4y2Ba
xgydF4y2Ba
)gydF4y2Ba
)gydF4y2Ba
。gydF4y2Ba
概率分布的方差:gydF4y2Ba
(24)gydF4y2Ba
概率gydF4y2Ba
分布gydF4y2Ba
函数gydF4y2Ba
=gydF4y2Ba
PgydF4y2Ba
(gydF4y2Ba
一个gydF4y2Ba
≤gydF4y2Ba
xgydF4y2Ba
≤gydF4y2Ba
bgydF4y2Ba
]gydF4y2Ba
=gydF4y2Ba
∫gydF4y2Ba
一个gydF4y2Ba
bgydF4y2Ba
fgydF4y2Ba
(gydF4y2Ba
xgydF4y2Ba
)gydF4y2Ba
dgydF4y2Ba
xgydF4y2Ba
方差gydF4y2Ba
的gydF4y2Ba
分布gydF4y2Ba
=gydF4y2Ba
方差gydF4y2Ba
(gydF4y2Ba
PgydF4y2Ba
)gydF4y2Ba
。gydF4y2Ba
的自相关:gydF4y2Ba
(25)gydF4y2Ba
自相关gydF4y2Ba
函数gydF4y2Ba
=gydF4y2Ba
RgydF4y2Ba
(gydF4y2Ba
年代gydF4y2Ba
,gydF4y2Ba
tgydF4y2Ba
)gydF4y2Ba
=gydF4y2Ba
EgydF4y2Ba
(gydF4y2Ba
(gydF4y2Ba
xgydF4y2Ba
tgydF4y2Ba
- - - - - -gydF4y2Ba
μgydF4y2Ba
)gydF4y2Ba
*gydF4y2Ba
(gydF4y2Ba
xgydF4y2Ba
tgydF4y2Ba
+gydF4y2Ba
rgydF4y2Ba
- - - - - -gydF4y2Ba
μgydF4y2Ba
)gydF4y2Ba
]gydF4y2Ba
σgydF4y2Ba
tgydF4y2Ba
σgydF4y2Ba
年代gydF4y2Ba
,gydF4y2Ba
在哪里gydF4y2Ba
年代gydF4y2Ba
和gydF4y2Ba
tgydF4y2Ba
在时间序列不同的时间,gydF4y2Ba
μgydF4y2Ba
的意思是gydF4y2Ba
XgydF4y2Ba
,gydF4y2Ba
σgydF4y2Ba
的标准偏差gydF4y2Ba
XgydF4y2Ba
,“gydF4y2Ba
EgydF4y2Ba
“是期望值算子:gydF4y2Ba
(26)gydF4y2Ba
总和gydF4y2Ba
的gydF4y2Ba
自相关gydF4y2Ba
=gydF4y2Ba
总和gydF4y2Ba
(gydF4y2Ba
RgydF4y2Ba
(gydF4y2Ba
年代gydF4y2Ba
,gydF4y2Ba
tgydF4y2Ba
)gydF4y2Ba
)gydF4y2Ba
。gydF4y2Ba
自协方差的意思是:gydF4y2Ba
(27)gydF4y2Ba
CgydF4y2Ba
(gydF4y2Ba
年代gydF4y2Ba
,gydF4y2Ba
tgydF4y2Ba
)gydF4y2Ba
=gydF4y2Ba
EgydF4y2Ba
(gydF4y2Ba
(gydF4y2Ba
xgydF4y2Ba
tgydF4y2Ba
- - - - - -gydF4y2Ba
μgydF4y2Ba
tgydF4y2Ba
)gydF4y2Ba
*gydF4y2Ba
(gydF4y2Ba
xgydF4y2Ba
年代gydF4y2Ba
- - - - - -gydF4y2Ba
μgydF4y2Ba
年代gydF4y2Ba
)gydF4y2Ba
]gydF4y2Ba
,gydF4y2Ba
在哪里gydF4y2Ba
年代gydF4y2Ba
和gydF4y2Ba
tgydF4y2Ba
在时间序列不同的时间,gydF4y2Ba
μgydF4y2Ba
的意思是gydF4y2Ba
XgydF4y2Ba
,“gydF4y2Ba
EgydF4y2Ba
“是期望值算子:gydF4y2Ba
(28)gydF4y2Ba
的意思是gydF4y2Ba
的gydF4y2Ba
自相关gydF4y2Ba
=gydF4y2Ba
的意思是gydF4y2Ba
(gydF4y2Ba
CgydF4y2Ba
(gydF4y2Ba
年代gydF4y2Ba
,gydF4y2Ba
tgydF4y2Ba
)gydF4y2Ba
)gydF4y2Ba
。gydF4y2Ba
特征提取后执行,总特性集将包含小波变换步骤;6系数近似详细(5 + 1)* 7(每个系数特性)= 42(每个频带特性列组件)。总共有5个乐队中提取组件,因此,42(每个频带特性组件)* 5(不同的乐队组件)= 210 (EEG)的总特征。这些计算功能然后发送到机器学习阶段分类、训练和测试。gydF4y2Ba
3.4。机器学习和分类模型gydF4y2Ba
在这个应用程序的结果在各种信号处理后获得心理上的信号是一个大型的特性集。由于数据收集系统和控制环境,从各自的部分的信号提取的特征可分为以下两个假定类别:“关注”和“nonattention。“因此监督学习方法用于本研究开发分类启发式。gydF4y2Ba
三个不同的机器学习算法实现了对这次试验和测试。这些都是如下。gydF4y2Ba
3.4.1。通过回归分类gydF4y2Ba
连续变量有不同的模型来预测或分类变量的连续预测和/或分类因素影响如一般线性模型(glm)和一般回归模型(grm)。回归问题是那些尝试预测一个连续变量的值从一个或多个连续的和/或分类预测变量(gydF4y2Ba
28gydF4y2Ba,gydF4y2Ba
38gydF4y2Ba,gydF4y2Ba
39gydF4y2Ba]。这是一种非参数方法意味着没有分布假设的数据而漠视它是已知或假定数据遵循特定的线性模型二项或泊松等。在回归分类器,将决策树是由基于变量最好区分目标的类别分类标签变量。这里的决定将由基于回归树。回归树中的每个节点分为两个子节点。随着回归树一定停止规则应用于阻止树木生长。gydF4y2Ba
更普遍而言,通过树构建算法分析的目的是确定一组假设的逻辑(分裂)条件允许准确预测或分类的病例。树分类技术,如果应用得当,产生准确的预测或预测分类基于一些合理假设的条件。回归树分类器的优势在许多的另一种技术是,他们生产简单输出分类结果。这种简单性不仅是有用的为目的的快速分类的新观测,但是也可以提供一个更简单的“模型”解释了为什么观察进行分类或预测一个特定的方式。计算分类和回归树的过程可以描述为包括四个基本步骤:指定的标准预测准确性,选择分裂,决定何时停止分裂,并选择适当的树。gydF4y2Ba
3.4.2。C4.5分类方法gydF4y2Ba
C4.5也是decision-tree-based分类算法,由昆兰(gydF4y2Ba
39gydF4y2Ba,gydF4y2Ba
40gydF4y2Ba]。它开发了基于ID3机器学习算法的基本原理(gydF4y2Ba
41gydF4y2Ba]。的决策树C4.5计算输入数据形成一个基于分治策略。在C4.5树中的每个节点与一组相关联的病例。每次分配权重处理未知的属性值。首先整个训练集开始作为根的权重分配给所有1.0例。从这里树计算每个属性的信息增益的训练集。对离散属性信息增益相对于分裂不同值的情况下在每一个节点。选择最高的属性信息增益作为测试节点。后由递归分割的分而治之的策略属性在每个节点形式的孩子节点根据每个节点的属性的信息增益。C4.5已用于多个应用程序在医疗信息gydF4y2Ba
42gydF4y2Ba,gydF4y2Ba
43gydF4y2Ba]。gydF4y2Ba
3.4.3。随机森林gydF4y2Ba
Breiman发达随机森林分类方法基本上就是一个分类器,由多个决策树(gydF4y2Ba
44gydF4y2Ba]。这是一个非常准确的分类器显示巨大的成功与多个数据集。这是特别有用超大数据集和数据库与数据挖掘。与其他两个提到的基于树的随机森林分类器使用多个树或森林开发决策和分类。虽然在这项研究中它被用来开发基于监督数据模型,随机森林也可以用于无监督数据学习(gydF4y2Ba
45gydF4y2Ba,gydF4y2Ba
46gydF4y2Ba]。随机森林也很受欢迎在biosignal和生物医学应用(gydF4y2Ba
46gydF4y2Ba]。gydF4y2Ba
上述机器学习方法都是已知的神经网络等方法在生理和医学应用gydF4y2Ba
47gydF4y2Ba]。此外,神经网络等方法,在分析使用统计学习理论,证明是容易过度拟合的问题(gydF4y2Ba
48gydF4y2Ba- - - - - -gydF4y2Ba
50gydF4y2Ba),因此进一步鼓励使用上面描述的方法,特别是当数据或对象的数量用于训练和测试是有限的。gydF4y2Ba
在机器学习的步骤中,所提到的三个分类器是独立的提取特征上实现心电图和脑电图和每一个分类器的结果进行了比较。这是较早开发的基于一个设置在初始阶段的这个实验。这个实验ECG信号来自21个主题和脑电图信号来自12个科目已被收集起来。gydF4y2Ba