【参考资料】
【1】《差分方程及其应用》
【2】http://www.docin.com/p-651898262.html
【3】http://www.docin.com/p-590864759.html
设t取离散的等间隔整数值,有 y t = f ( t ) y_t=f(t) yt=f(t),这里 y t y_t yt的取值是一个序列。以二阶差分为例:
y t + 2 + a ( t ) y t + 1 + b ( t ) y t = f ( t ) y_{t+2} + a(t)y_{t+1} + b(t)y_t = f(t) yt+2+a(t)yt+1+b(t)yt=f(t)
注:这里 a ( t ) b ( t ) f ( t ) a(t) \quad b(t) \quad f(t) a(t)b(t)f(t)都是自变量为t的已知函数。如果f(t)不为零,则称其为二阶非齐次线性差分方程。
若函数 y 1 ( t ) y_1(t) y1(t)和 y 2 ( t ) y_2(t) y2(t)是二阶齐次线性差分方程的线性无关特解,则 y C ( t ) = C 1 y 1 ( t ) + C 2 y 2 ( t ) y_C(t)=C_1y_1(t) + C_2y_2(t) yC(t)=C1y1(t)+C2y2(t)是该方程的通解,其中 C 1 C 2 C_1 \quad C_2 C1C2是任意常数。
若函数 y ∗ ( t ) y^*(t) y∗(t)是二阶线性非齐次方程的一个特解, y C ( t ) y_C(t) yC(t)是对应齐次方程的一个通解,则该差分方程的通解为: y t = Y C ( t ) + y ∗ ( t ) y_t = Y_C(t) + y^*(t) yt=YC(t)+y∗(t)
有一阶常系数差分方程 y t + 1 + a y t = f ( t ) y_{t+1} + ay_t = f(t) yt+1+ayt=f(t)
第一步:我们有对应的齐次方程 y t + 1 + a y t = 0 y_{t+1} + ay_t = 0 yt+1+ayt=0,由于 y t + 1 y_{t+1} yt+1是 y t y_t yt的常数倍,因此令 y t = λ t y_t = \lambda^t yt=λt
第二步:带入得到 λ t + 1 + a λ t = 0 \lambda^{t+1} + a \lambda^t = 0 λt+1+aλt=0,得到解 λ = − a \lambda = -a λ=−a
第三步:我们得到齐次方程的通解 y t = C ( − a ) t y_t = C(-a)^t yt=C(−a)t,其中C为任意常数
例题
n阶常系数线性差分方程的一般形式:
y t + n + a 1 y t + n − 1 + . . . + a n − 1 y t + 1 + a n y t = f ( t ) y_{t+n} + a_1y_{t+n-1} + ... + a_{n-1}y_{t+1} + a_ny_t=f(t) yt+n+a1yt+n−1+...+an−1yt+1+anyt=f(t)
第一步:将 y t = λ t y_t = \lambda^t yt=λt带入齐次方程,得到特征方程:
λ n + a 1 λ n − 1 + . . . + a n − 1 λ + a n = 0 \lambda^n + a_1\lambda^{n-1} + ... + a_{n-1}\lambda + a_n = 0 λn+a1λn−1+...+an−1λ+an=0
第二步: 计算通解,用两种可能:
{ r t c o s ω t , t r t c o s ω t , . . . . , , t k − 1 r t c o s ω t r t s i n ω t , t r t s i n ω t , . . . . , , t k − 1 r t s i n ω t \begin{cases} r^tcos \omega t, t r^tcos \omega t, ...., , t^{k-1} r^tcos \omega t \\ r^tsin \omega t, t r^tsin \omega t, ...., , t^{k-1} r^tsin \omega t \end{cases} {rtcosωt,trtcosωt,....,,tk−1rtcosωtrtsinωt,trtsinωt,....,,tk−1rtsinωt
其中 r = a 2 + b 2 r = \sqrt{a^2 + b^2} r=a2+b2, t a n ω = b a tan \omega = \dfrac{b}{a} tanω=ab
第三步: 采用待定系数法计算特解
将 y t ˉ = B t n − 1 \bar{y_t}=Bt^{n-1} ytˉ=Btn−1,包括 B B t B t 2 B t n − 1 B Bt Bt^2 Bt^{n-1} BBtBt2Btn−1带入原始计算
第四步: 组合通解和特解得到最终解