首页 一阶常系数线性差分方程通解求法
文章
取消

一阶常系数线性差分方程通解求法

最近遇到要求解此类差分方程的问题,查阅了相关资料,进行了完善并记录下来

求一阶常系数齐次线性差分方程的通解

一阶常系数齐次线性差分方程的一般形式为 $y_{n+1}-ay_n=0,(a \neq 0)$

迭代法

给定初始值为 $y_0$ ,则 $y_1=ay_0, y_2=ay_1=a^2y_0, y_3=ay_2=a(a^2y_0)=a^3y_0, \dots , y_n=a^ny_0$

其中初始值 $y_0$ 为常数,令 $y_0=C$ , 则通解可表示为 $Y_n=Ca^n$

当存在某一个 $y_x$ 已知时,将其代入通解,可以求得 $C$

特征根法

将原方程变形 $y_{n+1}-ay_n=0,(a \neq 0) \iff y_{n+1}-y_n+(1-a)y_n=0 \iff \Delta y_n+(1-a)y_n=0,(a \neq 0)$

根据 $\Delta \lambda^n=(\lambda-1)^n$ 可以看出 $y_n$ 的形式一定为某一指数函数

设 $y_n=\lambda^n(\lambda \neq 0)$ ,代入原方程得 $\lambda^{n+1}-a\lambda^n=0$ ,即 $\lambda-a=0 \iff \lambda=a$

于是 $y_n=a^n$ 是原方程的一个解,从而 $y_n=Ca^n$ 是原方程的通解

举例

【例1】求 $y_{n+1}-y_n=0$ 的通解

【解】特征方程为 $\lambda-1=0$ ,解得特征根为 $\lambda=1$ ,所以原方程的通解为 $Y_n=C$

【例2】求 $y_{n+1}-2y_n=0$ 的通解

【解】特征方程为 $\lambda-2=0$ ,解得特征根为 $\lambda=2$ ,所以原方程的通解为 $Y_n=C\cdot2^n$

【例3】已知 $y_0=1$ ,求 $y_{n+1}+y_n=0$ 的通解

【解】特征方程为 $\lambda+1=0$ ,解得特征根为 $\lambda=-1$ ,所以原方程的通解为 $Y_n=C(-1)^n$

将 $y_0=1$ 代入,得到 $1=C(-1)^0 \iff C=1$ ,所以原方程的通解为 $Y_n=(-1)^n$

求一阶常系数非齐次线性差分方程的通解

一阶常系数非齐次线性差分方程的一般形式为 $y_{n+1}-ay_n=f(n),(a \neq 0)$

当 $f(n)=0$ 时,方程为 $y_{n+1}-ay_n=0$ ,称它为原方程对应的齐次方程

一阶常系数非齐次线性差分方程的通解为对应的齐次方程通解 $Y_n$ 与原方程的特解 $y^_n$ 之和,即 $y_n=Y_n+y^_n$

当 $f(n)$ 为某些特殊类型的函数时,采用待定系数法求其特解 $y^*_n$ 较为方便

右端函数为m阶多项式类型

原方程变形为 $\Delta y_n+(1-a)y_n=f(n),(a \neq 0)$

由于 $f(n)$ 为多项式,因此 $y^*_n$ 也应该是多项式

当 $a\neq1$ 时,令 $y^*_n=\theta_0 n^m+\theta_1 n^{m-1}+\dots+\theta_m$

当 $a=1$ 时,令 $y^*_n=n(\theta_0 n^m+\theta_1 n^{m-1}+\dots+\theta_m)$

举例

【例1】求 $y_{n+1}-y_n=n^2$ 的通解

【解】对应的齐次方程为 $y_{n+1}-y_n=0$ ,特征方程为 $\lambda-1=0$ ,特征根为 $\lambda=1$ ,齐次方程的通解为 $Y_n=C$

设原方程的特解为 $y^*_n=an^3+bn^2+cn$ ,代入原方程得 $a(n+1)^3+b(n+1)^2+c(n+1)-an^3-bn^2-cn=n^2$

原方程要恒成立,用待定系数法得到 $a=\frac{1}{3}, b=\frac{1}{2}, c=\frac{1}{6}$

所以原方程的通解为 $y_n=\frac{1}{3}n^3+\frac{1}{2}n^2+\frac{1}{6}n+C$

右端函数为指数函数与m阶多项式相乘

设原方程为 $y_{n+1}-ay_n=\mu^nP_m(n),(a \neq 0)$

当 $\mu=0,1$ 时,属于上面一种情况

当 $\mu \neq 0,1$ 时,设 $y_n=\mu^n \cdot z_n$

代入原方程得 $\mu^{n+1}z{n+1}-a\mu^nz_n=\mu^nP_m(n)$

消去 $\mu^n$ ,得 $\mu z_{n+1}-az_n=P_m(n)$ ,就成为了上面一种类型,于是 $y^_n=\mu^n \cdot z^_n$

参考资料

本文由作者按照 CC BY 4.0 进行授权

博客园为代码块添加复制功能

Excel公式初级教程