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

期刊菜单

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

文章导航

  • ●Abstract
  • ●Full-Text PDF
  • ●Full-Text HTML
  • ●Full-Text ePUB
  • ●Linked References
  • ●How to Cite this Article
Advances in Applied Mathematics 应用数学进展, 2013, 2, 10-14
http://dx.doi.org/10.12677/aam.2013.21002 Published Online February 2013 (http://www.hanspub.org/journal/aam.html)
Stability and Optimal Tax for a Kind of Leslie-Gower
Model with Prey Refuge*
Fengy ing Wei, Yuting Guo
College of Mathematics and Computer Science, Fuzhou University, Fuzhou
Email: weifengying@fzu.edu.cn, guoyutingstil l@163.com
Received: Nov. 24th, 2012; revised: Nov. 30th, 2012; accepted: Dec. 19th, 2012
Abstract: Stability and the optimal tax of Leslie-Gower model with prey refuge is considered in this paper.
By constructing Lyapunov functio n, sufficient conditions of the globally asymptotical stab ility of the positive
equilibrium point for the model are obtained. Moreover, the optimal taxation policy of the model is derived b y
means of Pontryagin maximum principle.
Keywords: Leslie-Gower Model; Equilibrium Point; Globally Asymptotic Stability; Optimal Tax a tion
一类食饵具有避难所的 Leslie-Gower 模型的
稳定性及最优税收*
魏凤英,郭瑜婷
福州大学数学与计算机科学学院,福州
Email: weifengying@fzu.edu.cn, guoyutingstil l@163.com
收稿日期:2012 年11 月24 日;修回日期:2012 年11 月30 日;录用日期:2012 年12 月19 日
摘 要:研究了一类食饵具有避难所的 Leslie-Gower 模型。通过构造合适的 Lyapunov 函数,得到了正
平衡点全局渐近稳定的充分性条件;根据 Pontryagin最大值原理,得到了模型的最优税收策略。
关键词:Leslie-Gower 模型;平衡点;全局渐近稳定;最优税收
1. 引言
Leslie-Gower 捕食者–食饵模型是数学生态模型研究中一类重要的模型[1-9]。众所周知,生物资源并非取之
不尽的,因此,对其进行最优收获策略方面的研究也就十分必要。目前,通过控制 Leslie-Gower 模型的税收,
进而控制其收获的模型还不多。最近,在陈凤德[1]的工作基础上,李有文[2]考虑了一类具有食饵避难的
Leslie-Gower 最优税收模型:

1
1
a mxy 
11
d
d
xrbxx
t ,

2
2
d
d1
ay
yryEqy
x





tm

,

qyc

d
d
EEp
t




ht Eqy

, (1)
并得到了系统(1)的正平衡点的全局渐近稳定的充分条件和其最优税收策略。由于文[2]的收获方式 是
线性的,对于实际的生态系统而言,非线性收获方式更贴近现实,故本文借鉴文[3]的收获函数 Eqy
ht aE by


,
建立如下模型:
*资助信息:国家自然科学基金(No:11201075)和福建省自然科学基金(No:2010J01005)共同资助。
Copyright © 2013 Hanspub
10
魏凤英,郭瑜婷  一类食饵具有避难所的 Leslie-Gower 模型的稳定性及最优税收

1
1
a mxy 
11
d
d
xrbxx
t ,

2
2
d
d1
ay
r
tm



yEqy
y
xaEby






,d
d
pqy
EEc
taEby








00 0
, (2)
初值条件:
 
0,0

,0
x
xyyE E,




,
x
tyt12
,rr

111 2
,, ,aa q
为食饵种群和捕食者种群的密度, 分别为食饵种群和
捕食者种群的内禀增长率,为食饵的避难率, ab 为正常数, 表示捕食者种群收获系数,

0mm ,b

表示总投入的刚性伸缩参数, 表示对 的收获努力量, 为单位资源 的价格,E ypy

表示单位资源 的 y
c税收, 为收获单位资源 的成本,则
y


pqy
Ec
aE by







为收获者的纯经济收入(见[4])。
2. 平衡点的分析
系统(2)有三个平衡点,分别是1
r

11
,0,0Pb


,


211
,,0Pxy,


3,,
P
xyE


21
Px

。
1) 满足:

1
,,0y

2
, 0
1
ay
mx
11 1 2
10rbxamy r

 
,解得:




12
21
1
122 112
1
0,
11
rr
ar
xy
ar mabar m

 1
22
2 1
0
m
ab




3,,
。
xyE
 
满足: 2)
P
 


2
11 12
10, 0,
1
ay Eq
rbxamy rm xaEbyaE

  
 0
pq
y
c
by


, (3)
解得:





2
11
2
1 11
am
rD
a m





1
1
121
1
1
xr
bab D

 ,



11
21
rD
yab
2
11
1
1
m
Da m


,11
2
21 11
1
1
pqbcrDm
EacabD am

 

。





(H1): ,
20qar
2
bc bc
pp
qar


q
;(H2):20qar

,bc
pq

。


0y



3,,PxyE当条件(H1)或者(H2)成立时,有 ,,则0,x

 0E




是系统(2)的唯一正平衡点,其中
12
Dr
pqbc
ap



。
3. 稳定性分析
系统(2)的Jacobi矩阵为:
 
 
 



111 1
2
22
2
22
2
21 1
21
1
()
0
rbxamya mx
ayay aqE
Jr
mx
mx aE byaE
bq p
aq pE
aE byaE


 



22
2
2
22
0
bqy
by
yc
by


























。
1) 平衡点 1
11
,0,0
r
Pb



的特征方程为:





0c

12
rr


,于是 1
11
,0,0
r
Pb



是鞍点。
Copyright © 2013 Hanspub 11
魏凤英,郭瑜婷  一类食饵具有避难所的 Leslie-Gower 模型的稳定性及最优税收
2) 平衡点 点的特征方程为:

211
,,0Pxy12
3
0

32


,其中

22
12 1
11
aa yK
xx
12 1
111 2211211 23112
,,
aa y
bx rKbxrbx rKbxrK
 








,
其中 pq
K
c





 123
0, 0
b,当 时,有0K0,


 0K,由 Routh-Hurwitz判别法[6]知:当

时,

,,0Pxy

3,,PxyE
 12
3
=0



是局部渐近稳定的。
211
3) 平衡点 点的特征方程为: ,其中
32



2
yBC



12
11 21 1
,aa
bx ADAbxbx ADx

  
,

2
12
31 1
bxBC
 


aa y
bx ADD
x


 

 ,

2
1
ay
Amx



 
2
bqy E
aE by



,






22
22
,
bq yaq pE
BC
aE byaE by




 
2
,
bc
Dc
qp




 








0A123
0, 0
,
当 ,(H1)或者(H2)成立时,有0,




3,,yE
 

3,,Pxy
,由 Routh-Hurwitz 判别法[6]知,Px 是局部渐
近稳定的。
定理 1 若正平衡点

E




存在,则 3,,PxyE


在第一象限
A
yxP




的上半平面内是全局渐近稳
定的,其中
2
m
PaaEby


1bqE。
证明:取正定函数

 
1
, ,lnln2
ln
x
xy E
xyEdy yydEEEVxy E
 
 


x

 





 ,,d为
12
d任意正常
数,由(3)式得到


1
1a my

,

1
r



pqy
c

aE by

1
bx

2
2
r1
ay
Eq
mx
aE by



 
,




,
于是,



 










12
2
11
1
12
2
2
11
1
d1 ddd
11 1
ddd d
1
1
1
VxxyyEE
dd
txxtytE t
bxxdbqyyEy
am
xx yy
xx aE byaE
dayyE yEy
yydp aqEE
mx
xaE byaE
bxxd bqEy
am
xx yy
xx

 








 
 
 
 


 






 


Ey
by
by





 



 










2
12
1
2
2
2
1
y
aE byaE by
yy EEday
dbqyx xyy
mxx
aEbyaE by
dpaqyEEyy EE
d paqE
aE byaEbyaE byaE



 



 

 






 

 
2
12
1
dayy
mx
by







Copyright © 2013 Hanspub
12
魏凤英,郭瑜婷  一类食饵具有避难所的 Leslie-Gower 模型的稳定性及最优税收
选取正数


2
1
1
2
1,
amx
d
ay


1
2
bd y
d
a pE




,则
 







222
11 122
d
d(1
bxxdbqEyydayydp
V
tx mx
aE byaEbyaE by

 
 

 


2
)
aqyEE
aE by



,
当


2
1E m
aE byaaE

 
bq x
by

时,则有 d0
d
V
t


3,,PxyE


在第一象限
A
yxP 的上半平 。因此,正平衡点
面内是全局渐近稳定的。证毕。
4. 最优税收政策
本节的目的是通过控制税收使得系统的社会总收入最大,社会收入的贴现值为:
0
tpq
ed
y
J
cEt
by





aE


,

et

为贴现因子,我们的目的是利用 Pontrjagin 最大值原理找到最优税收政策

为资本贴现率,其中

t

,
使得
J
同时满足系统(2)和控制约束条件min max


时的值最大。该控制问题的 Hamilton 函数是:


   
111 1
2
22 3
e1
1
tpqy
HcEtrbxxa
aE by
ay Eqy
tr ytE
m xaEby










mxy
pqy
c
aE by








 










 
123
,,ttt

min


,
其中 是伴随变量[9]。假设Hamilton 函数的最优解不发生在




和max

时,我们有


30
y
tq
yE
H
aE b





30t


,从而 。构造辅助方程:
  

2
2
22
1
ay
tmx









111
11
d21
d
tHtr bxamy
tE


  

, (4)


 

2
22
112 2
2
d2
e1
d1
t
tay
HapqE tamxt r
ty m
aE by


2
2
aqE
x
aE by





 











, (5)

 

22
2
bqy
aE by
32
2
de
d
t
tHbpqy ct
tE aE by









 









30t



, (6)
由 及(6)得到:

2
caEby
y
2et
tp
bq








 







3,,PxyE
 
, (7)
为了得到最优解,在正平衡点处整理方程(4)得到


2
et
Pt P
111
d
d
t
t




 

,
其中 ,
11
Pbx




22
2
aE by
x






2
21
abpq yc
Pbq m


,解得 2
1
1
et
P
tP




。

Copyright © 2013 Hanspub 13
魏凤英,郭瑜婷  一类食饵具有避难所的 Leslie-Gower 模型的稳定性及最优税收
Copyright © 2013 Hanspub
14
同理得到:


2
et
Qt Q
212
d
d
t
t



,
其中


2
12
2
1
ay
Qr mx


 


2
aqE
aE by








,


21
221
1Pamx
apqE
QP
aE by









,解得
2
1
et
Q
tQ
2






, (8)
将(7)代入(8 )得到:

2
2
1
et
Q
y
Q
2etca
Eb
tp
bq y















Px ,,
 , (9)
将的值

,,yE
 
3yE
 
代入(9)得到关于

的方程,令
x
x


为(9)的解。再将




代入 ,,yE
 
,,
则得到最优
平衡解
x
xyyEE



,,
 。
5. 结论
通过对一类具有避难所的Leslie-Gower 模型的最优税收研究,本文得到了该系统的平衡点及其局部渐近稳
定的充分性条件。利用构造合适的 Lyapunov 函数,进一步得到正平衡点的全局渐近稳定的充分条件。进而由
Pontryagin 最大值原理得到最优平衡解
x
xyyEE


 。我们的研究结果说明:通过控制税收来控制收获能
达到可持续发展和经济效益最大化的双赢结果,有利于资源的合理开发。
参考文献 (References)
[1] F. D. Chen. On a Leslie—Gower predator—Prey model incorporating a prey refuge. Nonlinear Analysis: Real World Applications, 2009, 10:
2905-2908.
[2] 李有文, 杨洪娴, 田广立, 张运丽. 具有食饵避难的 Leslie-Gower最优税收模型分析[J]. 数学实践与认识, 2011, 41: 167-171.
[3] 贾春莹, 高建国. 改进了的捕食者–食饵系统的征税模型分析[J]. 数学的实践与认识, 2010, 40: 237-242.
[4] T. K. Kar. Conservation of a fishery through optimal taxation: A dynamic reaction model. Communications in Nonlinear Science and Numeri-
cal Simulation, 2005, 10: 121-131.
[5] 张睿. 具有比率依赖型功能反应函数的食饵–捕食者征税模型分析[J]. 生物数学学报, 2008, 23: 484-488.
[6] 马知恩, 周义仓. 常微分方程定性与稳定性方法[M]. 北京: 科学出版社, 2001.
[7] 陈兰荪, 宋新宇, 陆征一. 数学生态学模型与研究方法[M]. 成都: 四川科学技术出版社, 2006.
[8] 马知恩. 种群生态学的数学建模与研究 [M]. 合肥: 安徽教育出版社,1996.
[9] C. W. Clark. Mathematical bioeconomics: The optimal management of renewable resources. New York: John Wiley , 1976: 88-120.

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