没有合适的资源?快使用搜索试试~ 我知道了~
资源推荐
资源详情
资源评论
一 类 干 扰 血 液 模 型 的 Hopf 分 支
倡
李 娜, 陈斯养
(陕西师范大学数学与信息科学学院,陕西 西安 710062)
摘 要: 研究了一类具有离散时滞和干扰的血液模型的 Hopy 分支周期解。 利用函数的单调性,分支
理论及周期函数正交性等方法得到了该模型正平衡态存在唯一的充要条件,分支周期解存在条件和近
似表达式,举出实例且运用 Matlab 绘出了血液模型数值解的拟合图,并分析了参数对周期解的周期,振
幅及正平衡态的影响。
关键词: 干扰血液模型;周期解;时滞;Hopf 分支
中图分类号: O175畅7 文献标识码: A 文章编号: 1007 -9793(2009)05 -0022 -07
1 模型引入
在生态学模型中考虑时滞对生物种群增长的影响极其重要
[6] [12]
,研究时滞产生 Hopf 分支周期解
问题具有更重要的实际意义
[10 -11]
。 1997 年,Mackey and Glass
[1]
提出了一个血液细胞新陈代谢的数学
模型
dp(t)
dt
=
βθ
n
θ
n
+p
n
(t -τ)
-δp(t)
其中 p(t) 表示血液循环过程中成熟干细胞的密度;τ是骨髓中未成熟干细胞由产生到成为成熟干细胞
所需时间:β,θ,n,δ 都是正参数,利用线性化技巧讨论了该模型正平衡态的局部渐近稳定性。对该模型
做变换 p(t) =θx(t),α=
β
θ
可得
dx(t)
dt
=
α
1 +x
n
(t -τ)
-δx(t)
[10]讨论了正平衡态的全局吸引性和全体解关于正平衡态的振动性。 本文在[2 -3][5]的基础上讨论
如下模型
dx(t)
dt
=-ax(t) +
cx
m
(t -τ) +d
m
x(t)
1 +bx
n
(t -τ)
(1畅1)
其中 x(t),τ的生态意义如上,a,b,n 为正参数,c,d,m 为非负实参数。
注:文献[5]中 δ =0 的情况为系统(1畅1) 中 m =1 的特例;文献[10] 的模型为系统(1畅1) 中 m =
0 的特例。
基于生态意义,假设系统(1畅1)满足初始条件
x(s) =φ(s) s ∈ [ -τ,0] φ∈ C([ -τ,0],R
+
) φ(0) >0 (1畅2)
本文共分如下三部分:首先证明了正平衡态的稳定性与局部 Hopf 分支;其次给出近似分支周期解
的近似表达式;最后举出实例,用 Matlab 绘图并对参数对周期解周期,振幅及正平衡态的影响作出分
第 29 卷第 5 期
2009 年 9 月
云南师范大学学报
Journal of Yunnan Normal University
Vol.29 No.5
Sep.2009
倡
收稿日期:2009 -03 -12
基金项目:国家自然科学基金资助项目(60671063),国家自然科学基金资助项目(10871122)畅
作者简介:李娜(1984 -),女,山东省德州市人,硕士研究生,主要从事生态数学方面研究 畅
通讯作者:陈斯养,教授 畅
析。
2 正平衡态的稳定性与局部 Hopf 分支
引理 2畅1 模型(1畅1)满足初始条件(1畅2)的解都是正的。
证明 假设存在 t
1
>0,使得在[0,t
1
) 上,有 x(t) >0,而 x(t
1
) =0,则有
dx(t)
dt
t =t
1
<0 而
dx(t)
dt
t =t
1
=
cx
m
(t -τ)
1 +bx
n
(t -τ)
>0
与之矛盾。 故不存在 t ∈ [0,t
1
) 使得 x(t) =0 的,又 x(t) 关于 t 连续可微,所以 x(t) >0 在其定义域
上恒成立,引理得证。
引理 2畅2 模型(1畅1)满足初始条件(1畅2)的解都是全局存在的。
证明 假设存在一个解不是全局存在的,则存在 t
1
,使得lim
t
→
t
-
1
x(t) =+
∞
,设 t
2
( <t
1
) 为第一个满足
x(t
2
) =δ =
a
c +d
1 -m
的时刻,则
dx(t
2
)
dt
辰 0 而
dx(t
2
)
dt
=-ax(t
2
) +
cx
m
(t
2
-τ) +dx
m
(t
2
)
1 +bx
n
(t
2
-τ)
=-aδ +
cx
m
(t
2
-τ) +dδ
m
1 +bx
m
(t
2
-τ)
<-aδ +
cδ
m
+dδ
m
1 +bx
n
(t
2
-τ)
<-aδ +cδ
m
+dδ
m
=0
矛盾,故引理得证。
定理 2畅1 若 0 尘 m <1,则系统(1畅1) 有唯一的正平衡态 x
倡
。
证明 令方程(1畅1)右端等于零,正平衡态满足方程 (c +d)x
m
=ax(1 +bx
n
) 即(c +d)x
m -1
=a(1
+bx
n
)
令:g(x) =(c +d)x
m -1
,h(x) =a(1 +bx
n
),对任意的 x >0 且 0 尘 m <1,有 g′(x) =(m -1)(c
+d)x
m -2
<0,h′(x) =abnx
n -1
>0,即 g(x) 在区间(0, +
∞
) 上严格单调递减,h(x) 严格单调递增,且
lim
t
→
0
+
g(x) =+
∞
,h(0) =a <+
∞
。故系统(1畅1) 有唯一的正平衡态 x
倡
。证毕。
注:当 m =1 时,若系统(1畅1) 满足 c +d >a,则系统(1畅1) 存在唯一的正平衡态。同[5] 中 δ =0
的情况
对系统(1畅1) 作变换 x(t) =x
倡
e
y( t)
,则 y(t) 满足
dy(t)
dt
=
c(x
倡
)
m -1
e
my( t -τ) -y( t)
+d(x
倡
)
m -1
e
( m -1) y( t)
1 +b(x
倡
)
n
e
ny( t -τ)
-a (2畅1)
(2畅1)式右端在 y(t) =0,y(t -τ) =0 处泰勒展开得
dy(t)
dt
+Ay(t) +By(t -τ) =f(y(t),y(t -τ)) (2畅2)
其中
f(y(t),y(t -τ)) =Cy
2
(t) +Dy(t)y(t -τ) +Ey
2
(t -τ) +ο
3
(y(t),y(t -τ))
A =
cx
m -1
倡
+(1 -m)dx
m -1
倡
1 +bx
n
倡
,B =
abnx
n
倡
-mcx
m -1
倡
1 +bx
n
倡
,C =
cx
m -1
倡
+(m -1)
2
dx
m -1
倡
2(1 +bx
n
倡
)
,
·32·
第 5 期 李 娜,等: 一类干扰血液模型的 Hopf 分支
剩余6页未读,继续阅读
资源评论
weixin_38734269
- 粉丝: 3
- 资源: 930
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
最新资源
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功