CN102727223B - 基于多生理参数pca融合的脑力负荷测量方法 - Google Patents
基于多生理参数pca融合的脑力负荷测量方法 Download PDFInfo
- Publication number
- CN102727223B CN102727223B CN201210201512.5A CN201210201512A CN102727223B CN 102727223 B CN102727223 B CN 102727223B CN 201210201512 A CN201210201512 A CN 201210201512A CN 102727223 B CN102727223 B CN 102727223B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- matrix
- cov
- variables
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 230000003340 mental effect Effects 0.000 title claims abstract description 86
- 238000000034 method Methods 0.000 title claims abstract description 31
- 238000000513 principal component analysis Methods 0.000 title claims abstract description 28
- 210000001747 pupil Anatomy 0.000 claims abstract description 27
- 238000005516 engineering process Methods 0.000 claims abstract description 14
- 238000004364 calculation method Methods 0.000 claims abstract description 12
- 239000011159 matrix material Substances 0.000 claims description 61
- 230000004927 fusion Effects 0.000 claims description 42
- 239000013598 vector Substances 0.000 claims description 27
- 230000003595 spectral effect Effects 0.000 claims description 21
- 238000001514 detection method Methods 0.000 claims description 20
- 238000009826 distribution Methods 0.000 claims description 17
- 238000000605 extraction Methods 0.000 claims description 14
- 230000008859 change Effects 0.000 claims description 13
- 238000005259 measurement Methods 0.000 claims description 12
- 238000003909 pattern recognition Methods 0.000 claims description 10
- 238000007619 statistical method Methods 0.000 claims description 4
- 230000003247 decreasing effect Effects 0.000 claims description 3
- 238000002474 experimental method Methods 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 238000013519 translation Methods 0.000 claims description 3
- 230000017105 transposition Effects 0.000 claims description 3
- 238000000354 decomposition reaction Methods 0.000 claims description 2
- 239000002184 metal Substances 0.000 abstract 2
- 238000004519 manufacturing process Methods 0.000 abstract 1
- 210000004556 brain Anatomy 0.000 description 9
- 238000000691 measurement method Methods 0.000 description 6
- 230000008901 benefit Effects 0.000 description 5
- 238000011160 research Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 230000029058 respiratory gaseous exchange Effects 0.000 description 3
- 230000007547 defect Effects 0.000 description 2
- 206010016256 fatigue Diseases 0.000 description 2
- 238000012847 principal component analysis method Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 241001364929 Havel River virus Species 0.000 description 1
- 208000019914 Mental Fatigue Diseases 0.000 description 1
- 208000006550 Mydriasis Diseases 0.000 description 1
- 238000012356 Product development Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 210000000467 autonomic pathway Anatomy 0.000 description 1
- 230000036772 blood pressure Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 230000008451 emotion Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000002595 magnetic resonance imaging Methods 0.000 description 1
- 210000003205 muscle Anatomy 0.000 description 1
- 210000005036 nerve Anatomy 0.000 description 1
- 210000002589 oculomotor nerve Anatomy 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000006461 physiological response Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000002207 retinal effect Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 210000005070 sphincter Anatomy 0.000 description 1
- 230000035900 sweating Effects 0.000 description 1
- 230000002889 sympathetic effect Effects 0.000 description 1
Landscapes
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
Abstract
本发明涉及医疗器械领域。为有效地提高脑力负荷检测***准确性和简便性,为达到上述目的,本发明采取的技术方案是,基于多生理参数PCA融合的脑力负荷测量方法,包括如下步骤:测量心率变异性HRV、瞳孔直径、皮肤电阻SR三个生理参数,利用PCA技术得出三个参数的权重系数,根据参数融合计算公式计算脑力负荷的参数融合分值MWS,MWS是mentalworkload?score的缩写,MWS等于各参数与其权重之积的和,并将MWS作为脑力负荷的测量指标。
Description
技术领域
本发明涉及医疗器械领域,具体讲,涉及应用于医疗器械领域的基于多生理参数PCA融合的脑力负荷测量方法。
背景技术
脑力负荷,指作业人员为达到业绩标准而付出的注意力大小,其涉及到完成某项任务时的工作要求、时间压力、作业人员的能力和努力程度,以及任务不顺利时的挫折感等。低脑力负荷会使人感到厌倦,注意力不集中;高强度的脑力负荷会影响作业者工作绩效,如此长时间的脑力负荷会使人产生脑力疲劳,而疲劳会使人反应迟钝,灵活性低,容易被无关信息干扰,出错率增加,重复犯错。
因为脑力负荷是一个多维的概念,所以脑力负荷的测量方法是多种多样的。目前,脑力负荷测量的方法归纳起来主要有如下4类:(1)主观测量法:是以劳动者对作业或***功能的成绩判断为基础建立的一些心理学方法,如主观劳动负荷测量技术和作业负荷指数等;(2)主任务测量法:也叫工作成绩测量,通过操作者完成作业或***功能的成绩评价劳动负荷,如完成作业的负荷量,作业速度,时间和成绩,错误率等;(3)次任务测量法;(4)生理学测量法:通过作业者对***或作业需要的生理反应进行评价,如心率及其变异性、呼吸、眨眼频率、瞳孔直径、皮肤电阻、眼电图、脑事件相关电位、脑地形图、脑磁图、磁共振成像、正电子发射扫描等。由于生理参数的客观、实时性,因此生理学测量法是目前测量脑力负荷的常用方法。
在当前的研究中,学者们多采用单一参数如EEG、fNIRS、ERP等作为检测脑力负荷的指标,但由于现实工作中任务的复杂性,单一测量指标不能全面反映脑力负荷的变化。研究表明,不同的生理测量技术能从不同角度表现出对脑力负荷的敏感性。任何单一的生理指标对脑力负荷的测量都是片面的,只有多种生理指标的综合运用才能全面反映脑力负荷的变化。另外,由于某些生理参数如EEG、ERP等采集相对复杂,在实际应用中受到限制。因此基于多参数的,简单便捷的脑力负荷检测方法还在不断的探索之中。
发明内容
本发明旨在克服现有技术的不足,有效地提高脑力负荷检测***准确性和简便性,为达到上述目的,本发明采取的技术方案是,基于多生理参数PCA融合的脑力负荷测量方法,包括如下步骤:测量心率变异性HRV、瞳孔直径、皮肤电阻SR三个生理参数,利用PCA技术得出三个参数的权重系数,根据参数融合计算公式计算脑力负荷的参数融合分值MWS,MWS是mental workload score的缩写,MWS等于各参数与其权重之积的和,并将MWS作为脑力负荷的测量指标。
测量心率变异性HRV、瞳孔直径、皮肤电阻SR三个生理参数是:
1)心率变异性HRV提取:采用傅立叶变换FFT计算HRV的低频频率LF、高频频率HF,将LF/HF即低频与高频的比值记为P1,将总功率TP的功率谱密度记为P2,其中LF频率变化在0.04~0.15HZ,HF的频率变化在0.15~0.40HZ,TP表示频率在0.00~0.40HZ的总频率变化值;
2)瞳孔直径特征提取:采用AR模型计算瞳孔直径的功率谱密度记为P3;
3)SR特征提取:采用傅里叶变换FFT,计算皮肤电阻在0.03~0.5HZ的功率谱密度记为P4。利用PCA技术得出三个参数的权重系数是,将提取的四个特征变量低频与高频的比值P1、总功率TP的功率谱密度P2、瞳孔直径的功率谱密度P3和皮肤电阻的功率谱密度P4,这四个特征变量构成一个四维空间,样本数为n,对在该四维空间下的所测样本进行变换:设其原始变量的坐标系为P1a、P2a、、P3a、、P4a,在对原始坐标经过坐标平移、尺度伸缩、旋转变换后,得到一组新的、相互正交的坐标轴v1、v2、v3、v4,根据原始变量在新坐标系上投影值的方差来确定这四个特征变量的权重系数W1、W2、W3、W4:
首先对该4个特征变量进行n次观测得到的观测数据可用下面的矩阵表示:
用主成分分析求取特征变量的权重系数W1、W2、W3、W4的步骤如下:
(1)对原始数据矩阵P进行标准化处理,以消除其量纲、数量级上的差异,使其具有可比性,然后,用矩阵中的元素减去所在列的均值,然后除以所在列的标准差,使得原始数据矩阵P中每个特征变量变为均值为0,方差为1,得到矩阵Y:
Y=[yij]n×4,i=1,2,…,n(2-2)
其中, pij、pij表示观测矩阵P中的元素,表示观测矩阵P所在列的均值;
(2)对标准化后的矩阵Y求协方差矩阵,Y的每一列对应一个变量的n个测量值,任意两列之间可以计算两变量间的协方差,得到协方差矩阵
其中k=1,2,3,4,和分别为Y中第k列和第m列的均值,当k=m时,
(3)特征分解:计算协方差矩阵Z的特征值和特征向量,由式(2-6)
|Z-λI|=0(2-6)
求出协方差矩阵Z的4个特征值并将其按照由大到小排列,λ1≥λ2≥λ3≥λ4,特征值对应的特征向量分别为U1,U2,U3,U4,则协方差矩阵Z可以写成下式:
Z=UΛUT(2-7)
其中,Λ-Z的特征值按照由大到小所组成的对角阵,T表示转置,λ表示协方差矩阵Z的特征值,I表示与Z相应的单位矩阵,Λ表示由Z的特征值组成的对角矩阵;
U-Z的特征向量按列组成的正交阵,它构成了新的矢量空间,作为新变量即主成分的坐标轴,又称为载荷轴;特征值表示新变量即主成分方差的大小;得到的特征向量的方差比前一个特征向量的更小,也就是依次递减;特征向量相互正交,即不相关;
(4)求主成分得分——新的变量值
Fn×m=Yn×4U4×m(2-8)
矩阵Fn×m的每一行相当于原数据矩阵的所有行即原始变量构成的向量在主成分坐标轴即载荷轴上的投影,这些新的投影构成的向量就是主成分得分向量;
由以上步骤求得这四个特征变量P1、P2、P3、P4对应的权重系数W1、W2、W3和W4;
参数融合计算公式具体为:
MWS=W1P1+W2P2+W3P3+W4P4,将该参数融合值作为其检测指标,在后续的模式识别中作为输入特征向量。
采用模糊模式识别对脑力负荷的大小进行判别,具体包括如下步骤:
(1)模糊特征提取由参数融合计算公式得到的参数融合值MWS,根据一定的模糊化规则把该参数融合值MWS分成多个模糊变量,使每个模糊变量表达参数融合值MWS的某一局部特性,用模糊变量进行模式识别;
(2)隶属函数建立在模糊特征有效抽取的基础上,采用模糊分布的方法,通过对上述模糊变量统计分析描出大致曲线,选择一个与给出的几种分布最接近的一个,再根据实验确定模糊变量的实际值,从而确定隶属函数,模糊分布包括有矩形与半矩形分布、梯形与半梯形分布、正态分布;
(3)根据最大隶属原则对结果进行判别:设由参数融合值MWS组成的集合X为所要识别的对象全体,Ai属于F(X),F(X)表示X上的模糊子集的全体,i=1,2,3,表示3个模糊模式即“低度脑力负荷”、“中度脑力负荷”和“高度脑力负荷”;对于X中任一元素x,设x对于模式Ai的隶属度为μAk(x),要识别它属于哪一个模式,可按下列原则作判断,即若μAk(x)=max{μA1(x),μA2(x),μA3(x)},则认为x相对属于Ak所代表的那一类,从而判断出脑力负荷程度。
本发明的技术特点及效果:
本发明采用三参数测量。解决了脑力负荷检测指标测量复杂,不方便的问题,并且基于PCA的多参数融合,更适合脑力负荷多维的特点,可以更准确,敏感的反应其大小变化,可有效地提高脑力负荷检测***的准确性和简便性,给实际应用带来便捷,并且可应用于多种工作场景。
附图说明
图1基于PCA多参数融合脑力负荷检测流程图。
具体实施方式
提出了由主成分分析(PCA)来确定脑力负荷测量中所记录的多生理参数的权重系数,根据定义的计算公式将多参数测量值融合成一个数值作为脑力负荷的检测指标的方法。其技术流程是:测量心率变异性(HRV)、瞳孔直径、皮肤电阻(SR)三个生理参数,利用PCA技术得出三个参数的权重系数,根据参数融合计算公式MWS(mental workload score)=W1P1+W2P2+W3P3+W4P4计算脑力负荷的参数融合分值,并将该数值作为脑力负荷的测量指标。该方法属于人机交互与人因工程领域的脑力负荷检测新方法。
本发明的主旨是提出一种基于多生理参数PCA融合的脑力负荷测量方法,通过主成分分析方法计算心率变异性(HRV)、皮肤电阻(SR)及瞳孔直径的权重系数,继而将该三参数融合的数值作为脑力负荷的检测指标,从而准确、客观、简便的检测脑力负荷。该项发明可有效地提高脑力负荷检测***准确性和简便性,并获得可观的社会效益和经济效益。
基于PCA多参数融合脑力负荷检测流程图如图1所示。脑力负荷检测的整个流程为:首先采集心率变异性(HRV)、皮肤电阻(SR)、瞳孔直径三个生理参数,然后分别对各参数进行频域分析,提取特征,形成参数矩阵,利用PCA确定各参数特征的权重系数,将权重系数与参数特征融合加成至脑力负荷的特征向量,作为脑力负荷的检测指标,并将其作为后续模式识别的输入,从而对脑力负荷的大小进行识别判断。
1 信号采集
心率变异(Heart rate variability,HRV)是指窦性心率在一定时间内周期性改变的现象。HRV主要受自主神经调节,同时亦受呼吸、血压、皮肤温度、肾素-血管紧张素等多种因素影响。研究表明,心率变异对许多不同难度的操纵敏感,是评价脑力负荷大小的很好指标。HRV典型的频谱可有三个峰,大致位于0.04Hz以下,0.05-0.15Hz和大于0.15Hz,分别称为极低频(VLF)、低频(LF)和高频(HF)峰,在脑力负荷评价中,HF(高频功率)的有效性较差;而LF(低频功率)和LF/HF(低高频功率之比)两个指标的有效性较好。
瞳孔是位于人眼虹膜中央的圆形缺损,直径为2~6mm。瞳孔相当于眼光学***中的孔径光阑,它可以通过放大和缩小来调节进入眼内光线的量,从而影响视网膜像差大小。瞳孔大小由动眼神经支配的瞳孔***和交感神经支配的瞳孔开大肌共同控制,它们彼此在中枢紧密联系并相互拮抗。研究指出,瞳孔面积的变化能够评价脑力负荷的不同:随着工作难度水平的增加,瞳孔缩小;而随着紧张程度的增加,瞳孔扩大;当到达一定疲劳程度之后,瞳孔开始缩小。早期就有很多研究利用这一参数检测脑力负荷的变化,均取得了较理想的效果。
人体的皮肤对电流和电压会呈现一定的电阻,我们称之为皮肤电阻。人体的皮肤电阻受呼吸、情绪、出汗等多种条件的影响。研究表明,皮肤电阻的大小会随着脑力负荷的变化呈现明显的变化,也是测量脑力负荷的另一敏感指标。
在本技术中,由于该三个参数测量简单、方便,测量方法成熟,并且能敏感的反应脑力负荷的大小变化,在实际应用中更具有实时性与优越性。
2 特征提取
4)HRV特征提取:采用FFT计算HRV的LF(低频频率)、HF(高频频率),将(LF/HF)即低频与高频的比值记为P1,将TP(总功率)的功率谱密度记为P2,其中LF频率变化在0.04~0.15HZ,HF的频率变化在0.15~0.40HZ,TP表示频率在0.00~0.40HZ的总频率变化值。
5)瞳孔直径特征提取:人的瞳孔直径一般在0.2~0.8mm之间,采用AR模型计算瞳孔直径的功率谱密度,将此作为该参数的特征,记为P3。
6)SR特征提取:采用傅里叶变换,计算皮肤电阻在0.03~0.5HZ的功率谱密度,作为该参数的特征,记为P4。
3 PCA算法
主成分分析(Principal Component Analysis,PCA)是一种常用的基于变量协方差矩阵对信息进行处理、压缩和抽提的有效方法。主成分分析的目的是压缩变量个数,用较少的变量去解释原始数据中的大部分变量,剔除冗余信息。即将许多相关性很高的变量转化成个数较少、能解释大部分原始数据方差且彼此互相独立的几个新变量,也就是所谓的主成分。这样就可以消除原始变量间存在的共线性,克服由此造成的运算不稳定、矩阵病态等问题。
主成分分析可以看作是坐标变换。在本发明中,我们共提取了四个特征变量,它们分别是低频与高频的比值P1、TP(总功率)的功率谱密度P2、瞳孔直径的功率谱密度P3和皮肤电阻的功率谱密度P4,这四个变量构成一个四维空间。我们对在该四维空间下的所测样本(样本数为n)进行变换,设其原始变量的坐标系为P1a、P2a、P3a、P4a,在对原始坐标经过坐标平移、尺度伸缩、旋转等变换后,得到一组新的、相互正交的坐标轴v1、v2、v3、v4,根据原始变量在新坐标系上投影值的方差来确定这四个变量的权重系数W1、W2、W3、W4。
在本发明中,首先对该4个变量进行n次观测得到的观测数据可用下面的矩阵表示
用主成分分析求取四个变量的权重系数W1、W2、W3、W4的步骤如下:
(1)对原始数据矩阵P进行标准化处理,矩阵中的元素减去所在列的均值,然后除以所在列的标准差,使得每个变量变为均值为0,方差为1,得到矩阵Y:
Y=[yij]n×4,i=1,2,…,n(2-2)
其中, pij、pij表示观测矩阵P中的元素,表示观测矩阵P所在列的均值;
在本发明中确定的四个变量,有不同的量纲、不同的数量级,由于不同的量纲会引起个别变量取值的分散程度差异较大,这时总体方差受方差较大的变量的控制。所以不同量纲、不同数量级的数据不能放在一起直接进行比较,也不能直接用于多元统计分析,需要对这四个变量的数值进行标准化处理,以消除其量纲、数量级上的差异,使其具有可比性。
(2)对标准化后的矩阵Y求协方差矩阵,Y的每一列对应一个变量的n个测量值,任意两列之间可以计算两变量间的协方差,得到协方差矩阵
其中k=1,2,3,4,和分别为Y中第k列和第m列的均值,当k=m时,
(3)特征分解。计算协方差矩阵Z的特征值和特征向量,由式(2-6)
|Z-λI|=0(2-6)
求出4个特征值并将其按照由大到小排列,λ1≥λ2≥λ3≥λ4,特征值对应的特征向量分别为U1,U2,U3,U4,则协方差矩阵Z可以写成下式:
Z=UΛUT(2-7)
其中,Λ-Z的特征值按照由大到小所组成的对角阵,T表示转置,λ表示协方差矩阵Z的特征值,I表示与Z相应的单位矩阵,Λ表示由Z的特征值组成的对角矩阵;
U-Z的特征向量按列组成的正交阵,它构成了新的矢量空间,作为新变量(主成分)的坐标轴,又称为载荷轴。特征值表示新变量(主成分)方差的大小。得到的特征向量的方差比前一个特征向量的更小,也就是依次递减;特征向量相互正交,即不相关。
(4)求主成分得分——新的变量值
Fn×m=Yn×4U4×m(2-8)
F阵的每一行相当于原数据矩阵的所有行(即原始变量构成的向量)在主成分坐标轴(载荷轴)上的投影,这些新的投影构成的向量就是主成分得分向量。
由以上步骤我们就求得这四个变量P1、P2、P3、P4对应的权重系数W1、W2、W3和W4。
该方法的特点:主成分是原变量的线性组合;各个主成分之间互不相关;主成分按照方差从大到小依次排列,第一主成分对应最大的方差(特征值);每个主成分的均值为0,其方差为协方差阵对应的特征值;不同的主成分轴(载荷轴)之间相互正交;原来有4个变量,我们选取4个主成分,这4个主成分的变化可以完全反映原来全部4个变量的变化。
该方法的优点:它能找到表现原始数据阵最重要的变量的组合;通过表示最大的方差能有效地直观反映样本之间的关系;能从这几个主成分的得分来近似反映原始的数据阵的信息,从而更准确的反映脑力负荷的变化。
4 参数融合计算公式
我们将脑力负荷的参数融合值记为MWS(Mental Workload Score),定义计算公式为MWS=W1P1+W2P2+W3P3+W4P4。由所测得的变量值P1、P2、P3、P4和PCA确定的权重系数W1、W2、W3、W4计算得脑力负荷的参数融合值,将该参数融合值作为其检测指标,在后续的模式识别中作为输入特征向量。
5 模式识别
由于脑力负荷的概念有一定的模糊性,对脑力负荷大小的判定没有一个明确的判别边界对其进行精确地分类,所以在该发明中我们采用模糊模式识别对脑力负荷的大小进行判别。
模糊模式识别是运用模糊数学的理论和方法将待识别类、对象作为模糊集或其元素,然后对这些模糊集或元素进行分类。在该技术中,经以上特征融合得到融合值之后,我们对其进行模式识别步骤如下:
(1)模糊特征提取由参数融合计算公式得到四个生理参数的参数融合值MWS,根据一定的模糊化规则把该参数融合值分成多个模糊变量,使每个模糊变量表达原变量的某一局部特性,用这些新的特征代替原来的特征进行模式识别,这样,每个模糊特征的取值实际上是一个新的连续变量,它们表示的不再是具体的数值,而是对于特征程度或状况的描述,在该技术中,我们将脑力负荷分为三个模糊特征:“低度脑力负荷”、“中度脑力负荷”和“高度脑力负荷”。
(2)隶属函数建立在模糊特征有效抽取的基础上,识别***的成功将依赖于建立恰当的描述目标类别模式和输入客体特性的隶属函数,以及设计有效的模糊分类器。我们采用模糊分布的方法,通过上述模糊特征统计分析描出大致曲线,选择一个与给出的几种分布最接近的一个,再根据实验确定实际的参数,从而确定隶属函数。常用的F分布有矩形与半矩形分布、梯形与半梯形分布、正态分布等。
(3)根据最大隶属原则对结果进行判别:设由参数融合值MWS组成的集合X为所要识别的对象全体,Ai属于F(X),F(X)表示X上的模糊子集的全体,i=1,2,3,表示3个模糊模式即“低度脑力负荷”、“中度脑力负荷”和“高度脑力负荷”。对于X中任一元素x,设x对于模式Ai的隶属度为μAk(x),要识别它属于哪一个模式,可按下列原则作判断,即若μAk(x)=max{μA1(x),μA2(x),μA3(x)},则认为x相对属于Ak所代表的那一类,从而判断出脑力负荷程度。
有益效果
本技术解决了脑力负荷检测指标测量复杂,不方便的问题,并且基于PCA的多参数融合,更适合脑力负荷多维的特点,可以更准确,敏感的反应其大小变化,可有效地提高脑力负荷检测***的准确性和简便性,给实际应用带来便捷,并且可应用于多种工作场景。
本发明的主旨是提出一种基于PCA的多参数融合脑力负荷测量方法,通过主成分分析方法计算心率变异性(HRV)、皮肤电阻(SR)及瞳孔直径的权重系数,继而将这三参数融合的数值作为脑力负荷的检测指标,从而准确、客观、简便的检测脑力负荷的大小。该项发明可有效地提高脑力负荷检测***准确性和简便性,并获得可观的社会效益和经济效益。最佳实施方案拟采用专利转让、技术合作或产品开发。由于该技术操作简单,敏感性强,基于该技术开发的产品可应用于多种场景,如汽车驾驶、航天员操作及特定工作环境等。
Claims (2)
1.一种基于多生理参数PCA融合的脑力负荷测量方法,其特征是,包括如下步骤:测量心率变异性HRV、瞳孔直径、皮肤电阻SR三个生理参数,利用PCA技术得出三个参数的权重系数,根据参数融合计算公式计算脑力负荷的参数融合分值MWS,MWS是mental workloadscore的缩写,MWS等于各参数与其权重之积的和,并将MWS作为脑力负荷的测量指标;
其中,测量心率变异性HRV、瞳孔直径、皮肤电阻SR三个生理参数是:
1)心率变异性HRV提取:采用傅立叶变换FFT计算HRV的低频频率LF、高频频率HF,将LF/HF即低频与高频的比值记为P1,将总功率TP的功率谱密度记为P2,其中LF频率变化在0.04~0.15HZ,HF的频率变化在0.15~0.40HZ,TP表示频率在0.00~0.40HZ的总频率变化值;
2)瞳孔直径特征提取:采用AR模型计算瞳孔直径的功率谱密度记为P3;
3)SR特征提取:采用傅里叶变换FFT,计算皮肤电阻在0.03~0.5HZ的功率谱密度记为P4;
其中,利用PCA技术得出三个参数的权重系数是,将提取的四个特征变量低频与高频的比值P1、总功率TP的功率谱密度P2、瞳孔直径的功率谱密度P3和皮肤电阻的功率谱密度P4,这四个特征变量构成一个四维空间,样本数为n,对在该四维空间下的所测样本进行变换:设其原始变量的坐标系为P1a、P2a、P3a、P4a,在对原始坐标经过坐标平移、尺度伸缩、旋转变换后,得到一组新的、相互正交的坐标轴v1、v2、v3、v4,根据原始变量在新坐标系上投影值的方差来确定这四个特征变量的权重系数W1、W2、W3、W4:
首先对该4个特征变量进行n次观测得到的观测数据可用下面的矩阵表示:
用主成分分析求取特征变量的权重系数W1、W2、W3、W4的步骤如下:
(1)对原始数据矩阵P进行标准化处理,以消除其量纲、数量级上的差异,使其具有可比性,然后,用矩阵中的元素减去所在列的均值,然后除以所在列的标准差,使得原始数据矩阵P中每个特征变量变为均值为0,方差为1,得到矩阵Y:
Y=[yij]n×4,i=1,2,…,n (2-2)
其中, pij表示观测矩阵P中的元素,表示观测矩阵P所在列的均值;
(2)对标准化后的矩阵Y求协方差矩阵,Y的每一列对应一个变量的n个测量值,任意两列之间可以计算两变量间的协方差,得到协方差矩阵
其中k=1,2,3,4,和分别为Y中第k列和第m列的均值,当k=m时,
(3)特征分解:计算协方差矩阵Z的特征值和特征向量,由式(2-6)
|Z-λI|=0 (2-6)
求出协方差矩阵Z的4个特征值并将其按照由大到小排列,λ1≥λ2≥λ3≥λ4,特征值对应的特征向量分别为U1,U2,U3,U4,则协方差矩阵Z可以写成下式:
Z=UΛUT (2-7)
其中,Λ表示Z的特征值按照由大到小所组成的对角阵,T表示转置,λ表示协方差矩阵Z的特征值,I表示与Z相应的单位矩阵;
U表示Z的特征向量按列组成的正交阵,它构成了新的矢量空间,作为新变量即主成分的坐标轴,又称为载荷轴;特征值表示新变量即主成分方差的大小;得到的特征向量的方差比前一个特征向量的更小,也就是依次递减;特征向量相互正交,即不相关;
(4)求主成分得分——新的变量值
Fn×m=Yn×4U4×m (2-8)
矩阵Fn×m的每一行相当于原数据矩阵的所有行即原始变量构成的向量在主成分坐标轴即载荷轴上的投影,这些新的投影构成的向量就是主成分得分向量;
由以上步骤求得这四个特征变量P1、P2、P3、P4对应的权重系数W1、W2、W3和W4;
参数融合计算公式具体为:
MWS=W1P1+W2P2+W3P3+W4P4,将该参数融合值作为其检测指标,在后续的模式识别中作为输入特征向量。
2.如权利要求1所述的基于多生理参数PCA融合的脑力负荷测量方法,其特征是,采用模糊模式识别对脑力负荷的大小进行判别,具体包括如下步骤:
(1)模糊特征提取由参数融合计算公式得到的参数融合值MWS,根据一定的模糊化规则把该参数融合值MWS分成多个模糊变量,使每个模糊变量表达参数融合值MWS的某一局部特性,用模糊变量进行模式识别;
(2)隶属函数建立在模糊特征有效抽取的基础上,采用模糊分布的方法,通过对上述模糊变量统计分析描出大致曲线,选择一个与给出的几种分布最接近的一个,再根据实验确定模糊变量的实际值,从而确定隶属函数,模糊分布包括有矩形与半矩形分布、梯形与半梯形分布、正态分布;
(3)根据最大隶属原则对结果进行判别:设由参数融合值MWS组成的集合X为所要识别的对象全体,Ai属于F(X),F(X)表示X上的模糊子集的全体,i=1,2,3,表示3个模糊模式即“低度脑力负荷”、“中度脑力负荷”和“高度脑力负荷”;对于X中任一元素x,设x对于模式Ai的隶属度为μAk(x),要识别它属于哪一个模式,可按下列原则作判断,即若μAk(x)=max{μA1(x),μA2(x),μA3(x)},则认为x相对属于Ak所代表的那一类,从而判断出脑力负荷程度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210201512.5A CN102727223B (zh) | 2012-06-18 | 2012-06-18 | 基于多生理参数pca融合的脑力负荷测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210201512.5A CN102727223B (zh) | 2012-06-18 | 2012-06-18 | 基于多生理参数pca融合的脑力负荷测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102727223A CN102727223A (zh) | 2012-10-17 |
CN102727223B true CN102727223B (zh) | 2015-04-29 |
Family
ID=46984040
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210201512.5A Active CN102727223B (zh) | 2012-06-18 | 2012-06-18 | 基于多生理参数pca融合的脑力负荷测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102727223B (zh) |
Families Citing this family (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103610447B (zh) * | 2013-12-04 | 2015-12-09 | 天津大学 | 一种基于前额脑电信号的脑力负荷在线检测方法 |
CN103784119B (zh) * | 2014-01-24 | 2015-07-08 | 北京航空航天大学 | 基于飞行模拟器的飞行员脑力负荷测定*** |
CN105310703B (zh) * | 2014-07-02 | 2018-01-19 | 北京邮电大学 | 一种基于用户瞳孔直径数据获取主观满意度的方法 |
CN104523237B (zh) * | 2015-01-08 | 2016-08-24 | 中国人民解放军装甲兵工程学院 | 一种移动舱室乘员工作项目调整方法 |
WO2016189711A1 (ja) * | 2015-05-27 | 2016-12-01 | 糧三 齋藤 | 携帯端末用ストレス評価プログラム及び該プログラムを備えた携帯端末 |
CN104983434B (zh) * | 2015-06-17 | 2018-05-01 | 重庆邮电大学 | 基于层次分析法的多参数心理压力评估方法及装置 |
EP4026487A1 (en) | 2016-10-26 | 2022-07-13 | Telefonaktiebolaget LM Ericsson (publ) | Identifying sensory inputs affecting working memory load of an individual |
CN108451526A (zh) * | 2017-02-17 | 2018-08-28 | 祥明大学校产学协力团 | 使用瞳孔反应来检测频域心脏信息的方法及*** |
CN107788970B (zh) * | 2017-11-15 | 2020-08-07 | 中国航天员科研训练中心 | 一种基于多生理参数方法融合的脑力负荷评估方法 |
CN109394225A (zh) * | 2018-09-04 | 2019-03-01 | 国家电网公司华东分部 | 一种多特征融合的人因风险实时检测方法及装置 |
CN111568423B (zh) * | 2019-02-15 | 2023-03-03 | 南京昇翠人工智能有限公司 | 用户呼吸的共振频率估计值的量测方法及装置 |
CN112686495A (zh) * | 2020-12-03 | 2021-04-20 | 中广核工程有限公司 | 核电厂操作员工作负荷评估方法、***、介质及电子设备 |
CN113740277A (zh) * | 2021-10-15 | 2021-12-03 | 北方民族大学 | 一种基于光谱多组分分析的环境安全分析方法 |
CN117643461B (zh) * | 2024-01-30 | 2024-04-02 | 吉林大学 | 基于人工智能的心率智能监测***及方法 |
-
2012
- 2012-06-18 CN CN201210201512.5A patent/CN102727223B/zh active Active
Non-Patent Citations (3)
Title |
---|
D DE WAARD.THE MEASUREMENT OF DRIVERS'MENTAL WORKLOAD.《THE MEASUREMENT OF DRIVERS'MENTAL WORKLOAD》.V STUDIECENTRUM,2006, * |
double1125.主成分分析华东理工版.《百度文库》.2011, * |
S Miyake.Multivariate workload evaluation combining physiological and subjective mureeass.《International journal of psychophysiology》.2001,第40卷(第3期), * |
Also Published As
Publication number | Publication date |
---|---|
CN102727223A (zh) | 2012-10-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102727223B (zh) | 基于多生理参数pca融合的脑力负荷测量方法 | |
Zeng et al. | A LightGBM‐based EEG analysis method for driver mental states classification | |
KR102221264B1 (ko) | 인간 감정 인식을 위한 딥 생리적 정서 네트워크를 이용한 인간 감정 추정 방법 및 그 시스템 | |
Wang et al. | A k-nearest-neighbor classifier with heart rate variability feature-based transformation algorithm for driving stress recognition | |
US20200367800A1 (en) | Method for identifying driving fatigue based on cnn-lstm deep learning model | |
Liu et al. | EEG-based estimation of mental fatigue by using KPCA–HMM and complexity parameters | |
US8862581B2 (en) | Method and system for concentration detection | |
EP0699413B1 (en) | Apparatus and method for analyzing information relating to physical and mental condition | |
CA2962083A1 (en) | System and method for detecting invisible human emotion | |
Ahmadi et al. | Automated detection of driver fatigue from electroencephalography through wavelet-based connectivity | |
Das et al. | Unsupervised approach for measurement of cognitive load using EEG signals | |
Li et al. | A novel method for the evaluation of fashion product design based on data mining | |
CN114533086A (zh) | 一种基于空域特征时频变换的运动想象脑电解码方法 | |
KR20150029969A (ko) | 뇌파 측정을 이용한 감성상태 분류 방법 | |
Mustafa et al. | The analysis of eeg spectrogram image for brainwave balancing application using ann | |
CN114693170A (zh) | 一种基于多源参数融合的飞行员跳伞训练心理评估方法 | |
Shangguan et al. | Feature extraction of EEG signals based on functional data analysis and its application to recognition of driver fatigue state | |
CN113349780A (zh) | 一种情绪设计对在线学习认知负荷影响的评估方法 | |
CN112185493A (zh) | 人格偏好诊断装置及基于该装置的项目推荐*** | |
Maaoui et al. | Unsupervised stress detection from remote physiological signal | |
Rezaee et al. | EEG-based driving fatigue recognition using hybrid deep transfer learning approach | |
Lopes et al. | Ensemble deep neural network for automatic classification of EEG independent components | |
KR20080107961A (ko) | 오감 정보 처리기법 및 뇌파를 활용한 사용자 적응형임상진단/치료시스템 | |
Fan et al. | EEG-TransMTL: A transformer-based multi-task learning network for thermal comfort evaluation of railway passenger from EEG | |
CN117883082A (zh) | 一种异常情绪识别方法、***、设备及介质 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |