数学建模（一）对变化进行建模及其解

杨利霞

数学建模（一）对变化进行建模及其解
1.1用差分方程对变化进行建模
5 D7 O& l8 D% ]5 }# B# l人们常常用数学（使用函数或方程）来描述世界中的某种现象：数学模型，就是对现实世界现象的理想化，但从不是完全精确的表示，任何数学模型都有其局限性，模型没有最好只有最适合，适合的模型可以提供有价值的结论和结果，本次学习内容将主要为：对变化进行建模。
数学模型：人们在对对象进行建模时，常常是为了预测未来某个时刻变量的值，例如模型的对象可能是人口、房价、传染病人数等，数学模型常常能帮助人们更好地了解对象地变化趋势以及内在属性，对于未来地分析与应对具有更好地辅助性。


对变化进行建模：
3 c3 E  X! U( h* c% g未来值=现在值+变化；
9 O; u0 |6 N1 |' D& L1 S6 a变化=未来值-变化；
因此我们想要预测未来值就需要先研究变化。通过一段时间中实际问题所产生的数据，将该数据画成图形，通过对图形的研究，我们常常可以识别出能够抓住这种变化趋势的模式。如果这中行为是在离散时间段上发生的，那么就主要采用差分方程的方法进行建模；若该行为是在连续时间上发生的，就可以采用微分方程进行建模。
4 U8 z: b7 ^9 s( h6 d- m
* f! H3 U8 @9 H+ \9 R0 k
. B4 o- I5 O, s% u& V) s1.1用差分方程对变化进行建模
" }- z/ W5 f+ ~) E& k/ w在本节中主要通过差分方程方法以了解变化发生的原因和形式，分析不同的条件对行为的影响并且去预测未来值。
定义：

, O+ [9 t5 h4 j" i4 ~% x6 \' \% z
例1：储蓄存款类问题
% R$ ^8 \" L# [' S+ |7 w2 h假设本金为1000元的储蓄存单在月利率1%的条件下进行增值，下面的数值表示在该利率下每月的价值：
A={1000，1010，1020.10，1030.30，}其一阶差分为：
* C4 {  p2 o7 \* ~5 Q1 X
! c) k1 e0 c+ Z. b
5 b5 z' z: f- ^注：一阶差分表示在一个时间周期里数列的变化，在该例子中表示的就是所得利息。在对发生在离散时间段上的变化的建模是有很大帮助的，在此例子中，从本月到下一个月变化的仅仅是所得的利息，如果n是月数而an是n个月后存储的总价值，那么每个月价值的变化由第n个分差：


2 y# D. D% ~. ^# i这其中an是n个月后利息累计总值，上述方程可以表示为无穷个代数方程，称为动力系统*，动力系统可以描述从一个周期到下一个周期的变化，如果已知该序列中的某一项，就可以通过差分方程算出下一项。
如果再对该例子进行变化，在上述例子中，如果每月从账户中提取50元，那么一个周期里的变化就应该是该周期里利息减去存款：
! X7 a- W5 q+ Y' s( s& C/ ^
7 \) t- G8 E, O. w
/ m2 {- e* O4 L5 l* Z; q' Y在大多数例子中，用这种数学方式描述变化不会像这种例子一样准确，常常需要画出变化，观察模式，然后用数学术语来描述变化：变化=Δan=某个函数f；或者是变化=Δan=an+1-an=f;用这种方式进行变化建模，就需要决定或者近似决定表示该变化的函数f。
举例：抵押贷款
  W- o- {' d! @小明贷款80000以月利率1%，每月还款880.87元共还20年的方式买了一套房，已经还了72个月，问还欠多少贷款？
: A. E0 F. x& F- V2 \解：分析欠款随着利息而增加而又随着每月还款而减少，所以可以给出下列动力系统模型：
  D, w: e) O9 F% _) c3 h4 g设b(n)为n个月后的欠款总额。则有;
4 H9 ?  ~1 c$ Fb(n+1)=b(n)+0.01b(n)-880.87,b(0)=80000;
, l1 q) \3 E( v所以b(1)=80000+0.0180000-880.87=79919.13;
" E; i% h, o2 [! o: A3 O# o所以b(2)=79919.13+0.0179919.13-880.87=79837.45;
$ p$ o, |+ w& ]. X, R( r% M) s从而得到b(72)=71523.11元，所以还了72个月之后还需要还款71523.11元。
总结：
序列：定义域为全体非负整数集合上的一个函数，其值域为实数的一个子集；
动力系统：序列各项之间的一种关系；
数值解：满足该动力系统的一张数值表。
! |" D" h4 \4 h( x( \

1.2用差分方程近似描述变化
2 ?/ O: w7 ?9 z) d+ w) i; Q) f5 J在我们构建变化模型时，需要区别出这些变化是在离散时间（例如存款利息）上发生的还是在连续时间（物品温度变化）上发生的。
9 W0 f& S0 ]( f, U$ Z% e差分方程：表示在离散时间区间的变化，下面将考察自离散时间区间上的数据近似描述连续变化。
& Z5 M# _; ~' G, X+ p2 M例题1：酵母菌培养物的增长，下图的数据表示的是在一个受限制的区域内，随时间增长而发生的酵母菌培养的数量变化：
) |0 C: M/ x- K2 H! ~2 W, T

从图中可以发现起初酵母菌的增长方式很快，但到达该容器的极限时，就增长很慢了，所以可以假设该容器内的容量为665，该值并不是一个准确值只是在该值附近，变化几乎不变了。通过分析我们可以用以下数学模型：
: L( G9 S9 I- T8 z9 E△P(n)=P(n+1)-P(n)=k(665-P(n))*P(n):
因此我们可以以P(n+1)-P(n)为因变量，以(665-P(n))*P(n)为自变量构成一个线性函数：
  @9 d! S" Q8 W* |; {& V
5 a; S1 i0 `. \/ G- M8 c9 c# `4 }6 G( w2 w
可以求出k=0.00082,所以可以得到动力系统方程：P(n+1)-P(n)=0.00082(665-P(n))*P(n)，这是二次非线性的模型。
' n2 [- u* c! ?
# [* S" H0 t; l9 c- K& A/ d+ f
8 S% a* ~- w" I1.3动力系统的解法
建立某些动力系统解法，从一个初值开始，并对后面的值迭代充分多次，以确定相关模式，动力系统预测的行为是由该系统的数学结构表征的，有时行为极大的变化只是由于动力系统很小的初值引起的，同时还要考虑很小的比例常数变化会造成极大的预测值的变动。
方法：猜测法：
步骤：
1.观察模式：
0 z8 X9 R2 |" \2.猜测动力系统解的形式：
3.用代入法测试该猜测：
  \0 b7 H; M: n# c9 b4.接收或者拒绝该猜测取决于在代入和代数运算后结果是否满足该动力系统。
5 V( \/ ^/ k2 V) f3 k* A# [例如1： 有一款存单，开始存有1000元，然后按每月利率1％给利息，不再存取钱：
解;可以给出如下动力方程：b(n+1)=b(n)+0.01b(n),也就是b(n+1)=1.01b(n),从而b(n)=1.01(k)1000;
例如2： 污水处理中，每小时去除掉剩余污水的12％的污水，问多长时间可以把污水量减少一半？
解：可以给出如下动力方城：b(n+1)=b(n)-0.12b(n)=0.88b(n),所以b(n)=0.88(k)b(0),所以0.5b(0)=0.88(k)b(0);
% K/ \. C8 R4 W  o5 z! P所以0.88(k)=0.5,所以用k=ln(0.5)/ln(0.88),所以k=5.42；
总结：
) p6 l0 ^( a. R# i' N1 }以上两种例题都是可以写成：**b(n+1)=rb(n)**的形式，对于不同的r将会有不同的变化趋势：
' K7 \- n, F5 _8 l/ E1 P同时还有另一种形式：**b(n+1)=rb(n)+b**,其中r和b为常数；
定义：对于b(0)=d,任意k=a,a+1,有b(k)=c，则称c为在初始情景b(0)=d下的平衡点。
**例如3：**当初始存入多少金额时，按照每月1％的利率同时每月取1000元可以使得初始金额不减少。
) B, {; s5 m3 \! _' ]**解：**可以列出如下差分方程：b(n+1)=1.01b(n)-1000;可以算出当b(0)=100000时，金额不会减少。
0 \# A# P. T) k1 b6 S求平衡点并对其进行分类：
" B3 O9 U8 f5 d& z5 w对于b(n+1)=rb(n)+b,当其存在平衡点时，会存在如下关系：b(n)=rb(n)+b，则有：b(n)=b/(1-r),(r不为1和0，r和b为常数)；
6 {0 E# d0 r" N) \2 f定义：对于动力方程：b(n)=r*b(n)+b，其平衡点为：b(n)=b/(1-r)，如果r=1而b=0，则每个点都是平衡点，若r=1,b不为0则不存在平衡点。当|r|<1,时有稳定的平衡点，当|r|>1时有不稳定平衡点，当r=1时，没有平衡点，是一条直线。
6 Y8 M  [: b- p% s. y- R4 ]同时对于上述解我们可以发现结果为：b(n)=r(k)*c+b/(1-r),其中c为初始值，b/(1-r)为平衡值。
8 W0 n0 b7 x9 p3 I3 D**例如4：**对于例3问题中，问初始存多少钱，可以在20年后刚好用完？
% G( e$ \* R) e: G9 P7 w+ c解：b(240)=1.01(240)c+100000=0,可以得到c=90819.42.所以如果初始存90819.42，那么在20年后能刚好取完。
1 ]. f: a4 C; L& J  S+ }1.4差分方程组
上面一节中，平衡点时因变量的一种取值，一旦达到了平衡点，系统就不会再发生变化，而本节则是先求出平衡点，然后对平衡点附近的起始值进行探究，从一个靠近平衡点的初始值开始，变化并观察其：
3 m% x$ d( X1 _* ~7 _  M1.系统是否仍然靠近该平衡点；
2.趋近该平衡点；
9 t! W& h" E3 U6 J0 p3.不再靠近平衡点；
- R; ^: o, G, N: A从而探究其是否有周期性的行为或者振荡行为；或者动力系统对于初始条件或者比例常数的微小变化敏感吗？
  C4 `9 c. o  C9 }: w6 ]例如5：一家汽车租赁公司在A与B设立了两个公司，都可以租还车，从A地开车可以到B地还，也可以从B地开车在A地还，公司需要知道应该在两地各投放多少辆车，可以很好的满足使用。已知每天在A地租的车有60％还到A地，有40％还到B地，而在B地租的车有70％还到了B地，30％还到了A地。
2 P+ O' v; Q3 b! H2 H( W
; N" m; ?9 S" s6 h) z  ^/ f, f% q. ^4 K
1 K' o. @" S2 e' T! G# j9 B8 F若初始投入车辆不是3000和4000，但是总数仍然是7000，则有：
% \$ B3 \8 [# \( d8 v

可以发现一周内都是会区域平衡点，说明平衡点时稳定的而且对初始值是不敏感的，所以我们可以较为确定的倾向于总数3/7的车放在A地，总数4/7的车放在B地。
5 l# s: R% F8 K) ~例如6：
# e# x2 f* k$ F; a) x3 a& u竞争猎兽模型： 斑点猫头鹰和隼，这两个种群在没有其他种群存在时，每个都可以单独的无限增长，但是如果将二者放在一个区间内，则二者将会存在竞争关系，其增长率会发生变化，假设这种增长率的减少和两个种群之间的相互之间的关系成正比。

" q# F+ Z" M1 m
5 l  J  q8 T7 s; `
0 U/ n8 T4 I# U- Y  }
对初始条件的敏感性和长期行为：可以发现当(150,200)时不改变，但是当在该平衡点附近时，都将会有一种将灭绝，也就是该模型对初始条件是及其敏感的。
8 ]: T0 h! ]7 e& y* \( B# |例如7：
离散流行病模型:
考虑一种流行病毒会在人群中传染，我们把人群分为三种：易感染者、已感染者、移除者，同时考虑一下前提条件：
1.每个人都是易感染者S（都能感染这种流感），已感染者I（当前已经感染这种疾病、并且具有传播性），移除者R（得过这种病毒已经不会再感染的人）
2.一开始，每个人都是S或者I；
3.这种流感的持续时间为（1+2/3）周，在这期间是容易传染的；
% H! b+ N- r. Z- o/ q, u+ A, S& \4.该模型按照时间周期周记；
) }# G& i/ Y' z4 Z. Y$ b$ Y1 ^7 F
4 N) g9 Y8 e- J+ G: ]# W
+ p" q4 n. [( T3 R1 k6 U
2 ^' \+ X3 A% w9 x' T. T- c
————————————————
版权声明：本文为CSDN博主「两个鼻孔的猪」的原创文章，遵循CC 4.0 BY-SA版权协议，转载请附上原文出处链接及本声明。
原文链接：https://blog.csdn.net/DOUBLE121PIG/article/details/93375834
. f' m& M' D, I1 M/ q- W7 ]
$ V- D2 O9 ~% W! T9 V1 V! v