1. 引言
近年来,用非线性发展方程刻画模型已引起了人们的广泛关注,并且成功地解决了一些模型的性态演化问题。如在文献 [1],利用CRE方法研究了两类非线性方程得到了新的相互作用解,并且描述了非线性发展方程在数学物理中的重要性,而在合适的假设下利用微分方程刻画肿瘤细胞的增长过程以及影响因素,也为肿瘤治疗提供了新的思路和方法 [2]。
在许多专家学者的研究中,通过采用指数函数法、齐次平衡法、双曲正切法和F-展开法来求解非线性偏微分方程,得到了丰富的精确解。在文献 [3] [4] 中Rotenberg利用Fokker-Planck逼近给出了一个描述细胞增长的线性扩散型方程。本文在上述文献的启示下研究一个非线性的偏微分方程 [5]
(1)
并运用F-展开法构造方程(1)的一个精确解,其中
是时间,
是成熟速度,
是成熟度。
本文的结构如下:在第2节,我们概述F-展开法的主要思路,然后利用Mathematica软件绘制图形,得到方程的精确解并且描述了方程模型。第3节对本文内容进行了简单的总结。
2. 求解非线性偏微分方程
2.1. 方法概述
求解非线性偏微分方程的F-展开法可以看作是齐次平衡原则的新应用,为了简单起见,本文以两个自变量的情形为例来说明F-展开法的主要思路 [6] [7]。
给定非线性偏微分方程:
(2)
P为其变元的多项式,其中含有非线性项和高阶偏导数项。
第一步:为了将偏微分转化为微分方程,我们令
(3)
其中
为待定常数。将(3)式代入方程(2),则(2)式变为常微分方程
(4)
第二步:用齐次平衡法原则确定N。
第三步:设
为
的洛朗级数形式
(5)
其中
是待定常数,
满足Jacobi椭圆方程
(6)
其中
是常数。具体的解的分类可见 [8]
第四步:将式(5)代入方程(4),利用方程(6)可将方程(4)的左端变成
的多项式,合并
多项式。
第五步:令
的各次幂项的系数为零,得
的含有
的代数方程组,利用符号计算软件Mathematica解出
,用
表示的
。
第六步:将
的结果代入式子(5),得方程的行波解的一般形式,再由
与
关系代入
的值得出由雅各比椭圆函数表示的解。
2.2. 具体步骤
本节借助以上思路来构造非线性方程(1)的精确解。考虑如下形式的变换
(7)
其中
或者
。将变换(7)代入非线性方程(1),得到关于
的常微分方程
(8)
这里考虑一个特殊情形令
,得到方程
(9)
利用F-展开法的思路,利用方程(9)中的非线性项与最高阶导数项
与
,齐次平衡后就可确定
,则方程(8)的解具有形式
(10)
将(10)式代入方程(9)中,同时利用(6)式,合并F的同次幂项的系数,则方程(9)的左端
归结为
的8次多项式:
令
的各次幂项的系数为零,得到关于
代数方程组,解此代数方程组可得一组解:
特别的,利用方程(6)的特解,我们可得
将
与上面的解代入(10)式中可得
(11)
最后将(11)式代入(7)式中可得到方程(1)的精确解为
(12)
和
(13)
2.3. 图形分析
对于精确解(12)当分别取
,
时,利用Mathematica软件 [9] 绘制解随时间的演化性态。
1)
在
时解的演化过程见图1,此时有一个波峰。
2)
在
时解的演化过程见图2,此时有两个波峰。
3)
在
时解的演化过程见图3,此时有三个波峰。
Figure 1. When
, the solution evolves with v and
图1. 当
时解随
的演化过程
Figure 2. When
, the solution evolves with v and
图2. 当
时解随
的演化过程
Figure 3. When
, the solution evolves with v and
图3. 当
时解随
的演化过程
对精确解(13)当分别取
,
时利用Mathematica软件绘制解随时间演化的图形。
1)
在
时解的演化过程见图4,此时只有一个波峰。
2)
在
时解的演化过程如图5,此时只有一个波峰。
3)
在
时解的演化过程见图6,此时有两个波峰。
Figure 4. When
, the solution evolves with v and
图4. 当
时解随
的演化过程
综上:当
时,由解(12)随时间的的演化图形知,随着时间的增大,波峰个数变为一个、
两个、三个。
当
时,由解(13)随时间的演化图形知,随着时间的增大,波峰个数逐渐变为一个、一个、
两个,而造成这两种不同的原因可能与参数m的取值有关。
Figure 5. When
, the solution evolves with v and
图5. 当
时解随
的演化过程
Figure 6. When
,the solution evolves with v and
图6. 当
时解随
的演化过程
3. 结论
本文利用F-展开法 [10],并借助Mathematica软件,研究了描述肿瘤细胞增长模型的Fokker-Planck方程,得到了该方程的两个精确解,通过参数取值绘制解的图形,直观简洁地展示了解的演化性态。通过设定不同的待定函数,采用相同的思路还可用该方法求得Fokker-Planck方程更多的精确解。
致谢
感谢中央民族大学本科生研究训练计划(URTP)——BEIJ2020110002的支持和张智勇导师的指导和帮助。