CN113143275A - 一种样本与特征质量联合量化评估的脑电疲劳检测方法 - Google Patents
一种样本与特征质量联合量化评估的脑电疲劳检测方法 Download PDFInfo
- Publication number
- CN113143275A CN113143275A CN202110317792.5A CN202110317792A CN113143275A CN 113143275 A CN113143275 A CN 113143275A CN 202110317792 A CN202110317792 A CN 202110317792A CN 113143275 A CN113143275 A CN 113143275A
- Authority
- CN
- China
- Prior art keywords
- formula
- sample
- electroencephalogram
- equation
- electroencephalogram data
- 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.)
- Granted
Links
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/16—Devices for psychotechnics; Testing reaction times ; Devices for evaluating the psychological state
- A61B5/18—Devices for psychotechnics; Testing reaction times ; Devices for evaluating the psychological state for vehicle drivers or machine operators
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7203—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7225—Details of analog processing, e.g. isolation amplifier, gain or sensitivity adjustment, filtering, baseline or drift compensation
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
- A61B5/7264—Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Signal Processing (AREA)
- Public Health (AREA)
- Surgery (AREA)
- Veterinary Medicine (AREA)
- Artificial Intelligence (AREA)
- Psychiatry (AREA)
- General Health & Medical Sciences (AREA)
- Animal Behavior & Ethology (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Physiology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Child & Adolescent Psychology (AREA)
- Fuzzy Systems (AREA)
- Evolutionary Computation (AREA)
- Power Engineering (AREA)
- Mathematical Physics (AREA)
- Developmental Disabilities (AREA)
- Educational Technology (AREA)
- Hospice & Palliative Care (AREA)
- Psychology (AREA)
- Social Psychology (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
本发明提供一种样本与特征质量联合量化评估的脑电疲劳检测方法。本发明步骤如下:1、多个被试者分别在模拟驾驶***下进行脑电数据采集。2、对步骤1所得的所有脑电数据进行预处理和特征提取。3、建立机器学习模型以实现样本与特征质量联合量化评估的脑电疲劳检测。4、求得衡量样本质量的描述因子v、特征的描述因子θ。5、对新的被试者脑电的数据进行疲劳回归预测。本发明将v和θ嵌入最小二乘模型后,获得的衡量样本质量与特征的权重描述因子提供了执行脑电数据样本选择和特征选择的有效工具,对质量较好的样本和特征赋予更高的权重,能够根据脑电数据准确地获取被测者的疲劳情况。
Description
技术领域
本发明属于脑电信号处理技术领域,具体涉及一种样本与特征质量联合量化评估的脑电疲劳检测方法。
背景技术
随着我国交通事业的日益发展,研究更为实用客观的驾驶疲劳检测对于提高主动交通安全具有重要的意义。综合国内外的疲劳检测方法,主要分为主观评价和客观检测两种。主观评价方法主要是通过记录主观调查问卷如皮尔逊疲劳量表、斯坦福睡眠尺度表等,判断是否处于疲劳状态。客观检测领域主要分为:基于车辆行为特征的检测、基于驾驶员特征行为特征的检测、基于生理电特征的检测。其中,基于客观检测方法中的脑电信号作为中枢神经信号活动的表征,具有检测准确率高的特点,是疲劳判定的金标准。但目前脑电采集过程当中由于受到跨时段采集和电极帽位置的影响,导致采集到的脑电数据质量缺乏可靠性。
脑电信号作为非稳态信号,样本中含有较多的噪声,如果能在学习过程中对每个脑电数据刻画样本质量和样本特征的描述因子,实现对脑电数据质量进行区分,选择有利于模型训练的特征,从而提升机器学习模型性能,保证模型具有较好的鲁棒性。
发明内容
本发明的目的在于提供一种样本与特征质量联合量化评估的脑电疲劳检测方法。通过该方法可以获得衡量样本质量好坏的描述因子v与衡量样本特征的描述因子Θ,嵌入最小二乘模型并计算得到量化因子v以刻画样本的质量以及通过自学习得到样本特征的描述因子Θ,避免因大噪声脑电样本对模型鲁棒性的负面影响,使用得到的描述因子来进行样本选择可以得到更高的疲劳检测精度。
本发明的具体步骤如下:
步骤1、对处于逐渐疲劳的过程中的多个被试者进行脑电数据采集。
步骤2、对步骤1所得的所有脑电数据进行预处理和特征提取。处理后的每一组数据均作为一个样本矩阵X。每个样本矩阵X均对应一个标签向量y;标签向量y对应受试者的疲劳程度。
步骤3、建立用于脑电疲劳检测的机器学习模型。
3-1.建立嵌入有描述因子v和θ目标函数如式(1)所示。
式(1)中,描述因子θ根据描述矩阵Θ来确定,其具体关系为Θjj为描述矩阵Θ的第j个对角元素;θj为描述因子θ的第j个元素;w、b分别为权重、偏差;Θ为描述矩阵;xi为样本矩阵X的第i列元素;yi为标签y的第i个元素;vi为描述因子v的第i个元素,描述样本矩阵中各样本的质量;表示2范数的平方计算;γ表示正则项参数;f(λ,v)表示正则项函数。
3-2.建立正则项函数f(λ,v)的表达式如式(2)所示:
式(2)中,n为样本矩阵中的样本个数;λ为正则项函数f(λ,v)中的参数。
由拉格朗日乘子法可求得θ的最优解如(4)式所示:
据此进一步将目标函数式(1)改写成式(5)。
作为优选,步骤2中,通过让被试者用模拟驾驶环境在一条笔直的道路持续模拟驾驶的方式来使得被试者逐渐疲劳。
作为优选,步骤2中所述的标签向量y通过受试者在脑电数据采集过程中眼睛闭合的时长来确定。所述的标签y以一个0-1之间的值来表示。
对式(9)中的b求导并且令导数为0,得到b的更新规则为式(10):
b=(UTU)-1(UTU)(y-XTw) (10)
作为优选,步骤2中预处理的过程如下:
2-1.将脑电数据下采样至200Hz后并对其进行带通滤波至1-50Hz范围;根据5频段法,将其划分为Delta、Theta、Alpha、Beta和Gamma五个频段
2-2.针对这5个频段的脑电数据分别进行时间窗口为4秒且不重叠的短时傅立叶变换,提取微分熵特征h(X)如式(14)所示:
h(X)=-∫xf(x)ln(f(x))dx (14)
式(14)中,X为输入的样本矩阵,c为输入的样本矩阵中的元素;f(x)为概率密度函数;
经过更新后的微分熵特征h(X)如式(15)所示;
式(15)中,σ为概率密度函数的标准差;μ为概率密度函数的期望。
作为优选,所述的脑电数据采集采用17导联,选用5个频段;5个频段分别为1-4Hz、4-8Hz、8-14Hz、14-31Hz、31-50Hz。
本发明具有的有益效果是:
1、本发明获得的衡量样本质量与特征的描述因子为非稳态的脑电数据的样本选择提供了有效工具,通过数学模型对质量较好的样本赋予更高的权重,能够根据脑电数据准确地获取被测者的疲劳情况。
2、本发明针对脑电信号噪声较多的特点,设置了用于刻画样本质量和样本特征的描述因子,其中v是通过对损失函数的度量计算得到,当损失函数越大时,由公式(7)可以得到越小的υi,从而使得噪声样本点在损失函数中所占的比重较小,达到去噪的效果;θ的值是通过自学习得到,由公式(4)可知θi的大小只与相关,当中的值越大时,由公式(4)可知θi的值也越大,而的值由公式(1)可知刚好是对样本矩阵的特征进行刻画,这个过程从数学理论的角度解释了为什么θi可以对样本矩阵进行特征选择,从而对不同脑电数据的质量和特征进行区分,以此来提高机器学习模型对脑电信号的识别效果。
3、脑电数据的采集过程中是由多个电极帽获得,该样本数据受到实验时间与导联位置的影响,每一个导联代表一个特征维度。本发明通过对样本进行特征选择,可以求得对于模型训练更有利的导联位置,即对样本进行特征选择,从而在训练过程中提高模型的鲁棒性与准确性。
附图说明
图1为本发明的流程示意图;
图2为本发明中分别对应样本质量和特征重要程度的两个描述因子的嵌入过程示意图。
具体实施方式
以下结合附图对本发明作进一步说明。
本发明解决在疲劳检测中脑电信号的重要性特征挖掘这一重要问题基于如下的出发点:我们认为在疲劳检测中,脑电信号作为非稳态信号,样本中含有较多的噪声,如果能在学习过程中对每个样本质量进行刻画与每个样本特征维度的选择,以此挑选出有利于模型训练的样本与特征,将获得一个鲁棒性较好的模型。因此,我们可以选择质量较好的样本与特征进行学习,这对提升疲劳检测的准确率具有重要意义。
如图1和2所示,一种样本与特征质量联合量化评估的脑电疲劳检测方法,具体步骤如下:
步骤1、利用模拟驾驶***(本实施例中使用包含明显诱发受试者疲劳的模拟驾驶场景)采集受试者的诱发疲劳状态脑电数据。
对N个受试者在相同的模拟驾驶环境下进行M次脑电数据采集,得到N·M组脑电数据,每一组数据的数据量均为d*n,其中,d是每组数据的维度,n是单次采集获得的与时间有关的脑电数据样本个数。一组数据包括一次采集中获得的多个时刻的脑电数据。每一组数据均作为一个样本矩阵X。每个样本矩阵X均对应一个标签y;标签y对应受试者的疲劳程度,其通过受试者在脑电数据采集过程中眼睛闭合的时长占采集总时长的百分比来确定。眼睛闭合的时长越长表示受试者越疲劳。
统一受试者在测试任务中对应的被测任务场景,每个受试者均完成所有的被测任务场景的测试。在本实施例中,被测任务场景为笔直而单调的道路,实验时间为午饭后进行,以使受试者更容易疲劳。这些不同批次的实验所做的内容是相同的,但是他们在不同日期所处的状态可能是不同的;这反应了脑电的本质特征不会随着时间或实验场景的改变而出现较大改变。这些时间的不同,或者实验时的场景的不同即为差异化部分。
步骤2、对步骤1所得的所有脑电数据进行预处理和特征提取。本发明基于17导联、5频段(Delta(1-4Hz),Theta(4-8Hz),Alpha(8-14Hz),Beta(14-31Hz)和Gamma(31-50Hz)),并提取微分熵特征来进行。在实际应用中,导联数目取决于数据采集时候受试者佩戴的脑电帽;频段的划分也遵循具有生理学意义的5频段划分;脑电信号最常用的特征为功率谱密度和微分熵。人的脑电信号是非常微弱的,这就意味着脑电信号易受干扰,采集的结果很难直接做实验,这就对脑电信号预处理提出了要求:
预处理的过程如下:
2-1.将脑电数据下采样至200Hz后并对其进行带通滤波至1-50Hz范围。根据5频段法,将其划分为Delta、Theta、Alpha、Beta和Gamma五个频段
2-2.将这5个频段的脑电数据分别作为样本矩阵,分别进行时间窗口为4秒且不重叠的短时傅立叶变换,提取微分熵特征。微分熵特征h(X)定义为:
h(X)=-∫xf(x)ln(f(x))dx (14)
式(14)中,X为输入的样本矩阵(即某一个频段的脑电数据),x为输入的样本矩阵中的元素;f(x)为概率密度函数。对于遵循高斯分布的样本矩阵X,其微分熵特征h(X)可以计算为式(15)所示:
式(15)中,σ为概率密度函数的标准差;μ为概率密度函数的期望。
步骤3、建立机器学习模型以实现基于样本与特征质量联合量化评估的脑电疲劳检测。衡量样本质量与特征的描述因子v和θ具有非负与归一化特性;描述因子v的第i个元素vi用以刻画第i个样本的质量;描述因子θ的第i个样本元素θi用以刻画样本的第i维特征的重要程度。那么,将vi和θi嵌入到最小二乘学习模型中,便得到如式(16)所示的样本与特征质量联合量化评估的脑电疲劳检测模型。
3-1.将描述因子v和θ嵌入到最小二乘模型中,得到“样本与特征质量联合量化评估的脑电疲劳检测方法”的目标函数如式(5)所示:
目标函数式(5)是在初始目标函数式(17)的基础上增加一个正则项函数得到;
从式(17)中可以看出,当vi全为0时,目标函数式(17)最小,不符合要求,因此给目标函数加一个正则项,更新目标函数为式(5)。
3-2.建立正则项函数f(λ,v)的表达式如式(2)所示:
式(2)中,n为一个样本矩阵中的样本个数;υi表示任意一个衡量样本质量的描述因子,υi的大小由计算得到,当样本的质量越好时,vi的大小越接近1,样本质量越差则越接近0,每一个样本都有一个描述因子vi;λ为正则项函数f(λ,v)中的参数。
通过固定v、w、b,式(1)可以写成如式(3)所示:
由拉格朗日乘子法可求得θ的最优解如(4)式所示:
式(1)可以被写成
对式(9)中的b求导并且令导数为0,我们可以得到b的更新规则为式(10):
b=(UTU)-1(UTU)(y-XTw) (10)
Claims (6)
1.一种样本与特征质量联合量化评估的脑电疲劳检测方法,其特征在于:步骤1、对处于逐渐疲劳的过程中的多个被试者进行脑电数据采集;
步骤2、对步骤1所得的所有脑电数据进行预处理和特征提取;处理后的每一组数据均作为一个样本矩阵X;每个样本矩阵X均对应一个标签向量y;标签向量y对应受试者的疲劳程度;
步骤3、建立用于脑电疲劳检测的机器学习模型;
3-1.建立嵌入有描述因子v和θ目标函数如式(1)所示;
式(1)中,描述因子θ根据描述矩阵Θ来确定,其具体关系为Θjj为描述矩阵Θ的第j个对角元素;θj为描述因子θ的第j个元素;w、b分别为权重、偏差;Θ为描述矩阵;xi为样本矩阵X的第i列元素;yi为标签y的第i个元素;vi为描述因子v的第i个元素,描述样本矩阵中各样本的质量;表示2范数的平方计算;γ表示正则项参数;f(λ,v)表示正则项函数;
3-2.建立正则项函数f(λ,v)的表达式如式(2)所示:
式(2)中,n为样本矩阵中的样本个数;λ为正则项函数f(λ,v)中的参数;
由拉格朗日乘子法可求得θ的最优解如(4)式所示:
据此进一步将目标函数式(1)改写成式(5);
2.根据权利要求1所述的一种样本与特征质量联合量化评估的脑电疲劳检测方法,其特征在于:步骤2中,通过让被试者用模拟驾驶环境在一条笔直的道路持续模拟驾驶的方式来使得被试者逐渐疲劳。
3.根据权利要求1所述的一种样本与特征质量联合量化评估的脑电疲劳检测方法,其特征在于:步骤2中所述的标签向量y通过受试者在脑电数据采集过程中眼睛闭合的时长来确定;所述的标签y以一个0-1之间的值来表示。
对式(9)中的b求导并且令导数为0,得到b的更新规则为式(10):
b=(UTU)-1(UTU)(y-XTw) (10)
5.根据权利要求1所述的一种样本与特征质量联合量化评估的脑电疲劳检测方法,其特征在于:步骤2中预处理的过程如下:
2-1.将脑电数据下采样至200Hz后并对其进行带通滤波至1-50Hz范围;根据5频段法,将其划分为Delta、Theta、Alpha、Beta和Gamma五个频段
2-2.针对这5个频段的脑电数据分别进行时间窗口为4秒且不重叠的短时傅立叶变换,提取微分熵特征h(X)如式(14)所示:
h(X)=-∫xf(x)ln(f(x))dx (14)
式(14)中,X为输入的样本矩阵,x为输入的样本矩阵中的元素;f(x)为概率密度函数;
经过更新后的微分熵特征h(X)如式(15)所示;
式(15)中,σ为概率密度函数的标准差;μ为概率密度函数的期望。
6.根据权利要求1所述的一种样本与特征质量联合量化评估的脑电疲劳检测方法,其特征在于:所述的脑电数据采集采用17导联,选用5个频段;5个频段分别为1-4Hz、4-8Hz、8-14Hz、14-31Hz、31-50Hz。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110317792.5A CN113143275B (zh) | 2021-03-24 | 2021-03-24 | 一种样本与特征质量联合量化评估的脑电疲劳检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110317792.5A CN113143275B (zh) | 2021-03-24 | 2021-03-24 | 一种样本与特征质量联合量化评估的脑电疲劳检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113143275A true CN113143275A (zh) | 2021-07-23 |
CN113143275B CN113143275B (zh) | 2022-04-26 |
Family
ID=76884692
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110317792.5A Active CN113143275B (zh) | 2021-03-24 | 2021-03-24 | 一种样本与特征质量联合量化评估的脑电疲劳检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113143275B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116746931A (zh) * | 2023-06-15 | 2023-09-15 | 中南大学 | 一种基于脑电的增量式驾驶员不良状态检测方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5846208A (en) * | 1996-09-04 | 1998-12-08 | Siemens Aktiengesellschaft | Method and apparatus for the evaluation of EEG data |
CN101155548A (zh) * | 2005-03-04 | 2008-04-02 | 曼提斯库拉Ehf.公司 | 评估神经病症的方法和*** |
CN105957047A (zh) * | 2016-05-06 | 2016-09-21 | 中国科学院自动化研究所 | 有监督的多模态脑影像融合方法 |
CN111067513A (zh) * | 2019-12-11 | 2020-04-28 | 杭州电子科技大学 | 一种特征权重自学习的睡眠质量检测关键脑区判定方法 |
-
2021
- 2021-03-24 CN CN202110317792.5A patent/CN113143275B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5846208A (en) * | 1996-09-04 | 1998-12-08 | Siemens Aktiengesellschaft | Method and apparatus for the evaluation of EEG data |
CN101155548A (zh) * | 2005-03-04 | 2008-04-02 | 曼提斯库拉Ehf.公司 | 评估神经病症的方法和*** |
CN105957047A (zh) * | 2016-05-06 | 2016-09-21 | 中国科学院自动化研究所 | 有监督的多模态脑影像融合方法 |
CN111067513A (zh) * | 2019-12-11 | 2020-04-28 | 杭州电子科技大学 | 一种特征权重自学习的睡眠质量检测关键脑区判定方法 |
Non-Patent Citations (2)
Title |
---|
WEI-LONG ZHENG: "《Identifying Stable Patterns over Time for Emotion Recognition from EEG》", 《 IEEE TRANSACTIONS ON AFFECTIVE COMPUTING》 * |
WEI-LONG ZHENG: "《Investigating Critical Frequency Bands and Channels for EEG-Based Emotion Recognition with Deep Neural Networks》", 《 IEEE TRANSACTIONS ON AUTONOMOUS MENTAL DEVELOPMENT 》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116746931A (zh) * | 2023-06-15 | 2023-09-15 | 中南大学 | 一种基于脑电的增量式驾驶员不良状态检测方法 |
CN116746931B (zh) * | 2023-06-15 | 2024-03-19 | 中南大学 | 一种基于脑电的增量式驾驶员不良状态检测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113143275B (zh) | 2022-04-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113627518B (zh) | 利用迁移学习实现神经网络脑电情感识别模型的方法 | |
US20200367800A1 (en) | Method for identifying driving fatigue based on cnn-lstm deep learning model | |
CN108960182B (zh) | 一种基于深度学习的p300事件相关电位分类识别方法 | |
CN110070105B (zh) | 基于元学习实例快速筛选的脑电情绪识别方法、*** | |
CN108186011A (zh) | 房颤检测方法、装置及可读存储介质 | |
CN106108894A (zh) | 一种提高情绪识别模型时间鲁棒性的情绪脑电识别方法 | |
CN113052113B (zh) | 一种基于紧凑型卷积神经网络的抑郁症识别方法及*** | |
CN102727223B (zh) | 基于多生理参数pca融合的脑力负荷测量方法 | |
CN100333688C (zh) | 评估脑皮层中神经元损伤程度的装置 | |
WO2024098649A1 (zh) | 一种基于生理唤醒识别的街道绿化品质检测方法 | |
CN113662560B (zh) | 发作间期痫样放电的检测方法、存储介质和设备 | |
CN112426160A (zh) | 一种心电信号类型识别方法及装置 | |
CN112773378B (zh) | 一种特征权重自适应学习的脑电情绪识别方法 | |
CN114970599A (zh) | 注意缺陷关联脑电信号的识别方法、识别装置、存储介质 | |
CN113576481B (zh) | 一种脑力负荷评估方法、装置、设备、介质 | |
CN111067513B (zh) | 一种特征权重自学习的睡眠质量检测关键脑区判定方法 | |
CN111568412A (zh) | 一种利用脑电信号重建视觉图像的方法及装置 | |
CN113988123A (zh) | 一种基于自加权增量rvfl网络的脑电疲劳预测方法 | |
CN113143275B (zh) | 一种样本与特征质量联合量化评估的脑电疲劳检测方法 | |
CN114237046A (zh) | 一种基于sift数据特征提取算法与bp神经网络模型的局部放电模式识别方法 | |
CN115713246A (zh) | 一种关于虚拟场景的多模态人机交互的性能评估方法 | |
CN112438741A (zh) | 一种基于脑电特征迁移学习的驾驶状态检测方法及*** | |
CN112244767A (zh) | 一种基于脑电的认知功能评估方法 | |
CN114305452B (zh) | 一种基于脑电和领域适配的跨任务认知负荷识别方法 | |
CN115736920A (zh) | 基于双模态融合的抑郁状态识别方法及*** |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |