设为首页 加入收藏 期刊导航 网站地图
  • 首页
  • 期刊
    • 数学与物理
    • 地球与环境
    • 信息通讯
    • 经济与管理
    • 生命科学
    • 工程技术
    • 医药卫生
    • 人文社科
    • 化学与材料
  • 会议
  • 合作
  • 新闻
  • 我们
  • 招聘
  • 千人智库
  • 我要投搞
  • 办刊

期刊菜单

  • ●领域
  • ●编委
  • ●投稿须知
  • ●最新文章
  • ●检索
  • ●投稿

文章导航

  • ●Abstract
  • ●Full-Text PDF
  • ●Full-Text HTML
  • ●Full-Text ePUB
  • ●Linked References
  • ●How to Cite this Article
Pure Mathematics 理论数学, 2019, 9(8), 890-907
Published Online October 2019 in Hans. http://www.hanspub.org/journal/pm
https://doi.org/10.12677/pm.2019.98116
文章引用: 卢旸, 王倩兰, 毕波. 非自治时滞捕食–食饵模型正周期解的存在性[J]. 理论数学, 2019, 9(8): 890-907.
DOI: 10.12677/pm.2019.98116
The Existence of Positive Periodic Solutions
for the Non-Autonomous Time-Delayed
Predator-Prey Model
Yang Lu
1
, Lan Qian Wang
2
, Bo Bi
3
Department of Applied Mathematics, College of Mathematics and Statistics, Northeast Petroleum University,
Daqing Heilongjiang
Received: Sep. 20
th
, 2019; accepted: Oct. 8
th
, 2019; published: Oct. 15
th
, 2019
Abstract
In this paper, a time-varying predator-prey model of Beddington-DeAngelis functional response
with stage-structured is established. By defining two normal numbers
R
∗
and
R
∗
, we obtain suf-
ficient conditions for the persistence or extinction of system solutions. If
R1
∗
>
, the solution of
the system is uniformly persistent; if
R
1
∗
<
, the predator will be extinct. In addition, according to
the continuous theorem of degree theory, when the system solution is uniformly persistent, it is
concluded that the system has at least one positive periodic solution. Numerical simulation veri-
fies and complements the results of qualitative theoretical analysis, and concludes that the short
maturity period of immature predators is beneficial to the persistent survival of prey and preda-
tor populations.
Keywords
Predator-Prey Model, Non-Autonomous, Time Delay, Uniform Persistence,
Positive Periodic Solution
非自治时滞捕食–食饵模型正周期解的
存在性
卢
旸
,王倩兰,毕 波
东北石油大学数学与统计学院应用数学系,黑龙江 大庆
收稿日期:2019年9月20日;录用日期:2019年10月8日;发布日期:2019年10月15日
卢旸 等
DOI:
10.12677/pm.2019.98116 891
理论数学
摘 要
本文建立了一个捕食者具有阶段结构的Beddington-DeAngelis型功能反应的时变捕食–食饵模型,通过
定义两个正常数
R
∗
和
R
∗
得到了系统解持久或灭绝的充分条件。若
R
1
∗
>
,则系统的解是一致持久的,
若
R1
∗
<
,则捕食者种群灭绝。此外根据度理论的连续性定理,当系统解一致持久时,得到了系统至少
存在一个正周期解。数值模拟验证并补充了定性理论分析的结果,得出幼年捕食者较短的成熟期有利于
食饵和捕食者种群的持久生存。
关键词
捕食–食饵模型,非自治,时滞,一致持久性,正周期解
Copyright © 2019 by author(s) and Hans Publishers Inc.
This work is licensed under the Creative Commons Attribution International License (CC BY).
http://creativecommons.org/licenses/by/4.0/
1. 引言
近年来,具有阶段结构的捕食–食饵模型得到了生物学学者的广泛关注[1]-[6],一般来说,幼年捕食
者从出生到成熟有一个时间差,如果幼年捕食者和成年捕食者用一个固定的年龄来区分,那么这个时间
差称为成熟期,Liu 和Yuan [7]在功能反应项里考虑了成熟期时滞
τ
;S. Liu 和E. Beretta [8]在自己的模型
中也考虑了成熟期时滞
τ
对模型动力学行为的影响,其模型如下:
()()
()
()()
()()
()
()()
()()
()
()
()()
()()
()()
()()
()
12
12
12
12
1
1
e
1
1
e
1
i
i
d
j
d
jj
xtbxtyt
xtrxt
Kkxtkyt
nbxtyt
ytdyt
kxtkyt
nbxtyt
yt
kxtkyt
nbxtyt
dyt
kxtkt
τ
τ
ττ
ττ
ττ
ττ
−
−


′
=−−


++




−−
′
=−

+−+−



′
=

++



−−
−−

+−+−

(1)
其中,
()
xt
为食饵的密度,
()
yt
和
()
j
yt
分别表示成年捕食者和幼年捕食者的密度;n 为捕食者的出生率,
2
k
为捕食者之间干扰效应;b为捕食者对食饵的捕获率大小;
1
k
为捕食者处理食物的时间;假设幼年捕
食者的死亡率为
j
d
,
τ
为幼年捕食者的成熟期。S. Liu和E. Beretta [8]得出如果系统是持久的,那么捕食
者之间足够大的干扰不仅能使系统更加稳定,还可以增加成年捕食者的出生率。
在现实世界中,恒定的自然环境是很少的,任何生物或环境都会受到季节周期变化的影响。正如
Cushing [9]指出,有必要考虑例如,天气、食物供应、收获季节、交配习惯、狩猎等季节性因素的影响。
Shi [10]的模型中的全部参数都假定为周期系数,利用度连续理论讨论了正周期解的存在性。在本文中我
们就食饵具有Logistec 增长率的情形同时针对系统参数为周期时变情形做了定性的分析。
Open Access
卢旸 等
DOI:
10.12677/pm.2019.98116 892
理论数学
2. 模型的建立
()()()
()
()
()()()
()()()()
()
()()()()
()
()()()()
()()
()
()()()()
()()()()
()()
()()()()
()
()()()()
12
d
12
d
1212
1
1
e
1
e
11
t
j
t
t
j
t
dss
dss
jjj
xtbtxtyt
x txtrt
Ktktxtktyt
ntbtxtyt
y tdtyt
ktxtktyt
ntbtxtytntbtxtyt
ytdtyt
ktxtktytktxtktyt
τ
τ
ττττ
ττττ
ττττ
ττττ
−
−
−
−
∫
∫


′
=−−


++





−−−−
′
=−

+−−+−−


−−−−

′
=−−

+++−−+−−

(2)
其中
(
)
xt
和
()
yt
分别为t时刻食饵和成年捕食者的密度,
()
j
yt
为t时刻未成年捕食者的密度。模型(2)
中函数
()()()()()()()()
12
,,,,,,,
j
rtKtbtktktntdtdt
的意义与S. Liu和E. Beretta [8]的模型中的参数意义一
致。
3. 准备工作
首先,为叙述并证明主要结果,在此给出系统(2)的一些假设和符号说明。系统(2)的基本假设为:
(A
1
)
()
()
(
)
()
()
()
(
)(
)
12
,,,,,,,
j
rtKtbtktktntdtdt
为时间t 的连续有界正周期函数。
(A
2
) 若函数
(
)
gt
为
[
)
0,+∞
上的连续有界正周期函数,则
()
0
1
d
T
ggtt
T
=
∫
,
[]
()
0,
min
L
tT
ggt
∈
=
,
[]
(
)
0,
max
M
tT
ggt
∈
=
系统(2)的初值条件为:
()
()
()()()
()
()()(
)()
()()()()()(
)
()
()(
)
0
d
0
12
123
123
e
0d,
1
,,,
0,
0,0,00.
j
s
duu
j
j
nsxsys
ybts
k txsktys
xyy
τ
θφθθφθθφθ
τθ
φφφ
−
−
∫


=

++


===

−≤≤


>

∫
(3)
其中对所有的
[]
,0
θτ
∈−
,有
()
()()()
()
T
123
,,
φθφθφθφθ
=
使得
()
()
01,2,3
i
i
φθ
≥=
且C 表示从
[]
,0
τ
−
到
3
R
上的连续映射的全体构成的Banach 空间
[]
()
3
,,CR
τθ
−
元素
φ
在C 上的范数为:
()()()
{}
123
0
max,,
τθ
φφθφθφθ
−≤≤
=
.
4. 持久性和灭绝性
4.1. 持久性
定义
()
*
1
e
1
M
j
d
LLM
MMM
nbk
R
dkK
τ
−
=
+
,
()
*
1
e
1
L
j
d
MML
LLL
nbk
R
dkK
τ
−
=
+
定理1 假设具有初始条件(3)的系统(2)满足
*
1R>
,那么系统(2)是持久的。更具体地说,我们得到了
如下的结果:
()
1
lim inf
t
xtm
→+∞
≥
,
()
2
lim inf
t
ytm
→+∞
≥
卢旸 等
DOI:
10.12677/pm.2019.98116 893
理论数学
其中
()
2
1
2
LLLM
LL
Krkb
m
rk
−
=
,
()
()
2
exp
M
mdT
ε
ετ
=−+
证明:我们将给出以下几个命题来完成这个定理的证明,证明方法类似于文献[11],[12]中的证明方
法。
定义1 如果存在正常数
,
δ
∆
,
0
δ
<<∆
,使得对于所有初值为正的(2)的解满足:
()(
)
{}
minliminf,liminf
tt
xtyt
δ
→∞→∞
≥
,
()()
{}
maxliminf,liminf
tt
xtyt
→∞→∞
≤∆
则系统(2)的解是持久的。
命题1 具有初始条件(3)的系统(2)的解,对所有的
0t≥
,均为正。
证明:设
()()
()
,xtyt
为具有初始条件(3)的系统(2)的解,首先考虑
()
[]
,0,yttt
∈
时解的正性。由系统
(2)的第二个方程可知
()
()()()()
()
()()()()
()()
()()
d
12
e
1
t
j
t
dss
ntbtxtyt
y tdtyt
ktxtktyt
dtyt
τ
ττττ
ττττ
−
−
∫
−−−−
′
=−
+−−+−−
≥−
(4)
由于
(
)(
)
,0
xy
θθ
≥
连续的依赖于
0
τθ
−≤≤
,并且
(
)
(
)
0,00
xy>
。从而,由比较原理可得对于
[
]
0,
tt
∈
,
(
)
0yt
>
时
()
()
(
)
0
d
0e
t
dss
yty
−
∫
≥
(5)
由系统(2)的第一个方程,对于
[]
0,tt∈
,可得:
()()
()( )
()
()()
()()()()
0
12
d
1
0e0
t
rsxsbsys
s
Ksktxsktys
xtx

−−

++


∫
≥>
(6)
由初始条件(3),对
()
0,t
τ
∈
,可以重新改写
()
j
yt
如下:
()
(
)(
)
()
()
()
()(
)
()
()
()
()
()()
(
)
(
)(
)
()
()
(
)(
)
()
()
()
()()
()
()
d
12
dd
0
0
1212
e
d
1
ee
dd
11
0
t
j
s
tt
jj
ss
duu
t
j
t
duuduu
t
t
nsbsxsys
yts
ksxsksys
nsbsxsysnsbsxsys
ss
ksxsksysksxsksys
τ
τ
−
−
−−
−
∫
∫∫
=
++
=+
++++
>
∫
∫∫
(7)
因此可得基于
[]
,
ττ
−
区间上的
()()()
,,
j
xtytyt
的正性,对于区间
[]
()
,2,,,1,nnnN
ττττ
+∈



,可以采取
同样的方式证明。因此,具有初始条件(3)的系统(2)的解对于任何时刻
0t≥
都是正的。
引理1 [13] 考虑如下方程:
()()()()
2
xtaxtbxtcxt
τ
′
=−−−
其中
,,,abc
τ
为正常数,对任意
[
]
,0
t
τ
∈−
,
()
0xt>
。可得
(1) 如果
ab>
,那么
()
lim
t
ac
xt
c
→+∞
−
=
;
(2) 如果
ab<
,那么
()
lim0
t
xt
→+∞
=
。
命题 2 具有初始条件(3)的系统(2)的解是一致有界的。
卢旸 等
DOI:
10.12677/pm.2019.98116 894
理论数学
证明 设
()()()
()
,ztxtyt=
是具有初始条件(3)的系统(2)的任意一个解,由系统(2)的第一个方程,可
得
()()()
()
()
1
xt
x txtrt
Kt

′
≤−



从而
()
lim
L
t
xtK
→∞
≤
。
定义
()()()()
j
txtytyt
η
=++
,计算
()
t
η
沿着系统(2)的解轨线的导数,可得:
()(
)
()
()
()
()
(
)()
()()()
()
()()()()
()()
(
)(
)
()()()()
()
(
)(
)()
12
12
2
2
1
1
1
jj
LMM
MLL
jj
ML
xtbtxtyt
txtrt
Ktktxtktyt
ntbtxtyt
dtytdtyt
ktxtktyt
rxnb
rxtdytdtyt
Kk
η

′
=−−

++


+−−
++

≤++−−


(8)
对正常数
{}
(
)
min,
LL
j
dd
µµ
<
,由式(8)可得
()()
()
()
2
2
L
MM
M
LM
rxt
nb
ttrxt
kK
ηµηµ

′
+≤++−


因此,由
(
)
lim
L
t
xtK
→∞
≤
,可得对于任意时刻t,都存在正常数
,BT
使得:
()()
()
22
2
:
MLMLMM
LL
Kkrknb
ttB
kr
µ
ηµη
++
′
+<=
从而
()()
limsup0e
t
t
BB
t
µ
ηη
µµ
−
→∞

<+−


(9)
因此,
()()
,xtyt
是一致有界的。结合(7)容易得到
()
j
yt
也是一致有界的。
命题3 系统(2)的任意解
()()()
()
,,
j
xtytyt
满足:
()
()
2
1
2
liminf:0
LLLM
LL
t
Krkb
xtm
kr
→∞
−
≥=>
(10)
证明 由命题2,对于
1
0
tT>>
,存在正常数M,使得
()
ytM≤
。因此,当
1
tT
τ
>+
时,由系统(2)
的第一个方程可得:
()()()
()
()
()
()()
(
)()()()
()
()(
)
(
)
()
12
2
2
1
1
1
M
L
LL
M
M
L
LL
xtbtxtyt
x txtrt
Ktktxtktyt
xtb Mxt
xtr
KkM
rxt
b
xtr
kK

′
=−−

++



≥−−



=−−



从而,
()
()
2
2
liminf
LLLM
LL
t
Krkb
xt
kr
→∞
−
≥
。
卢旸 等
DOI:
10.12677/pm.2019.98116 895
理论数学
命题4 假设
()
*
1
e
1
1
M
j
d
LLM
MMM
nbk
R
dkK
τ
−
=>
+
,当
()
()
*1
2
11
0
MM
M
RkK
k
ε
−+
<<
,
0T
ε
>
时,对于具有初始条件
(3)的系统(2)的任意解
()()()
()
,,
j
xtytyt
,可得:
()
()
()
2
lim infexp:
M
t
ytdTm
ε
ετ
→+∞
≥−+=
(11)
证明 已知
()
*
1
e
1
1
M
j
d
LLM
MMM
nbk
R
dkK
τ
−
=>
+
,则当
0
ε
→
时
()
()
()
2
1
M
M
L
M
b
rt
k
K
rt
Kt
ε
ε

−

+

→




定义
()
()
()
2
1
:
M
M
L
b
rt
k
x
rt
Kt
ε
ε
ε
∆

−

+

−=




则
()
12
e
1
1
M
j
d
LL
MM M
bbx
dkxk
τ
ε
−
∆
∆
>
++
。
考虑以下函数
()
Vt
()
()
(
)(
)()()
()
(
)(
)(
)
()
d
12
e
d
1
s
j
s
d
t
t
nsbsxsys
Vtyts
ksxsksys
τ
θθ
τ
+
−
−
∫
=+
++
∫
(12)
计算
(
)
Vt
沿系统(2)的解轨线的导数
()()()
()
()(
)()
()
()
()()()
(
)
(
)()
()
d
12
12
e
1
e
1
t
j
t
M
j
dss
d
LL
M
MM
ntbtxtyt
Vtdtyt
ktxtktyt
bbxt
dyt
kxtkyt
τ
τ
+
−
−
∫
′
=−+
++


≥−

++


(13)
可知对任意的
()
11
0t tt≥≥
,
()
yt
ε
≤
是不可能成立的。假如
()
yt
ε
>
,可得当
1
tt
τ
≥+
时,
()()()
()
()
()()()
()()()()
()()
()()
()
12
2
1
1
1
M
L
xtbtxtyt
x txtrt
Ktktxtktyt
rtxt
b
xtrt
Kt
k
ε
ε

′
=−−

++



≥−−

+


则
()
()
()
()
2
1
lim
M
M
L
t
b
rt
k
xt
rt
Kt
ε
ε
→∞

−

+

≥




卢旸 等
DOI:
10.12677/pm.2019.98116 896
理论数学
若
()
()
()
1
e
1
M
j
d
LLM
MMM
nbK
dkK
τ
ε
ε
−
−
+−
中的
ε
为足够小的正常数,且满足
()
()
*1
2
11
0
MM
M
RkK
k
ε
−+
<<
,则必存在一个
0T
ε
>
,使得对于所有的
tT
ε
>
,
()
xtx
∆
>
成立。因此当
12
:ttTt
ε
τ
≥++=
时,对于任意的
2
tt≥
,
()
()
()
()
2
1
:
M
M
L
b
rt
k
xtx
rt
Kt
ε
ε
ε
∆

−

+

≥−=




(14)
成立。
由上述不等式,可得:
()
()
()
12
e
1
1
M
j
d
LL
M
MM M
nbx
Vtdyt
dkxk
τ
ε
−
∆
∆


′
≥−

++


.
令
[]
()
*2
,0
minyyt
θτ
τθ
∈−
=++
.
接下来将证明对于任意的
2
tt
≥
,
()
*
yty≥
成立。如果假设不成立,则存在
0T>
使得对于任意的
22
tttT
τ
≤≤++
,有:
()
*
yty≥
,
()
2*
ytTy
τ
++=
,
()
2
0ytT
τ
′
++≤
成立。另一方面,由系统(2)的第二个方程,当
2
ttT
τ
=++
时,
(
)
()(
)
()()
(
)
(
)()()
()
()()
d
12
*
12
e
1
e
0
1
t
j
t
M
j
dss
d
LL
M
MM
ntbtxtyt
y tdtyt
ktxtktyt
nbx
dy
kxk
τ
τ
ττττ
ττττ
ε
−
−
−
∆
∆
∫
−−−−
′
=−
+−−+−−


≥−>

++

显然是矛盾的。从而对于任意的
2
tt≥
,
()
*
yty≥
成立。因此对于任意的
2
tt≥
,可得
()
()
*
12
e
10
1
M
j
d
LL
M
MM M
nbx
Vtdy
dkxk
τ
ε
−
∆
∆


′
≥−>

++


即当
t→+∞
时,
()
Vt→+∞
。由命题2,
()
Vt
是有界的,这显然是矛盾的。接下来,我们将从下面两种
情形分别进行讨论:
(1) 对于任意充分大的时间t,
()
yt
ε
≥
;
(2) 对于任意充分大的时间t,
()
yt
是震荡的。
最终,将得到当t 充分大时,
()
()
()
exp
M
ytdT
ε
ετ
≥−+
成立。显然,我们只需要考虑情形(2),令
1
t
和
2
t
为充分大的时间,且满足当
()
12
,ttt∈
时:
()()
12
ytyt
ε
==
,
()
yt
ε
<
若
21
ttT
ε
τ
−≤+
,因为
()()
()
y tdtyt
′
≥−
,并且
()
1
yt
ε
=
表明对于
[]
12
,ttt∈
,有
()
()
()
exp
M
ytdT
ε
ετ
≥−+
成立。若
21
ttT
ε
τ
−>+
,则对于
[]
11
,tttT
ε
τ
∈++
有
()
()
()
exp
M
ytdT
ε
ετ
≥−+
成立。由(14),可知,对任意
卢旸 等
DOI:
10.12677/pm.2019.98116 897
理论数学
的
[]
12
,
ttTt
ε
τ
∈++
,有
()
xtx
∆
≥
成立。因此,对任意的
[
]
12
,
ttTt
ε
τ
∈++
,有
()
()
()
exp
M
ytdT
ε
ετ
≥−+
成
立。
如果对任意时间
*
11
,
tttTT
ε
τ

∈+++

,存在
0
T
∗
使得
()
()
()
exp
M
ytdT
ε
ετ
≥−+
()
()
()
*
1
exp
M
ytTTdT
εε
τετ
+++=−+
()
*
1
0ytTT
ε
τ
′
+++≤
由系统(2)的第二个方程,当
*
1
ttTT
ε
τ
=+++
时,可得:
(
)
(
)
()
()
(
)
(
)
(
)(
)(
)
()
()
(
)
()
()
d
12
12
e
1
e
exp0
1
t
j
t
M
j
dss
d
LL
MM
MM
ntbtxtyt
y tdtyt
ktxtktyt
nbx
ddT
kxk
τ
τ
ε
ττττ
ττττ
ετ
ε
−
−
−
∆
∆
∫
−−−−
′
=−
+−−+−−


≥−−+>

++

矛盾,因此对所有的
[]
12
,ttt∈
,
()
()
()
exp
M
ytdT
ε
ετ
≥−+
成立。
由系统(2)的第三个方程及命题1,可得:
(
)
()(
)()
()
()
()(
)()()
d
12
12
3
1122
e
d
1
e
:0
1
t
j
s
M
j
duu
t
j
t
d
LL
MM
nsbsxsys
yts
ksxsksys
nbmm
m
kmkm
τ
τ
−
−
−
∫
=
++
≥=>
++
∫
从而,结合命题2~4,系统(2)是持久的。
4.2. 灭绝性
定理2 若
*
1R<
,则对具有初始条件(3)的系统(2)的任意正解
()()
()
,
j
ytyt
满足
()
lim0
t
yt
→+∞
=
,
()
lim0
j
t
yt
→+∞
=
。
证明 由
*
1R<
,存在充分小的
0
ε
>
,使得
()
()
1
e
1
1
L
j
d
MM
LLL
nb
dkk
τ
ε
−
<
++
,
由
()
lim
L
t
xtK
→+∞
≤
,存在时间
0T>
,使得当
tT>
时,
()
(
)()()()
()
()()()()
()()
()
()
()
()()
d
12
1
e
1
e
1
t
j
t
M
j
dss
d
MML
L
LL
ntbtxtyt
y tdtyt
ktxtktyt
nbKt
dtyt
kK
τ
τ
ττττ
ττττ
ετ
ε
−
−
−
∫
−−−−
′
=−
+−−+−−
+−
≤−
++
利用
*
1R<
,可得
()
0yt
′
=
,
()
lim0
t
yt
→+∞
′
=
则
()
lim0
t
yt
→+∞
=
,
卢旸 等
DOI:
10.12677/pm.2019.98116 898
理论数学
从而,容易得到
()
lim0
j
t
yt
→+∞
=
。
注1 当系统(2)的系数
()()()
()()()()()
01
,,,,,,,
j
rtktbtntdtdtktkt
和
()
2
kt
均为常数时,则系统(2)就退
化为以下的自治系统:
()()
()
()()
()()
()
()()
()
()()
()
()
()()
()()
()
()()
()
()()
12
12
1212
1
1
e
1
e
11
j
j
d
d
jjj
xtbxtyt
x txtr
Kkxtkyt
nbxtyt
ytdyt
kxtkyt
nbxtytnbxtyt
ytdyt
kxtkytkxtkyt
τ
τ
ττ
ττ
ττ
ττ
−
−


′

=−−

++




−−

′
=−

+−+−



−−

′
=−−

+++−+−

(15)
则:
()
()
*
*0
1
e
1
j
d
nbK
RRR
dkK
τ
−
===
+
由定理1 不难发现,当
0
1R<
时,本文对Liu 和Beretta [14]中的主要结果进行了改进和推广。
5. 周期解的存在性
为了获得系统(2)的正周期解的存在性,引入连续性引理如下:
首先,取实巴拿赫空间及线性变换
DomLXY
⊂→
,取连续映射
:NXY→
。如果
dimKerL=
codimImL <+∞
且
ImL
在Y 上是闭的,则这个映射L称为指标为零的Fredholm映射,如果L为指标为
零的Fredholm映射且存在连续投影
:PXX→
和
:QYY→
使得
ImPKerL=
,
()
ImImLKerQIQ==−
,
则
()
:Im
DomLKerP
LIPXL
∩
−→
可逆,设
P
L
的逆映射为
P
K
。设
Ω
为X中的有界开集,如果
()
QNΩ
有界
且
()
:
p
kIQNX−Ω→
是紧的,则称N在
Ω
上是紧的。由于
ImQ
与
KerL
同构,所以存在一个同构
:ImJQKerL→
。
引理2 [15] 设
XΩ∈
是一个有界开集,L 是指标为零的Fredholm映射,N 在
Ω
上是紧的,假设
(a) 对任意的
()
0,1
λ
∈
,
xDomL∈∂Ω∩
,
LxNx
λ
≠
;
(b) 对任意的
xKerL∈∂Ω∩
,
0QXN≠
;
(c)
()
deg,,00JQNKerLΩ∩≠
。
则
LxNx=
在
DomL∩Ω
内至少有一个解。
接下来证明系统(2)的正周期解的存在性,方法类似于文献[10]和[16]的方法。
定理3 假设:
(H
1
) 引理2 成立;
(H
2
)
()()()()()()()()
01
,,,,,,,
j
rtktbtntdtdtktkt
和
()
2
kt
均为周期为
0
ω
>
的连续正周期函数;
(H
3
)
()
2
ˆ
ˆˆ
exp2
b
rr
k
ω

−>



;
(H
4
)
*
1R>
,
2
21
MMMM
KdkK
ω
+<
。
则具有初始条件(3)的系统(2)至少存在一个正周期解。
证明:首先考虑具有初始条件:
卢旸 等
DOI:
10.12677/pm.2019.98116 899
理论数学
(
)
()
(
)(
)(
)
()
(
)(
)
()
()
()
()
()()()
0
d
0
12
e
0d
1
,0,0
0,0,00
j
s
duu
j
j
nsxsys
ybts
k txsktys
xy
xyy
τ
θθτθ
−
−
∫


=

++

≥−≤≤


>

∫
在连
续
(16)
的子系统
()
(
)(
)
(
)
(
)
()
(
)(
)
(
)
(
)
()
()
()
()
()
(
)(
)
(
)
(
)(
)
()
(
)
(
)
(
)
12
d
12
1
1
e
1
t
j
t
dss
xtbtxtyt
x txtrt
Ktktxtktyt
ntbtxtyt
y tdtyt
ktxtktyt
τ
ττττ
ττττ
−
−
∫


′
=−−


++





−−−−
′
=−

+−−+−−

(17)
令
()()()()
12
ln,lnutxtutyt==


(18)
将(18)代入(17)可得:
()()
()
()
()
()
()
()
()
()
()
()
()
()
()()()
()
()
()
()
()
2
1
12
122
12
1
12
d
2
12
e
e
1
1ee
ee
1ee
t
j
t
ut
ut
utut
dss
ututut
utut
bt
u trt
Kt
ktkt
ntbt
utdt
ktkt
τ
ττ
ττ
ττ
ττ
−
−
−−+−
−−
∫


′
=−−


++






−−
′
=−

+−+−


(19)
从而,系统(19)有一个
ω
周期解。
接下来首先证明系统(19)至少有一个
ω
周期解。为了将引理2 应用于(19),首先定义
()()
()
()
()
{}
T
2
12
,,:(),1, 2
ii
XYututCRRututi
ω
==∈+==
和
()()
()
[]
()
[
]
()
T
1212
0,0,
,maxmax
tt
utututut
ωω
∈∈
=+
其中
.
表示欧几里德范数,则容易得到X 和Y 是范数为
.
的巴拿赫空间。
令
:LDomLXX
∩→
,
()()
()
()()
()
TT
1212
,,Lutututut
′′
=
。其中
()()
()
()
{}
T
12
12
,,DomLututCRR=∈
且
:NXX→
。
()
(
)
()
(
)
()
()
(
)
()
(
)
(
)
()
()
()
()
()
(
)
()
()
()
()()
()
()
(
)
(
)
()
2
1
12
122
12
1
12
1
d
2
2
12
e
e
1
1ee
ee
1ee
t
j
t
ut
ut
utut
dss
ututut
utut
bt
u trt
Kt
ktkt
Nut
Nut
ntbt
utdt
ktkt
τ
ττ
ττ
ττ
ττ
−
−
−−+−
−−
∫


′
=−−


++




=




−−

′
=−

+−+−


.
定义
(
)
(
)
(
)
()
()
()
1
1
11
1
22
2
1
d
1
d
ut t
PuQu
PuQu
utt
ω
ω
ω
ω



==






∫
∫
,
1
2
u
XY
u

∈=


易知
卢旸 等
DOI:
10.12677/pm.2019.98116 900
理论数学
(
)()()
{
}
2
121212
,,,KerLuuXuuhhR
=∈=∈
且
()()
{}
11
1212
Im(,)d0,d0LuuY uttutt
ωω
=∈==
∫∫
是Y 中的闭集,
dimcodimIm2
KerLL
==
,且P 和Q 为连续投影
ImPKerL
=
,
()
ImImKerQLIQ==−
,
由此可知,L 是指标为零的Fredholm 映射。同时
P
L
的逆映射为
p
k
,则
:Im
p
kLDomLKerP→∩
,且
()
()
()
()
()
()
1
11
00
1
1
2
22
00
1
ddd
1
ddd
t
P
t
P
ussusst
ku
ku
ussusst
ω
ω
ω
ω
ω
ω

−


=




−


∫∫∫
∫∫∫
则由
:QNXY→
及
()
:
P
kIQN XY−→
可得
()
()
(
)
()
()
()
()
()
()
()
()
()
()()()
()
()
()
()
()
2
1
12
122
12
0
12
d
0
12
e
1e
1d
1ee
ee
1
d
1ee
t
j
t
ut
ut
utut
dss
ututut
utut
bt
rtt
Kt
ktkt
QNX
ntbt
dtt
ktkt
τ
ω
ττ
ω
ττ
ω
ττ
ω
ττ
−
−
−−+−
−−
∫




−−



++





=


−−


−


+−+−


∫
∫
(
)()
()
()
0000
11
dddd
2
tt
P
t
kIQNxNxssNxsstNxss
ωω
ωω

−=−−−


∫∫∫∫
显然
QN
和
()
P
kIQN−
是连续的。
为应用引理2,需要寻找适当的开的有界
Ω
子集。相应的算子方程为
(
)
,0,1LxNx
λλ
=∈
,则:
()
()
()
()
(
)
()
()
()
()
()
()
(
)
()
()
()()()
()
()
()
()
(
)
2
1
12
122
12
12
1
d
2
12
e
e
1
1ee
ee
1ee
t
j
t
ut
ut
utut
dss
ututut
utut
bt
rt
Kt
ktkt
ut
ut
ntbt
dt
ktkt
τ
ττ
ττ
λ
ττ
λ
ττ
−
−
−−+−
−−
∫




−−


++




′


=


′

−−



−


+−+−


(20)
假设
()()
()
T
12
,ututX∈
是(20)当常数
(
)
0,1
λ
∈
时的解,则沿着区间
[]
0,
ω
积分可得:
()
()
()
()
()
()
()
()
()
()
(
)
()
()
()
()()()
()
()
()
()
12
12
122
12
0
12
0
d
0
0
12
ee
d
1ee
d
ee
d
d
1ee
t
j
t
utut
utut
dss
ututut
utut
rtbt
t
Kt
ktkt
rtt
ntbt
dtt
t
ktkt
τ
ω
ω
ω
ττ
ω
ττ
ττ
ττ
−
−
−−+−
−−
∫



−


++




=



−−






+−+−


∫
∫
∫
∫
(21)
由
()()
()
T
12
,ututX∈
,则存在
[]
,0,
ii
ξηω
∈
,使得
()
[]
()
0,
min
iii
t
uut
ω
ξ
∈
=
,
()
[
]
()
0,
max
iii
t
uut
ω
η
∈
=
,
1,2i =
由(21)的第一个方程,可得
卢旸 等
DOI:
10.12677/pm.2019.98116 901
理论数学
()
(
)
(
)
()
(
)
()
()
(
)
()
12
12
0
12
ee
ˆ
d
1ee
utut
utut
rtbt
rt
Kt
ktkt
ω
ω

=−

++

∫
此外
()
()
()
()
11
11
0
e
ˆ
ˆ
de
u
u
rt
r
rt
Ktk
ξ
ω
ξ
ωω


>=




∫
可得
()
0
ˆ
d2uttr
ω
ω
′
≤
∫
, (22)
()
111
ˆ
ln:
ˆ
r
ul
r
k
ξ



≤=






, (23)
由式(22)和(23)可得
()
()(
)
111111
0
ˆ
d2:
utuutt lrH
ω
ξω
′
≤+≤+=
∫
, (24)
由式(21)及式(H
3
)的第一个方程可得
(
)
(
)
(
)
(
)
()
()
11
11
0
22
ˆ
e
ˆ
ˆ
de
u
u
rtbt
br
rt
KtktkK
η
ω
η
ωωω



≤+=+







∫
,
从而
()
2
111
ˆ
ˆ
ln:
ˆ
b
r
k
uL
r
K
η


−





≥=







. (25)
由式(22),(25)以及(H
3
)可得
()()()
2
11112
0
ˆ
ˆ
ˆ
dln2:
ˆ
b
r
k
utuuttrH
r
K
ω
ηω


−





′
≥−≥−=







∫
, (26)
从而,结合式(24)可得
[]
()
{}
1121
0,
maxmax,:
t
utHHB
ω
∈
≤=
.
由(21)的第二个式子,可得
()
()
()
()
()()
()
(
)
()
()
()
()
122
2
12
d
00
12
ee
de
1ee
t
j
t
dss
ututut
ut
utut
ntbt
tdt
ktkt
τ
ττ
ωω
ττ
ττ
ττ
−
−
−−+−
−−
∫

−−

=

+−+−

∫∫
, (27)
和
()
()
()
()
()
()()()
()
()
()
()
122
2
12
d
00
12
ee
edd
1ee
t
j
t
dss
ututut
ut
utut
ntbt
dttt
ktkt
τ
ττ
ωω
ττ
ττ
ττ
−
−
−−+−
−−
∫

−−

≤

+−+−

∫∫
,
卢旸 等
DOI:
10.12677/pm.2019.98116 902
理论数学
从而
()()
21
00
2
e
eded
L
j
d
MM
utut
L
L
nb
dtt
k
τ
ωω
−
≤
∫∫
.
结合式(24)可得:
()
1
221
2
e
ln:
L
j
Hd
MM
L
nb
ul
kd
τ
ξ
−
′
≤=
, (28)
再由式(20)可得
()
1
0
d2
M
uttd
ω
ω
′
≤
∫
, (29)
从而,由式(28)和(29)可得
()()()
1
22221
0
2
e
dln2:
L
j
Hd
MM
M
L
nb
utuuttdH
kd
τ
ω
ξω
−
′′
≤≤+=
∫
. (30)
由(27),可得
()
()
(
)
()()()()
()()
12211
2
112 2
d
00
12
ee
edd
1ee
t
M
j
j
t
dss
utututud
LL
ut
uu
MM
nb
dttt
kk
τ
ττξτ
ωω
ξη
−
−
−+−++−
∫


≥

++

∫∫
,
连同式(26)和(H
4
),可得
()
2
1
221
22
e
1
lne:
M
j
d
LLM M
H
MMM
nbdk
uL
dkk
τ
η
−


−


′
≥−=




. (31)
由式(29),(31)和(H
4
)可推得
()()()
221212
0
d2:
M
utuuttLdH
ω
ηω
′′′
≥−≥−=
∫
, (32)
结合(30),可得
[]
()
{}
2122
0,
maxmax,:
t
utHHB
ω
∈
≤=
. (33)
显然
1
B
、
2
B
与
λ
无关,定义
012
BBBB=++
,这里
0
B
足够大使得代数方程组的唯一解
()
T
**
,
uv
()
0
12
d
0
12
ˆ
1e
ed0
1ee
1e
d
1ee
t
j
t
u
u
uv
udss
uv
rb
rt
K
kk
nb
td
kk
τ
ω
ω
ω
ω
−
−
∫



−−=




++








=


++



∫
∫
(34)
满足
()
T
****
,uvuvM=+<
令
()()
()
()
{}
T
T
1212
,:,ut u tXuuMΩ=∈<
满足定理1 条件(H
1
)。当
()()
()
T
2
12
,ututKerLR∈∂Ω∩=∂Ω∩
,则
()
T
12
,uu
是满足
12
uuM+<
的
2
R
中的一个常数向量,则
卢旸 等
DOI:
10.12677/pm.2019.98116 903
理论数学
()
()
()
1
1
12
1
12
0
12
1
d
2
0
12
ˆ
1
ed
1ee
0
0
1e
d
1ee
t
j
t
u
u
uu
udss
uu
rbe
rt
K
kk
QNu
QNu
nb
td
kk
τ
ω
ω
ω
ω
−
−
∫



−−




++





=≠








−


++



∫
∫
则定理3 的条件(H
2
)是满足的。
最后,将证明定理3 的条件(H
3
)是满足的,定义
[]
:0,1domLX
ψ
×→
如下
()
()
()
()
1
1
12
1
12
d
0
12
12
d
0
12
ˆ
e
1e
d
1ee
,,
1e
d
1ee
t
j
t
t
j
t
u
udss
uu
udss
uu
r
r
nb
K
t
kk
uu
nb
t
d
kk
τ
τ
ω
ω
µ
ω
ψµ
µ
ω
−
−
−
−
∫
∫





−






−




++
=+











++



∫
∫
这里
[
]
0,1
µ
∈
,
(
)
()
(
)
T
2
12
,
ututR
∈∂Ω∩
,
(
)
T
12
,uu
是
2
R
中的一个常向量,满足
()()
(
)
T
12
,ututB
=
.
接下来证明
()()
()
T
12
,ututKerL∈∂Ω∩
,
(
)
12
,,0uu
ψµ
≠
。假设结论不真,具有范数
()
()
()
T
12
,ututB=
的常向量
()
T
12
,uu
满足
(
)
12
,,0
uu
ψµ
=
,
()
0
12
d
0
12
ˆ
e
ed0
1ee
1e
d0.
1ee
t
j
t
u
u
uv
udss
uv
rb
rt
K
kk
nb
td
kk
τ
ω
ω
µ
µ
ω
−
−
∫



−−=




++








−=


++



∫
∫
,
结合(23)~(26)以及(28),(30),(31)~(32),可得
()()
()
T
12
,
ututB<
,
这与
()()
(
)
T
12
,ututB=
相矛盾,利用拓扑度理论,取
()()
TT
1212
:Im,,,JIQKerLuuuu=→→
,
直接计算可得
()
()
()()
()
()()
()
()
()
1
1
12
T
T
12
T
12
d
T
0
12
deg,,0,0
deg,,1,,0,0
deg,,0,,0,0
ˆ
1e
dege,d,,0,0
1ee
0
t
j
t
udss
u
uu
JQNKerL
uuKerL
uuKerL
rnb
rtKerL
K
kk
τ
ω
ψ
ψ
ω
−
−
∫
Ω∩
=Ω∩
=Ω∩







=−Ω∩




++




≠
∫
即,定理3 中的条件(H
3
)是满足的,其中
()
T
**
,uv
是(34)的唯一解。
最后,容易证明集合
()
{}
p
kIQNxx−∈Ω
是等度连续且一致有界的,通过使用Arzela-Ascoli定理,
可得,
()
{}
p
kIQNxx−∈Ω
是紧的。因此,N 是L-紧的。
卢旸 等
DOI:
10.12677/pm.2019.98116 904
理论数学
现已证明
Ω
满足定理3 所有的条件,从而,(19)至少有一个
ω
周期解。系统(16)至少有一个正的
ω
周
期解。
令
()()
()
T
**
,xtyt
是系统(16)的一个正
ω
周期解,可以证明
()
()()
()
()()
()()()()
d
**
0
**
12
e
d
1
t
j
t
dss
j
nsbsxsys
yts
ksxsksys
τ
τ
−
−
∗
−
∫


=

++

∫
也是
ω
周期解。从而
()()
()
T
**
,xtyt
为具有初始条件(3)的系统(2)的一个正
ω
周期解。
6. 数值模拟
6.1. 当
R1
∗
>
时,系统(2)的解的持久性以及正周期解的存在性
影响种群持久性的因素很多,其中一个重要因素就是种群成熟期的长短,接下来将针对幼年捕食者
的不同成熟期进行数值模拟:
(1) 选择第一组参数
()()
2sin2rtt=+π
,
()
2Kt
=
,
()
10bt=
,
()
1
1kt=
,
()
2
1kt=
,
()()
32sin2
ntt
=+π
,
()
(
)
32sin2dtt
=+π
,
()
0.1
j
dt
=
。当
2
τ
=
,
*
9.82481R=>
;
6
τ
=
,
*
5.33831R=>
时,由定理1,系
统(2)的解是持久的;由定理1,系统(2)至少存在一个正周期解(见图1 和图2)。
(2) 选择第二组参数
()()
2sin2rtt=+π
,
()
1000Kt=
,
()()
2sin2btt=+π
,
()
1
1kt=
,
()
2
1
kt=
,
(
)(
)
32sin2
ntt
=+π
,
()
(
)
32sin2dtt
=+π
,
()
0.1
j
dt=
。当
6
τ
=
,
*
1R>
时,由定理1,系统(2)的解是
持久的;由定理1,系统(2)至少存在一个正周期解(见图3)。
Figure 1. Basic behavior of solutions of system (3) with
2
τ
=
,
*
9.82481R=>
图1. 系统(2)在
2
τ
=
,
*
9.82481R=>
时的解的基本性态
6.2. 当
1R
∗
<
时,系统(2)中捕食者的灭绝性
选择参数
2
τ
=
,
()()
3sinrtt=+
,
(
)
2Kt=
,
()
1bt=
,
()
1
1kt=
,
()
2
0.5kt=
,
()()
1sinntt=+
,
()()
4sindtt=+
,时
()
1
j
dt=
,
*
0.06011R=<
,捕食者种群灭绝(见图4)。
6.3.
R
∗
与
R
∗
对幼年捕食者成熟期时滞
τ
的敏感度分析
除时滞τ外,当其它参数取值同图1 情形时,
*
R
与
*
R
对幼年捕食者成熟期时滞τ的敏感度分析如图5。
卢旸 等
DOI:
10.12677/pm.2019.98116 905
理论数学
Figure 2. Basic behavior of solutions of system (2) with
6
τ
=
,
*
5.33831
R=>
图2. 系统(2)在
6
τ
=
,
*
5.33831
R=>
时的解的基本性态
Figure 3. Basic behavior of solutions of system (2) with
6
τ
=
,
*
1
R>
图3. 系统(2)在
6
τ
=
,
*
1R>
时的解的基本性态
卢旸 等
DOI:
10.12677/pm.2019.98116 906
理论数学
Figure 4. Basic behavior of solutions of system (2) with
2
τ
=
,
*
0.06011R=<
图4. 系统(2)在
2
τ
=
,
*
0.06011R=<
时的解的基本性态
Figure 5. The relationships between
*
*
,RR
and parameter
τ
, respectively
图5.
*
*
,RR
与
τ
的关系
7. 结论
本文提出了一个基于时变系数和阶段结构的Beddington-DeAngelis型的时滞捕食–食饵模型。通过
应用微分方程比较原理和度理论的连续性定理,对模型进行了定性的理论分析。
在定性理论分析的同时,通过定义两个正常数
*
R
和
*
R
,得到当
*
1R>
时,系统(2)是持久的(见定理1
和图1~3);当
*
1R<
时,捕食者是灭绝的(见定理2 和图4)。
当
*
1R>
时,系统是持久的,本文证明了食饵和成年捕食者存在正的下界(见定理1)。同时,基于食
饵和成年捕食者的正的下界,得到了时变参数下食饵和成年捕食者的种群数量随
*
R
和
*
R
的变化关系。本
文中系统(2)持久性的充分条件弱于文献[17]中相应定理的条件。当系统(2)所有的参数都为常数时,
*
R
和
*
R
就还原为相应自治系统的基本再生数(见注1)。当系统(2)的解持久时,利用度理论中的连续性定理,
得到了系统(2)至少存在一个正
ω
周期解(见图1~3)。
卢旸 等
DOI:
10.12677/pm.2019.98116 907
理论数学
数值模拟中的敏感度分析显示了较短的幼年捕食者成熟期对系统解的持久性是有益的,反之较长的
成熟期会导致捕食者种群的灭绝(见图5)。
值得注意的是系统(2)的动力学行为对于
*
*
1
RR
<<
的情形仍然是不清楚的,本文就此情形并未给出
定性的理论分析证明以及数值模拟。关于模型(2)中成年捕食者方程中含有密度制约项,以及模型(2)存在
多少个正周期解的问题都将留给未来去研究解决。
致 谢
在此特别感谢哈工大的刘胜强教授对本文前期工作的建设性意见,以及东北石油大学校引导性基金
(
项目代码:15011030604)的支持。
参考文献
[1] Aiello, W. and Freeman, H.I. (1990) A Time-Delay Model of Single-Species Growth with Stage-Structure. Mathemat-
ical BioSciences, 101, 139-153.
https://doi.org/10.1016/0025-5564(90)90019-U
[2] Al-omari, J. and Gourley, S. (2003) Stability and Traveling Fronts in Lotka-Volterra Competition Models with Stage
Structure. SIAM Journal on Applied Mathematics, 63, 2063-2086. https://www.jstor.org/stable/4096076
https://doi.org/10.1137/S0036139902416500
[3] Liu, S., Chen, L., Luo, G. and Jiang, J. (2002) Asymptotic Behavior of Competitive Lotka-Volterra System with Stage
Structure. Journal of Mathematical Analysis and Applications, 271, 124-138.
https://doi.org/10.1016/S0022-247X(02)00103-8
[4] Liu, S., Chen, L. and Agarwal, R. (2002) Recent Progress on Stage-Structured Population Dynamics. Mathematical
and Computer Modelling, 36, 1319-1360. https://doi.org/10.1016/S0895-7177(02)00279-0
[5] Liu, S., Chen, L. and Luo, G. (2002) Extinction and Permanence in Competitive Stage Structure System with Time
Delays. Nonlinear Analysis, 51, 1317-1361. https://doi.org/10.1016/S0362-546X(01)00901-4
[6] Ou, L., Lou, Y., Jiang, Y. and Li, Y. (2003) The Asymptotic Behaviors of a Stage Structured Autonomous Preda-
tor-Prey System with Time Delay. Journal of Mathematical Analysis and Applications, 283, 534-548.
https://doi.org/10.1016/S0022-247X(03)00283-X
[7] Liu, Z. and Yuan, R. (2004) Stability and Bifurcation in a Delayed Predator-Prey System with Beddington-DeAngelis
Functional Response. Journal of Mathematical Analysis and Applications, 296, 521-537.
https://doi.org/10.1016/j.jmaa.2004.04.051
[8] Liu, S. and Edoardo, B. (2010) Predator-Prey Model of Beddington-DeAngelis Type with Maturation and Gestation
Delays. Nonlinear Analysis, 11, 4072-4091. https://doi.org/10.1016/j.nonrwa.2010.03.013
[9] Cushing, J.M. (1997) Periodic Time-Dependent Predator-Prey System. SIAM Journal on Applied Mathematics, 32,
82-95.
https://doi.org/10.1137/0132006
[10] Yang, S.J. and Shi, B. (2008) Periodic Solution for a Three-Stage-Structured Predator-Prey System with Time Delay.
Journal of Mathematical Analysis and Applications, 341, 287-294.
https://doi.org/10.1016/j.jmaa.2007.10.025
[11] Zhang, T. and Teng, Z. (2009) Permanence and Extinction for a Nonautonomous SIRS Epidemic Model with Time
Delay. Applied Mathematical Modelling, 33, 1058-1071. https://doi.org/10.1016/j.apm.2007.12.020
[12] Wang, X., Liu, S. and Song, X. (2013) Dynamic of a Nonautonomous HIV-1 Infection Model with Delays. Interna-
tional of Biomathematics, 6, Article ID: 1350030. https://doi.org/10.1142/S1793524513500307
[13] Song, X. and Chen, L. (2001) Optimal Harvesting and Stability for a Two-Species Competitive System with
Stage-Structure. Mathematical Biosciences, 170, 173-186. https://doi.org/10.1016/S0025-5564(00)00068-7
[14] Liu, S. and Edoardo, B. (2006) A Stage-Structured Predator-Prey Model of Beddington-DeAngelis Type. SIAM Jour-
nal on Applied Mathematics, 4, 1001-1129. https://doi.org/10.1137/050630003
[15] Gaines, R.E. and Mawhin, J.L. (1977) Coincidence Degree and Nonlinear Differential Equations. Springer-Verlag,
Berlin. https://doi.org/10.1007/BFb0089537
[16] Zhang, Z. and Luo, J. (2010) Multiple Periodic Solutions of a Delayed Predator-Prey System with Stage Structure for
the Predator. Nonlinear Analysis, 11, 4109-4120. https://doi.org/10.1016/j.nonrwa.2010.03.015
[17] Xu, R., Chaplain, M.A.J. and Davidson, F.A. (2005) Permanence and Periodicity of a Delayed Ratio-Dependent Pre-
dator-Prey Model with Stage Structure. Journal of Mathematical Analysis and Applications, 303, 602-621.
https://doi.org/10.1016/j.jmaa.2004.08.062

版权所有:汉斯出版社 (Hans Publishers) Copyright © 2012 Hans Publishers Inc. All rights reserved.