数学建模理论自制笔记2：差分方程及其模型

最新推荐文章于 2023-09-15 09:55:40 发布

绿镇

最新推荐文章于 2023-09-15 09:55:40 发布

阅读量1.6k

点赞数

分类专栏：数学建模文章标签：数学

本文链接：https://blog.csdn.net/qq_42914565/article/details/124829147

版权

数学建模专栏收录该内容

6 篇文章 1 订阅

订阅专栏

1、差分方程基础概念：

差分：这里的差分常指向前差分，即对于数列 $\left \{x_n \right \}$ ，差分算子 $\Delta :\Delta x_n=x_{n+1}-x_n$ 为在 $n$ 处的向前差分；向后差分即是指： $\Delta x_n=x_{n}-x_{n-1}$ ；
在 $n$ 处的 $k$ 阶差分： $\Delta ^{k}x_n=\Delta \left ( ^{k-1}\left ( x_n \right ) \right )$ ；其中在 $n$ 处的2阶差分为 $\Delta ^{2}x_n=\Delta \left ( \Delta x_n \right )$ ，反映的是增量的增量；
差分方程：由 $x_n$ 以及差分 $\Delta x_n,\Delta ^{2}x_n,\cdots$ 所构成的方程，比如 $k$ 阶差分方程形式为 $\Delta ^{k}x_n=f\left ( n,x_n,\Delta x_n,\cdots ,\Delta ^{k-1}x_n \right )$ ，其中这个方程可以化为 $x_{n+k}=F\left ( n,x_n,x_{n+1},\cdots ,x_{n+k-1} \right )$ ，而这种形式是我们经常用的差分形式，线性差分方程的一般形式可以写成： $g(n)=p_{k(n)}x_{n+k}+p_{k-1(n)}x_{n+k-1}+\cdots +p_{0(n)}x_{n}=\sum_{i=0}^{k}p_{i(n)}x_{n+i}$ ；
齐次线性差分方程：线性差分方程式子里 $g(n)\equiv 0$ ；
常系数线性差分方程：线性差分方程式子里 $p_{i(n)}$ 均为常数 $c_i$ ；
差分方程的解：使 $k$ 阶差分方程对所有的 $n$ 都成立的 $x_n$ ；
差分方程的平衡解/平衡点：若常数 $\overline{x}$ 是 $x_n$ 的解，那 $x_n=\overline{x}$ 也就是差分方程的平衡解；
差分算子的性质：
- $\Delta \left ( \alpha x_n+\beta y_n \right )=\alpha \Delta \left ( x_n \right )+\beta \Delta y_n$ ；
- $\Delta \left ( x_ny_n \right )=y_{n+1}\Delta x_n+x_n\Delta y_n=x_{n+1}y_{n+1}-x_ny_n$ ；
- $\Delta \left ( \frac{x_n}{y_n} \right )=\frac{x_{n+1}}{y_{n+1}}-\frac{x_n}{y_n}=\frac{y_n\Delta x_n-x_n\Delta y_n}{y_{n+1}y_n}$ ；

2、简单常系数线性微分方程的解：

一阶常系数线性微分方程：；
1. 当 $g(n)=0$ 时：易知通解 $x_n=x_0a^n$ ；
2. 当 $g(n)=cb^n$ 时：有特解 $x=\left\{\begin{matrix} kb^n & b\neq a\\ knb^n & b=a \end{matrix}\right.$ （其中 $k$ 为待定系数，对应的方程可以求出该值），那么通解为 $x_n=x_0a^n+\left\{\begin{matrix} kb^n & b\neq a\\ knb^n & b=a \end{matrix}\right.$ （即1中的通解+2中的特解）；比如 $x_{n+1}-2x_n=2^n$ 的通解为： $x_n=\left ( x_0 +\frac{1}{2}n\right )2^n$ （其中求出 $k=\frac{1}{2}$ ）；
3. 当 $g(n)=\sum b_mn^m$ 时：有特解 $x=\left\{\begin{matrix} \sum B_mn^m & a\neq 1 \\ n\left ( \sum B_mn^m \right ) & a=1 \end{matrix}\right.$ （其中 $B_m$ 为待定系数，对应的方程可以求出该值），那么通解为 $x_n=x_0a^n+\left\{\begin{matrix} \sum B_mn^m & a\neq 1 \\ n\left ( \sum B_mn^m \right ) & a=1 \end{matrix}\right.$ （即1中的通解+2中的特解）；比如 $x_{n+1}-2x_n=n+1$ 的通解为： $x_n=x_02^n-n-2$ （其中求出 $B_1=-1,B_0=-2$ ）；
二阶常系数线性微分方程： $x_{n+2}+ax_{n+1}+bx_n=g(n)$ ；
1. 当 $g(n)=0$ 时：得到对应的特征方程为 $\lambda ^2+a\lambda +b=0$ ，若存在实根 $\lambda _1,\lambda _2$ ，则方程的特解和实根的情况相对应： $x_n=\left\{\begin{matrix} c_1\lambda _{1}^{n}+c_2\lambda _{2}^{n} & \Delta =a^2-4b> 0\\ \left (c_1+c_2n \right )\lambda _{1}^{n} & \Delta =a^2-4b=0 \end{matrix}\right.$ ；比如 $x_{n+2}+3x_{n+1}+2x_n=0$ 的通解为： $x_n=(-1)^n\left ( c_1+2^nc_2 \right )$ ；
2. 当 $g(n)=cd^n$ 时：有特解 $x=\left\{\begin{matrix} kd^n & \lambda _1\neq d,\lambda _2\neq d\\ knd^n & \lambda _1\neq \lambda _2=d \textup{ or } \lambda _2\neq \lambda _1=d\\ kn^2d^n & \lambda _1=\lambda _2=d \end{matrix}\right.$ （其中 $k$ 为待定系数，对应的方程可以求出该值），那么通解为 $x_n=\left\{\begin{matrix} c_1\lambda _{1}^{n}+c_2\lambda _{2}^{n} & \Delta =a^2-4b> 0\\ \left (c_1+c_2n \right )\lambda _{1}^{n} & \Delta =a^2-4b=0 \end{matrix}\right.+\left\{\begin{matrix} kd^n & \lambda _1\neq d,\lambda _2\neq d\\ knd^n & \lambda _1\neq \lambda _2=d \textup{ or } \lambda _2\neq \lambda _1=d\\ kn^2d^n & \lambda _1=\lambda _2=d \end{matrix}\right.$ （即1中的通解+2中的特解）；比如 $x_{n+2}+3x_{n+1}+2x_n=2^n$ 的通解为： $x_n=(-1)^n\left ( c_1+2^nc_2 \right )+\frac{1}{3}2^{n-2}$ （其中求出 $k=\frac{1}{12}$ ）；
3. 当 $g(n)=\sum d_mn^m$ 时：有特解 $x=\left\{\begin{matrix} \sum D_mn^m & a+b\neq 0 \\ n\left ( \sum D_mn^m \right ) & a+b=0,a\neq -2\\ n^2\left ( \sum D_mn^m \right ) & a=-2,b=1 \end{matrix}\right.$ （其中 $D_m$ 为待定系数，对应的方程可以求出该值），那么通解为 $x_n=\left\{\begin{matrix} c_1\lambda _{1}^{n}+c_2\lambda _{2}^{n} & \Delta =a^2-4b> 0\\ \left (c_1+c_2n \right )\lambda _{1}^{n} & \Delta =a^2-4b=0 \end{matrix}\right.+\left\{\begin{matrix} \sum D_mn^m & a+b\neq 0 \\ n\left ( \sum D_mn^m \right ) & a+b=0,a\neq -2\\ n^2\left ( \sum D_mn^m \right ) & a=-2,b=1 \end{matrix}\right.$ ；（即1中的通解+3中的特解）；比如 $x_{n+2}+3x_{n+1}+2x_n=n+1$ 的通解为： $x_n=(-1)^n\left ( c_1+2^nc_2 \right )+\frac{1}{6}n+\frac{1}{36}$ （其中求出 $D_1=\frac{1}{6},D_0=\frac{1}{36}$ ）；

3、差分方程建模的步骤：

设定好实际问题中的未知函数，按照已知的相关学科的规律来建立相邻的自变量的未知函数取值间的依赖关系，从而建立差分方程模型；
对上述建立的差分方程模型，若能直接求解的则求出其解，若不能直接求解的或者直接求解比较困难的，则用定性的方法讨论其解的性质；
对解得的模型结果与实际情形加以对照，进行讨论。

4、差分方程建模模型举例：

银行存款与利率模型：假如你在银行开设了一个1000元的存款账户，银行的年利率为7%。用 $a_n$ 表示 $n$ 年后你账户上的存款额，那么下面的数列就是你每年的存款额：。
- 设 $r$ 为年利率，由于 $a_{n+1}=(1+r)a_n$ ，因此存款问题的数学模型是： $a_0=1000,r=7%,a_{n+1}=(1+r)a_n,n=0,1,2,\cdots$ ；
- 很明显该方程为一阶齐次常系数线性微分方程，解得 $a_n=1000\cdot 1.07^n$ ；
家庭教育基金模型：从1994年开始，我国逐步实行了大学收费制度。为了保障子女将来的教育费用，小张夫妇从他们的儿子出生时开始，每年向银行存入 $x$ 元作为家庭教育基金。若银行的年利率为 $r$ ，试写出第 $n$ 年后教育基金总额的表达式。预计当子女18岁入大学时所需的费用为100000元，按年利率3%计算，小张夫妇每年应向银行存入多少元？
- 设 $n$ 年后教育基金总额为 $a_n$ ，每年向银行存入 $x$ 元，依据复利率计算公式，得到家庭教育基金的数学模型为： $a_0=x,a_{n+1}=(1+r)a_n+x,n=0,1,2,\cdots$ ；
- 该方程为一阶常系数线性微分方程，解得 $a_n=\frac{x[(1+r)^{n+1}-1]}{r}$ ；这里我们需要求出 $x$ ，则化为 $x=\frac{ra_n}{(1+r)^{n+1}-1}$ ；代入 $a_{18}=10^5,r=0.03,n=18$ ，则 $x=3981.39$ ；
抵押贷款模型：小李夫妇要购买二居室住房一套，共需30万元。他们已经筹集10万元，另外20万元申请抵押贷款。若贷款月利率为0.6%，还贷期限为20年，问小李夫妇每月要还多少钱？
- 设贷款额为 $a_0$ ，每月还贷额为 $x$ ，月利率为 $r$ ，第 $n$ 个月后的欠款额为 $a_n$ ，则 $a_0=2\times 10^5,\left\{\begin{matrix} a_1=(1+r)a_0-x\\ a_2=(1+r)a_1-x\\ \cdots \cdots \\ a_{n}=(1+r)a_{n-1}-x \end{matrix}\right. ,n=1,2,\cdots$ ；
- 该方程为一阶常系数线性微分方程，解得 $a_n=a_0(1+r)^n-\frac{x[(1+r)^{n+1}-1]}{r}$ ；若在第 $N$ 个月还清贷款，令 $a_N=0$ ，则 $x=\frac{a_0r(1+r)^N}{(1+r)^N-1}$ ；代入 $a_0=2\times 10^5,r=0.006,N=20\times 12=240$ ，则 $x=1574.7$ 。