技术领域
本发明属于脑电信号处理技术领域，尤其是涉及一种基于多重分形消除 趋势波动分析的脑电测谎方法。
背景技术
测谎技术作为司法审讯的一种辅助手段，对刑事案件事实的侦查和证 明发挥着重要的作用，现已得到广泛应用。如何准确有效地识别谎言，对 于司法机关工作人员来说就显得尤为重要，研究谎言的检测方法具有重大 的现实意义。1912年，美国研制的多导测谎仪主要通过呼吸、脉搏、血压、 皮肤电压等生理变化指标来判别谎言，由于上述参数易受主观因素的影 响，使得检测结果难以保证准确无误，误报和漏报率高。1989年，Rosenfeld 提出事件相关电位(ERP)的测谎改进方法，其原理是利用受试者接受与 案情相关的刺激时产生的P300诱发电位的脑电波变化来判定谎情。
事件相关电位(ERP，event-relatedpotential)是一种特殊的脑诱 发电位，通过有意地赋予刺激以特殊的心理意义，利用多个或多样的刺激 所引起的脑电位，它反映了认知过程中大脑的神经电生理的变化，也被称 为认知电位，也就是指当人们对某课题进行认知加工时，从头颅表面记录 到的脑电位，也称为ERP脑电位，对应的脑电波信号为ERP脑电波信号。
事件相关电位(ERP)测谎研究中主要涉及的诱发电位是P300和N400， 其中P300诱发电位的脑电波信号为250ms～700ms的正向波，N400诱发 电位的脑电波信号为潜伏期为200ms～500ms的负向波。P300诱发电位的 脑电波信号的典型特征是潜伏期、波幅和波面积。由于P300诱发电位是 大脑在“无意识”过程中产生的，不会出现“把无当有”的情况，较科学 有效地帮助警方侦破案件。随后，众多研究者对ERP测谎技术提出多种改 进，如幅值差(即BAD方法)、线性相关系数(即BCD方法)、频谱估计 值(即功率谱估计)、小波系数、脑网络结构等，上述处理方法主要以脑 电信号的时域、频域特征参数为评价标准展开研究，如BAD方法和BCD方 法这两种经典的信号分析处理方法，都利用单通道的P300诱发电位的脑 电波信号幅值信息区分诚实和欺骗性反应。赵敏等提出的多通道幅值和受 试者的反应时间相结合的方法以及DeyS提出的小波特征提取同样能有效 地反映出欺骗和诚实之间的显著差异，并与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-relatedpotential) 是一种特殊的脑诱发电位，通过有意地赋予刺激以特殊的心理意义，利用 多个或多样的刺激所引起的脑电位，它反映了认知过程中大脑的神经电生 理的变化，也被称为认知电位，也就是指当人们对某课题进行认知加工时， 从头颅表面记录到的脑电位，也称为ERP脑电位，对应的脑电信号为ERP 脑电信号。其中，赋予刺激指的是对被测试者赋予一次刺激。此处，对被 测试者赋予刺激指的是对被测试者进行一次测谎提问。因而，步骤二中所 述的对被测试者赋予刺激指的是引起事件相关电位的刺激，具体是对被测 试者进行一次提问。
其中，步骤302中所述的多重分形消除趋势波动分析模块为多重分形 消除趋势波动分析法(multifractaldetrendedfluctuationanalysis， MF-DFA)模块。
多重分形消除趋势波动分析法(multifractaldetrended fluctuationanalysis，MF-DFA)是2002年Kantelhardt等[20]在消除趋 势波动分析法(DFA)的基础上提出的，详见《Multifractaldetrended fluctuationanalysisofnon-stationarytimeseries》，作者 KantelhardtJW,ZschiegnerSA,Koscielny-BundeE,etal.，2002 年发表于PhysicaA-statisticalMechanics&ItsApplications。
本实施例中，所述处理器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脑电电极， 投入成本低，经济实用。
以上所述，仅是本发明的较佳实施例，并非对本发明作任何限制，凡是 根据本发明技术实质对以上实施例所作的任何简单修改、变更以及等效结构 变化，均仍属于本发明技术方案的保护范围内。