差分方程基础知识备注

【参考资料】
【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为任意常数

第四步:根据f(t)求非齐次方程的特解,如下:
差分方程基础知识备注_第1张图片

例题

差分方程基础知识备注_第2张图片


n阶常系数线性差分方程求解

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+n1+...+an1yt+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λn1+...+an1λ+an=0

第二步: 计算通解,用两种可能:

  1. 有实特征根 λ \lambda λ,其重复次数为m,则 λ t , t λ t , . . . , t m − 1 λ t \lambda^t, t \lambda^t, ..., t^{m-1} \lambda^t λt,tλt,...,tm1λt是其m个线性无关特征解
  2. 有共轭复根 λ = a ± b i \lambda = a \pm bi λ=a±bi,其重复次数为k,则

{ 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,....,,tk1rtcosωtrtsinωt,trtsinωt,....,,tk1rtsinω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ˉ=Btn1,包括 B B t B t 2 B t n − 1 B Bt Bt^2 Bt^{n-1} BBtBt2Btn1带入原始计算

第四步: 组合通解和特解得到最终解


例题:
差分方程基础知识备注_第3张图片


差分方程基础知识备注_第4张图片


你可能感兴趣的:(数学)