一种脑电波中疼痛成分的提取和表达方式的制作方法

文档序号:11267521阅读:618来源:国知局
一种脑电波中疼痛成分的提取和表达方式的制造方法与工艺

本发明涉及一种疼痛成分的脑电波特征指标的提取方法和表达方式,可以是一体机或架构在无线物联网平台,可应用于居家、医院一体化的疼痛测量方法的实现。



背景技术:

当今时代脑科学是人们重点关注的领先科学。脑科学中,脑状态的测量技术发展是其中的重要内容,也是脑科学应用研究的目标之一。在临床医学领域,直到目前,并没有一种能够实时客观定量度量疼痛、镇痛、焦虑、紧张、瞻妄、遗忘、舒适等等脑状态的测量方法应用于实践。由于这种实时客观定量测量技术的缺失,临床疼痛治疗中,人们对于疼痛的度量只能依赖于疼痛的主观量表完成,主观量表是人们用主观的自我判断,以0-10个刻度,形容自己的疼痛程度。这样的方法充满了不确定性和不可重复性,致使有关疼痛的治疗效果无法准确度量,镇痛程度(疼痛阈值)也无法测量,对于药物的使用也无法实现客观和适度。特别是疼痛带来的焦虑、紧张、谵妄、遗忘更加无法度量和施加相应的处理。其次,互联网技术和通讯技术是引领科学技术进步的支撑力量之一,结合临床医疗服务,在互联网平台上以无线通讯手段实现医疗信息集成交换,对于慢性疼痛患者,比如癌痛病人,在家庭环境下的镇痛治疗并不充分,没有实时客观定量的疼痛指标,往往使医生无法根据患者的客观疼痛程度施加药物治疗,患者往往是在极度疼痛中去世。对于急性疼痛患者,如术后患者的疼痛治疗,更加需要包含疼痛深度在内的脑状态定量数据作为医生治疗疼痛的客观依据,互联网技术的发展,为医生实时了解居家化或病床上的患者的疼痛、镇痛、焦虑、紧张、瞻妄、遗忘、舒适程度提供了一种可以被应用的技术手段。

基于互联网通讯的急、慢性疼痛治疗、监护信息管理平台的建立,可以对于疼痛治疗和管理方法的科学研究提供强有力的支持手段。特别是对于患者的疼痛、镇痛、焦虑、紧张、瞻妄、遗忘、舒适程度的定量数据的掌握,对于分析药物疗效、总结治疗方法、调整治疗方案、研究分析疼痛机理、收集大样本长时间临床疼痛相关信息等方面有重要的实用价值。



技术实现要素:

本发明的目的是提供一种脑电波中疼痛成分的提取和表达方式,以解决通过采集脑电波信号,提取出其中的疼痛、镇痛、焦虑、紧张、谵妄、遗忘、舒适程度的特征成分的方法和表达方式,并形成疼痛测量设备,实现临床医疗服务和居家化健康管理中疼痛管理的应用的技术问题。

为了实现上述发明目的,本发明所采用的技术方法如下:

一种脑电波中疼痛成分的提取和表达方式,包括在各种环境状态下,人体由于疼痛引起的大脑脑电波中关于疼痛、镇痛、焦虑、紧张、瞻妄、遗忘、舒适的特征成分的提取方法和表达方式;利用脑电波作为原始信号处理单元,离散化采集到计算机中,施加数学计算分析,将脑电波分解成低频、中频、高频、超高频单元,提取其中和疼痛、镇痛、焦虑、紧张、瞻妄、遗忘、舒适相关的特征成分,经过归一化处理,表达为0-100的无量纲数据,实时反应人的疼痛、镇痛、焦虑、紧张、瞻妄、遗忘、舒适的实时客观定量深度;利用脑电传感器,佩戴于头部的前额部和耳部部位,无创采集脑部多个部位的电位信号,信号经由前置放大器和多级整合电路、以及模数转换器进入计算单元,信号带宽包含脑电波中的超过30hz以上的超高频成分,脑电采集导联设置为最少包含左右脑前额叶的两个导联或以上,经由前置放大、模数转换,进入单片计算机系统,再通过数据加密、压缩后,经由有线方式直接传送至计算机单元,或经由无线通讯控制单元的tcp/ip协议包装并发送到互联网数据服务器,计算机或服务器中的软件系统对接收到的脑电信号实时处理计算、显示、存储、转发,其中计算部分采用小波分析、波形识别算法,分解脑电波中的各种相关于频率和时域的有效成分和伪差成分,并组合成为新的数据流,再经由模糊识别分析、拟合回归分析、相关性分析、多维谱分析,提取出疼痛、镇痛、焦虑、紧张、瞻妄、遗忘、舒适的特征指标,经过归一化处理,获得从0-100的多个定量数据,实时客观定量反映人体的疼痛、镇痛、焦虑、紧张、瞻妄、遗忘、舒适的程度,得到疼痛深度、镇痛深度、焦虑深度、紧张深度、瞻妄深度、遗忘深度、舒适深度指数。

脑电信号采集单元具备有线及无线通讯能力,可以和计算机组合成为一体的疼痛深度测量系统,也可以通过互联网平台和服务器组合成为云计算模式的疼痛深度测量系统,形成对疼痛测量数据的远程医疗服务共享平台,形成疼痛测量的大数据收集新模式,

利用232接口标准、3g,4g、蓝牙、wifi通讯技术,实现近远程、多环境、实时的脑电采集传输终端,对于采集放大的脑电波模拟信号,单片计算机对其实现模数转换的离散化处理,将多路脑电波离散信号实时压缩、加密、打包,实现动态存储队列管理,识别通讯线路,自动切换通讯模式,控制通讯发送,实时数据包或者直接通过232通讯接口传送到计算机中或经由无线互联网平台,直接发送到云计算服务器中,经由计算机或云计算服务器的计算处理,形成的疼痛深度指数、镇痛深度指数、焦虑深度指数、紧张深度指数、谵妄深度指数、遗忘深度指数、舒适深度指数显示在计算机终端或移动手机端。

疼痛数值伴随镇痛、焦虑、紧张、谵妄、遗忘、舒适的心理精神状态的定量测量,自脑电波中同步提取出镇痛、焦虑、紧张、瞻妄、遗忘、舒适的定量表达,于疼痛深度的定量表达一起,完成疼痛测量的过程,达到疼痛测量的完整性和精准性。

在小波处理计算中,可以分解出这两种生物电信号序列,作为信号强度和干扰信号的标记特征使用;采用小波公式:

对脑电向量组

fi(x)=[x1x2x3…xm-2xm-1xm]

i:脑电波导联数量,m:向量元素数量,x:脑电波数据

实时计算处理,以多尺度滤波器组算法,分解出各尺度窗口下的小波基函数

(wf(2^j,fj(x)))j∈z

(wf(2^0,fj(x))),(wf(2^1,f0(x))),,,(wf(2^n,fn(x)))

j:小波基数量

由小波基函数和各尺度数据的反变换,获得一组时域重构函数如下:

fi,j(x)=∑wf(2^j,x)*x2^j(x)

fi,1(x),fi,2(x),,,fi,n(x);重构函数,x2^1(x),x2^2(x),,,x2^n(x),各尺度小波数据点;n:阶数;j:尺度数量,i:导联数;

各个重构函数代表了不同尺度下的脑电波,眼动电波和肌电波的表现;尺度也对应信号的频率成分,分布于脑电图的常规节律和高频节律;对于分解的重构函数的各个序列数据,采用模式识别算法中的波形识别算法提取出数据的特征点:

tj(x)∈z;

t:特征值向量,j:尺度,x:离散数据,代表波形特点;z:时域空间;

向量tj(x)中的特征数据包含特异点、幅度、变异、斜率、面积、自相关、互相关等计算结果,计算来自于基本算法:

数据序列:y(t,i)=(f(j)-f(j-1))/δt

j:离散数据下标,i:尺度;

获取序列y(t,i)中最大数值,可以获得特征指数之一,其中的正负反相点,就是特异点,特异点的数量由t值的大小表示;

对f(x)序列数据施加迭代微分算法:

d(j,i)=∑(f(j+i)-f(j+i-1)/(δt+i))

i:δt的增量,从1..n,j:数值序号;

对于矩阵d(j,i)中的各个向量,将向量中的数据点排序以及相加,选取每个向量中的最大值和相加和,作为斜率和幅度;

针对重构函数中的各个尺度向量fi(x)=[x1x2x3…xm-2xm-1xm],将处理得到的结果向量y(t,i)、d(j,i)生成模式类ω1,ω2,…,ωc,然后利用模式之间的距离函数计算各个重构函数的距离;得到重构函数的变异、自相关、互相关数值;

计算重构函数的积分:

si(t)=∫f(j,i)*δt

求得各个函数的面积数值;

tj(x)向量表达了重构函数的波形特点和规律;对于采集的脑电波中的眼动电信号和额肌电信号,分布于特定尺度的重构函数中,针对这些重构函数,首先,求取其中的一阶导数:

di(x)=(fi(x)-fi(x+m))/δx

x:横坐标,δx:横坐标增量;

对di(x)排序,得到极大值和极小值,设定阈值a,可以求取

di(x)的正负极性变异点,得到一组极值点向量:

mi(j);高点;

mi(j);低点;

j:极值点的数量;

对相关重构函数采用积分算法:

e=∫fi(t)^2*δt

e1=∫fj(t)^2*δt

得到肌电和眼动的功率值大小;针对上述结果数据,经过归一组合算法:

sq=exp(a*abs(mi-mi)+b*e);

a,b:加权系数,由信号表达范围决定;

可以得到指示干扰的信号强度的定量表达sq,作为显示结果之一;

重构函数中消除了眼动和肌电的脑电波重构成分,采用功率谱算法:

f(x):时域脑电波离散数据,x(w):频谱功率大小;

可以获得脑电波中功率谱的各个成分,包括αβδθ波段的数值,sef,mef等边缘频率,优势频率等数值:

fj={α,β,δ,θ,,,sef,mef};

结合波形识别所得到的特征向量,组合成为一组涵盖了时域和频域,以及非线性复杂度结果的数据向量组:

gi(x)={tj,fj};

i:尺度数量;

数据向量gi(x)组,作为脑电波的一次处理结果,被命名为脑电波一级处理的元数据组,元数据组中包含了数百个脑电波中的规律性的特征成分,可以被作为二次计算的基础数据;

提取各尺度重构函数元数据中的某些特征数据,计算其中各尺度相对应的数据的比率、各尺度数据的斜率、变化率以及积分面积:

y(x,t)={gi(x,t)-gi(x,t-1)/δt}∈z;

zi,j(x)={gi(x1)/gj(x1)}∈z;

s(x,t)={∫g(x,t)*δt}∈z;

d(x,t)={abs(g(x,t1)-g(x,t2))}∈z;

t:计算时间点,i,j:尺度序列,x:特征向量下标;z:计算空间大小;

对于上述各个数据序列,经过数据加权,得到一系列数据计算公式:

ei={a,b,c,d,,,n}*{yi,zi,si,di};

i:数据指标序号;

对ei数据施加归一化计算:

疼痛指数=(esp(e0))×100

镇痛指数=(esp(e1))×100

紧张指数=(esp(e2))×100

谵妄指数=(esp(e3))×100

焦虑指数=(esp(e4))×100

遗忘指数=(esp(e5))×100

舒适指数=(esp(e5))×100;

得到脑电波中的疼痛、镇痛、焦虑、紧张、瞻妄、遗忘、舒适的实时客观定量特征指标。

本发明的优点:

本发明通过实时采集脑电波信号,利用离散信号处理的数学计算算法,提取出脑电波中有关疼痛、镇痛、焦虑、紧张、瞻妄、遗忘、舒适感觉的特征成分的方法,并表达成为0-100的无量纲数据的表达方式。利用电生理的信号传感技术采集脑电波信号,并对信号施加多种数学计算算法,找到并提取其中的疼痛、镇痛、焦虑、紧张、瞻妄、遗忘、舒适成分的特征指标,经过归一化处理,得到测量以上特征指标的定量数值,形成实时疼痛测量的技术手段和应用设备。本发明对于疼痛治疗和管理方法的科学研究提供强有力的支持手段,对于分析药物疗效、总结治疗方法、调整治疗方案、研究分析疼痛机理、收集大样本长时间临床疼痛相关信息等方面有重要的实用价值。

附图说明

图1为本发明的脑电信号采集电路。

图2为本发明的移动信息采集传输设备各部分之间的连接关系。

具体实施方式

本发明的原理及结构参见图1、2所示。一种脑电波中疼痛成分的提取和表达方式,包括在各种环境状态下,人体由于疼痛引起的大脑脑电波中关于疼痛、镇痛、焦虑、紧张、瞻妄、遗忘、舒适的特征成分的提取方法和表达方式;利用脑电波作为原始信号处理单元,离散化采集到计算机中,施加多种数学计算算法,将脑电波分解成低频、中频、高频、超高频单元,提取其中和疼痛、镇痛、焦虑、紧张、瞻妄、遗忘、舒适相关的特征成分,经过归一化处理,表达为0-100的无量纲数据,实时反应人的疼痛、镇痛、焦虑、紧张、瞻妄、遗忘、舒适的实时客观定量深度;

利用脑电传感器,佩戴于头部的前额部和耳部部位,无创采集脑部多个部位的电位信号,信号经由前置放大器和多级整合电路、以及模数转换器进入计算单元,信号带宽包含脑电波中的超高频成分(超过30hz以上的脑电波成分),脑电采集导联设置为最少包含左右脑前额叶的两个导联或以上,经由前置放大、模数转换,进入单片计算机系统,再通过数据加密、压缩后,经由有线方式直接传送至计算机单元,或经由无线通讯控制单元的tcp/ip协议包装并发送到互联网数据服务器,计算机或服务器中的专用软件系统对接收到的脑电信号实时处理计算、显示、存储、转发,其中计算部分采用小波分析、波形识别算法,分解脑电波中的各种相关于频率和时域的有效成分和伪差成分,并组合成为新的数据流,再经由模糊识别分析、拟合回归分析、相关性分析、多维谱分析等多种算法,提取出疼痛、镇痛、焦虑、紧张、瞻妄、遗忘、舒适的特征指标,经过归一化处理,获得从0-100的多个定量数据,实时客观定量反映人体的疼痛、镇痛、焦虑、紧张、瞻妄、遗忘、舒适的程度,得到疼痛深度、镇痛深度、焦虑深度、紧张深度、瞻妄深度、遗忘深度、舒适深度指数。

脑电信号采集单元具备有线及无线通讯能力,可以和计算机组合成为一体的疼痛深度测量系统,也可以通过互联网平台和服务器组合成为云计算模式的疼痛深度测量系统,形成对疼痛测量数据的远程医疗服务共享平台,形成疼痛测量的大数据收集新模式,利用232接口标准、3g,4g、蓝牙、wifi通讯技术,实现近远程、多环境、实时的脑电采集传输终端,对于采集放大的脑电波模拟信号,单片计算机对其实现模数转换的离散化处理,将多路脑电波离散信号实时压缩、加密、打包,实现动态存储队列管理,识别通讯线路,自动切换通讯模式,控制通讯发送,实时数据包或者直接传送到计算机中(利用232通讯接口)或经由无线互联网平台,直接发送到云计算服务器中,经由计算机或云计算服务器的计算处理,形成的疼痛深度指数、镇痛深度指数、焦虑深度指数、紧张深度指数、谵妄深度指数、遗忘深度指数、舒适深度指数显示在计算机终端或移动手机端。

疼痛数值伴随镇痛、焦虑、紧张、谵妄、遗忘、舒适的心理精神状态的定量测量,自脑电波中同步提取出镇痛、焦虑、紧张、瞻妄、遗忘、舒适的定量表达,于疼痛深度的定量表达一起,完成疼痛测量的过程,达到疼痛测量的完整性和精准性;疼痛是一种脑状态的表现,镇痛、焦虑、紧张、谵妄、遗忘、舒适也是脑状态的表现,疼痛可以引起人们的焦虑、紧张和谵妄、遗忘;焦虑、紧张又具有加强疼痛程度的作用,特别是在紧张和焦虑状态,人们承受疼痛的阈值(镇痛)将会随之下降;焦虑、紧张、谵妄、遗忘、舒适的程度可以伴随疼痛的发生,也可以发生在没有疼痛的状态下。

本发明涉及制造一个模块或一台设备,移动信息采集传输终端以及相关移动通讯终端、固定终端、计算和数据管理软件等,包含多导脑电波的采集处理部分和无创电极、通讯控制模块以及计算机、服务器的疼痛、镇痛、焦虑、紧张、瞻妄、遗忘、舒适程度的脑电波特征成分提取方法和归一化表达方式。

移动信息采集传输设备包含了前置放大电路、单片机控制电路、动态数据链路缓存电路、组合式无线通讯接入控制电路、电源电路等部分。

通过病人身上的脑电波感受器电极(生物电信号传感器),将患者的多导脑电波信号采集传送至前置放大电路的滤波、噪音控制、放大输入部分,经相应通道分别送入终端的单片计算机控制部分的模、数转换电路中。单片计算机控制部分获得数字化的脑电波(包括超过30hz以上的高频脑电波后)数据,经过加密、压缩的数据处理过程,得到处理后的数据流,送入动态数据链路缓存队列中。动态数据链路缓存队列是变化的数据存储和输出结构。根据网络状态的不同,数据在存储区域内的结构不同。由单片机控制电路在获得网络状态变化的中断事件触发后,控制采集到的数据的不同排列组合。在写指令控制下,写入存储队列数据。队列中的数据在单片控制电路的读指令控制下,经数据端口输出到蓝牙、wifi和3g、4g互联网络通讯模块。互联网络通讯模块完成对网络的自动拨号、网络状态识别、tcp模式信号调制、分包、输出的功能。采集的多导脑电数据可以直接以有线、无线方式送入到计算机中,并同时通过互联网平台自动传输至云计算服务器中。

数据中心服务器管理所有的移动信息采集传输终端发送的数据,每个终端设定唯一的地址编码,由机器号、服务器号加服务器的固定ip地址组成:

终端机器地址=服务器号范围内的唯一地址号+服务器号+组网服务器的固定ip地址

组网内最大机器数量=65535

移动信息采集传输终端对于转换成数字信号的数据实施加密化处理和压缩。整合的数据流进入链路存储队列中。数据窗口为7-15个bety

l流窗口=m+addr+asyn+data

m为模块传送的加密标记,addr为机器地址,asyn为同步,data为加密数据,包括数值型数据和波形数据

终端数据发送和状态接受、命令接受是系统实用化的根本保证。也是系统能够安全使用的关键环节之一。当利用近距离的计算机作为处理平台时,终端数据直接通过计算机上配置的wifi或蓝牙通讯端口接收数据,并由计算机的操作系统和疼痛测量应用软件完成数据计算、显示、存储的所有工作。同时终端将原始脑电波数据发送到互联网平台,设定的云处理服务器链中,数据被配有专用webserver应用数据流控制管理软件的相对应服务器接受、计算、存储、显示,并以相同的加密格式将结果数据发送至互联网平台,供相关人员在不同地理位置、不同的终端上共享。由于gprs或3g、4g通讯技术在移动中并不是一个完全稳定的无线通讯平台。诸多环节可以影响数据传送的效果。对于数值型数据的传送,可以使用udp包。采用多次发送的方式,保证服务器可以接受到一次有效数据实体。对于波形图像的传输,要求连续性和排他性。波形图像的传送必须采用tcp包。利用单片计算机的流控制能力,实时处理数据发送的状态和内部缓存大小。自动识别和自我修正网络通讯的实时进程,进而控制数据流的方向和网络呼叫拨号。设置了多达六分钟的数据链路动态队列区。保证数据的完整和容错能力的提高。

数据在链路存储队列中等待入网。网络状态触发电路触发一个网络状态事件,计算控制单元的事件处理流程根据事件性质,分别进入不同的线程处理单元。在正常状态下,取得链路数据,以tcp协议分包,通过无线网络向网络中设置的网址和端口发送协议包。

tcp协议包+=(save_point++)*real_data→save_point++

save_point为数据地址指针,real_data为实时数据。

动态数据链路管理包含了内存数据队列管理,数据加密管理,数据压缩管理。涵盖了当前数据发送的动态相互协调的功能。

系统的数据传输,采用公网和无线模式,数据的加密采用数据结构随数据内容动态改变的方法,改变规则其一:是根据各导联数值的动态大小,自动排列各导联在数据流中的位置。其二,采用16位二进制码加权算法,对数据流实时加权计算,生成加密的数据流发送至互联网平台。其三,数据流不包含使用者的个人身份信息,对应信息在中心工作站的使用登录系统中绑定。

针对移动应用,各个移动信息采集传输终端的实时数据首先传输到终端地址指定的云计算服务器链中的某一台服务器中,服务器运行专用的云计算处理实时软件包、数据转发管理软件包、用户身份同步管理软件包。服务器完成数据的计算、身份匹配、原始和结果数据存储、转发、命令生成等等智能化信息处理监控软件。服务器云计算实时提取疼痛、镇痛、焦虑、紧张、瞻妄、遗忘、舒适等等脑状态的定量表达,分布于不同地点的计算机、手机、移动平板等等智能终端负责数据的实时显示和报警。服务器需要连接公共internet互联网,要求在路由器上配置固定ip地址。服务器可以和医院内网实现互联,建立内外网的vpn通道,在多重防火墙保护下使计算结果信息进入医院的内网,连接医院内网的计算机(包括台式计算机、手机、平板电脑等等)都可以作为数据显示终端使用,由使用人员的权限控制进入本系统的软件服务。云计算服务器链路中配置一台用户身份管理服务器,作为使用者(包括医生、关注者、患者)的身份登记工作站使用,通常此服务器设置于医院的内网中。其上的软件将数据采集传输终端地址与使用者身份绑定,并分配权限。所有使用过程的疼痛、镇痛、焦虑、紧张、瞻妄、遗忘、舒适深度的测量结果数据和使用者的时间、身份等等所有信息统一存储于其相应的数据库中保存。

对于医疗实时监护的非移动应用场合,为了满足及时性的数据显示表达的需求,采用数据实时上传到互联网平台的上述数据处理模式,同时,终端通过其上配置的wifi或蓝牙通讯模块和近距离的计算机(包括台式计算机、平板等等)的wifi、蓝牙通讯接口实现数据交换。计算机上的数据计算本地版软件对数据实施计算处理,获得的结果显示在计算机屏幕上,完成医疗用监护仪的传统功能模式的应用。同时,云计算结果同步的传送到医务人员的手机上,实现随身手机监护仪的应用模式。

医务人员的个人电脑、手机可以通过internet互联网在相应权限下直接获取不同患者的疼痛诊疗信息,个人电脑上要安装专用软件包,手机需要下载专用的app应用软件。权限管理要划入医院整体管理权限范围。同时,所有操作记录将发回云计算服务器,作为追踪、挖掘、分析的原始数据保存。医患之间和一线二线医生之间的实时在线互动由连接在服务器端的计算机或个人移动通讯终端之间实现,交流信息通过服务器的双向转发,达到实时双向通讯的目的。实时互动功能由系统软件功能实现。

数据采集传输终端所采用的无线通讯控制接入模块,可以采用手机代替,手机通过其上配置的蓝牙功能,在疼痛测量应用app软件的支持下,接受终端发送的加密数据流。并经过显示、呼叫、转发的过程,将数据实时同步发送至互联网平台。手机和终端之间设置蓝牙通讯,配对,手机app中采用socket线程,初始化蓝牙通道,进入接收流程,对于接收到的数据信号,解码,并实时显示在手机屏幕上。同时对于解码的数据,提取出其中的云计算服务器地址和端口地址,将接受的原始数据不加任何处理,向相应的云计算服务器发送,app利用另一个socket线程接受相应服务器发送的命令数据和结果数据,数据的解密计算格式相同于终端发送的数据计算格式。对于结果数据,app将数据实时显示在手机的固定屏幕位置,供使用者或医生实时观察疼痛深度等身体状态的变化。疼痛测量应用软件的下载需要用户设定相应的真实身份信息,利用身份证号作为唯一用户识别码,在下载的同时,在用户管理服务器中建立使用者的个人档案,其中包含患者和使用的医生的身份信息,如果有数据保密和归属要求,服务器架构可以设置在医院的内网上,再通过vpn加密防护通道进入互联网平台。一旦生成了用户信息后,手机app软件就可以成为数据的转发和结果显示终端了。下载疼痛应用软件的手机数量不限。

对于脑电波的处理,以采样频率为1400/s、采样时间窗口为1.25s、采样精度为10位bite的离散化处理,采用波形识别、小波分析算法,计算出大脑的不同成分序列和各序列的时域特征、多尺度复杂度特征。并采用功率谱算法,提取出脑电波中的节律成分和各节律波段的功率大小,得到一级处理计算的计算结果序列,本序列对应脑电波中的数百个独立或弱相关的规律特征。对于额部双侧采集的脑电波中,包含有眼动电成分和额肌电成分。在小波处理计算中,可以分解出这两种生物电信号序列,作为信号强度和干扰信号的标记特征使用。采用小波公式:

对脑电向量组

fi(x)=[x1x2x3…xm-2xm-1xm]

i:脑电波导联数量,m:向量元素数量,x:脑电波数据

实时计算处理,以多尺度滤波器组算法,分解出各尺度窗口下的小波基函数

(wf(2^j,fj(x)))j∈z

(wf(2^0,fj(x))),(wf(2^1,f0(x))),,,(wf(2^n,fn(x)))

j:小波基数量

由小波基函数和各尺度数据的反变换,获得一组时域重构函数如下:

fi,j(x)=∑wf(2^j,x)*x2^j(x)

fi,1(x),fi,2(x),,,fi,n(x);重构函数,x2^1(x),x2^2(x),,,x2^n(x),各尺度小波数据点。n:阶数;j:尺度数量,i:导联数。

各个重构函数代表了不同尺度下的脑电波,眼动电波和肌电波的表现。尺度也对应信号的频率成分,分布于脑电图的常规节律和高频节律。对于分解的重构函数的各个序列数据,采用模式识别算法中的波形识别算法提取出数据的特征点:

tj(x)∈z;

t:特征值向量,j:尺度,x:离散数据,代表波形特点。z:时域空间。

向量tj(x)中的特征数据包含特异点、幅度、变异、斜率、面积、自相关、互相关等计算结果,计算来自于基本算法:

数据序列:y(t,i)=(f(j)-f(j-1))/δt

j:离散数据下标,i:尺度。

获取序列y(t,i)中最大数值,可以获得特征指数之一,其中的正负反相点,就是特异点,特异点的数量由t值的大小表示。

对f(x)序列数据施加迭代微分算法:

d(j,i)=∑(f(j+i)-f(j+i-1)/(δt+i))

i:δt的增量,从1..n,j:数值序号。

对于矩阵d(j,i)中的各个向量,将向量中的数据点排序以及相加,选取每个向量中的最大值和相加和,作为斜率和幅度。

针对重构函数中的各个尺度向量fi(x)=[x1x2x3…xm-2xm-1xm],将处理得到的结果向量y(t,i)、d(j,i)生成模式类ω1,ω2,…,ωc,然后利用模式之间的距离函数计算各个重构函数的距离。得到重构函数的变异、自相关、互相关数值。

计算重构函数的积分:

si(t)=∫f(j,i)*δt

求得各个函数的面积数值。

tj(x)向量表达了重构函数的波形特点和规律。对于采集的脑电波中的眼动电信号和额肌电信号,分布于特定尺度的重构函数中,针对这些重构函数,首先,求取其中的一阶导数:

di(x)=(fi(x)-fi(x+m))/δx

x:横坐标,δx:横坐标增量。

对di(x)排序,得到极大值和极小值,设定阈值a,可以求取

di(x)的正负极性变异点,得到一组极值点向量:

mi(j);高点。

mi(j);低点。

j:极值点的数量。

对相关重构函数采用积分算法:

e=∫fi(t)^2*δt

e1=∫fj(t)^2*δt

得到肌电和眼动的功率值大小。针对上述结果数据,经过归一组合算法:

sq=exp(a*abs(mi-mi)+b*e);

a,b:加权系数,由信号表达范围决定。

可以得到指示干扰的信号强度的定量表达sq,作为显示结果之一。

重构函数中消除了眼动和肌电的脑电波重构成分,采用功率谱算法:

f(x):时域脑电波离散数据,x(w):频谱功率大小。

可以获得脑电波中功率谱的各个成分,包括αβδθ波段的数值,sef,mef等边缘频率,优势频率等数值:

fj={α,β,δ,θ,,,sef,mef};

结合波形识别所得到的特征向量,组合成为一组涵盖了时域和频域,以及非线性复杂度结果的数据向量组:

gi(x)={tj,fj};

i:尺度数量。

数据向量gi(x)组,作为脑电波的一次处理结果,被命名为脑电波一级处理的元数据组,元数据组中包含了数百个脑电波中的规律性的特征成分,可以被作为二次计算的基础数据。

提取各尺度重构函数元数据中的某些特征数据,计算其中各尺度相对应的数据的比率、各尺度数据的斜率、变化率以及积分面积:

y(x,t)={gi(x,t)-gi(x,t-1)/δt}∈z;

zi,j(x)={gi(x1)/gj(x1)}∈z;

s(x,t)={∫g(x,t)*δt}∈z;

d(x,t)={abs(g(x,t1)-g(x,t2))}∈z;

t:计算时间点,i,j:尺度序列,x:特征向量下标。z:计算空间大小。

对于上述各个数据序列,经过数据加权,得到一系列数据计算公式:

ei={a,b,c,d,,,n}*{yi,zi,si,di};

i:数据指标序号。

对ei数据施加归一化计算:

疼痛指数=(esp(e0))×100

镇痛指数=(esp(e1))×100

紧张指数=(esp(e2))×100

谵妄指数=(esp(e3))×100

焦虑指数=(esp(e4))×100

遗忘指数=(esp(e5))×100

舒适指数=(esp(e5))×100;

得到脑电波中的疼痛、镇痛、焦虑、紧张、瞻妄、遗忘、舒适的实时客观定量特征指标。

当前第1页1 2 
网友询问留言 已有0条留言
  • 还没有人留言评论。精彩留言会获得点赞!
1