CN105615879B - 基于多重分形消除趋势波动分析的脑电测谎方法 - Google Patents

基于多重分形消除趋势波动分析的脑电测谎方法 Download PDF

Info

Publication number
CN105615879B
CN105615879B CN201610207819.4A CN201610207819A CN105615879B CN 105615879 B CN105615879 B CN 105615879B CN 201610207819 A CN201610207819 A CN 201610207819A CN 105615879 B CN105615879 B CN 105615879B
Authority
CN
China
Prior art keywords
eeg signals
lie
fractal
testee
eeg
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
Application number
CN201610207819.4A
Other languages
English (en)
Other versions
CN105615879A (zh
Inventor
艾玲梅
陈慧君
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shaanxi Normal University
Original Assignee
Shaanxi Normal University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shaanxi Normal University filed Critical Shaanxi Normal University
Priority to CN201610207819.4A priority Critical patent/CN105615879B/zh
Publication of CN105615879A publication Critical patent/CN105615879A/zh
Application granted granted Critical
Publication of CN105615879B publication Critical patent/CN105615879B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316—Modalities, i.e. specific diagnostic methods
    • A61B5/369—Electroencephalography [EEG]
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316—Modalities, i.e. specific diagnostic methods
    • A61B5/369—Electroencephalography [EEG]
    • A61B5/377—Electroencephalography [EEG] using evoked responses
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/68—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient
    • A61B5/6801—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient specially adapted to be attached to or worn on the body surface
    • A61B5/6813—Specially adapted to be attached to a specific body part
    • A61B5/6814—Head
    • 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
    • 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/7221—Determining signal validity, reliability or quality
    • 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
    • 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/7271—Specific aspects of physiological measurement analysis
    • A61B5/7282—Event detection, e.g. detecting unique waveforms indicative of a medical condition
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
    • G16Z99/00—Subject matter not provided for in other main groups of this subclass
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
    • G06F2218/02—Preprocessing
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
    • G06F2218/08—Feature extraction
    • G06F2218/10—Feature extraction by analysing the shape of a waveform, e.g. extracting parameters relating to peaks
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
    • G06F2218/12—Classification; Matching

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Surgery (AREA)
  • General Health & Medical Sciences (AREA)
  • Biophysics (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Physics & Mathematics (AREA)
  • Animal Behavior & Ethology (AREA)
  • Pathology (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Psychiatry (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physiology (AREA)
  • Signal Processing (AREA)
  • Psychology (AREA)
  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)

Abstract

本发明公开了一种基于多重分形消除趋势波动分析的脑电测谎方法，该方法包括以下步骤：一、脑电测谎判断阈值设定：采用与处理器连接的参数设置单元对脑电信号提取装置的脑电测谎判断阈值进行设定；二、脑电信号采集：采用Pz电极对被测试者头部顶点的脑电信号进行实时提取，并通过脑电信号采集设备按照预先设定的采样频率f对Pz电极提取的脑电信号同步进行采集，且将所采集的脑电信号同步传送至处理器；三、脑电信号分析处理，过程如下：脑电信号接收与同步存储、脑电信号特征提取、说谎与否判断和说谎与否判断结果输出。本发明方法步骤简单、设计合理且实现方便、使用效果好，能简单、快速对说谎状态进行准确检测。

Description

基于多重分形消除趋势波动分析的脑电测谎方法
技术领域
本发明属于脑电信号处理技术领域，尤其是涉及一种基于多重分形消除趋势波动分析的脑电测谎方法。
背景技术
测谎技术作为司法审讯的一种辅助手段，对刑事案件事实的侦查和证明发挥着重要的作用，现已得到广泛应用。如何准确有效地识别谎言，对于司法机关工作人员来说就显得尤为重要，研究谎言的检测方法具有重大的现实意义。1912年，美国研制的多导测谎仪主要通过呼吸、脉搏、血压、皮肤电压等生理变化指标来判别谎言，由于上述参数易受主观因素的影响，使得检测结果难以保证准确无误，误报和漏报率高。1989年，Rosenfeld提出事件相关电位(ERP)的测谎改进方法，其原理是利用受试者接受与案情相关的刺激时产生的P300诱发电位的脑电波变化来判定谎情。
事件相关电位(ERP，event-related potential)是一种特殊的脑诱发电位，通过有意地赋予刺激以特殊的心理意义，利用多个或多样的刺激所引起的脑电位，它反映了认知过程中大脑的神经电生理的变化，也被称为认知电位，也就是指当人们对某课题进行认知加工时，从头颅表面记录到的脑电位，也称为ERP脑电位，对应的脑电波信号为ERP脑电波信号。
事件相关电位(ERP)测谎研究中主要涉及的诱发电位是P300和N400，其中P300诱发电位的脑电波信号为250ms～700ms的正向波，N400诱发电位的脑电波信号为潜伏期为200ms～500ms的负向波。P300诱发电位的脑电波信号的典型特征是潜伏期、波幅和波面积。由于P300诱发电位是大脑在“无意识”过程中产生的，不会出现“把无当有”的情况，较科学有效地帮助警方侦破案件。随后，众多研究者对ERP测谎技术提出多种改进，如幅值差(即BAD方法)、线性相关系数(即BCD方法)、频谱估计值(即功率谱估计)、小波系数、脑网络结构等，上述处理方法主要以脑电信号的时域、频域特征参数为评价标准展开研究，如BAD方法和BCD方法这两种经典的信号分析处理方法，都利用单通道的P300诱发电位的脑电波信号幅值信息区分诚实和欺骗性反应。赵敏等提出的多通道幅值和受试者的反应时间相结合的方法以及Dey S提出的小波特征提取同样能有效地反映出欺骗和诚实之间的显著差异，并与BAD方法和BCD方法对比得到了较好的结果。但上述处理方法虽然提高了单次响应下诱发脑电信号的信噪比，并且直观地反映出了诚实和欺骗行为之间的差异，但是没有考虑到脑电信号的非线性和混沌特性，造成部分信息损失。
发明内容
本发明所要解决的技术问题在于针对上述现有技术中的不足，提供一种基于多重分形消除趋势波动分析的脑电测谎方法，其方法步骤简单、设计合理且实现方便、使用效果好，能简单、快速对说谎状态进行准确检测。
为解决上述技术问题，本发明采用的技术方案是：一种基于多重分形消除趋势波动分析的脑电测谎方法，其特征在于：该方法包括以下步骤：
步骤一、脑电测谎判断阈值设定：采用与处理器连接的参数设置单元对脑电信号提取装置的脑电测谎判断阈值进行设定；
所述脑电信号检测装置包括Pz电极，所述Pz电极为按照国际标准“10/20”电极安放方法布设在被测试者头部顶点的脑电电极；
所述脑电信号提取装置的脑电测谎判断阈值包括所述Pz电极的脑电测谎判断阈值，所述Pz电极的脑电测谎判断阈值包括Hurst指数判断阈值HPz和多重分形谱宽度判断阈值WPz，其中HPz＝1.0142～1.0967，WPz＝0.894～1.11；
步骤二、脑电信号采集：测谎过程中，对被测试者赋予刺激前的时间段t1内和对被测试者赋予刺激后的时间段t2内，均采用所述Pz电极对被测试者头部顶点的脑电信号进行实时提取，并通过脑电信号采集设备按照预先设定的采样频率f对所述Pz电极提取的脑电信号同步进行采集，且将所采集的脑电信号同步传送至处理器；
所述Pz电极提取的脑电信号记作Pz脑电信号，所述Pz脑电信号为因赋予刺激诱发的诱发脑电信号且其诱发电位为事件相关电位，所述事件相关电位为P300；
本步骤中，采样时间为T且T＝t1+t2；其中，t1＝400ms～600ms；t2＝1200ms～1800ms；
步骤三、脑电信号分析处理，过程如下：
步骤301、脑电信号接收与同步存储：所述处理器将所接收到的脑电信号xPz(i)同步存储至存储器内；
所述存储器与处理器连接；
所述脑电信号xPz(i)包括N个采样周期内采集的Pz脑电信号，其中i为正整数且i＝1、2、3、…、N，f的单位为Hz，T的单位为ms；
步骤302、脑电信号特征提取：所述处理器调用多重分形消除趋势波动分析模块计算得出脑电信号xPz(i)的Hurst指数HPz’和多重分形谱宽度WPz’；
步骤303、说谎与否判断：根据步骤302中得出的Hurst指数HPz’和多重分形谱宽度WPz’，并结合步骤一中设定的Hurst指数判断阈值HPz和多重分形谱宽度判断阈值WPz，所述处理器调用差值比较模块对此时被测试者说谎与否进行判断：当HPz’＞HPz且WPz’＜WPz时，判断为此时被测试者没有说谎；否则，判断为此时被测试者说谎；
步骤304、说谎与否判断结果输出：所述处理器将步骤303中作出的说谎与否判断结果输出。
上述基于多重分形消除趋势波动分析的脑电测谎方法，其特征是：步骤302中进行脑电信号特征提取之前，还需对当前分析处理周期内所采集脑电信号xPz(i)进行预处理；
对当前分析处理周期内所采集脑电信号xPz(i)进行预处理时，先对脑电信号xPz(i)进行去除眼电伪迹处理，再对去除眼电伪迹处理后的脑电信号xPz(i)进行滤波处理。
上述基于多重分形消除趋势波动分析的脑电测谎方法，其特征是：步骤一中进行脑电测谎判断阈值设定时，过程如下：
步骤101、按照常规的ERP实验方法，对被测试者随机赋予三种刺激，三种刺激分别为探测刺激、无关刺激和靶刺激；三种刺激的测试次数均为多次，其中无关刺激出现的概率为70％，探测刺激和靶刺激出现的概率均为15％；
对被测试者赋予任一种刺激时，均按照步骤二中所述的方法进行脑电信号采集；
步骤102、将赋予多次探测刺激时脑电信号采集设备采集的多个脑电信号xPz(i)进行叠加平均，得到探测刺激下的脑电信号XPz(i)，再按照步骤301至步骤302中所述的方法计算得出脑电信号XPz(i)的Hurst指数HPzM和多重分形谱宽度WPzM；
同时，将赋予多次无关刺激时脑电信号采集设备采集的多个脑电信号xPz(i)进行叠加平均，得到无关刺激下的脑电信号X’Pz(i)，再按照步骤301至步骤302中所述的方法计算得出脑电信号X’Pz(i)的Hurst指数HPzm和多重分形谱宽度WPzm；
步骤103、根据公式和公式计算得出Hurst指数判断阈值HPz和多重分形谱宽度判断阈值WPz。
上述基于多重分形消除趋势波动分析的脑电测谎方法，其特征是：步骤二中所述的t1＝500ms，t2＝1500ms。
上述基于多重分形消除趋势波动分析的脑电测谎方法，其特征是：步骤一中所述脑电信号检测装置还包括Fz电极；所述Fz电极为按照国际标准“10/20”电极安放方法布设在被测试者头部中央点的脑电电极；
所述脑电信号提取装置的脑电测谎判断阈值还包括所述Fz电极的脑电测谎判断阈值，所述Fz电极的脑电测谎判断阈值包括Hurst指数判断阈值HFz和多重分形谱宽度判断阈值WFz，其中HFz＝0.9035～0.9985，WFz＝0.6972～0.8377；
步骤二中进行脑电信号采集时，对被测试者赋予刺激前的时间段t1内和对被测试者赋予刺激后的时间段t2内，均采用所述Fz电极对被测试者头部中央点的脑电信号进行实时提取，并通过脑电信号采集设备按照预先设定的采样频率f对所述Fz电极提取的脑电信号同步进行采集且将所采集的脑电信号同步传送至处理器；
所述Fz电极提取的脑电信号记作Fz脑电信号，所述Fz脑电信号为因赋予刺激诱发的诱发脑电信号且其诱发电位为事件相关电位，所述事件相关电位为P300；
步骤三中进行脑电信号分析处理，过程如下：
步骤301中进行脑电信号接收与同步存储时，所述处理器还需将所接收到的脑电信号xFz(i)同步存储至存储器内；
所述脑电信号xFz(i)包括N个采样周期内采集的Fz脑电信号；
步骤302中进行脑电信号特征提取时，所述处理器还需调用多重分形消除趋势波动分析模块计算得出脑电信号xFz(i)的Hurst指数HFz’和多重分形谱宽度WFz’；
步骤303中得出的说谎与否判断结果为根据脑电信号xFz(i)判断得出的说谎与否判断结果；
步骤303中完成说谎与否判断后，还需根据步骤302中得出的Hurst指数HFz’和多重分形谱宽度WFz’，并结合步骤一中设定的Hurst指数判断阈值HFz和多重分形谱宽度判断阈值WFz，所述处理器调用差值比较模块对此时被测试者说谎与否进行判断：当HFz’＞HFz且WFz’＜WFz时，判断为此时被测试者没有说谎；否则，判断为此时被测试者说谎；此时，所得出的说谎与否判断结果为根据脑电信号xFz(i)判断得出的说谎与否判断结果；
步骤304中进行说谎与否判断结果输出之前，还需判断根据脑电信号xPz(i)判断得出的说谎与否判断结果和根据脑电信号xFz(i)判断得出的说谎与否判断结果是否一致：当根据脑电信号xPz(i)判断得出的说谎与否判断结果和根据脑电信号xFz(i)判断得出的说谎与否判断结果一致时，进入步骤304；否则，所述处理器输出本次测谎无效。
上述基于多重分形消除趋势波动分析的脑电测谎方法，其特征是：步骤一中对所述Fz电极的脑电测谎判断阈值进行设定时，过程如下：
步骤A1、按照常规的ERP实验方法，对被测试者随机赋予三种刺激，三种刺激分别为探测刺激、无关刺激和靶刺激；三种刺激的测试次数均为多次，其中无关刺激出现的概率为70％，探测刺激和靶刺激出现的概率均为15％；
对被测试者赋予任一种刺激时，在被测试者赋予刺激前的时间段t1内和对被测试者赋予刺激后的时间段t2内，均采用所述Fz电极对被测试者头部中央点的脑电信号进行实时提取，并通过脑电信号采集设备按照预先设定的采样频率f对所述Fz电极提取的脑电信号同步进行采集且将所采集的脑电信号同步传送至处理器；
步骤A2、将赋予多次探测刺激时脑电信号采集设备采集的多个脑电信号xFz(i)进行叠加平均，得到探测刺激下的脑电信号XFz(i)，再按照步骤301至步骤302中所述的方法计算得出脑电信号XFz(i)的Hurst指数HFzM和多重分形谱宽度WFzM；
同时，将赋予多次无关刺激时脑电信号采集设备采集的多个脑电信号xFz(i)进行叠加平均，得到无关刺激下的脑电信号X’Fz(i)，再按照步骤301至步骤302中所述的方法计算得出脑电信号X’Fz(i)的Hurst指数HFzm和多重分形谱宽度WFzm；
步骤A3、根据公式和公式计算得出Hurst指数判断阈值HFz和多重分形谱宽度判断阈值WFz。
上述基于多重分形消除趋势波动分析的脑电测谎方法，其特征是：步骤一中所述脑电信号检测装置还包括辅助测试用脑电电极；所述辅助测试用脑电电极的数量至少为两个且其数量不大于3个；所述辅助测试用脑电电极为按照国际标准“10/20”电极安放方法布设在FC1、FC2或Fz点的脑电电极；
步骤一中所述脑电信号提取装置的脑电测谎判断阈值还包括各辅助测试用脑电电极的脑电测谎判断阈值，所述辅助测试用脑电电极的脑电测谎判断阈值包括Hurst指数判断阈值HAs和多重分形谱宽度判断阈值WAs；
布设FC1点的脑电电极的Hurst指数判断阈值HAs＝0.9587～1.0242且其多重分形谱宽度判断阈值WAs＝0.7247～0.8522；
布设FC2点的脑电电极的Hurst指数判断阈值HAs＝0.9497～1.0392且其多重分形谱宽度判断阈值WAs＝0.6947～0.8082；
布设Fz点的脑电电极的Hurst指数判断阈值HAs＝0.9035～0.9985且其多重分形谱宽度判断阈值WAs＝0.6972～0.8377；
步骤二中进行脑电信号采集时，对被测试者赋予刺激前的时间段t1内和对被测试者赋予刺激后的时间段t2内，均采用各辅助测试用脑电电极对被测试者头部对应部位的脑电信号进行实时提取，并通过脑电信号采集设备按照预先设定的采样频率f对各辅助测试用脑电电极提取的脑电信号同步进行采集且将所采集的脑电信号同步传送至处理器；
所述辅助测试用脑电电极提取的脑电信号记作As脑电信号，所述As脑电信号为因赋予刺激诱发的诱发脑电信号且其诱发电位为事件相关电位，所述事件相关电位为P300；
步骤三中进行脑电信号分析处理，过程如下：
步骤301中进行脑电信号接收与同步存储时，所述处理器还需将所接收到的各脑电信号xAs(i)同步存储至存储器内；
每个所述脑电信号xAs(i)包括N个采样周期内同一个所述辅助测试用脑电电极采集的As脑电信号；
步骤302中进行脑电信号特征提取时，所述处理器还需调用基于多重分形消除趋势波动分析模块分别计算得出各脑电信号xAs(i)的Hurst指数HAs’和多重分形谱宽度WAs’；
步骤303中得出的说谎与否判断结果为根据脑电信号xPz(i)判断得出的说谎与否判断结果；
步骤303中完成说谎与否判断后，还需根据各脑电信号xAs(i)分别进行说谎与否判断；
其中，对任一个脑电信号xAs(i)进行说谎与否判断时，根据步骤302中得出的该脑电信号xAs(i)的Hurst指数HAs’和多重分形谱宽度WAs’，所述处理器调用差值比较模块对此时被测试者说谎与否进行判断：当该脑电信号xAs(i)的Hurst指数HAs’大于步骤一中设定的与该脑电信号xAs(i)对应的辅助测试用脑电电极的Hurst指数判断阈值HAs且该脑电信号xAs(i)的多重分形谱宽度WAs’小于步骤一中设定的与该脑电信号xAs(i)对应的辅助测试用脑电电极的多重分形谱宽度判断阈值WAs时，判断为此时被测试者没有说谎；否则，判断为此时被测试者说谎；此时，所得出的说谎与否判断结果为根据脑电信号xAs(i)判断得出的说谎与否判断结果；
步骤304中进行说谎与否判断结果输出之前，还需判断根据脑电信号xPz(i)判断得出的说谎与否判断结果和根据各脑电信号xAs(i)判断得出的说谎与否判断结果是否一致：当根据脑电信号xPz(i)判断得出的说谎与否判断结果和根据各脑电信号xAs(i)判断得出的说谎与否判断结果均一致时，进入步骤304；否则，所述处理器输出本次测谎无效。
上述基于多重分形消除趋势波动分析的脑电测谎方法，其特征是：步骤一中所述脑电信号提取装置为32导联脑电电极或64导联脑电电极。
本发明与现有技术相比具有以下优点：
1、所采用的脑电测谎系统结构简单、设计合理且接线方便，投入成本较低。
2、所采用的特征提取方法设计合理，所提取的特征参数包括Hurst指数和多重分形谱宽度，其中Hurst指数为长程幂函数相关性的一个定量化指标，多重分形谱宽度的大小能反映时间序列信号分布的复杂程度和不均匀程度，将Hurst指数和多重分形谱宽度相结合能全面、综合反映脑电信号的长程幂律性和复杂性，能真实且全面反映所分析脑电信号的特征，能克服传统特征提取方法中仅重点考虑脑电信号的线性特征造成部分关键信息损失的缺陷，从而确保测谎检测结果的准确性。
3、所采用的脑电测谎方法步骤简单、设计合理且实现方便，主要包括脑电测谎判断阈值设定、脑电信号采集和脑电信号分析处理三个步骤，由数据处理设备自动完成，同步性好，在脑电信号采集的同时能及时、快速得出测谎结果。
4、脑电测谎判断阈值设定方法步骤简单、实现方便且所设定脑电测谎判断阈值准确性好，通过常规的ERP实验方法进行有效利用，能直接得出与被测试者相适应的脑电测谎判断阈值，该脑电测谎判断阈值具有适应性，根据各被测试者的脑电信号分析处理结果进行确定，因而能准确反映个体差异，使得测谎结果更加准确，针对性更强。
5、脑电测谎方法灵活，既能根据一个单一脑电电极(具体是Pz电极)所提取脑电信号的分析处理结果直接得出测谎结果，数据处理过程简单，所选用的Pz电极合理，该Pz电极对应的分析处理结果能准确、明显反映被测试者的说谎状态；并且，本发明在Pz电极的基础上，也能结合一个或多个辅助测试用脑电电极所提取脑电信号的分析处理结果进行综合判断后得出测谎结果，从而避免单一电极分析处理结果的片面性和不准确性，进一步提高测谎结果的准确性。
6、使用效果好且实用价值高，推广应用前景广泛，本发明考虑到脑电信号的非线性和混沌特性，能有效保证测谎结果的准确性。
综上所述，本发明方法步骤简单、设计合理且实现方便、使用效果好，能简单、快速对说谎状态进行准确检测。。
下面通过附图和实施例，对本发明的技术方案做进一步的详细描述。
附图说明
图1为本发明所采用脑电测谎系统的电路原理框图。
图2为采用本发明进行脑电测谎时的方法流程框图。
图3为本发明32导联脑电电极的布设位置示意图。
图4为本发明64导联脑电电极的布设位置示意图。
图5为本发明ERP实验时探测刺激和无关刺激下所提取Hurst指数的对比示意图。
图6为本发明ERP实验时探测刺激和无关刺激下所提取多重分形谱宽度的对比示意图。
附图标记说明:
1—处理器； 2—参数设置单元； 3—脑电信号采集设备；
4—存储器； 5—报警单元； 6—显示器；
7—计时电路； 8—脑电信号提取装置； 9—脑电信号放大器。
具体实施方式
如图2所示的一种基于多重分形消除趋势波动分析的脑电测谎方法，包括以下步骤：
步骤一、脑电测谎判断阈值设定：采用与处理器1连接的参数设置单元2对脑电信号提取装置的脑电测谎判断阈值进行设定；
所述脑电信号检测装置包括Pz电极，所述Pz电极为按照国际标准“10/20”电极安放方法布设在被测试者头部顶点的脑电电极；
所述脑电信号提取装置的脑电测谎判断阈值包括所述Pz电极的脑电测谎判断阈值，所述Pz电极的脑电测谎判断阈值包括Hurst指数判断阈值HPz和多重分形谱宽度判断阈值WPz，其中HPz和WPz均为常数，HPz＝1.0142～1.0967，WPz＝0.894～1.11；
步骤二、脑电信号采集：测谎过程中，对被测试者赋予刺激前的时间段t1内和对被测试者赋予刺激后的时间段t2内，均采用所述Pz电极对被测试者头部顶点的脑电信号进行实时提取，并通过脑电信号采集设备3按照预先设定的采样频率f对所述Pz电极提取的脑电信号同步进行采集，且将所采集的脑电信号同步传送至处理器1；
所述Pz电极提取的脑电信号记作Pz脑电信号，所述Pz脑电信号为因赋予刺激诱发的诱发脑电信号且其诱发电位为事件相关电位，所述事件相关电位为P300；
本步骤中，采样时间为T且T＝t1+t2；其中，t1＝400ms～600ms；t2＝1200ms～1800ms；
步骤三、脑电信号分析处理，过程如下：
步骤301、脑电信号接收与同步存储：所述处理器1将所接收到的脑电信号xPz(i)同步存储至存储器4内；
所述存储器4与处理器1连接；
所述脑电信号xPz(i)包括N个采样周期内采集的Pz脑电信号，其中i为正整数且i＝1、2、3、…、N，f的单位为Hz，T的单位为ms；
步骤302、脑电信号特征提取：所述处理器1调用多重分形消除趋势波动分析模块计算得出脑电信号xPz(i)的Hurst指数HPz’和多重分形谱宽度WPz’；
步骤303、说谎与否判断：根据步骤302中得出的Hurst指数HPz’和多重分形谱宽度WPz’，并结合步骤一中设定的Hurst指数判断阈值HPz和多重分形谱宽度判断阈值WPz，所述处理器1调用差值比较模块对此时被测试者说谎与否进行判断：当HPz’＞HPz且WPz’＜WPz时，判断为此时被测试者没有说谎；否则，判断为此时被测试者说谎；
步骤304、说谎与否判断结果输出：所述处理器1将步骤303中作出的说谎与否判断结果输出。
根据本领域公知常识，事件相关电位(ERP，event-related potential)是一种特殊的脑诱发电位，通过有意地赋予刺激以特殊的心理意义，利用多个或多样的刺激所引起的脑电位，它反映了认知过程中大脑的神经电生理的变化，也被称为认知电位，也就是指当人们对某课题进行认知加工时，从头颅表面记录到的脑电位，也称为ERP脑电位，对应的脑电信号为ERP脑电信号。其中，赋予刺激指的是对被测试者赋予一次刺激。此处，对被测试者赋予刺激指的是对被测试者进行一次测谎提问。因而，步骤二中所述的对被测试者赋予刺激指的是引起事件相关电位的刺激，具体是对被测试者进行一次提问。
其中，步骤302中所述的多重分形消除趋势波动分析模块为多重分形消除趋势波动分析法(multifractal detrended fluctuation analysis，MF-DFA)模块。
多重分形消除趋势波动分析法(multifractal detrended fluctuationanalysis，MF-DFA)是2002年Kantelhardt等[20]在消除趋势波动分析法(DFA)的基础上提出的，详见《Multifractal detrended fluctuation analysis of non-stationary timeseries》，作者Kantelhardt J W,Zschiegner S A,Koscielny-Bunde E,et al.，2002年发表于Physica A-statistical Mechanics&Its Applications。
本实施例中，所述处理器1与报警单元5连接，所述报警单元5由处理器1进行控制。并且，所述处理器1分别与显示器6和计时电路7连接。
同时，所述脑电信号提取装置8与脑电信号采集设备3之间接有脑电信号放大器9，通过脑电信号放大器9对脑电信号提取装置8所提取的脑电信号进行放大处理。
本实施例中，所述报警单元5为声光报警装置。如图1所示，所述处理器1、参数设置单元2、存储器4、脑电信号放大器9、脑电信号采集设备3、报警单元5、显示器6和计时电路7组成脑电测谎系统。
本实施例中，步骤302中进行脑电信号特征提取之前，还需对当前分析处理周期内所采集脑电信号xPz(i)进行预处理；
对当前分析处理周期内所采集脑电信号xPz(i)进行预处理时，先对脑电信号xPz(i)进行去除眼电伪迹处理，再对去除眼电伪迹处理后的脑电信号xPz(i)进行滤波处理。
本实施例中，步骤304中进行说谎与否判断结果输出时，通过所述报警单元5同步进行报警提示。
本实施例中，步骤一中进行脑电测谎判断阈值设定时，过程如下：
步骤101、按照常规的ERP实验方法，对被测试者随机赋予三种刺激，三种刺激分别为探测刺激、无关刺激和靶刺激；三种刺激的测试次数均为多次，其中无关刺激出现的概率为70％，探测刺激和靶刺激出现的概率均为15％；
对被测试者赋予任一种刺激时，均按照步骤二中所述的方法进行脑电信号采集；
步骤102、将赋予多次探测刺激时脑电信号采集设备3采集的多个脑电信号xPz(i)进行叠加平均，得到探测刺激下的脑电信号XPz(i)，再按照步骤301至步骤302中所述的方法计算得出脑电信号XPz(i)的Hurst指数HPzM和多重分形谱宽度WPzM；
同时，将赋予多次无关刺激时脑电信号采集设备3采集的多个脑电信号xPz(i)进行叠加平均，得到无关刺激下的脑电信号X’Pz(i)，再按照步骤301至步骤302中所述的方法计算得出脑电信号X’Pz(i)的Hurst指数HPzm和多重分形谱宽度WPzm，详见图5和图6；
步骤103、根据公式和公式计算得出Hurst指数判断阈值HPz和多重分形谱宽度判断阈值WPz。
步骤102中所述的“叠加平均”指的是“求和后取平均值”。
其中，ERP实验方法为事件相关电位实验方法，并且对被测试者随机赋予三种刺激，按照常规的新异刺激Oddball范式进行刺激。其中，Oddball范式是诱发P300等与刺激概率有关的ERP成分时常用的经典实验范式。
研究显示，自传体信息更容易诱发出相对明显的事件相关电位P300波即诱发电位为P300的脑电信号。此时，ERP实验中分别采用被测试者的家乡和无关城市作为测试信息，以“你知道XX吗？”的形式出现，刺激序列的编排和出现依据典型的新异刺激Oddball范式，包括三种刺激信息。其中被试的家乡作为探测刺激，即要求故意隐瞒并否认的项；另选5个对被试无特殊意义的城市名称作为无关刺激，问题“你知道北京吗？”作为靶刺激，要求被测试者对无关刺激和靶刺激做诚实回答。
上述七个刺激随机出现160次，其中无关刺激112次即出现概率为70％，探测刺激和靶刺激均为24次即出现概率为15％，每个刺激显示1.5s，相邻两个刺激的时间间隔为2s。被测试者按照实验要求，对随机出现的刺激问题判断后进行回答。每一次赋予刺激过程中，均按照步骤二中所述的方法进行脑电信号采集。
本实施例中，步骤二中所述的t1＝500ms，t2＝1500ms。
实际使用时，可根据具体需要，对t1和t2的取值大小进行相应调整。
其中，步骤二中所述的“在被测试者赋予刺激前的时间段t1内和对被测试者赋予刺激后的时间段t2内”，指的是从赋予刺激前500ms至赋予刺激后1500ms这2000ms时间段内进行脑电信号采集。
本实施例中，步骤二中进行脑电信号采集时，采样频率f为1kHz。
本实施例中，步骤一中所述脑电信号检测装置还包括Fz电极；所述Fz电极为按照国际标准“10/20”电极安放方法布设在被测试者头部中央点的脑电电极；
所述脑电信号提取装置的脑电测谎判断阈值还包括所述Fz电极的脑电测谎判断阈值，所述Fz电极的脑电测谎判断阈值包括Hurst指数判断阈值HFz和多重分形谱宽度判断阈值WFz，其中HFz和WFz均为常数，HFz＝0.9035～0.9985，WFz＝0.6972～0.8377；
步骤二中进行脑电信号采集时，对被测试者赋予刺激前的时间段t1内和对被测试者赋予刺激后的时间段t2内，均采用所述Fz电极对被测试者头部中央点的脑电信号进行实时提取，并通过脑电信号采集设备3按照预先设定的采样频率f对所述Fz电极提取的脑电信号同步进行采集且将所采集的脑电信号同步传送至处理器1；
所述Fz电极提取的脑电信号记作Fz脑电信号，所述Fz脑电信号为因赋予刺激诱发的诱发脑电信号且其诱发电位为事件相关电位，所述事件相关电位为P300；
步骤三中进行脑电信号分析处理，过程如下：
步骤301中进行脑电信号接收与同步存储时，所述处理器1还需将所接收到的脑电信号xFz(i)同步存储至存储器4内；
所述脑电信号xFz(i)包括N个采样周期内采集的Fz脑电信号；
步骤302中进行脑电信号特征提取时，所述处理器1还需调用多重分形消除趋势波动分析模块计算得出脑电信号xFz(i)的Hurst指数HFz’和多重分形谱宽度WFz’；
步骤303中得出的说谎与否判断结果为根据脑电信号xFz(i)判断得出的说谎与否判断结果；
步骤303中完成说谎与否判断后，还需根据步骤302中得出的Hurst指数HFz’和多重分形谱宽度WFz’，并结合步骤一中设定的Hurst指数判断阈值HFz和多重分形谱宽度判断阈值WFz，所述处理器1调用差值比较模块对此时被测试者说谎与否进行判断：当HFz’＞HFz且WFz’＜WFz时，判断为此时被测试者没有说谎；否则，判断为此时被测试者说谎；此时，所得出的说谎与否判断结果为根据脑电信号xFz(i)判断得出的说谎与否判断结果；
步骤304中进行说谎与否判断结果输出之前，还需判断根据脑电信号xPz(i)判断得出的说谎与否判断结果和根据脑电信号xFz(i)判断得出的说谎与否判断结果是否一致：当根据脑电信号xPz(i)判断得出的说谎与否判断结果和根据脑电信号xFz(i)判断得出的说谎与否判断结果一致时，进入步骤304；否则，所述处理器1输出本次测谎无效。
本实施例中，步骤一中对所述Fz电极的脑电测谎判断阈值进行设定时，过程如下：
步骤A1、按照常规的ERP实验方法，对被测试者随机赋予三种刺激，三种刺激分别为探测刺激、无关刺激和靶刺激；三种刺激的测试次数均为多次，其中无关刺激出现的概率为70％，探测刺激和靶刺激出现的概率均为15％；
对被测试者赋予任一种刺激时，在被测试者赋予刺激前的时间段t1内和对被测试者赋予刺激后的时间段t2内，均采用所述Fz电极对被测试者头部中央点的脑电信号进行实时提取，并通过脑电信号采集设备3按照预先设定的采样频率f对所述Fz电极提取的脑电信号同步进行采集且将所采集的脑电信号同步传送至处理器1；
步骤A2、将赋予多次探测刺激时脑电信号采集设备3采集的多个脑电信号xFz(i)进行叠加平均，得到探测刺激下的脑电信号XFz(i)，再按照步骤301至步骤302中所述的方法计算得出脑电信号XFz(i)的Hurst指数HFzM和多重分形谱宽度WFzM；
同时，将赋予多次无关刺激时脑电信号采集设备3采集的多个脑电信号xFz(i)进行叠加平均，得到无关刺激下的脑电信号X’Fz(i)，再按照步骤301至步骤302中所述的方法计算得出脑电信号X’Fz(i)的Hurst指数HFzm和多重分形谱宽度WFzm，详见图5和图6；
步骤A3、根据公式和公式计算得出Hurst指数判断阈值HFz和多重分形谱宽度判断阈值WFz。
步骤A2中所述的“叠加平均”指的是“求和后取平均值”。
其中，所述Fz电极的脑电测谎判断阈值的设定方法与所述Pz电极的脑电测谎判断阈值的设定方法相同。
实际使用时，步骤一中所述脑电信号检测装置还包括辅助测试用脑电电极；所述辅助测试用脑电电极的数量至少为两个且其数量不大于3个；所述辅助测试用脑电电极为按照国际标准“10/20”电极安放方法布设在FC1、FC2或Fz点的脑电电极；
步骤一中所述脑电信号提取装置的脑电测谎判断阈值还包括各辅助测试用脑电电极的脑电测谎判断阈值，所述辅助测试用脑电电极的脑电测谎判断阈值包括Hurst指数判断阈值HAs和多重分形谱宽度判断阈值WAs；各辅助测试用脑电电极的Hurst指数判断阈值HAs和多重分形谱宽度判断阈值WAs均为常数；
布设FC1点的脑电电极的Hurst指数判断阈值HAs＝0.9587～1.0242且其多重分形谱宽度判断阈值WAs＝0.7247～0.8522；
布设FC2点的脑电电极的Hurst指数判断阈值HAs＝0.9497～1.0392且其多重分形谱宽度判断阈值WAs＝0.6947～0.8082；
布设Fz点的脑电电极的Hurst指数判断阈值HAs＝0.9035～0.9985且其多重分形谱宽度判断阈值WAs＝0.6972～0.8377；
步骤二中进行脑电信号采集时，对被测试者赋予刺激前的时间段t1内和对被测试者赋予刺激后的时间段t2内，均采用各辅助测试用脑电电极对被测试者头部对应部位的脑电信号进行实时提取，并通过脑电信号采集设备3按照预先设定的采样频率f对各辅助测试用脑电电极提取的脑电信号同步进行采集且将所采集的脑电信号同步传送至处理器1；
所述辅助测试用脑电电极提取的脑电信号记作As脑电信号，所述As脑电信号为因赋予刺激诱发的诱发脑电信号且其诱发电位为事件相关电位，所述事件相关电位为P300；
步骤三中进行脑电信号分析处理，过程如下：
步骤301中进行脑电信号接收与同步存储时，所述处理器1还需将所接收到的各脑电信号xAs(i)同步存储至存储器4内；
每个所述脑电信号xAs(i)包括N个采样周期内同一个所述辅助测试用脑电电极采集的As脑电信号；
步骤302中进行脑电信号特征提取时，所述处理器1还需调用基于多重分形消除趋势波动分析模块分别计算得出各脑电信号xAs(i)的Hurst指数HAs’和多重分形谱宽度WAs’；
步骤303中得出的说谎与否判断结果为根据脑电信号xPz(i)判断得出的说谎与否判断结果；
步骤303中完成说谎与否判断后，还需根据各脑电信号xAs(i)分别进行说谎与否判断；
其中，对任一个脑电信号xAs(i)进行说谎与否判断时，根据步骤302中得出的该脑电信号xAs(i)的Hurst指数HAs’和多重分形谱宽度WAs’，所述处理器1调用差值比较模块对此时被测试者说谎与否进行判断：当该脑电信号xAs(i)的Hurst指数HAs’大于步骤一中设定的与该脑电信号xAs(i)对应的辅助测试用脑电电极的Hurst指数判断阈值HAs且该脑电信号xAs(i)的多重分形谱宽度WAs’小于步骤一中设定的与该脑电信号xAs(i)对应的辅助测试用脑电电极的多重分形谱宽度判断阈值WAs时，判断为此时被测试者没有说谎；否则，判断为此时被测试者说谎；此时，所得出的说谎与否判断结果为根据脑电信号xAs(i)判断得出的说谎与否判断结果；
步骤304中进行说谎与否判断结果输出之前，还需判断根据脑电信号xPz(i)判断得出的说谎与否判断结果和根据各脑电信号xAs(i)判断得出的说谎与否判断结果是否一致：当根据脑电信号xPz(i)判断得出的说谎与否判断结果和根据各脑电信号xAs(i)判断得出的说谎与否判断结果均一致时，进入步骤304；否则，所述处理器1输出本次测谎无效。
实际使用时，通过将根据辅助测试用脑电电极提取脑电信号的说谎与否判断结果与根据脑电信号xPz(i)判断得出的说谎与否判断结果进行一致性判断，弥补仅根据脑电信号xPz(i)判断得出的说谎与否判断结果的片面性，确保说谎与否判断结果的准确性。
并且，所述辅助测试用脑电电极的脑电测谎判断阈值的设定方法与所述Pz电极的脑电测谎判断阈值的设定方法相同。
对任一个所述辅助测试用脑电电极的脑电测谎判断阈值进行设定时，过程如下：
步骤B1、按照常规的ERP实验方法，对被测试者随机赋予三种刺激，三种刺激分别为探测刺激、无关刺激和靶刺激；三种刺激的测试次数均为多次，其中无关刺激出现的概率为70％，探测刺激和靶刺激出现的概率均为15％；
对被测试者赋予任一种刺激时，在被测试者赋予刺激前的时间段t1内和对被测试者赋予刺激后的时间段t2内，均采用所述As电极对被测试者头部中央点的脑电信号进行实时提取，并通过脑电信号采集设备3按照预先设定的采样频率f对所述As电极提取的脑电信号同步进行采集且将所采集的脑电信号同步传送至处理器1；
步骤B2、将赋予多次探测刺激时脑电信号采集设备3采集的多个脑电信号xAs(i)进行叠加平均，得到探测刺激下的脑电信号XAs(i)，再按照步骤301至步骤302中所述的方法计算得出脑电信号XAs(i)的Hurst指数HAsM和多重分形谱宽度WAsM；
同时，将赋予多次无关刺激时脑电信号采集设备3采集的多个脑电信号xAs(i)进行叠加平均，得到无关刺激下的脑电信号X’As(i)，再按照步骤301至步骤302中所述的方法计算得出脑电信号X’As(i)的Hurst指数HAsm和多重分形谱宽度WAsm，详见图5和图6；
步骤B3、根据公式和公式计算得出Hurst指数判断阈值HAs和多重分形谱宽度判断阈值WAs。
步骤B2中所述的“叠加平均”指的是“求和后取平均值”。
其中，所述Fz电极为按照国际标准“10/20”电极安放方法布设在被测试者头部额中点的脑电电极；
本实施例中，步骤一中所述脑电信号提取装置8为32导联脑电电极，32导联脑电电极中各脑电电极的布设位置详见图3。
所述32导联脑电电极中的32个脑电电极分别布设在被测试者头部的FP1、FP2、F7、F8、F3、F4、Fz、T7、T8、C3、C4、Cz、P3、P4、Pz、P7、P8、O1、Oz、O2、FC1、FC2、CP1、CP2、FC5、FC6、CP5、CP6、M1、M2、HEO和VEO点。其中，各导联的阻抗均小于5kΩ。
其中，64导联脑电电极中各脑电电极的布设位置详见图4，64导联脑电电极中的64个脑电电极分别布设在被测试者头部的FP1、FPz、FP2、AF3、AF4、F7、F5、F3、F1、Fz、F2、F4、F6、F8、FT7、FC5、FC3、FC1、FCz、FC2、FC4、FC6、FT8、T7、C5、C3、C1、Cz、C2、C4、C6、T8、TP7、CP5、CP3、CP1、CPz、CP2、CP4、CP6、TP8、P7、P5、P3、P1、Pz、P2、P4、P6、P8、PO7、PO5、PO3、POz、PO4、PO6、PO8、CB1、O1、Oz、O2、CB2、HEO(即“HEOG”)和VEO(即“VEOG”)点。
图3和图4中，“VEOG”表示垂直眼电，“HEOG”表示水平眼电。M1和M2点为双侧乳突位置。
实际使用时，所述脑电信号提取装置8也可以仅采用Pz脑电电极，投入成本低，经济实用。
以上所述，仅是本发明的较佳实施例，并非对本发明作任何限制，凡是根据本发明技术实质对以上实施例所作的任何简单修改、变更以及等效结构变化，均仍属于本发明技术方案的保护范围内。

Claims (8)

1.一种基于多重分形消除趋势波动分析的脑电测谎方法，其特征在于：该方法包括以下步骤：
步骤一、脑电测谎判断阈值设定：采用与处理器(1)连接的参数设置单元(2)对脑电信号提取装置的脑电测谎判断阈值进行设定；
所述脑电信号提取装置包括Pz电极，所述Pz电极为按照国际标准“10/20”电极安放方法布设在被测试者头部顶点的脑电电极；
所述脑电信号提取装置的脑电测谎判断阈值包括所述Pz电极的脑电测谎判断阈值，所述Pz电极的脑电测谎判断阈值包括Hurst指数判断阈值HPz和多重分形谱宽度判断阈值WPz，其中HPz和WPz均为常数，HPz＝1.0142～1.0967，WPz＝0.894～1.11；
步骤二、脑电信号采集：测谎过程中，对被测试者赋予刺激前的时间段t1内和对被测试者赋予刺激后的时间段t2内，均采用所述Pz电极对被测试者头部顶点的脑电信号进行实时提取，并通过脑电信号采集设备(3)按照预先设定的采样频率f对所述Pz电极提取的脑电信号同步进行采集，且将所采集的脑电信号同步传送至处理器(1)；
所述Pz电极提取的脑电信号记作Pz脑电信号，所述Pz脑电信号为因赋予刺激诱发的诱发脑电信号且其诱发电位为事件相关电位，所述事件相关电位为P300；
本步骤中，采样时间为T且T＝t1+t2；其中，t1＝400ms～600ms；t2＝1200ms～1800ms；
步骤三、脑电信号分析处理，过程如下：
步骤301、脑电信号接收与同步存储：所述处理器(1)将所接收到的脑电信号xPz(i)同步存储至存储器(4)内；
所述存储器(4)与处理器(1)连接；
所述脑电信号xPz(i)包括N个采样周期内采集的Pz脑电信号，其中i为正整数且i＝1、2、3、…、N，f的单位为Hz，T的单位为ms；
步骤302、脑电信号特征提取：所述处理器(1)调用多重分形消除趋势波动分析模块计算得出脑电信号xPz(i)的Hurst指数HPz’和多重分形谱宽度WPz’；
步骤303、说谎与否判断：根据步骤302中得出的Hurst指数HPz’和多重分形谱宽度WPz’，并结合步骤一中设定的Hurst指数判断阈值HPz和多重分形谱宽度判断阈值WPz，所述处理器(1)调用差值比较模块对此时被测试者说谎与否进行判断：当HPz’＞HPz且WPz’＜WPz时，判断为此时被测试者没有说谎；否则，判断为此时被测试者说谎；
步骤304、说谎与否判断结果输出：所述处理器(1)将步骤303中作出的说谎与否判断结果输出。
2.按照权利要求1所述的基于多重分形消除趋势波动分析的脑电测谎方法，其特征在于：步骤302中进行脑电信号特征提取之前，还需对脑电信号xPz(i)进行预处理；
对脑电信号xPz(i)进行预处理时，先对脑电信号xPz(i)进行去除眼电伪迹处理，再对去除眼电伪迹处理后的脑电信号xPz(i)进行滤波处理。
3.按照权利要求1或2所述的基于多重分形消除趋势波动分析的脑电测谎方法，其特征在于：步骤一中进行脑电测谎判断阈值设定时，过程如下：
步骤101、按照常规的ERP实验方法，对被测试者随机赋予三种刺激，三种刺激分别为探测刺激、无关刺激和靶刺激；三种刺激的测试次数均为多次，其中无关刺激出现的概率为70％，探测刺激和靶刺激出现的概率均为15％；
对被测试者赋予任一种刺激时，均按照步骤二中所述的方法进行脑电信号采集；
步骤102、将赋予多次探测刺激时脑电信号采集设备(3)采集的多个脑电信号xPz(i)进行叠加平均，得到探测刺激下的脑电信号XPz(i)，再按照步骤301至步骤302中所述的方法计算得出脑电信号XPz(i)的Hurst指数HPzM和多重分形谱宽度WPzM；
同时，将赋予多次无关刺激时脑电信号采集设备(3)采集的多个脑电信号xPz(i)进行叠加平均，得到无关刺激下的脑电信号X’Pz(i)，再按照步骤301至步骤302中所述的方法计算得出脑电信号X’Pz(i)的Hurst指数HPzm和多重分形谱宽度WPzm；
步骤103、根据公式和公式计算得出Hurst指数判断阈值HPz和多重分形谱宽度判断阈值WPz。
4.按照权利要求1或2所述的基于多重分形消除趋势波动分析的脑电测谎方法，其特征在于：步骤二中所述的t1＝500ms，t2＝1500ms。
5.按照权利要求1或2所述的基于多重分形消除趋势波动分析的脑电测谎方法，其特征在于：步骤一中所述脑电信号提取装置还包括Fz电极；所述Fz电极为按照国际标准“10/20”电极安放方法布设在被测试者头部中央点的脑电电极；
所述脑电信号提取装置的脑电测谎判断阈值还包括所述Fz电极的脑电测谎判断阈值，所述Fz电极的脑电测谎判断阈值包括Hurst指数判断阈值HFz和多重分形谱宽度判断阈值WFz，其中HFz＝0.9035～0.9985，WFz＝0.6972～0.8377；
步骤二中进行脑电信号采集时，对被测试者赋予刺激前的时间段t1内和对被测试者赋予刺激后的时间段t2内，均采用所述Fz电极对被测试者头部中央点的脑电信号进行实时提取，并通过脑电信号采集设备(3)按照预先设定的采样频率f对所述Fz电极提取的脑电信号同步进行采集且将所采集的脑电信号同步传送至处理器(1)；
所述Fz电极提取的脑电信号记作Fz脑电信号，所述Fz脑电信号为因赋予刺激诱发的诱发脑电信号且其诱发电位为事件相关电位，所述事件相关电位为P300；
步骤三中进行脑电信号分析处理，过程如下：
步骤301中进行脑电信号接收与同步存储时，所述处理器(1)还需将所接收到的脑电信号xFz(i)同步存储至存储器(4)内；
所述脑电信号xFz(i)包括N个采样周期内采集的Fz脑电信号；
步骤302中进行脑电信号特征提取时，所述处理器(1)还需调用多重分形消除趋势波动分析模块计算得出脑电信号xFz(i)的Hurst指数HFz’和多重分形谱宽度WFz’；
步骤303中得出的说谎与否判断结果为根据脑电信号xFz(i)判断得出的说谎与否判断结果；
步骤303中完成根据脑电信号xPz(i)判断得出的说谎与否判断后，还需根据步骤302中得出的Hurst指数HFz’和多重分形谱宽度WFz’，并结合步骤一中设定的Hurst指数判断阈值HFz和多重分形谱宽度判断阈值WFz，所述处理器(1)调用差值比较模块对此时被测试者说谎与否进行判断：当HFz’＞HFz且WFz’＜WFz时，判断为此时被测试者没有说谎；否则，判断为此时被测试者说谎；此时，所得出的说谎与否判断结果为根据脑电信号xFz(i)判断得出的说谎与否判断结果；
步骤304中进行说谎与否判断结果输出之前，还需判断根据脑电信号xPz(i)判断得出的说谎与否判断结果和根据脑电信号xFz(i)判断得出的说谎与否判断结果是否一致：当根据脑电信号xPz(i)判断得出的说谎与否判断结果和根据脑电信号xFz(i)判断得出的说谎与否判断结果一致时，进入步骤304；否则，所述处理器(1)输出本次测谎无效。
6.按照权利要求5所述的基于多重分形消除趋势波动分析的脑电测谎方法，其特征在于：步骤一中对所述Fz电极的脑电测谎判断阈值进行设定时，过程如下：
步骤A1、按照常规的ERP实验方法，对被测试者随机赋予三种刺激，三种刺激分别为探测刺激、无关刺激和靶刺激；三种刺激的测试次数均为多次，其中无关刺激出现的概率为70％，探测刺激和靶刺激出现的概率均为15％；
对被测试者赋予任一种刺激时，在被测试者赋予刺激前的时间段t1内和对被测试者赋予刺激后的时间段t2内，均采用所述Fz电极对被测试者头部中央点的脑电信号进行实时提取，并通过脑电信号采集设备(3)按照预先设定的采样频率f对所述Fz电极提取的脑电信号同步进行采集且将所采集的脑电信号同步传送至处理器(1)；
步骤A2、将赋予多次探测刺激时脑电信号采集设备(3)采集的多个脑电信号xFz(i)进行叠加平均，得到探测刺激下的脑电信号XFz(i)，再按照步骤301至步骤302中所述的方法计算得出脑电信号XFz(i)的Hurst指数HFzM和多重分形谱宽度WFzM；
同时，将赋予多次无关刺激时脑电信号采集设备(3)采集的多个脑电信号xFz(i)进行叠加平均，得到无关刺激下的脑电信号X’Fz(i)，再按照步骤301至步骤302中所述的方法计算得出脑电信号X’Fz(i)的Hurst指数HFzm和多重分形谱宽度WFzm；
步骤A3、根据公式和公式计算得出Hurst指数判断阈值HFz和多重分形谱宽度判断阈值WFz。
7.按照权利要求1或2所述的基于多重分形消除趋势波动分析的脑电测谎方法，其特征在于：步骤一中所述脑电信号提取装置还包括辅助测试用脑电电极；所述辅助测试用脑电电极的数量至少为两个且其数量不大于3个；所述辅助测试用脑电电极为按照国际标准“10/20”电极安放方法布设在FC1、FC2或Fz点的脑电电极；
步骤一中所述脑电信号提取装置的脑电测谎判断阈值还包括各辅助测试用脑电电极的脑电测谎判断阈值，所述辅助测试用脑电电极的脑电测谎判断阈值包括Hurst指数判断阈值HAs和多重分形谱宽度判断阈值WAs；
布设FC1点的脑电电极的Hurst指数判断阈值HAs＝0.9587～1.0242且其多重分形谱宽度判断阈值WAs＝0.7247～0.8522；
布设FC2点的脑电电极的Hurst指数判断阈值HAs＝0.9497～1.0392且其多重分形谱宽度判断阈值WAs＝0.6947～0.8082；
布设Fz点的脑电电极的Hurst指数判断阈值HAs＝0.9035～0.9985且其多重分形谱宽度判断阈值WAs＝0.6972～0.8377；
步骤二中进行脑电信号采集时，对被测试者赋予刺激前的时间段t1内和对被测试者赋予刺激后的时间段t2内，均采用各辅助测试用脑电电极对被测试者头部对应部位的脑电信号进行实时提取，并通过脑电信号采集设备(3)按照预先设定的采样频率f对各辅助测试用脑电电极提取的脑电信号同步进行采集且将所采集的脑电信号同步传送至处理器(1)；
所述辅助测试用脑电电极提取的脑电信号记作As脑电信号，所述As脑电信号为因赋予刺激诱发的诱发脑电信号且其诱发电位为事件相关电位，所述事件相关电位为P300；
步骤三中进行脑电信号分析处理，过程如下：
步骤301中进行脑电信号接收与同步存储时，所述处理器(1)还需将所接收到的各脑电信号xAs(i)同步存储至存储器(4)内；
每个所述脑电信号xAs(i)包括N个采样周期内同一个所述辅助测试用脑电电极采集的As脑电信号；
步骤302中进行脑电信号特征提取时，所述处理器(1)还需调用基于多重分形消除趋势波动分析模块分别计算得出各脑电信号xAs(i)的Hurst指数HAs’和多重分形谱宽度WAs’；
步骤303中得出的说谎与否判断结果为根据脑电信号xPz(i)判断得出的说谎与否判断结果；
步骤303中完成根据脑电信号xPz(i)判断得出的说谎与否判断后，还需根据各脑电信号xAs(i)分别进行说谎与否判断；
其中，对任一个脑电信号xAs(i)进行说谎与否判断时，根据步骤302中得出的该脑电信号xAs(i)的Hurst指数HAs’和多重分形谱宽度WAs’，所述处理器(1)调用差值比较模块对此时被测试者说谎与否进行判断：当该脑电信号xAs(i)的Hurst指数HAs’大于步骤一中设定的与该脑电信号xAs(i)对应的辅助测试用脑电电极的Hurst指数判断阈值HAs且该脑电信号xAs(i)的多重分形谱宽度WAs’小于步骤一中设定的与该脑电信号xAs(i)对应的辅助测试用脑电电极的多重分形谱宽度判断阈值WAs时，判断为此时被测试者没有说谎；否则，判断为此时被测试者说谎；此时，所得出的说谎与否判断结果为根据脑电信号xAs(i)判断得出的说谎与否判断结果；
步骤304中进行说谎与否判断结果输出之前，还需判断根据脑电信号xPz(i)判断得出的说谎与否判断结果和根据各脑电信号xAs(i)判断得出的说谎与否判断结果是否一致：当根据脑电信号xPz(i)判断得出的说谎与否判断结果和根据各脑电信号xAs(i)判断得出的说谎与否判断结果均一致时，进入步骤304；否则，所述处理器(1)输出本次测谎无效。
8.按照权利要求1或2所述的基于多重分形消除趋势波动分析的脑电测谎方法，其特征在于：步骤一中所述脑电信号提取装置为32导联脑电电极或64导联脑电电极。
CN201610207819.4A 2016-04-05 2016-04-05 基于多重分形消除趋势波动分析的脑电测谎方法 Active CN105615879B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610207819.4A CN105615879B (zh) 2016-04-05 2016-04-05 基于多重分形消除趋势波动分析的脑电测谎方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610207819.4A CN105615879B (zh) 2016-04-05 2016-04-05 基于多重分形消除趋势波动分析的脑电测谎方法

Publications (2)

Publication Number Publication Date
CN105615879A CN105615879A (zh) 2016-06-01
CN105615879B true CN105615879B (zh) 2018-07-17

Family

ID=56031543

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610207819.4A Active CN105615879B (zh) 2016-04-05 2016-04-05 基于多重分形消除趋势波动分析的脑电测谎方法

Country Status (1)

Country Link
CN (1) CN105615879B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108498106B (zh) * 2018-02-08 2021-10-29 陕西师范大学 基于多重分形消除趋势波动分析的cnv脑电测谎方法
CN109700458B (zh) * 2019-01-14 2021-09-24 广西医科大学第一附属医院 一种eeg脑功能网络构建方法、装置及存储介质
CN110192880A (zh) * 2019-05-24 2019-09-03 中南民族大学 基于多导脑电信号格兰杰因果的测谎方法
CN110192878A (zh) * 2019-05-24 2019-09-03 中南民族大学 基于多导脑电信号定向转移函数的测谎方法
CN111616702A (zh) * 2020-06-18 2020-09-04 北方工业大学 一种基于认知负荷增强的测谎分析系统
CN114533066B (zh) * 2022-04-28 2022-08-19 之江实验室 基于复合表情加工脑网络的社交焦虑评估方法和系统

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1883384A (zh) * 2006-06-22 2006-12-27 复旦大学 一种脑电信号伪迹成分的自动识别和去除方法
WO2010093007A1 (ja) * 2009-02-12 2010-08-19 国立大学法人長岡技術科学大学 感性状態判定装置
WO2011159592A1 (en) * 2010-06-15 2011-12-22 Flints Hills Scientific, Llc A systems approach to disease state, health, and comorbidity
WO2013049156A1 (en) * 2011-09-26 2013-04-04 President And Fellows Of Harvard College Quantitative methods and systems for neurological assessment
WO2014176286A1 (en) * 2013-04-22 2014-10-30 The Regents Of The University Of California Fractal index analysis of human electroencephalogram signals
CN105249963A (zh) * 2015-11-16 2016-01-20 陕西师范大学 基于样本熵的n400诱发电位测谎方法
CN105426822A (zh) * 2015-11-05 2016-03-23 郑州轻工业学院 基于双树复小波变换的非平稳信号多重分形特征提取方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1883384A (zh) * 2006-06-22 2006-12-27 复旦大学 一种脑电信号伪迹成分的自动识别和去除方法
WO2010093007A1 (ja) * 2009-02-12 2010-08-19 国立大学法人長岡技術科学大学 感性状態判定装置
WO2011159592A1 (en) * 2010-06-15 2011-12-22 Flints Hills Scientific, Llc A systems approach to disease state, health, and comorbidity
WO2013049156A1 (en) * 2011-09-26 2013-04-04 President And Fellows Of Harvard College Quantitative methods and systems for neurological assessment
WO2014176286A1 (en) * 2013-04-22 2014-10-30 The Regents Of The University Of California Fractal index analysis of human electroencephalogram signals
CN105426822A (zh) * 2015-11-05 2016-03-23 郑州轻工业学院 基于双树复小波变换的非平稳信号多重分形特征提取方法
CN105249963A (zh) * 2015-11-16 2016-01-20 陕西师范大学 基于样本熵的n400诱发电位测谎方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
Functional brain network and multichannel analysis for the P300-based brain computer interface system of lying detection;Hong Wang et al.;《Expert Systems with Applications》;20160701;117-128 *
基于P300幅值几何差和脑网络特征的测谎方法研究;常文文等;《仪器仪表学报》;20150430;822-829 *
基于脑电样本熵的测谎分析;高军峰等;《电子学报》;20170831;1836-1841 *
多导脑电复杂度特征的谎言测试研究;高军峰等;《电子科技大学学报》;20170731;636-640 *

Also Published As

Publication number Publication date
CN105615879A (zh) 2016-06-01

Similar Documents

Publication Publication Date Title
CN105615879B (zh) 基于多重分形消除趋势波动分析的脑电测谎方法
CN105852885B (zh) 基于改进超限学习机的n400诱发电位测谎方法
CN105249963B (zh) 基于样本熵的n400诱发电位测谎方法
US6356781B1 (en) Functional magnetic resonance imaging capable of detecting the occurrence of neuronal events with high temporal accuracy
CN108498106A (zh) 基于多重分形消除趋势波动分析的cnv脑电测谎方法
CN103399639B (zh) 基于ssvep和p300联合脑机接口方法及装置
US6804661B2 (en) Drug profiling apparatus and method
CN104173064B (zh) 基于心率变异分析的测谎方法及测谎装置
CN100411580C (zh) 基于人脑记忆特征客观检测的身份识别方法
CN101477405A (zh) 基于左右视野两个频率刺激的稳态视觉诱发脑机接口方法
US20050124863A1 (en) Drug profiling apparatus and method
CN108309329B (zh) 一种基于脑网络分析的cnv脑电测谎方法
CN103472922A (zh) 一种基于p300与ssvep混合式脑机接口的目的地选择系统
CN104571504B (zh) 一种基于想象动作的在线脑‑机接口方法
CN106874722A (zh) 一种基于心电信号的身份识别方法及其装置
CN105701336B (zh) 基于脑电数据的中医辨证分型系统、模型建立的方法和系统
CN110192874B (zh) 基于多导脑电信号裕度因子的测谎方法
CN108985029A (zh) 基于视觉刺激的脑电身份识别方法
CN105769221B (zh) 基于Lyapunov指数和关联维数的脑电测谎方法
CN111616702A (zh) 一种基于认知负荷增强的测谎分析系统
CN107550486A (zh) 一种利用去除qrs波形后的体表心电预测房颤发生的方法
CN110192875B (zh) 基于多导脑电信号波形因子的测谎方法
CN110613446A (zh) 一种信号处理方法及装置
CN107169434A (zh) 一种具备排他性的脑电身份识别方法
Jiao et al. Multilayer correlation maximization for frequency recognition in SSVEP brain-computer interface

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant