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

期刊菜单

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

文章导航

  • ●Abstract
  • ●Full-Text PDF
  • ●Full-Text HTML
  • ●Full-Text ePUB
  • ●Linked References
  • ●How to Cite this Article
Advances in Applied Mathematics 应用数学进展, 2014, 3, 35-40
http://dx.doi.org/10.12677/aam.2014.31006 Published Online February 2014 (http://www.hanspub.org/journal/aam.html)
A Qualitative Analysis of Model of Differential Equations
of Gambling
Liqun Zhang1,2, Benyan Ding1,2
1Department of Mathematics, Shandong Normal University, Jinan
2Department of Mathematics, Linyi University, Linyi
Email: liqun0911@163.com, dingbenyan2012@163.com
Received: Dec. 16th, 2013; revised: Jan. 14th, 2014; accepted: Jan. 28th, 2014
Copyright © 2014 Liqun Zhang, Benyan Ding. This is an open access article distributed under the Creative Commons Attribution Li-
cense, which permits un restricted use, distribution, and reprod uction in any medium, provid ed the original work is properly cited. In ac-
cordance of the Creative Commons Attribution License all Copyrights © 2014 are reserved for Hans and the owner of the intellectual
property Liqun Zhang, Benyan Ding. All Copyright © 2014 are guarded by law and by Hans as a guardian.
Abstract: This paper gave three models of differential equations for describing social phenomena of gam-
bling. Through qualitative analysis of differential equations, we can analyze the model of differential equa-
tions and the trend of gambling, and get the proper results.
Keywords: Gambling; Singular; Phase Diagram
赌博的微分方程模型及定性分析
张利群 1,2,丁本艳 1,2
1山东师范大学数学系,济南
2临沂大学理学院,临沂
Email: liqun0911@163.com, dingbenyan2012@163.com
收稿日期:2013 年12 月16 日;修回日期:2014年1月14 日;录用日期:2014年1月28 日
摘 要:本文建立三个描述赌博社会现象的数学模型,运用微分方程定性的基本理论来分析微分方程
模型和赌博流行的趋势,得到基本符合实际的效果。
关键词:赌博;奇点;相图
1. 引言
从二十世纪八十年代起,赌博之风死灰复燃。因赌博意识的渗透,高科技手段的运用,赌博之风蔓延异常
迅速,由此带来的社会危害也在日益加剧。而对于赌博的危害,一些人认识不足。有的人认为“赌博只是一种
娱乐而已,大多数人都可以享受赌博的乐趣而不会导致什么问题”。这种认识是极其错误的,大量事实表明,
赌博会给社会经济带来巨大的损失,给家庭和个人带来无比的灾难。为此本文建立三个描述赌博社会现象的数
学模型并通过分析该数学模型来了解赌博流行的趋势。
本文通过运用微分方程定性的基本理论,来分析微分方程模型和赌博流行的趋势,得到基本符合实际的效果。
2. 模型 1
首先借助克马克和麦肯德里克研究微生物传染病动力学时用过的模型。此模型的特点是:不考虑外界措施
等的客观因素,仅仅把赌博这种丑恶现象作为人们正常条件下常态,研究其流行的趋势及结局。
OPEN ACCESS 35
张利群,丁本艳 | 赌博的微分方程模型及定性分析
设
( )
xt
是赌博人员的人数,
( )
yt
是赌博高发人群的人数。建立微分方程模型[1]如下:
,xaxy bx
y axy
= −

= −



(2.1)
其中
,ab
是区间
( )
0,1
中的常数,各项的含义是:
1) (2.1)中的第一式
axy
项表示参与赌博人员与赌博高发人员接触造成的赌博人员人数增加的速度;
2) (2.1)中的第一式
bx
−
项表示由于某些原因(死亡)引起的赌博人员减少的速度;
3) (2.1)中的第二式
axy−
项表示参与赌博人员与赌博高发人员接触造成的赌博高发人员人数减少的速度;
4)
0x≥
,
0y≥
。
下面讨论模型(2.1),由问题的实际意义,只分析第一象限的情形
( )
0, 0xy
>>
。
( )
d10
d
xb
x
y ay
=−+ ≠
(2.2)
积分可得
( )
00 0
ln
by
xyxyy
ay
= +−+
(2.3)
由于(2.1)式的奇点为
y
轴上各点,第一象限无奇点,所以第一象限(2.1)无闭轨。
由(2.2) 可得当
b
ya
>
时,
x
关于
y
是单调递减的;当
b
y
a
<时,
x
关于
y
是单调递增的,又正半轴
x
是一条
轨线。
Figure 1. Phase diagram of model 1
图1. 模型 1的相图
由图 1可以看出,若赌博高发人员人数
( )
yt
不超过值
b
a
,那么当
t→ +∞
时,赌博人员人数
( )
xt
趋于零;
当赌博高发人员人数
( )
yt
超过值
b
a
时,赌博人员将会越来越多,并会发生一个时期的流行,但是赌博高发人员
人数总是越来越少;当赌博高发人员人数减到
b
a
时,赌博人员人数开始减少,且最终趋于零。由此克马克—麦
肯德里克模型是一个乐观的模型,赌博现象最终会消失。
3. 模型 2
模型 1没有考虑政府干预赌博的措施等其他因素,下面将考虑其他因素建立修正模型更客观地反应现象:
2
,xaxy bx
yaxycx ydx
= −

=−− +



(3.1)
其中
,,,abcd
是区间(0,1)中的常数。
OPEN ACCESS
36
张利群,丁本艳 | 赌博的微分方程模型及定性分析
1) (3.1)中
axy
,
bx
−
,
axy−
,
x
,
y
各项意义与(2.1)的意义相同;
2)
2
cx y−
表示政府采取措施后使得赌博高发人员人数减少的速度;政府采取措施的强度
c
与当时赌博高发
人员人数成正比,与当时赌博人员人数的平方成正比;
3)
dx
表示赌博人员人数增多时,赌博高发人员人数也会随之增多,与当时的赌博人员人数成正比。
接下来对(3.1)进行定性分析。
定理 3.1:系统(3.1)在第一象限无闭轨。
证明:取 Dulac 函数
() ()
1
,Bx yxy −
=
,则
( )()
( )
( )
( )
( )
22 12
,,DBx yaxybxBx yaxycx ycxdyxy
xy
δδ
− −−
∂∂
=−+− −++=−−
∂∂
因为
0, 0xy>>
又
0d>
0
δ
≥
,故
0D<
。由 Dulac 准则[2,3],系统(3.1)在第一象限内无闭轨。
注:定理 3.1 的现实意义是赌博的流行不会具有周期性,不用担心赌博在各种恶习会周而复始的出现。
定理 3.2: 1) 当
db>
时,
y
轴的非负半轴上的点是奇点,在第一象限内有唯一奇点
( )
( )
00
,,
ad b b
xy bc a
−

=

,
且
0, b
a



是鞍点。记
( )
( )
2
2
0 00000
42axcxaxdaycx y
∆=++− −
,当
0∆>
时,此种情况不存在;当
0
∆<
时,
( )
00
,xy
是稳定粗焦点,顺时针绕行;当
0
∆=
时,是稳定的退化结点;当
0
∆<
时,
( )
00
,xy
是稳定粗焦点,顺时针绕行。
2) 当
db≤
时,第一象限内无奇点,
y
轴的非负半轴上的点是奇点。
证明:第一步求奇点:
系统(3.1)的奇点应满足如下条件:
2
0,
0
axy bx
axy cxydx e
−=

−−+ +=

(3.2)
由(3.2)的第 一式 得
0x=
或
b
ya
=
;若把
0x=
带入(3.2)式的第二式,则有
y
轴的非负半轴上的点是奇点;若
把
b
ya
=
带入(3.2)式的第二式,则有
20
bb
ax cxdx
aa
− −+=
解得
( )
ad b
xbc
−
=
。
由于只在第一象限内考虑,故上式只能取正号,故要求
db>
,即得唯一奇点
( )
( )
00
,,
ad b b
xy bca
−

=

第二步分析奇点类型:
为了分析奇点
( )
00
,xy
的类型,取(3.1)在奇点
( )
00
,xy
点的线性近似系统
( )
( )
0
2
0000 0
,
2
xax y
ydaycxyxaxcxy
=


=−−− +




(3.3)
(3.3)的特征方程为
( )
( )
22
0 00000
20axcxaxdaycx y
λλ
+ +−−−=
(3.4)
OPEN ACCESS 37
张利群,丁本艳 | 赌博的微分方程模型及定性分析
其中
( )
0
ad b
xbc
−
=
,
0
b
ya
=
故
()
( )
2
2
00 00
20
adb
axdaycxybc
−
− −−=>
(3.5)
所以特征方程(3.4)式中左端的常数项为正,有两个非零实部的特征根:
( )( )
( )
2
22
1,2000000 00
142
2axcxaxcxaxdaycx y
λ

=− +±++−−


因此,
( )
( )
00
,,
ad b b
xy bca
−

=

不是中心, 1,2
λ
只能是两个负数或共轭复数,所以只能是结点或焦点,不是鞍点[2,3]。
在以
()
00
,
xy
为坐标原点经
xoy
平移后得到的坐标系中,线性近似系统(3.3)的第一式可知,上半平面
( )
xt
单
调上升,下半平面
( )
xt
单调下降,所以
( )
00
,xy
不会是双切结点或临界结点,由此可得
∆
不会是正的。
焦点是可能出现的,若取
12abcd== ==
,
00
1xy= =
则
()
( )
2
2
0 00000
42 10axcxaxdaycxy∆=++−−=− <
又由于
1
λ
和
2
λ
的实部均是负的,即
( )
2
00
10
2ax cx− +<
,则可知存在焦点且是稳定的粗焦点。
由于线性近似系统(3.3)的第二式中
0 00
20
daycx y−− <
,在 以
( )
00
,xy
为原点的坐标系中,右半
x
轴实轴轨线
向下,左半
x
轴实轴轨线向上,轨线顺时针绕行。
单切(退化)结点也有可能出现,例如取
abc= =
时,有
0
da
xa
−
=
,
0
1
b
ya
= =
,欲使
0∆=
只需
( )
( )
2
2
0 00000
4 20axcxaxdaycxy∆=++− −=
化简后有
22
4da=
,即
2da=
。
取
2da=
时,会出现退化结点,又特征根是负数即
( )
2
00
10
2ax cx− +<
,故是稳定的退化结点。
4. 模型 3
在模型 1和2的基础之上,我们既考虑政府打击力度的影响又考虑社会力量的影响后建立模型 3:
2
,xaxy bx
yaxy cxydx e
= −

=− −++



(4.1)
各项意义如下:
1) 此模型是在模型 2的基础上增加一项
e
,其他各项如前所述;
2)
e
表示易染赌博人群的受社会影响的自然增长率,且
0e≥
。
类似前边的讨论可有以下结论:
定理 4.1: 1) 若
0e>
,则系统(3.1)在第一象限有唯一奇点:
( )
( )()
2
2
00
4
,,
2
a deadbabceb
xy bc a

−+− −

=

记
( )
( )
2
2
0 00000
42axcxaxdaycxy∆=++− −
,当
0∆=
时,
( )
00
,xy
为稳定的退化结点;当
0∆<
时 ,则
( )
00
,xy
为顺时针绕行的稳定粗焦点;当
0∆>
时,此种情况不存在;
OPEN ACCESS
38
张利群,丁本艳 | 赌博的微分方程模型及定性分析
2) 若
0
e=
,结论如定理 3.1。
由此画出模型 2、3的相图为:
(a)
0, 0,edb= ∆<>
(b)
0, 0,edb= ∆=>
(c)
0,e db= ≤
(d)
0, 0e∆<>
(e)
0, 0e∆= >
Figure 2. Phase diagram of model 2, 3
图2. 模型 2、3的相图
从相图上可以得到以下结论:
1) 赌博高发人员人数
b
ya
=
是流行的阈值,
b
y
a
> 即图中的I区流行区域在这个区域参与赌博的人群人数
随时间的增加而单调递增,赌博发生流行。
2) 赌博高发人员人数
b
ya
<
时,即在
b
ya
=
的水平线下方,是赌博现象的流行缓解区域 II 与扑灭区 III;在
水平虚线表示出的扑灭区
t→+∞
时赌博人员人数趋于零而消失;在斜虚线表出的缓解区域,当
t→+∞
时赌博
人员人数与赌博高发人员的人数将趋于一个极限值,即
( )
0
xt x→
,
( )
0
yt y→
达到稳定平衡状态。
3) 当赌博高发人员人数超过阈值即
b
ya
>
,即使赌博人员人数减少也会发生流行;如果当时赌博人员人数
相当多,但赌博高发人员人数较少
b
ya
<
随着时间推移,赌博高发人员人数也要自然减少最后还是会消除赌博现
象。
4) 由于
( )()
2
2
0
4
2
a deadbabce
xbc
−+−−
=
表明
c
越大,即政府采取措施的强度越大
0
x
越小越易控制;反之
将失控。
5. 总结
优点:本文通过三个模型,逐个分析得出和实际情况较为符合的结论。此时提醒我们政府必须采取切实有
效的措施来控制赌博恶习的传播,使社会处在一个良性发展的空间里。
缺点:一方面没有考虑不易参与赌博的人群对赌博人员和赌博高发人员的影响,当然不易参与赌博人员的
OPEN ACCESS 39
张利群,丁本艳 | 赌博的微分方程模型及定性分析
增长速度受社会力量的干预和赌博高发人员的人数的影响。另一方面对总体人数也没有作一个明确的假定,只
是考虑了人群中两个比较重要的类,而不考虑不易参与赌博人群的影响。
致谢
本文是在金银来老师的指导和丁本艳同学的批评指正下完成的,在此论文发表之际向我的导师和朋友致以
深深地敬意和衷心的感谢,同时也感谢《应用数学进展》各位老师对本文提出的宝贵建议,祝愿《应用数学进
展》越办越好。
参考文献 (References)
[1] 王树禾 (1999) 微分方程模型与混沌. 中国科学技术大学出版社, 北京.
[2] 罗定军, 张祥, 董梅芳 (2001) 动力系统的定性与分支理论. 科学出版社, 北京.
[3] 张芷芬 (1997) 微分方程定性理论. 科学出版社, 北京.
OPEN ACCESS
40

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