差分方程的解析解

您所在的位置:网站首页 一阶差分方程的特解形式 差分方程的解析解

差分方程的解析解

2023-11-29 15:55| 来源: 网络整理| 查看: 265

差分方程的解析解

例题.求$y_t=a_0+a_1y_{t-1}+\epsilon_t$的解析解

当初始条件已知,即$y_0$已知时,迭代算法如下:

向前迭代算法

$$y_1=a_0+a_1y_0+\epsilon_1$$$$y_2=a_0+a_1y_1+\epsilon_2$$$$y_3=a_0+a_1y_2+\epsilon_3$$$$…$$$$y_t=a_0+a_1y_{t-1}+\epsilon_t$$

以上列举了该差分方程的所有项,我们要做的就是通过这些递推关系来逐步迭代,当迭代完所有项之后,最终可以得到目标解析解.

$$y_1=a_0+a_1y_0+\epsilon_1$$$$y_2=a_0+a_1y_1+\epsilon_2=a_0+a_1(a_0+a_1y_0+\epsilon_1)+\epsilon_2=$$$$a_0+a_0a_1+{a_1}^2y_0+a_1\epsilon_1+\epsilon_2=[a_0(1+a_1)]+[a_1^2y_0]+[a_1\epsilon_1+\epsilon_2]$$$$y_3=a_0+a_1y_2+\epsilon_3=a_0+a_1(a_0+a_0a_1+{a_1}^2y_0+a_1\epsilon_1+\epsilon_2)+\epsilon_3=$$$$[a_0+a_0a_1+a_0{a_1}^2]+[a_1^3y_0]+[a_1^2\epsilon_1+a_1\epsilon_2+\epsilon_3]=[a_0(1+a_1+a_1^2)]+[a_1^3y_0]+[a_1^2\epsilon_1+a_1\epsilon_2+\epsilon_3]$$$$…$$由数学归纳法,递推下去,则有:

$$y_t=[a_0(a_1+a_1^2+…+a_1^{t-1})]+[a_1^ty_0]+[a_1^{t-1}\epsilon_1+a_1^{t-2}\epsilon_2+…+a_1^0\epsilon_t]=$$$$a_0\sum_{i=0}^{t-1}a_1^i+a_1^ty_0+\sum_{i=0}^{t-1}a_1^i\epsilon_{t-i}$$上式便是目标解析解,因为$y_0$已知,即不含有任何未知项.

之所以称之为“向前迭代法”,是因为迭代的索引是从$y_1$开始,一步一步往前面跑去,逐步迭代得到到目标式$y_n$的.

你肯定会想到,既然有向前迭代,那么一定有向后迭代了,没错,向后迭代就是从最后一项$y_t$逐渐向$y_0$进行迭代.

向后迭代

同样是求解$y_t=a_0+a_1y_{t-1}+\epsilon_t$的解析解:$$y_t=a_0+a_1y_{t-1}+\epsilon_t$$$$y_{t-1}=a_0+a_1y_{t-2}+\epsilon_{t-1}$$$$y_{t-2}=a_0+a_1y_{t-3}+\epsilon_{t-2}$$$$…$$$$y_2=a_0+a_1y_1+\epsilon_2$$$$y_1=a_0+a_1y_0+\epsilon_1$$$$y_0=a_0+a_1y_{-1}+\epsilon_0$$开始迭代:(前面写了两次迭代过程,总共迭代$t-1$次)

$$y_t=a_0+a_1y_{t-1}+\epsilon_t=$$$$a_0+a_1(a_0+a_1y_{t-2}+\epsilon_{t-1})+\epsilon_t=$$$$a_0+a_1a_0+a_1^2y_{t-2}+a_1\epsilon_{t-1}+\epsilon_t=$$$$[a_0(1+a_1)]+[a_1^2y_{t-2}]+[a_1\epsilon_{t-1}+\epsilon_t]$$$$a_0+a_1a_0+a_1^2(a_0+a_1y_{t-3}+\epsilon_{t-2})+a_1\epsilon_{t-1}+\epsilon_t=$$$$[a_0+a_1a_0+a_1^2a_0]+[a_1^3y_{t-3}]+[a_1^2\epsilon_{t-2}+a_1\epsilon_{t-1}+\epsilon_t]=$$$$[a_0(1+a_1+a_1^2)]+[a_1^3y_{t-3}]+[a_1^2\epsilon_{t-2}+a_1\epsilon_{t-1}+\epsilon_t]=$$$$…$$$$=[a_0(1+a_1+a_1^2+…+a_1^{t-1})]+[a_1^ty_0]+[a_1^{t-1}\epsilon_{1}+a_1^{t-2}\epsilon_{2}+…+a_1^0\epsilon_t]=$$$$a_0\sum_{i=0}^{t-1}a_1^i+a_1^ty_0+\sum_{i=0}^{t-1}a_1^i\epsilon_{t-i}$$

而当初始条件未知,即$y_0$未知时,我们需要对上面得到的差分方程的解继续化:

$$y_t=a_0\sum_{i=0}^{t-1}a_1^i+a_1^ty_0+\sum_{i=0}^{t-1}a_1^i\epsilon_{t-i}=$$$$a_0\sum_{i=0}^{t-1}a_1^i+a_1^t(a_0+a_1y_{-1}+\epsilon_0)+\sum_{i=0}^{t-1}a_1^i\epsilon_{t-i}=$$$$a_0\sum_{i=0}^{t-1}a_1^i+(a_0a_1^t+a_1^{t+1}y_{-1}+a_1^t\epsilon_0)+\sum_{i=0}^{t-1}a_1^i\epsilon_{t-i}=$$$$a_0\sum_{i=0}^{t}a_1^i+a_1^{t+1}y_{-1}+\sum_{i=0}^{t}a_1^i\epsilon_{t-i}=$$$$…$$$$=a_0\sum_{i=0}^{(t-1)+m}a_1^i+a_1^{t+m}y_{-m}+\sum_{i=0}^{(t-1)+m}a_1^i\epsilon_{t-i}$$

由于初始条件未知,我们现在来讨论解是否存在,即解的敛散性:

当$|a_1|1$时,发散; 当$|a_1|=1$时,直接代回到最开始的式子,即得$$y_t=a_0t+y_0+\sum_{i=0}^{t-1}\epsilon_{t-i}=a_0t+y_0+\sum_{i=1}^{t}\epsilon_{i}$$

以上便是针对例题展开的求解过程

而对于一般的线性差分方程$$y_t=a_0+a_1y_{t-1}+a_2y_{t-2}…+a_ny_{t-n}+x(t)...(1)$$其齐次线性差分方程为$$y_t=a_1y_{t-1}+a_2y_{t-2}…+a_ny_{t-n}...(2)$$我们首先要明确,

对于齐次线性差分方程(2式),对任意常数$A$,如果$y_t^h$是其一个解,那么$Ay_t^h$也是其一个解. 而对于非齐次线性差分方程(1式),我们需要求其一个特解及其对应的齐次线性差分方程的一个通解.

我们接下来将介绍一种新的解法(不同于之前的迭代法),称之为备选(特征根法)解法,来求解非齐次线性差分方程,即(1)式.

步骤如下:

(1)建立齐次线性差分方程(2)式,求齐次方程的解$y_t^{h_1},y_t^{h_2},…,y_t^{h_n}$;(2)求出非齐次线性差分方程(1)式的一个特解$y_t^p$; 

(3)非齐次线性差分方程(1)式的通解即为其自身的一个特解加上其对应的齐次线性差分方程的一个通解,即$$y_t=y_t^p+A_1y_t^{h_1}+A_2y_t^{h_2}+…+A_ny_t^{h_n}$$

其中$A_1,A_2,…A_n$是任意常数;

(4)若已知初始条件$y_0,y_1,…$,则可以求出$A_1,A_2,…A_n$.(选)

其中第一步是最麻烦的,所以下面我们将针对第一步,即齐次线性差分方程的求解问题展开讨论

例题:解齐次线性差分方程$$y_t=a_1y_{t-1}+a_2y_{t-2}…+a_ny_{t-n}$$

一阶齐次差分方程$y_t=ay_{t-1}$,通解为$Aa^t$ ,可以用迭代法证明,用特征根法更为简单:

猜想其解的形式为$y_t=A\alpha^t$,代入上式,得$A\alpha^t=aA\alpha^{t-1}$,即$\alpha=a$,故通解为$y_t=aA\alpha^{t-1}=Aa^t$.

二阶齐次差分方程$y_t=a_1y_{t-1}+a_2y_{t-2}$ ,

利用特征根法,猜想其解的形式为$y_t=A\alpha^t$,代入上式,得:$$A\alpha^t=a_1A\alpha^{t-1}+a_2A\alpha^{t-2}$$化简,得:$$\alpha^2=a_1\alpha+a_2$$移项,得:$$\alpha^2-a_1\alpha-a_2=0$$从而解得特征根为$$\alpha_{1}=\frac{a_1{+}\sqrt{a_1^2+4a_2}}{2}$$$$\alpha_{2}=\frac{a_1{-}\sqrt{a_1^2+4a_2}}{2}$$其中$d=a_1^2+4a_2$ 

$d$有三种情况:$d>0,d=0,d0$时,$\alpha_1,\alpha_2$为互不相同的实根,此时的通解为$y_t=A_1\alpha_1^t+A_2\alpha_2^t$,其中$A_1,A_2$为任意常数.

当$|\alpha_1|$或者$|\alpha_2|$大于1时,发散;

当$|\alpha_1|$或者$|\alpha_2|$小于1时,收敛.

(2)当$d=a_1^2+4a_2=0$时,即$a_2=-\frac{a_1^2}{4}$时,$\alpha_1=\alpha_2$为重根,代回上面的特征根方程有$\alpha_1=\alpha_2=\frac{a_1}2$,此时$A\alpha_1^t=A{(\frac{a_1}2})^t$是方程的一个解(俩都一样,算一个),并且可以验证,$At{(\frac{a_1}2})^t$也是方程的一个解,所以此时方程的通解为$y_t=A_1{(\frac{a_1}2})^t+A_2t{(\frac{a_1}2})^t$,其中$A_1,A_2$为任意常数.

当$|\frac{a_1}2|1$时,发散.

(3)当$d=a_1^2+4a_2



【本文地址】


今日新闻


推荐新闻


CopyRight 2018-2019 办公设备维修网 版权所有 豫ICP备15022753号-3