基于内源性脑信号的闭环神经调控系统、方法和设备

1.本发明属于神经科学和电磁学领域，具体涉及了一种基于内源性脑信号的闭环神经调控系统、方法及设备。


背景技术：

2.脑神经调控技术是通过侵入性或非侵入性手段，利用光、磁、电、超声等物理或化学等外部技术手段改变脑部内源神经信号传递，从而引起脑功能变化的重要方法。神经调控不仅能在较短时间内引起神经元结构变化，而且特定设计的调控范式还可带来功能环路的改变，从而恢复或者提高神经突触连接可塑性，是解析神经元信号传导、研究神经环路、阐明脑功能与行为因果关系的得力工具，又是治疗临床神经系统疾病的重要手段。
3.然而现有的神经调控技术仍处于技术突破和发展阶段，由于刺激参数固定、刺激流程单一、刺激效果难以定量评估，目前应用于临床的神经调控技术常为开环控制并未实现患者个体化闭环的治疗，存在着适应症有限、治疗稳定性不佳以及个体差异明显等问题。无创神经调控的效果随着大脑状态的不同而不同，这就启示我们当前开环神经调控存在的缺点，很可能是因为没有结合个体的内源性脑信号来调控，神经调控的参数需结合给定时间的大脑功能活动进行动态控制。
4.脑功能活动包括神经元活动和局部能量代谢等多个过程，复杂的功能活动使得脑汇集了多个模态的信息，其中最为重要的是神经元的电活动和激活区域的血氧代谢变化，只有实现这两种信息的有效提取、分析和融合，才能将脑功能活动有机的联系起来。光电同步脑活动检测仪可以同时检测脑电信号和脑血氧信号从而从两个维度解析脑功能信号。截至目前，尚没有结合光电同步脑活动检测仪来实时评估脑状态，并以此来控制神经调控的参数，实现个体化闭环调控目的。


技术实现要素：

5.为了解决现有技术中的上述问题，即现有的神经调控技术无法结合光电同步脑活动实时对神经进行调控，本发明提供了一种基于内源性脑信号的闭环神经调控系统，包括脑信号采集模块、调控处理模块和调控参数优化模块；
6.所述脑信号采集模块，配置为通过脑信号采集设备实时获取多模态脑信号；所述多模态脑信号包括血氧水平信号fnirs和脑电信号eeg；
7.所述调控处理模块，包括静息阶段单元、调控阶段单元和评价阶段单元；
8.所述静息阶段单元，配置为基于所述多模态脑信号，进行解码获取静息脱氧血红蛋白浓度信息hbr、静息氧合血红蛋白浓度信息hbo、eeg各子频段的静息阶段频点和eeg各子频段的静息阶段幅值，并进行子频段修正生成个体化子频段；
9.所述调控阶段单元，配置为通过基于所述多模态脑信号，进行解码获取最新ns的脱氧血红蛋白浓度信息hbr和氧合血红蛋白浓度信息hbo，通过带通滤波和平均处理的方式获得调控脱氧血红蛋白浓度信息和调控氧合血红蛋白浓度信息；基于所述多模态脑信号截
取最新ms的脑电信号eeg，并对各子频段进行个体化子频段的带通滤波、快速傅里叶变换和希尔伯特变换计算各子频段的调控阶段幅值和调控阶段相位；n和m为非负数；
10.基于所述调控脱氧血红蛋白浓度信息和调控氧合血红蛋白浓度信息设置调控设备刺激的幅度控制信号，基于所述幅度控制信号和调控阶段相位生成控制信号的调控时刻，基于所述调控设备刺激的幅值和调控时刻生成调控信号，通过调控设备输出调控信号；
11.所述评价阶段单元，在调控阶段单元完成调控之后，基于所述多模态脑信号，进行解码获取评价脱氧血红蛋白浓度信息、评价氧合血红蛋白浓度信息、各子频段的评价阶段频点和各子频段的评价阶段幅值；
12.基于所述评价脱氧血红蛋白浓度信息、评价氧合血红蛋白浓度信息、静息脱氧血红蛋白浓度信息和静息氧合血红蛋白浓度信息计算血氧调控效果；
13.基于所述各子频段的评价阶段幅值和各子频段的静息阶段幅值计算脑电调控效果；
14.所述调控参数优化模块，配置为基于所述血氧调控效果和脑电调控效果调整调控参数。
15.在一些优选的实施方式中，所述静息阶段单元包括：
16.配置为对所述脱氧血红蛋白浓度信息hbr、氧合血红蛋白浓度信息hbo进行0.01hz
‑
0.2hz的带通滤波，并对4个通道脑血氧信息平均处理获得静息脱氧血红蛋白浓度信息hbr
rest
和静息氧合血红蛋白浓度信息hbo
rest
；
17.对所述脑电信号eeg进行0.5hz
‑
100hz的带通滤波以及50hz的陷波处理，根据脑电频率标准分为δ频段、θ频段、α频段和β频段，通过快速傅里叶变换计算出各子频段的静息阶段幅值和并计算各子频段幅值最高的静息阶段频点和对各子频段进行修正获得个体化的子频段对各子频段进行修正获得个体化的子频段和和
18.在一些优选的实施方式中，所述调控阶段单元包括：
19.配置为通过ns的矩形窗或汉宁窗截取血氧水平信号fnirs，获得最新ns的脱氧血红蛋白浓度信息hbr和氧合血红蛋白浓度信息hbo，基于所述最新最新ns的脱氧血红蛋白浓度信息hbr和氧合血红蛋白浓度信息hbo进行0.01hz
‑
0.2hz的带通滤波并对4个通道脑血氧信息平均处理获得调控脱氧血红蛋白浓度信息hbr
now
和调控氧合血红蛋白浓度信息hbo
now
；
20.通过ms的矩形窗或汉宁窗截取所述脑电信号eeg，获得最新ms的脑电信号，基于所述最新ms的脑电信号进行0.5hz
‑
100hz的带通滤波及50hz的陷波处理，接着进行100hz的带通滤波及50hz的陷波处理，接着进行和和的带通滤波，再在4个子频段中分别进行快速傅里叶变换和希尔伯特变换计算静息阶段频点和对应的调控阶段幅值和与调控阶段幅值相应的调控阶段相位为
和
21.基于相对于基线的脑血氧幅度控制信号、相对于基线的脑电幅度控制信号确定触发阈值，基于所述调控阶段相位获取相位控制信号进而获取调控时刻，根据所述触发阈值和调控时刻生成调控信号；
22.所述相对于基线的脑血氧幅度控制信号为所述相对于基线的脑血氧幅度控制信号为其中a和b为待调整参数，取值范围为[
‑
1,1]；
[0023]
所述相对于基线的脑电幅度控制信号为所述相对于基线的脑电幅度控制信号为其中c、d、e和f为待调整参数，取值范围为[
‑
1,1]；
[0024]
所述相位控制信号为为其中，g、h、i和j为预设的参数，取值范围为[
‑
1,1]；
[0025]
所述触发阈值为modu
thr
＝m*hb
modu
+n*eeg
modu
，调控时刻为time＝ang
modu
，其中m和n为待调整参数，取值范围为[
‑
1,1]；
[0026]
所述调控信号dev
modu
为：
[0027]
其中dev
base
为预设的调控设备刺激强度基准值，p为待调整参数，取值范围为[
‑
1,1]。
[0028]
在一些优选的实施例中，所述评价阶段单元包括：
[0029]
配置为对所述脱氧血红蛋白浓度信息hbr、氧合血红蛋白浓度信息hbo进行0.01hz
‑
0.2hz的带通滤波，并进行4个通道的脑血氧信息平均处理获得评价脱氧血红蛋白浓度信息hbo
post
和评价氧合血红蛋白浓度信息hbo
post
；
[0030]
对所述脑电信号eeg进行0.5hz
‑
100hz的带通滤波以及50hz的陷波处理，并进行100hz的带通滤波以及50hz的陷波处理，并进行和的带通滤波，再在4个子频段中分别进行快速傅里叶变换计算评价阶段频点和和对应的评价阶段幅值和
[0031]
血氧调控效果为：其中r和s为待调整参数，取值范围为[
‑
1,1]；
[0032]
脑电调控效果为：
其中u、v、w和x为待调整参数，取值范围为[
‑
1,1]。
[0033]
在一些优选的实施方式中，所述脑信号采集模块包括：
[0034]
通过光电同步脑活动检测仪neg同时收集多模态脑信号，光电同步脑活动检测仪的电极和光极按照neg
‑
32或neg
‑
8采集帽排布；fnirs光通道排布在电极周围，32或8通道帽子覆盖全脑或特定区域。
[0035]
在一些优选的实施方式中，所述调控设备包括经颅磁刺激设备tms或经颅电刺激设备tes。
[0036]
在一些优选的实施方式中，所述静息阶段单元和评价阶段单元均持续3分钟，使用者保持非任务状态下的闭眼或睁眼状态且调控设备处于关闭状态；所述调控阶段单元持续30分钟，调控设备处于开启状态，使用者保持睁眼状态或任务状态。
[0037]
本发明的另一方面，提出了一种基于内源性脑信号的闭环神经调控方法，所述方法包括：
[0038]
步骤s100，通过脑信号采集设备实时获取多模态脑信号；所述多模态脑信号包括血氧水平信号fnirs和脑电信号eeg；
[0039]
步骤s200，其具体包括步骤s210
‑
步骤s230；
[0040]
步骤s210，基于所述多模态脑信号，进行解码获取静息脱氧血红蛋白浓度信息hbr、静息氧合血红蛋白浓度信息hbo、eeg各子频段的静息阶段频点和eeg各子频段的静息阶段幅值，并进行子频段修正生成个体化子频段；
[0041]
步骤s220，通过基于所述多模态脑信号，进行解码获取最新ns的脱氧血红蛋白浓度信息hbr和氧合血红蛋白浓度信息hbo，通过带通滤波和平均处理的方式获得调控脱氧血红蛋白浓度信息和调控氧合血红蛋白浓度信息；基于所述多模态脑信号截取最新ms的脑电信号eeg，并对各子频段进行个体化子频段的带通滤波、快速傅里叶变换和希尔伯特变换计算各子频段的调控阶段幅值和调控阶段相位；n和m为非负数；
[0042]
基于所述调控脱氧血红蛋白浓度信息和调控氧合血红蛋白浓度信息设置调控设备刺激的幅度控制信号，基于所述幅度控制信号和调控阶段相位生成控制信号的调控时刻，基于所述调控设备刺激的幅值和调控时刻生成调控信号，通过调控设备输出调控信号；
[0043]
步骤s230，在调控阶段单元完成调控之后，基于所述多模态脑信号，进行解码获取评价脱氧血红蛋白浓度信息、评价氧合血红蛋白浓度信息、各子频段的评价阶段频点和各子频段的评价阶段幅值；
[0044]
基于所述评价脱氧血红蛋白浓度信息、评价氧合血红蛋白浓度信息、静息脱氧血红蛋白浓度信息和静息氧合血红蛋白浓度信息计算血氧调控效果；
[0045]
基于所述各子频段的评价阶段幅值和各子频段的静息阶段幅值计算脑电调控效果；
[0046]
步骤s300，基于所述血氧调控效果和脑电调控效果调整调控参数。
[0047]
本发明的第三方面，提出了一种电子设备，其特征在于，包括：至少一个处理器；以及与至少一个所述处理器通信连接的存储器；其中，所述存储器存储有可被所述处理器执行的指令，所述指令用于被所述处理器执行以实现上述的基于内源性脑信号的闭环神经调控方法。
[0048]
本发明的第四方面，提出了一种计算机可读存储介质，其特征在于，所述计算机可读存储介质存储有计算机指令，所述计算机指令用于被所述计算机执行以实现权利要求8所述的基于内源性脑信号的闭环神经调控方法。
[0049]
本发明的有益效果：
[0050]
(1)本发明提出的基于内源性脑信号的闭环神经调控系统，通过实时采集内源性的脑信号，同步调整外部调控设备的调控参数，实现了外源性刺激信号闭环动态地控制内源性脑信号。
[0051]
(2)本发明提出的基于内源性脑信号的闭环神经调控系统，通过光电同步脑活动检测仪采集多模态的内源性脑信号，融合高时间分辨率、高空间分辨率的采集信号，提高了脑部状态评估的精确度，使神经调控更精确。
[0052]
(3)本发明提出的基于内源性脑信号的闭环神经调控系统，过多模态脑信号实时定量的改变调控参数，方便易用，实现了一体化脑活动同步采集和个体化调控。
附图说明
[0053]
通过阅读参照以下附图所作的对非限制性实施例所作的详细描述，本申请的其它特征、目的和优点将会变得更明显：
[0054]
图1是本发明实施例的基于内源性脑信号闭环调控系统的结构框图；
[0055]
图2是本发明实施例的基于内源性脑信号闭环调控系统的原理示意图；
[0056]
图3是本发明实施例中运用的光电同步脑活动检测仪neg设备；
[0057]
图4是本发明实施例中neg设备的一种电极和光极排布；
[0058]
图5是本发明实施例中neg设备采集的多模态脑信号；
[0059]
图6是本发明实施例中的基于内源性脑信号闭环调控系统应用于经颅磁刺激设备tms的示意图；
具体实施方式
[0060]
下面结合附图和实施例对本申请作进一步的详细说明。可以理解的是，此处所描述的具体实施例仅用于解释相关发明，而非对该发明的限定。另外还需要说明的是，为了便于描述，附图中仅示出了与有关发明相关的部分。
[0061]
需要说明的是，在不冲突的情况下，本申请中的实施例及实施例中的特征可以相互组合。下面将参考附图并结合实施例来详细说明本申请。
[0062]
本发明提供一种基于内源性脑信号的闭环神经调控系统，本系统通过实时采集内源性的脑信号，同步调整外部调控设备的调控参数，实现了外源性刺激信号闭环动态地控制内源性脑信号。
[0063]
本发明的一种基于内源性脑信号的闭环神经调控系统，包括脑信号采集模块、调控处理模块和调控参数优化模块；
[0064]
所述脑信号采集模块，配置为通过脑信号采集设备实时获取多模态脑信号；所述多模态脑信号包括血氧水平信号fnirs和脑电信号eeg；
[0065]
所述调控处理模块，包括静息阶段单元、调控阶段单元和评价阶段单元；
[0066]
所述静息阶段单元，配置为基于所述多模态脑信号，进行解码获取静息脱氧血红
蛋白浓度信息hbr、静息氧合血红蛋白浓度信息hbo、eeg各子频段的静息阶段频点和eeg各子频段的静息阶段幅值，并进行子频段修正生成个体化子频段；
[0067]
所述调控阶段单元，配置为通过基于所述多模态脑信号，进行解码获取最新ns的脱氧血红蛋白浓度信息hbr和氧合血红蛋白浓度信息hbo，通过带通滤波和平均处理的方式获得调控脱氧血红蛋白浓度信息和调控氧合血红蛋白浓度信息；基于所述多模态脑信号截取最新ms的脑电信号eeg，对各子频段进行个体化子频段的带通滤波、快速傅里叶变换和希尔伯特变换计算各子频段的调控阶段幅值和调控阶段相位；n和m为非负数；
[0068]
基于所述调控脱氧血红蛋白浓度信息和调控氧合血红蛋白浓度信息设置调控设备刺激的幅度控制信号，基于所述幅度控制信号和调控阶段相位生成控制信号的调控时刻，基于所述调控设备刺激的幅值和调控时刻生成调控信号，通过调控设备输出调控信号；
[0069]
所述评价阶段单元，在调控阶段单元完成调控之后，基于所述多模态脑信号，进行解码获取评价脱氧血红蛋白浓度信息、评价氧合血红蛋白浓度信息、各子频段的评价阶段频点和各子频段的评价阶段幅值；
[0070]
基于所述评价脱氧血红蛋白浓度信息、评价氧合血红蛋白浓度信息、静息脱氧血红蛋白浓度信息和静息氧合血红蛋白浓度信息计算血氧调控效果；
[0071]
基于所述各子频段的评价阶段幅值和各子频段的静息阶段幅值计算脑电调控效果；
[0072]
所述调控参数优化模块，配置为基于所述血氧调控效果和脑电调控效果调整调控参数。
[0073]
为了更清晰地对本发明基于内源性脑信号的闭环神经调控系统进行说明，下面结合图1对本发明实施例中各功能模块展开详述。
[0074]
本发明第一实施例的基于内源性脑信号的闭环神经调控系统，包括步骤脑信号采集模块、调控处理模块和调控参数优化模块，各功能模块详细描述如下：
[0075]
本实施例将基于内源性脑信号的闭环神经调控方法应用于经颅磁刺激tms的运动区闭环调控；
[0076]
所述脑信号采集模块，配置为通过脑信号采集设备实时获取多模态脑信号；所述多模态脑信号包括血氧水平信号fnirs和脑电信号eeg；
[0077]
在本实施例中，脑信号采集模块包括：
[0078]
通过光电同步脑活动检测仪neg如图3所示，同时收集多模态脑信号，光电同步脑活动检测仪的电极和光极按照neg
‑
32或neg
‑
8采集帽排布，支持32或8个通道的fnirs信号、32或8个脑电通道排布，具体排布如图4所示，使用neg
‑
32的光电同步脑活动检测仪，该设备具有32通道的脑电和32通道的近红外，采集帽子的排布按照10
‑
20排布形式，其中用到的eeg通道为fc1、fc3、fc5、c1、c3、c5、cp1、cp3、cp5等9个通道、用到的近红外通道为两个光源s和两个探头d组成的4个通道；fnirs光通道排布在电极周围，32或8通道帽子覆盖全脑或特定区域。采集到的信号如图5所示。使用neg
‑
32的光电同步脑活动检测仪可直接采集到血氧水平信号fnirs和脑电信号eeg；
[0079]
所述调控处理模块，包括静息阶段单元、调控阶段单元和评价阶段单元；
[0080]
在本实施例中，所述静息阶段单元和评价阶段单元均持续3分钟，使用者保持非任务状态下的闭眼或睁眼状态且调控设备处于关闭状态；所述调控阶段单元持续30分钟，调
控设备处于开启状态，使用者保持睁眼状态或任务状态。
[0081]
所述静息阶段单元，配置为基于所述多模态脑信号，进行解码获取静息脱氧血红蛋白浓度信息hbr、静息氧合血红蛋白浓度信息hbo、eeg各子频段的静息阶段频点和eeg各子频段的静息阶段幅值，并进行子频段修正生成个体化子频段；
[0082]
在本实施例中，静息阶段具体包括：配置为对所述脱氧血红蛋白浓度信息hbr、氧合血红蛋白浓度信息hbo进行0.01hz
‑
0.2hz的带通滤波，并对4个通道脑血氧信息平均处理获得静息脱氧血红蛋白浓度信息hbr
rest
和静息氧合血红蛋白浓度信息hbo
rest
；本阶段的血氧信息的通道数可根据具体应用场景中的设备设置或用户需求进行调整，此处不做具体限定；
[0083]
在本实施例中，为突出运动区c3的脑电信号，可对9个电极采取拉普拉斯空间滤波处理；
[0084]
对所述脑电信号eeg进行0.5hz
‑
100hz的带通滤波以及50hz的陷波处理，根据脑电频率标准分为δ频段(0.5
‑
5hz)、θ频段(5
‑
8hz)、α频段(8
‑
12hz)和β频段(12
‑
30hz)，通过快速傅里叶变换计算出各子频段的静息阶段幅值和并计算各子频段幅值最高的静息阶段频点和对各子频段进行修正获得个体化的子频段段和所述个体化的子频段可适应不同个体的脑电信号，避免个体差异对本发明的效果造成影响。
[0085]
所述调控阶段单元，配置为通过基于所述多模态脑信号，进行解码获取最新ns的脱氧血红蛋白浓度信息hbr和氧合血红蛋白浓度信息hbo，通过带通滤波和平均处理的方式获得调控脱氧血红蛋白浓度信息和调控氧合血红蛋白浓度信息；基于所述多模态脑信号，截取最新ms的脑电信号eeg，并对各子频段进行个体化子频段的带通滤波、快速傅里叶变换和希尔伯特变换计算各子频段的调控阶段幅值和调控阶段相位；n和m为非负数，在本实施例中，优选的n取值为3，m取值为2；
[0086]
基于所述调控脱氧血红蛋白浓度信息和调控氧合血红蛋白浓度信息设置调控设备刺激的幅度控制信号，基于所述幅度控制信号和调控阶段相位生成控制信号的调控时刻，基于所述调控设备刺激的幅值和调控时刻生成调控信号，通过调控设备输出调控信号；
[0087]
在本实施例中，所述调控设备包括经颅磁刺激设备tms或经颅电刺激设备tes。
[0088]
在本实施例中，调控阶段单元包括：配置为通过ns的矩形窗或汉宁窗截取血氧水平信号fnirs，获得最新ns的脱氧血红蛋白浓度信息hbr和氧合血红蛋白浓度信息hbo，基于所述最新最新ns的脱氧血红蛋白浓度信息hbr和氧合血红蛋白浓度信息hbo进行0.01hz
‑
0.2hz的带通滤波并对4个通道脑血氧信息平均处理获得调控脱氧血红蛋白浓度信息hbr
now
和调控氧合血红蛋白浓度信息hbo
now
；
[0089]
通过ms的矩形窗或汉宁窗截取所述脑电信号eeg，获得最新ms的脑电信号，基于所述最新ms的脑电信号进行0.5hz
‑
100hz的带通滤波及50hz的陷波处理，接着进行
和和的带通滤波，再在4个子频段中分别进行快速傅里叶变换和希尔伯特变换计算静息阶段频点和对应的调控阶段幅值和与调控阶段幅值相应的调控阶段相位为与调控阶段幅值相应的调控阶段相位为和
[0090]
基于相对于基线的脑血氧幅度控制信号、相对于基线的脑电幅度控制信号确定触发阈值，基于所述调控阶段相位获取相位控制信号进而获取调控时刻，根据所述触发阈值和调控时刻生成调控信号；
[0091]
所述相对于基线的脑血氧幅度控制信号为所述相对于基线的脑血氧幅度控制信号为其中a和b为待调整参数，取值范围为[
‑
1,1]；在本实施例中，为突出关心运动区的激活，可优选的a取1，b取0,；
[0092]
所述相对于基线的脑电幅度控制信号为所述相对于基线的脑电幅度控制信号为其中c、d、e和f为待调整参数，取值范围为[
‑
1,1]；在本实施例中，由于大脑运动区α节律明显，优选地c、d、f取0，e取1；
[0093]
所述相位控制信号为为其中，g、h、i和j为预设的参数，取值范围为[
‑
1,1]；在本实施例中，由于在大脑运动区α节律明显，优选的g、h、j取0，i取1；
[0094]
所述触发阈值为modu
thr
＝m*hb
modu
+n*eeg
modu
，调控时刻为time＝ang
modu
，其中m和n为待调整参数，取值范围为[
‑
1,1]；在本实施例中，tms对eeg和fnirs信号均有影响，所以m和n均取1；
[0095]
所述调控信号为：
[0096]
其中dev
base
为预设的调控设备刺激强度基准值，p为待调整参数，取值范围为[
‑
1,1]；在本实施例中，tms的dev
base
按照临床常用的静息运动阈值rmt确定，p值按照具体的tms设备的情况而定，优选的取0.5；根据实际观测或治疗需要，ang
modu
可设定为其他角度或角度范围，调控设备可根据所述调控信号输出任意波形。本发明通生成的调控信号融合了脑电信号高时间分辨率的优点和血氧水平信号高空间分辨率的优点使神经调控更精确，能够适应更多的症状提高了神经调控的泛用性；由于存在个体差异，不同的使用者对脑血氧信号表现和脑电信号表现也存在差异，部分使用者脑血氧信号表现良好而脑电信号表现一般或脑血氧信号表现一般而脑电信号表现良好，通过本发明的融合两种模态信号的调控方法能够解决单一模态调控方法难以应对使用者个体差异的问题，即本发明并非将获取脑血氧信号进行调控和获取脑电信号进行调控的简单结合。另外，由于在实际的运用场景中，探测到的脑电信号容易受到干扰，本发明提出的神经调控方法由于其融合了脑血氧信号和脑电
信号提取的信息进行调控，一定程度上地减小了脑电信号受干扰对调控带来的影响，实现了持续的精确调控。
[0097]
所述评价阶段单元，在调控阶段单元完成调控之后，基于所述多模态脑信号，进行解码获取评价脱氧血红蛋白浓度信息、评价氧合血红蛋白浓度信息、各子频段的评价阶段频点和各子频段的评价阶段幅值；
[0098]
基于所述评价脱氧血红蛋白浓度信息、评价氧合血红蛋白浓度信息、静息脱氧血红蛋白浓度信息和静息氧合血红蛋白浓度信息计算血氧调控效果；
[0099]
基于所述各子频段的评价阶段幅值和各子频段的静息阶段幅值计算脑电调控效果；
[0100]
在本实施例中，评价阶段单元包括：配置为对所述脱氧血红蛋白浓度信息hbr、氧合血红蛋白浓度信息hbo进行0.01hz
‑
0.2hz的带通滤波，并进行4个通道的脑血氧信息平均处理获得评价脱氧血红蛋白浓度信息hbr
post
和评价氧合血红蛋白浓度信息hbo
post
；
[0101]
对所述脑电信号eeg进行0.5hz
‑
100hz的带通滤波以及50hz的陷波处理，并进行100hz的带通滤波以及50hz的陷波处理，并进行和的带通滤波，再在4个子频段中分别进行快速傅里叶变换计算评价阶段频点和和对应的评价阶段幅值和
[0102]
血氧调控效果为：其中r和s为待调整参数，取值范围为[
‑
1,1]；在本实施例中，由于关心运动区的激活，优选的r取1，s取0；
[0103]
脑电调控效果为：脑电调控效果为：其中u、v、w和x为待调整参数，取值范围为[
‑
1,1]在本实施例中，由于在大脑运动区α节律明显，优选的u、v、x取0，w取1。
[0104]
所述调控参数优化模块，配置为基于所述血氧调控效果和脑电调控效果调整调控参数；本发明在具体运用过程中，根据观测需要或是实际调节需要，人为地调控参数，以达到更好的控制效果。
[0105]
如图2所示，本发明通过广电同步脑活动检测仪对使用者进行实时脑信号采集，获得fnirs信号和eeg信号，对fnirs信号经过脑信号处理得hbo信号和hbr信号；同时对eeg信号提取子频段相位和子频段幅值；通过hbo信号、hbr信号eeg信号提取子频段相位和子频段幅值获取调控效果评价，基于调控效果评价进行调控参数优化进而实时调控刺激。
[0106]
基于内源性脑信号闭环调控系统应用于经颅磁刺激设备tms上如图6所示。
[0107]
需要说明的是，上述实施例提供的基于内源性脑信号的闭环神经调控的系统，仅以上述各功能模块的划分进行举例说明，在实际应用中，可以根据需要而将上述功能分配由不同的功能模块来完成，即将本发明实施例中的模块或者步骤再分解或者组合，例如，上述实施例的模块可以合并为一个模块，也可以进一步拆分成多个子模块，以完成以上描述
的全部或者部分功能。对于本发明实施例中涉及的模块、步骤的名称，仅仅是为了区分各个模块或者步骤，不视为对本发明的不当限定。
[0108]
本发明第二实施例的基于内源性脑信号的闭环神经调控方法，具体步骤包括：
[0109]
步骤s100，通过脑信号采集设备实时获取多模态脑信号；所述多模态脑信号包括血氧水平信号fnirs和脑电信号eeg；
[0110]
步骤s200，其具体包括步骤s210
‑
步骤s230；
[0111]
步骤s210，基于所述多模态脑信号，进行解码获取静息脱氧血红蛋白浓度信息hbr、静息氧合血红蛋白浓度信息hbo、eeg各子频段的静息阶段频点和eeg各子频段的静息阶段幅值，并进行子频段修正生成个体化子频段；
[0112]
步骤s220，通过基于所述多模态脑信号，进行解码获取最新ns的脱氧血红蛋白浓度信息hbr和氧合血红蛋白浓度信息hbo，通过带通滤波和平均处理的方式获得调控脱氧血红蛋白浓度信息和调控氧合血红蛋白浓度信息；基于所述多模态脑信号中的脑电信号eeg，对各子频段进行个体化子频段的带通滤波、快速傅里叶变换和希尔伯特变换计算各子频段的调控阶段幅值和调控阶段相位；
[0113]
基于所述调控脱氧血红蛋白浓度信息和调控氧合血红蛋白浓度信息设置调控设备刺激的幅度控制信号，基于所述幅度控制信号和调控阶段相位生成控制信号的调控时刻，基于所述调控设备刺激的幅值和调控时刻生成调控信号，通过调控设备输出调控信号；
[0114]
步骤s230，在调控阶段单元完成调控之后，基于所述多模态脑信号，进行解码获取评价脱氧血红蛋白浓度信息、评价氧合血红蛋白浓度信息、各子频段的评价阶段频点和各子频段的评价阶段幅值；
[0115]
基于所述评价脱氧血红蛋白浓度信息、评价氧合血红蛋白浓度信息、静息脱氧血红蛋白浓度信息和静息氧合血红蛋白浓度信息计算血氧调控效果；
[0116]
基于所述各子频段的评价阶段幅值和各子频段的静息阶段幅值计算脑电调控效果；
[0117]
步骤s300，基于所述血氧调控效果和脑电调控效果调整调控参数。
[0118]
所属技术领域的技术人员可以清楚地了解到，为描述的方便和简洁，上述描述的方法的具体工作过程及有关说明，可以参考上述系统实施例中的对应功能模块，在此不再赘述。
[0119]
本发明第三实施例的一种电子设备，其特征在于，包括：至少一个处理器；以及与至少一个所述处理器通信连接的存储器；其中，所述存储器存储有可被所述处理器执行的指令，所述指令用于被所述处理器执行以实现上述的基于内源性脑信号的闭环神经调控方法。
[0120]
本发明第四实施例的一种计算机可读存储介质，其特征在于，所述计算机可读存储介质存储有计算机指令，所述计算机指令用于被所述计算机执行以实现上述的基于内源性脑信号的闭环神经调控方法。
[0121]
所属技术领域的技术人员可以清楚地了解到，为描述的方便和简洁，上述描述的存储装置、处理装置的具体工作过程及有关说明，可以参考前述方法实施例中的对应过程，在此不再赘述。
[0122]
术语“第一”、“第二”等是用于区别类似的对象，而不是用于描述或表示特定的顺
序或先后次序。
[0123]
术语“包括”或者任何其它类似用语旨在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备/装置不仅包括那些要素，而且还包括没有明确列出的其它要素，或者还包括这些过程、方法、物品或者设备/装置所固有的要素。
[0124]
至此，已经结合附图所示的优选实施方式描述了本发明的技术方案，但是，本领域技术人员容易理解的是，本发明的保护范围显然不局限于这些具体实施方式。在不偏离本发明的原理的前提下，本领域技术人员可以对相关技术特征做出等同的更改或替换，这些更改或替换之后的技术方案都将落入本发明的保护范围之内。