一种基于图学习优化算法的动态脑功能网络学习方法

本发明涉及图机器学习领域，具体涉及一种基于图学习优化算法的动态脑功能网络学习方法。

背景技术：

1、动态脑网络研究是神经科学和认知科学领域的一个重要方向，它关注人脑在不同时间点、不同认知任务或不同认知状态下的神经网络变化。这种研究有助于更深入地理解脑网络在认知、学习、记忆和神经疾病等方面的功能，并且可以为精准医疗、脑机接口和认知增强等领域提供关键信息。

2、通过功能性磁共振成像(fmri)和脑电图(eeg)等脑成像技术可以记录脑活动的变化，可以帮助研究动态脑网络。fmri可以提供高空间分辨率的脑区活动信息，而eeg可以提供高时间分辨率的脑电活动数据。通过图机器学习以及网络分析方法对脑电数据进行分析可以帮助研究脑网络的拓扑结构和连接模式。这些方法可以揭示不同时间点或任务下神经网络之间的交互关系，以及脑节点之间的信息传递方式。

3、构建动态脑网络连接是对动态脑网络分析的基础也是第一步。传统的动态脑网络功能连接构建方法通常是基于滑动窗口获取fmri信号片段，再通过计算各脑区脑节点fmri信号的皮尔森系数来获得。这种方法受噪声影响较大且忽略了fmri信号动态变化的时间一致性。

技术实现思路

1、针对现有技术中存在的不足，本发明提供一种基于图学习优化算法的动态脑功能网络学习方法。

2、本发明方法可以从fmri信号中直接求取大脑的动态功能连接网络。动态脑功能连接网络能更加合理准确的反映出脑区之间的连接强度和脑网络的拓扑结构。

3、一种基于图学习优化算法的动态脑功能网络学习方法，包括步骤如下：

4、步骤(1)确定脑网络的基本数学模型；

5、步骤(2)预处理真实的神经影像数据；

6、步骤(3)构建求取动态脑功能连接网络的目标函数；

7、步骤(4)确定求取动态脑功能连接网络的求解方法。

8、步骤(1)具体方法如下：

9、用图论的方法表示脑网络是神经科学中的一个常见方法，可以将神经元和它们之间的连接关系以图的形式进行表达。具体表述为：

10、脑节点表示神经元：在图中，每个神经元或者脑区都被表示为一个脑节点(vertex)。每个脑节点代表一个神经元或神经元群体，这取决于所研究的尺度。

11、边表示连接：神经元之间的连接能够用边(edge)来表示。如果两个神经元或者神经元群体之间存在突触连接或功能性连接，就在相应的脑节点之间画一条边。边可以具有方向性，以表示连接的传递方向。

12、权重表示连接强度：边上附加权重，表示连接的强度或权重。在脑网络中，连接的强度与突触传递的效率或其他生理参数有关。

13、图的类型：根据问题的不同，使用不同类型的图。常见的图类型包括有向图(directed graph)和无向图(undirected graph)。有向图适用于描述单向的连接，而无向图则用于表示双向的连接或没有方向性的连接。

14、将大脑分成不同的脑区，用脑节点来表示大脑不同的脑区，脑节点集表示为：v＝{v1,v2,v3,...,vn}，以加权邻接矩阵的形式来表示脑功能连接网络，其中wij表示脑节点vi和vj之间的功能连接权重值，该值代表脑区之间功能连接的强度。

15、脑网络的每个脑节点上的信号为：x∈rn×m＝[x1,...,xn]t，其中x的每行xi∈rn代表脑节点vi上的脑电信号，本发明中选择的脑电信号是bold信号。

16、步骤(2)具体操作如下：

17、fmri是一种非侵入性的脑成像技术，通过测量血氧水平的变化来反映大脑活动。将fmri信号转化为血氧信号的步骤如下：fmri扫描会生成一系列三维体素(体积像素)，原始的fmri数据需要进行一些预处理步骤，例如头动校正、时间对准、去除噪声等。这些步骤有助于提高信号质量和准确性。

18、通过计算预处理后的数据中每个体素的信号强度，获得一个体素时间序列，其中每个时间点对应于一个三维体素的状态。将估计的血氧(bold)信号进行校正和标准化，以便在不同受试者之间进行比较和分析，校正处理后的血氧信号可以进行各种数据分析，如功能连接分析、脑区活动对比分析等，从而更深入地理解脑的功能。

19、利用滑动窗口的方法将处理后的血氧信号等分成多个相同尺寸的信号，具体操作为：设定尺寸为m×t的滑动窗口，其中m为roi脑区的数量，t为采样点的个数。设定窗口滑动步长为d个采样时间点。利用滑动窗口对bold信号进行数据分割，即可得到m×t×n的三维矩阵st∈rm×t。

20、步骤(3)具体操作如下：

21、对于高强度功能连接的脑节点对(vi,vj)，脑节点对上的bold信号呈现相似的特征。因为脑节点对上的信号能够视为一对向量，如式(1)所示，通过计算两个向量之间的欧氏距离的平方来衡量脑节点对上的bold信号相似性：

22、

23、其中，d(vi,vj)为成对距离，xi和xj分别表示脑节点vi和vj上的bold信号。计算所有脑节点对的成对距离并且乘脑节点对的功能连接权重值，如式(2)所示：

24、

25、其中x为滑动窗口所选的bold信号。

26、进一步地，将求解脑功能连接网络的拉普拉斯矩阵转化为求解如下最优化问题，如式(3)所示：

27、

28、由于式(3)优化所得的是拉普拉斯矩阵，从拉普拉斯矩阵中求解得到脑功能连接矩阵较为困难，因此改写式(3)，建立求解脑功能连接矩阵的目标函数,具体操作如下：

29、定义式(2)可等效为式(4):

30、

31、进一步地，如式(5)所示，求解脑功能连接网络转化为求解如下最优化问题：

32、

33、为了使每个脑节点至少与另一个脑节点保持连接状态，并且因为脑网络具有小世界属性，因此所求解的脑功能连接网络w应该是稀疏的，因此对目标函数式(5)增加约束项表示为式(6):

34、

35、其中1表示元素全为1的向量，即1＝[1,...,1]t,对数约束项log(w1)使w中脑节点的度为正数，但不会使w中的元素变为0，从而避免w出现平凡解。添加f范数约束项来控制图的稀疏性，参数α和β为约束性的权重参数，参数α作用为使w中的边保持较多的连接，参数β的作用为控制w的稀疏性，参数β的值越大，所求解的w越稀疏。

36、步骤(4)具体操作如下：

37、对于式(6)构建的目标函数，求解的具体步骤为：

38、对每个窗口的bold信号，通过式(6)求解每个窗口的脑功能连接网络w，即可得到动态脑功能连接网络由于所求解的是在wm空间，wm空间能够描述为：在采用具体的优化算法对式(6)的最优化问题进行求解时，使用双线性映射法将wm空间转化为wv空间，wv空间能够描述为：在wv空间优化会更加容易。例如：在wv空间中，可以不必处理功能连接矩阵w的对称性问题，从而降低优化难度。

39、因此利用空间转换，将式(6)改写为式(7):

40、

41、其中k是满足w1＝kw的线性算子。进一步地，定义：

42、f1(w)＝i{w≥0}+2wtz

43、f2(d)＝-α1tlog(d)

44、f3(w)＝β||w||2

45、其中i{}是指示函数，当满足括号中条件时，值为零，反之，值为无穷大，其中d＝kw∈rm。然后采用原始对偶算法对式(7)进行求解得到

46、本发明有益效果如下：

47、本发明首先对fmri数据处理得到bold信号，并使用滑动窗口将bold信号分成多个时间段；然后，设计基于图学习的优化算法，优化得到动态的大脑功能连接矩阵，即动态图学习。

48、动态大脑功能连接网络是一种衡量不同脑区之间协调性的方法，它研究脑区之间的信息传递和交互模式，并且可以应用于许多神经科学领域，研究意义包括但不局限于如下几点：

49、1.理解大脑的工作方式：动态脑功能连接网络可以帮助我们深入理解大脑的工作方式，包括不同区域之间的相互作用、信息处理和传递等。

50、2.探索神经系统疾病：研究发现，精神疾病的发病机制可能涉及到脑网络功能连接矩阵的改变。例如，精神分裂症患者可能存在前额叶-杏仁体通路的连接异常，抑郁症患者可能存在扣带回-杏仁体通路的连接减弱。因此本发明中所提出的动态图学习方法能够及时并且准确的表现大脑功能连接网络的演变过程，以及能够代表大脑拓扑结构的拉普拉斯矩阵，直接对诸多疾病的发病机制以及预防及治疗都有重要的作用。

51、总之，动态脑功能连接网络的研究对理解大脑的复杂性以及神经系统疾病的机制和治疗方法具有重要意义。