差分方程模型的稳定性分析.docx
- 文档编号:568833
- 上传时间:2022-10-11
- 格式:DOCX
- 页数:18
- 大小:366.65KB
差分方程模型的稳定性分析.docx
《差分方程模型的稳定性分析.docx》由会员分享,可在线阅读,更多相关《差分方程模型的稳定性分析.docx(18页珍藏版)》请在冰豆网上搜索。
差分方程模型的稳定性分析
摘要
微分方程是研究数学的一个重要分支,是本科期间我们必须掌握的基本知识,而本文我们研究的是一个递推关系式,也称差分方程。
它是一种离散化的微分方程,是利用描述客观事物的数量关系的一种重要的数学思想来建立模型的。
而利用差分方程建立模型解决问题的方法在生活中随处可见,比如在自由竞争市场经济中的蛛网模型是利用差分方程分析经济何时趋于稳定,又如金融问题中的养老保险也是利用差分方程来分析保险品种的实际投资价值。
而差分方程模型是描述客观世界中随离散时间变量演化规律的有力建模工具。
本文首先给出差分方程的定义以及求解过程并给出判断差分方程稳定性的判断方法,随后以同一环境下的羊群和草群的相互作用为模型分析其种群的数量变化过程,进而研究线性差分方程的稳定性,最后用一个实际模型来更好的说明差分方程的稳定性对解决实际问题有非常大的帮助。
关键字:
差分方程;差分方程模型;平衡点;稳定性
Abstract
Differenceequationisalsocalledrecursiveequation,itistodescribetherelationshipbetweenthenumberofobjectivethingsofakindofimportantmathematicalmodel.Andtheuseofthedifferentialequationmodelofthesolutioncanbefoundeverywhereinlife.Suchascobwebmodelinthefreemarketeconomyistousethedifferenceequationanalysiswhentheeconomicstability,andasthefinancialproblemofpensioninsurancebreeddifferenceequationisusedtoanalysistheactualinvestmentvalue.Thispapergivesthejudgethestabilityofdifferenceequationtojudgemethod,theninthesamegroupofsheepandgrassundertheenvironmentofinteractionanalysisforthemodelaprocess,thenumberofthepopulationchange,inturn,studythestabilityofthelineardifferenceequation.Intheend,onepracticalmodeltobetterexplainthestabilityofdifferenceequation.
Keywords:
Differenceequation;Differenceequationmodel;Balancepoint;Stability
引言
随着科学技术的不断发展,将数学思想融入实际生活解决社会问题变得非常普遍。
所以利用差分方程建立模型也显得至关重要。
在经济、社会、生态、医疗、网络、遗传学得某些数据都是按时、日、周、星期、月份、年等汇总和统计的,这时将时间离散化后建立差分方程模型更为方便,从而解决社会问题趋于稳定的状态,它是描述客观世界中随离散变量演变规律的一种重要的建立模型的方法,在现实生活中有很多问题都是借助差分方程模型来刻画并求解的,利用数学的思路与想法来研究实际问题,从而确保某个体系稳定运作的条件,进一步再结合其他条件分析,为客观体系的安全稳定运作提供理论上的保障,因此差分方程模型的稳定性分析是我们数学中研究的一个重要课题。
本文以同一空间下的羊群和草群的相互作用为模型分析这两物种的数量变化过程,进而研究线性差分方程的平衡点及其稳定性;最后根据差分方程的平衡点及其稳定性分析的相关理论解决实际问题。
我相信差分方程的稳定性相关理论将在未来更为应用普遍。
1、差分方程的定义及其分类
(1)差分算子:
定义1:
设
是定义在R上的函数,则
称
在
的差分,
称为差分算子,
称
在
的位移,E称为位移算子;用I表示恒等算子,即
这些算子都是线性算子,都是针对函数所定义的映射。
(2)差分方程:
定义2:
含有未知函数及未知函数差分的等式,我们称为差分方程,它的一般表达形式为:
由
(1)与
(2)的关系,可以将阶数为
的差分方程写为
或者
我们称
不显含
时的方程为自治差分方程。
形如
表示一阶差分方程;
表示n阶差分方程。
(2)差分方程的分类:
差分方程可以分为两大类:
其一为线性差分方程,它是指当
是
的线性函数时,称
为线性差分方程;也就是说
的次数都为
,其二为非线性差分方程,它是指当
是
的非线性函数时,称
为非线性差分方程。
显而易见,非线性差分方程求解比线性差分方程求解复杂,因此它的解的性态也比较难分析,本文我们只研究线性差分方程解的性态。
2.差分方程的求解与稳定性判断方法:
(1)差分方程的求解:
使得差分方程称为恒等式的序列称为差分方程的解。
满足方程及初始条件的序列称为初始值问题的解,形如
,
称为自治差分方程的初始问题;当
含有
时,
称为非自治差分方程的初始值问题。
那么,现在知道差分方程的解的定义,问题是如何求出一个差分方程的解呢?
这里我们给出普遍的解法----迭代法
定义3:
连续用变量的原值推算出新值的一种递推过程称为迭代法。
下面介绍一个具体的迭代过程:
类比常系数一阶微分方程的解法,我们可以容易求得常系数一阶差分方程的通解为:
式中
为任何常数。
现在将
代入通解中可得
,所以满足初始条件
的特解为
。
于是我们可得:
;
;
;
……
现在我们利用该方法来求解以下方程的初始值问题:
例1:
解:
其解序列的前几个为:
;
;
;
;
…
这个初始值问题解的一般形式是
。
那么此差分方程也满足其他初始条件的解,显然
和
都是此差分方程的解。
如果其方程满足初始值
,那么它的解的一般形式为
。
这里注意此差分方程的解当
时的极限:
当
时,有
。
例2
解:
将其转化为
其解序列的前几个利用迭代法可得:
;
;
;
;
…
(2).差分方程的平衡解稳定性判断方法:
定义4:
若有
,使
=
,则
为差分方程
的平衡点,
是差分方程
满足
如果对任意给定的正数
,有
>0,使得当
时,
对所有的
都成立,则称差分方程
的平衡解
是稳定的,否则,称为不稳定的。
我们也可以定理1.31(参考文献[4]32页)分析差分方程平衡解的稳定性。
其定理为;设
有连续的三阶导数,
为差分方程
的平衡解,则
时,
是渐进稳定的,
>
时,
是不稳定的;当
=
,
0时,
是不稳定的,当
=
,
=0,
>0时,
是不稳定的,当
=1,
=0,
<0时,
是稳定的;当
=
,
<0时,
是稳定的,当
>0,
是不稳定的。
对于阶数为
的线性差分方程平衡点的稳定性条件是它的特征解,也就是
次代数方程的解
均有
。
3.差分方程模型的应用:
3.1模型:
种群模型
3.11模型的引入与假设
在某个生态环境中,羊以草为食。
研究将羊群放入草场后羊和草两种群在同一环境下的种群数量变化。
草的生长遵循Logistic规律(当草群数量太大时,种群会发生生存竞争,草群的增长率受到环境最大容纳量等因素的影响,从而导致增长率的降低)每年固有增长率为0.7,最大密度为2800(密度单位),在草最茂盛时每只羊每年可吃掉1.2(密度单位)的草。
若没有草,羊群的年死亡率高达0.8,然而草的存在可使羊的死亡得以补偿,在草最茂盛的时候补偿率为1.1.在这种情况下,羊群和草群的种群数量将如何变化,羊会把草吃完从而导致羊的数量也减少还是两种群的数量趋于稳定呢?
这里我们以200只羊放入密度为1000和密度为2800的草场两种情况分析。
通常,以建立模型的用途为出发点,我们将连续对象离散化更为方便,从而我们采用差分方程建立模型来描述羊和草两物种的数量变化过程。
模型假设:
1.草场上除了羊群以外,没有其他以草为食的生物;
2.草独立生存且遵从logistic规律;
3.没有草的情况下羊一定会死;
4.假设每只羊每年的食草能力是草场密度的线性函数;
5.假设草对羊的补偿率也是草场密度的线性函数。
3.12线性差分方程模型的建立与求解
分析:
假设第
年草的密度为
,羊的数量为
,第
年草的密度为
,羊的数量为
。
记草的固有增长率为
,草的最大密度为
,羊独立生存时的年死亡率为
,草最茂盛时羊的吃草能力为
,草对鹿的年补偿作用为
。
建立差分方程模型:
草的增长差分模型为
(草独立生存时满足logistic增长规律)
但实际上羊对草增长有影响,草的数量会减少,则方程变为
(
)
(1)
羊的增长差分模型为
(羊独立生存时)
但实际上草的存在可以补偿羊的死亡率,则方程变为
(
)
(2)
此外,记初始状态草场的密度为
,初始状态羊的数量为
,各个参数值为
利用MATLAB软件分析计算该差分方程模型:
具体算法见附录
将密度为1000和密度为2800的草场上分别放200只羊时的两种情况如下(如图1):
由图中可以看到,蓝色曲线代表草场密度的初始值为1000时,两种群变化情况;而红色曲线则代表草场密度的初始值为2800时,两种群的变化情况。
经过仔细查看两种情况下曲线的变化规律,我们可以发现在大约25-50年时间后,两物种的数量将趋于平衡状态。
使用MATLAB软件可以计算出当
=(1900,550),意味着两种群数量的平衡点为(1900,550)。
3.13生态模型的平衡点及稳定性分析:
联立方程
(1)
(2)消去
得
对
递推不难得到
当
时
时,得到平衡点的稳定条件是
如果条件改变,结果也会不同。
即当
趋于无穷时,方程组的解会有不同情形。
那么现在我们以改变羊的数量初值,改变草场的最大密度
改变羊群独立生存时的死亡率。
改变这三个条件来研究这个生态模型两种群的稳定状态。
改变羊的数量初值有如下情况(草场初值取2000,羊群初值分别取10,100,500,2000):
由图2可以看到,从理论上看最终羊群与草群两种群数量的平衡值不受羊初始的数量影响。
然而,我们观察到,y0=2000的那条曲线(紫色曲线),在5-10区间降到了最低点,但这显而易见是不可能的,是不符合羊的繁殖客观现象的,因为羊的种群可以持续繁殖的最低数量是存在一定限制。
当种群数量不大于这个值时,在实际情况下,羊的种群就要灭绝。
改变草场的最大密度有如下情况:
如图4所示,如果草场密度的最大值N发生变化,则最终两种群数量的平衡点也会发生相应的变化。
结论:
如果N值增加,那么平衡点两种群的数量就增加;N值减少,相应的平衡点两种群的数量就减少。
改变羊群独立生存时的死亡率:
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 方程 模型 稳定性 分析
![提示](https://static.bdocx.com/images/bang_tan.gif)