本文研究了一类基于粒子集落刺激因子管理嗜中性粒子的血液病振荡数学模型。通常造血系统内部过程由激发分化机制和成熟机制两部分组成。造血干细胞分化得到了动物和人类必须的三类血液细胞:白细胞、红细胞和血小板。文中应用DDE-Biftool软件,数值模拟了造血系统中的长周期解的周期解分叉,得到了周期解的倍周期分叉。倍周期分叉后系统出现混沌解,应用反馈控制可以消除混沌或把系统稳定在期望的轨道上。 The oscillating phenomenon of a hematological cell model is investigated. Based on effective G-SCF administration, the mathematical blood cell model is described as DDEs with multi-delays. Com-monly, in blood cell models, the intrinsic mechanism is composed of triggering differentiating mechanism and maturation mechanism, and hematoplogical cell model differentiating into three types of necessary blood cells in the body. By applying numerical software DDE-Biftool, the bifur-cation of periodic solutions with long period are derived, and the period doubling bifurcation of periodic solutions is found at critical values. Chaotic solutions also appear afterwards period- doubling bifurcation, and the stabilization of dynamical chaos to the expected periodic solution is finished successfully by applying Pyragas feedback control method.
马苏奇
中国农业大学数学系,北京
收稿日期:2016年8月11日;录用日期:2016年8月25日;发布日期:2016年8月31日
本文研究了一类基于粒子集落刺激因子管理嗜中性粒子的血液病振荡数学模型。通常造血系统内部过程由激发分化机制和成熟机制两部分组成。造血干细胞分化得到了动物和人类必须的三类血液细胞:白细胞、红细胞和血小板。文中应用DDE-Biftool软件,数值模拟了造血系统中的长周期解的周期解分叉,得到了周期解的倍周期分叉。倍周期分叉后系统出现混沌解,应用反馈控制可以消除混沌或把系统稳定在期望的轨道上。
关键词 :造血系统,周期解,倍周期分叉,时滞
理解和研究造血系统对人类健康是非常必要的。造血干细胞分化能够得到动物和人类必须的三类血液细胞:白细胞、红细胞和血小板。而造血作用的失调会导致很多的血液学方面的疾病,如典型的例子有慢性的嗜中性白血球减少症。在病人中,一个或多个类型的血液循环细胞中出现明显的振荡。因而非线性动力学是研究血液病动力学的一个很好的工具。
通常造血系统由激发分化和增殖的机制以及成熟机制组成。从建模起,人们习惯于用Holling函数表达造血干细胞分化成为嗜中性粒子和自我更新的干细胞的数量。和组织细胞的相互作用中,细胞需要经历的成熟期一般为几天(2到6天),增殖放大的系数用指数形式表示。作为造血系统的一个极为明显的特征就是系统中引入了滞量,由于细胞的分化和成熟都需要一定的天数,所以系统是一多时滞模型,人们大都发展了包括造血干细胞房和嗜中性粒子细胞房在内的两房非线性反馈调控系统 [
实验里,在人类血液病——周期性嗜中性粒子(CN)中,以振荡周期为19~21天几乎不能检测到嗜中性粒子为特征。在癌症治疗中,医学上一般使用集粒落刺激因子(G-CSF)来管理CN,其一般原理是降低嗜中性粒子前体细胞的程序性凋亡率 [
本文应用DDE-Biftool数学软件,研究了两房血液病模型的周期解分叉。模型中发现三类典型的长周期解,利用DDE-Biftool软件,分别追踪得到周期从46天左右演化为186天的倍周期解分叉,和周期为79天和周期为159天的吸引子的连续延拓。倍周期解是通往混沌的路径之一。利用Pyragas方法研究了混沌的反馈控制,发现负反馈增益可以使系统稳定在不同的周期轨道上,并且增益的连续改变带来了混沌系统的倍周期分叉现象。正反馈增益使得系统稳定在期望的吸引态。
造血过程由激发分化机制和成熟机制两部分组成,具体地,造血干细胞经由细胞分裂分化成为全部血液细胞包括白细胞、红细胞和血小板,并能自我更新,形成了外周血液细胞的反馈调控环。基于G-CSF管理CN,考虑如下造血系统的两房模型
其中Holling函数
其生物学意义如下:干细胞以速率
选取
应用DDE-Biftool工具和DDE23,可以高效跟踪周期解随参数
图2(a)中1分支参数
让参数
混沌控制问题的研究已经被广泛关注,设计Pyragas延时控制回路,得到如下的反馈控制系统
图1. 系统(1)的三类长周期解,
图2. 周期解随参数
图3. 两种类型的吸引子。(a) 周期解分支1的全部吸引子;(b) 周期解分支2的全部吸引子
图4. 变化参数
选取不同的
在Pyragas的延时反馈控制方法中增加相控制器,采用相空间压缩法,
则可以将系统完全稳定到期望的轨道上,见图6(c),选取
图5. 负反馈增益将混沌控制在稳定的周期轨道。(a)
图6. 相空间压缩法将系统稳定在期望的周期轨道上。(a) 由(3)得到的Poincare截面上的控制效果;(b) 相应于截面(a),
造血系统中存在很丰富的周期动力学现象。在实验中,Harrie等人发现了嗜中性粒子的振荡周期可以分布在一个很广泛的范围内。本文研究了三类共存的长周期的周期解的连续分支,发现倍周期分叉是通往混沌的路径之一。利用Pyragas的延时反馈控制方法,负反馈增益可将系统稳定嵌入到倍周期分叉的不同周期轨道上。在Pyragas的延时反馈控制方法中增加相控制器,正反馈增益可将系统稳定到期望的吸引态。
马苏奇. 造血系统中的周期解分岔及反馈控制Bifurcation of Periodic Solution of a Hematopological System and Feedback Control[J]. 应用数学进展, 2016, 05(03): 487-493. http://dx.doi.org/10.12677/AAM.2016.53059