一种基于钙浓度的神经突触可塑性计算方法与流程

本发明涉及大脑模拟仿真领域，尤其是针对大脑模拟仿真中脉冲神经网络的神经突触可塑性的计算问题。

背景技术：

大脑中包含数以亿万计的神经元，这些神经元相互连接构成神经网络，神经元连接的部分称为突触。以人脑为例，每个神经元会与几千个至几万个神经元产生连接，因此，大脑中的连接数量十分庞大。感知和认知过程会改变突触连接的强弱，也就是权重。这种突触连接权重变化的现象称为突触可塑性，是大脑学习和记忆的主要方式。如何准确地计算这些突触权重的变化过程是进行大脑模拟仿真的关键。最接近现有技术的神经突触可塑性计算方法为graupner等人在calcium-basedplasticitymodelexplainssensitivityofsynapticchangestospikepattern,rate,anddendriticlocation一文中描述的方法，其问题在于，并没有考虑神经元内蛋白质浓度和树突棘状态对突触可塑性的影响，导致了计算不准确等问题。

技术实现要素：

本发明所要解决的技术问题是在大脑模拟中突触可塑性模型和计算问题。

为了解决这一问题，本发明采用的技术方案分为以下五个步骤：

1.钙离子浓度的计算：针对需要计算的突触，根据初始时刻t0的突触前神经元和突触后神经元的膜电位值，以及突触的初始连接权重w0，分别计算下一时刻t1的树突和树突棘中的钙离子浓度，其中w0是非负常数，在实际应用中，建议取值为1，原因是在初始时刻做归一化处理；

2.确定权重变化方向(flag)：根据树突棘中的钙离子浓度，通过和钙离子浓度阈值ca0s及ca1s对比得到权重变化的方向，其中ca0s和ca1s均为非负常数，且ca1s>ca0s，在实际应用中，根据经验，建议ca0s＝0.5mm,ca1s＝1mm；

3.计算突触状态标签(tag)：根据权重变化方向及突触状态变化常数αt，βt,ltd和βt,ltp，计算突触标签状态，其中αt，βt,ltd和βt,ltp均为大于零的常数，在实际应用中，根据经验，建议αt＝0.01,βt,ltd＝0.5,βt,ltp＝0.5；

4.计算可塑性相关蛋白质(prp)浓度：根据树突中的钙离子浓度计算可塑性相关蛋白质浓度；

5.计算突触权重的变化：根据突触状态标签tag和可塑性相关蛋白质prp浓度，计算突触权重的变化，得到t1时刻的新权重。

重复以上过程，计算模拟时间内的突触连接强度。

本发明的内容可以应用于构建类脑神经网络，完成类脑智能所需的学习和记忆过程的仿真，实现通用强人工智能，并应用于智能传媒、智能医疗、智能制造、智能教育和自动驾驶等方面。

根据本发明的第一个步骤，采用如下方法：

首先，设定突触前神经元和突触后神经元的初始膜电位值以及突触前神经元和突触后神经元之间的初始连接权重w0。其中突触前神经元和突触后神经元的初始膜电位值为小于零的常数，在实际应用中，建议设置为-70mv，原因是这是神经元的静息态膜电位，突触前神经元和突触后神经元之间的初始连接权重w0是非负常数，在实际应用中，建议取值为1，原因是在初始时刻做归一化处理。

其次，采用脉冲神经元模型对突触前神经元和突触后神经元进行建模，其中脉冲神经元模型来自于神经科学领域文献，模型的输入为各离子通道的电导、各离子通道的反转电位、输入电流和膜电容，以上输入数据均来自于神经科学领域文献，模型的输出为膜电位值。

再次，根据突触前神经元的脉冲神经元模型，利用欧拉法或龙格库塔法等数值计算方法计算t1时刻突触前神经元的膜电位。

然后，根据突触动力学模型对突触连接进行建模，其中突触动力学模型来自于神经科学领域文献，模型的输入为t1时刻突触前神经元膜电位，t1时刻突触后神经元膜电位，突触动力学时间常数以及突触受体的反转电位，突触动力学时间常数以及突触受体的反转电位均来自于神经科学领域参考文献，模型的输出为突触连接对突触后神经元产生的电流。

最后，根据突触动力学模型及突触后神经元的神经元模型，利用欧拉法或龙格库塔法等数值计算方法计算t1时刻突触后神经元的膜电位以及突触后神经元树突和树突棘中的钙离子浓度。突触后神经元树突中的钙离子源自位于细胞膜上的钙离子通道，突触后神经元树突棘中的钙离子源自位于树突棘上允许钙离子通过的突触受体。在t1时刻由钙离子通道进入神经元内的钙离子浓度[ca²⁺]ca-channel可以由包含该离子通道的脉冲神经元模型计算得出，在t1时刻由树突棘上允许钙离子通过的突触受体进入神经元内的钙离子浓度[ca²⁺]receptor可以由该突触的动力学模型计算得出。突触后神经元树突和树突棘中的钙离子可以相互扩散，记t1时刻树突棘和该树突棘所处的树突中的钙离子浓度分别为[ca²⁺]s和[ca²⁺]d。

根据本发明的第二个步骤，采用如下方法：

根据树突棘中的钙离子浓度[ca²⁺]s和钙离子浓度阈值ca0s及ca1s对比得到t1时刻该树突棘所处的突触的权重变化方向flag。其中ca0s和ca1s均为非负常数，且ca1s>ca0s，在实际应用中，根据经验，建议ca0s＝0.5mm,ca1s＝1mm。如果t1时刻树突棘中的钙离子浓度[ca²⁺]s低于阈值ca0s，则flag为0，该树突棘所处的突触不会发生连接权重变化；如果t1时刻树突棘中的钙离子浓度[ca²⁺]s处于阈值ca0s和阈值ca1s之间，则flag为-1，该树突棘所处的突触连接权重变化方向为减弱；如果t1时刻树突棘中的钙离子浓度[ca²⁺]s高于阈值ca1s，则flag为1，该树突棘所处的突触的连接权重变化方向为增加。

根据本发明的第三个方面，采用如下方法：

根据树突棘所处突触权重变化方向及突触状态变化常数αt，βt,ltd和βt,ltp得到t1时刻该树突棘所处的突触状态标签tag。如果tag为0，则该树突棘所处的突触不会发生连接权重变化；如果tag大于0，则该树突棘所处的突触进入连接权重可以增加的状态；如果tag小于0，则该树突棘所处的突触进入连接权重可以减弱的状态。tag的演化是一个与flag有关的动力学过程。

如果t1时刻flag为0，则其中αt是大于零的常数，在实际应用中，根据经验，建议αt＝0.01，若tag以t0时刻初始值为0开始演化，则tag关于时间的导数为0，tag不发生变化；

如果t1时刻flag为-1，则其中βt,ltd是大于零的常数，在实际应用中，根据经验，建议βt,ltd＝0.5，若tag以t0时刻初始值为0开始演化，则tag关于时间的导数为负，tag降低，直至降到tag关于时间的导数为0时，即时停止下降，并保持

如果t1时刻flag为1，则其中βt,ltp是大于零的常数，，在实际应用中，根据经验，建议βt,ltp＝0.5，若tag以t0时刻初始值为0开始演化，则tag关于时间的导数为正，tag升高，直至升至tag关于时间的导数为0时，即时停止上升，并保持

根据本发明的第四个步骤，采用如下方法：

根据树突中的钙离子浓度[ca²⁺]d和钙离子浓度阈值ca0d得到t1时刻树突中prp的浓度。prp浓度的演化是与[ca²⁺]d有关的动力学过程。其中ca0d为非负常数，在实际应用中，根据经验，建议ca0d＝1mm。树突中的钙离子可以引发生成prp的化学反应，单位时间内prp的生成率记为prprate。如果t1时刻树突中的钙离子浓度[ca²⁺]d低于阈值ca0d，则prprate为0或衰减为0，prp不发生变化；如果t1时刻树突中的钙离子浓度[ca²⁺]d高于阈值ca0d，则prprate按alpha函数的形式变化，prp作为prprate对时间的积分，先增加，后持平。

根据本发明的第五个步骤，采用如下方法：

首先，计算突触连接权重系数相关参数y，y的演化是与tag和prp浓度有关的动力学过程，其中τy是大于零的常数，在实际应用中，根据经验，τy＝1ms。由于prp浓度是非负的，若tag和prp浓度中有一个为0，则y对时间的导数为0，y不会发生变化；若tag为负，prp浓度为正，则y对时间的导数为负，y降低；若tag为正，prp浓度为正，则y对时间的导数为正，y升高。

其次，计算突触连接权重系数z，z与y之间的关系满足其中λ，zl，zh均为大于零的常数，zl代表权重系数z的下限，zh代表权重系数z的上限，在实际应用中，根据经验，λ＝1,zl＝0.5,zh＝2。

再次，将突触连接权重系数z乘以初始连接权重w0，得到t1时刻的连接权重w1＝z×w0。

最后，重复本发明的第一个步骤、第二个步骤、第三个步骤、第四个步骤和第五个步骤，计算整个模拟时间内神经元之间的突触连接权重变化。

本发明的有益效果是大大提高了大脑模拟中突触可塑性的计算精度，并可以应用于构建类脑神经网络，完成类脑智能所需的学习和记忆过程的仿真，实现通用强人工智能，并应用于智能传媒、智能医疗、智能制造、智能教育和自动驾驶等方面。

附图说明

图1为突触后神经元树突和树突棘示意图

图2为权重变化方向flag和突触态标签tag的计算流程图

图3为prp生成率和prp浓度的计算流程图

图4为突触连接权重系数相关参数y和突触连接权重系数z的计算流程图

具体实施方式

下面将结合本发明中的附图和具体实施方式对本发明的技术方式进行详细描述。

在钙离子浓度的计算中：

首先，设定突触前神经元和突触后神经元的初始膜电位值以及突触前神经元和突触后神经元之间的初始连接权重w0。其中突触前神经元和突触后神经元的初始膜电位值为小于零的常数，在实际应用中，建议设置为-70mv，原因是这是神经元的静息态膜电位，突触前神经元和突触后神经元之间的初始连接权重w0是非负常数，在实际应用中，建议取值为1，原因是在初始时刻做归一化处理。

其次，采用脉冲神经元模型对突触前神经元和突触后神经元进行建模，其中脉冲神经元模型来自于神经科学领域文献，模型的输入为各离子通道的电导、各离子通道的反转电位、输入电流和膜电容，以上输入数据均来自于神经科学领域文献，模型的输出为膜电位值。

再次，根据突触前神经元的脉冲神经元模型，利用欧拉法或龙格库塔法等数值计算方法计算t1时刻突触前神经元的膜电位。

然后，根据突触动力学模型对突触连接进行建模，其中突触动力学模型来自于神经科学领域文献，模型的输入为t1时刻突触前神经元膜电位，t1时刻突触后神经元膜电位，突触动力学时间常数以及突触受体的反转电位，模型的输出为突触连接对突触后神经元产生的电流。

最后，根据突触动力学模型及突触后神经元的神经元模型，利用欧拉法或龙格库塔法等数值计算方法计算t1时刻突触后神经元的膜电位以及突触后神经元树突和树突棘中的钙离子浓度。图1为突触后神经元树突和树突棘示意图，突触后神经元树突中的钙离子源自位于细胞膜上的钙离子通道，突触后神经元树突棘中的钙离子源自位于树突棘上允许钙离子通过的突触受体。在t1时刻由钙离子通道进入神经元内的钙离子浓度[ca²⁺]ca-channel可以由包含该离子通道的离子通道型神经元模型计算得出，在t1时刻由树突棘上允许钙离子通过的突触受体进入神经元内的钙离子浓度[ca²⁺]receptor可以由该受体的动力学模型计算得出。突触后神经元树突和树突棘中的钙离子可以相互扩散，记t1时刻树突棘和该树突棘所处的树突中的钙离子浓度分别为[ca²⁺]s和[ca²⁺]d。

在突触的权重变化方向flag和突触态标签tag的计算中：

首先，根据树突棘中的钙离子浓度[ca²⁺]s和钙离子浓度阈值ca0s及ca1s对比得到t1时刻该树突棘所处的突触的权重变化方向flag。其中ca0s和ca1s均为非负常数，且ca1s>ca0s，在实际应用中，根据经验，建议ca0s＝0.5mm,ca1s＝1mm。如果t1时刻树突棘中的钙离子浓度[ca²⁺]s低于阈值ca0s，则flag为0，该树突棘所处的突触不会发生连接权重变化；如果t1时刻树突棘中的钙离子浓度[ca²⁺]s处于阈值ca0s和阈值ca1s之间，则flag为-1，该树突棘所处的突触连接权重变化方向为减弱；如果t1时刻树突棘中的钙离子浓度[ca²⁺]s高于阈值ca1s，则flag为1，该树突棘所处的突触的连接权重变化方向为增加。

其次，根据树突棘所处突触权重变化方向及突触状态变化常数αt，βt,ltd和βt,ltp得到t1时刻该树突棘所处的突触状态标签tag。如果tag为0，则该树突棘所处的突触不会发生连接权重变化；如果tag大于0，则该树突棘所处的突触进入连接权重可以增加的状态；如果tag小于0，则该树突棘所处的突触进入连接权重可以减弱的状态。tag的演化是一个与flag有关的动力学过程。如果t1时刻flag为0，则其中αt是大于零的常数，在实际应用中，根据经验，建议αt＝0.01，若tag以t0时刻初始值为0开始演化，则tag关于时间的导数为0，tag不发生变化；如果t1时刻flag为-1，则其中βt,ltd是大于零的常数，在实际应用中，根据经验，建议βt,ltd＝0.5，若tag以t0时刻初始值为0开始演化，则tag关于时间的导数为负，tag降低，直至降到tag关于时间的导数为0时，即时停止下降，并保持如果t1时刻flag为1，则其中βt,ltp是大于零的常数，，在实际应用中，根据经验，建议βt,ltp＝0.5，若tag以t0时刻初始值为0开始演化，则tag关于时间的导数为正，tag升高，直至升至tag关于时间的导数为0时，即时停止上升，并保持

图2为权重变化方向flag和突触态标签tag的计算流程图。如果t1时刻树突棘中的钙离子浓度[ca²⁺]s低于阈值ca0s，则flag为0，若tag以t0时刻初始值为0开始演化，则tag不发生变化，仍然为0；如果树突棘中的钙离子浓度[ca²⁺]s处于阈值ca0s和阈值ca1s之间，则flag为-1，若tag以t0时刻初始值为0开始演化，则tag降低，直至降到时停止下降，并保持如果树突棘中的钙离子浓度[ca²⁺]s高于阈值ca1s，贼flag为1，若tag以t0时刻初始值为0开始演化，则tag升高，直至升至时停止上升，并保持在该计算流程示例中，ca0s＝0.5mm,ca1s＝1mm,αt＝0.01,βt,ltd＝0.5,βt,ltp＝0.5。

在prp生成率和prp的计算中：

根据树突中的钙离子浓度[ca²⁺]d和钙离子浓度阈值ca0d得到t1时刻树突中prp的浓度。prp浓度的演化是与[ca²⁺]d有关的动力学过程。其中ca0d为非负常数，在实际应用中，根据经验，建议ca0d＝1mm。树突中的钙离子可以引发生成prp的化学反应，单位时间内prp的生成率记为prprate。图3为prp生成率和prp的计算流程图。如果t1时刻树突中的钙离子浓度[ca²⁺]d低于阈值ca0d，若prprate以t0时刻初始值为0开始演化，则prprate不发生变化，仍然为0，prp是prprate对时间的积分，因此，prp的浓度为0；如果t1时刻树突中的钙离子浓度[ca²⁺]d高于阈值ca0d，若prprate以t0时刻初始值为0开始演化，则prprate上升，prp是prprate对时间的积分，因此，prp的浓度增加。在该计算流程示例中，ca0d＝1mm。

在突触连接权重系数相关参数y和突触连接权重系数z的计算中：

首先，计算突触连接权重系数相关参数y，y的演化是与tag和prp浓度有关的动力学过程，其中τy是大于零的常数，在实际应用中，根据经验，τy＝1ms。由于prp浓度是非负的，若tag和prp浓度中有一个为0，则y对时间的导数为0，y不会发生变化；若tag为负，prp浓度为正，则y对时间的导数为负，y降低；若tag为正，prp浓度为正，则y对时间的导数为正，y升高。

其次，计算突触连接权重系数z，z与y之间的关系满足其中λ，zl，zh均为大于零的常数，zl代表权重系数z的下限，zh代表权重系数z的上限，在实际应用中，根据经验，λ＝1,zl＝0.5,zh＝2。图4为突触连接权重系数相关参数y和突触连接权重系数z的计算流程图。由于prp浓度是非负的，若tag和prp浓度中有一个为0，则y不会发生变化；若tag为负，prp浓度为正，则y降低；若tag为正，prp浓度为正，则y升高。由y可以计算出z。在该计算流程示例中，τy＝1ms,λ＝1,zl＝0.5,zh＝2。

最后，将突触连接权重系数z乘以初始连接权重w0，得到t1时刻的连接权重w1＝z×w0。