差分方程齐次解单根例
差分方程方法

a1
k
k 1
a2
k 2
ak 0
称为差分方程(1)的特征方程,其特征方程的根 称为特征根。
3 2013年7月18日
1.常系数线性齐次差分方程
(1) 特征根为单根
设差分方程(1)有 k 个单特征根 1 , 2 ,, k , 则通解为
xn c c2 ck
9
如果仍用固定努力量的捕捞方式,该公司应采取怎样的策略 才能使总收获量最高。
18 2013年7月18日
四、案例:最优捕鱼策略问题
2. 模型的假设 (1)只考虑鱼的繁殖和捕捞的变化,不考虑鱼群迁 入与迁出; (2)各龄鱼在一年的任何时间都会发生自然死亡;
(3)所有鱼都在每年最后四个月内完成产卵孵化 的过程,成活的幼鱼在下一年初成为一龄鱼; (4)产卵发生于后四个月之初,产卵鱼的自然死 亡发生于产卵之后;
a x0 x1 xn1 b ,试求函数的导数值 f ( xk )(k 1,2,, n) 。
用差商代替微商,则有
f ( xk 1 ) f ( xk ) 向前差: f ( xk ) (k 1,2,, n) xk 1 xk f ( xk ) f ( xk 1 ) 向后差: f ( xk ) (k 1,2,, n) xk xk 1
!!无捕捞时 鱼群会无限 的增长吗??
由假设(1)和(2) :
dxi (t ) rxi (t ), i 1, 2,3, 4; dt k t k 1, k 0,1, 2,
各龄鱼都不会无限地增长!
21
No! I don’t know!
2013年7月18日
3、模型的建立与求解
差分方程方法

第四章差分方程方法在实际中,许多问题所研究的变量都是离散的形式,所建立的数学模型也是离散的,譬如,像政治、经济和社会等领域中的实际问题。
有些时候,即使所建立的数学模型是连续形式,例如像常见的微分方程模型、积分方程模型等等,但是,往往都需要用计算机求数值解。
这就需要将连续变量在一定条件下进行离散化,从而将连续型模型转化为离散型模型,因此,最后都归结为求解离散形式的差分方程解的问题。
关于差分方程理论和求解方法在数学建模和解决实际问题的过程中起着重要作用。
下面就不同类型的差分方程进行讨论。
所谓的差分方程是指:对于一个数列x n,把数列中的前n 1项x i i 0,1,2, n 关联起来所得到的方程。
4.1 常系数线性差分方程4.1.1 常系数线性齐次差分方程般形式为常系数线性齐次差分方程的一x n a1x n 1a2 x n 2a k x n k 0 (4.1)其中k 为差分方程的阶数,a i i 1,2, ,k为差分方程的系数,且a k 0 k n 。
对应的代数方程k k 1k2k a1k 1a2k 2a k0(4.2 )称为差分方程的(4.1)的特征方程,其特征方程的根称为特征根。
常系数线性齐次差分方程的解主要是由相应的特征根的不同情况有不同的形式。
下面分别就特征根为单根、重根和复根的情况给出差分方程解的形式。
1.特征根为单根设差分方程( 4.1)有k 个单特征根1, 2, 3, , k ,则差分方程( 4.1 )的通解为x n c1 1 c2 2c k k n,其中c1,c2,,c k 为任意常数,且当给定初始条件x i i0i 1,2, ,k (4.3)时,可以唯一确定一个特解。
2.特征根为重根设差分方程(4.1 )有|个相异的特征根1, 2, 3, , I 1 l k重数分别为lm1,m2 , ,m l且m i k 则差分方程(4.1 )的通解为i1k ,则差分方程的通解为为已知函数。
m i X ni 1 n C 1i n11m 2i 1 nQi n 2i 1mli5n同样的,由给定的初始条件3.特征根为复根4.3 )可以唯一确定一个特解。
差分方程(1)-基础知识省名师优质课赛课获奖课件市赛课一等奖课件

形如
yx+2 + ayx+1 + byx = f (x).
(10)
(其中 a , b 0, 且均为常数)旳方程, 称为二阶常系数线性 差分方程. 当 f (x) = 0 时, 即
yx+2 + ayx+1 + byx = 0
(11)
称为齐次差分方程; 当 f (x) 0时, 称为非齐次差分方程.
类似于二阶线性常微分方程, 二阶线性差分方程与 其有相同旳解旳构造. 故先求齐次方程(11)旳通解.
故所求通解为
yx
C1
C2 (2)x
10 3
x
2x2.
(2) f (x) = Cqx 设特解旳待定式为
y x Bq x (q不是特征根); y x Bxq x (q是特征方程单根); y x Bx2q x (q是二重特征根). 其中 B 为待定系数.
例11 求差分方程 yx+2 3yx+1 + 2yx = 2x旳一种特解.
x B0 B1x Bm xm (a 1) (6)
或
y x (B0 B1x Bm xm ) x (a 1) (7)
其中B0 , B1 , , Bm为待定系数.
例5 求差分方程 yx+1 2yx = 3x2 旳一种特解.
解 这里 a = 2, 设 y x B0 B1x B2 x2 , 代入差分方程, 得
解 相应旳齐次方程旳特征方程为
2 3 + 2 = 0.
方程旳根为
1 = 1, 2 = 2,
因为 q = 2 =2, 设特解为 y x Bx2x ,
代入原方程, 得
B(x+2)2x+23B(x+1)2x+1+2Bx2x = 2x,
线性差分方程

线性差分方程内容提要:1 齐次线性差分方程1-1 一阶齐次线性差分方程1-2 二阶齐次线性差分方程(容许复数解)1-3 二阶齐次线性差分方程(容许实数解)1-4 齐次线性差分方程2 线性差分方程3 例子本文主要参考文献.由于最近需要用到一些线性差分方程,所以这里做一个复习小结.注:由于阶数为 2 或者 2 以上,处理方法毫无区别,所以我们集中火力搞定 2 阶情形,一般情形则不加证明给出结果. 但不难由 2 阶情形照搬证明过去.1 齐次线性差分方程1-1 一阶齐次线性差分方程称如下形式的方程为序列 \{z_t, \ t\in \mathbb{Z} \} 的一阶齐次线性差分方程:z_t =a_1 z_{t-1} ,式中 a_1 为实数.\bullet 显然这个方程的解为z_t =C a_1^t . C 为任意实数.1-2 二阶齐次线性差分方程(容许复数解)称如下形式的方程为序列 \{z_t, \ t\in \mathbb{Z} \} 的二阶齐次线性差分方程:z_t =a_1 z_{t-1} + a_2 z_{t-2} ,式中 a_1, a_2 为实数.[特征方程与特征根] 我们把矩阵A={ \left[ \begin{array}{cc} a_1 & a_2 \\ 1 & 0\end{array} \right ]} 的特征多项式\lambda^{2}=a_{1}x+a_{2}称为齐次线性差分方程 z_t =a_1 z_{t-1} + a_2 z_{t-2} 的特征方程,而它的两个根\lambda_{1},\lambda_{2} (可能有重根)叫做特征根.[特解]z_{t}=\lambda_{i}^{t} ( i=1,2 ) 为方程的特解.[证明] 由\lambda_{i}^{2}=a_{1}\lambda_{i}+a_{2} ,两边同时乘以 \lambda_{i}^{t-2} ,得\lambda_{i}^{t}=a_{1}\lambda_{i}^{t-1}+a_{2}\lambda_{i}^{t-2}因此z_{t}=\lambda_{i}^{t} ( i=1,2 )满足原方程.1-2-1 不等特征根情形\bullet 如果 \lambda_{1}\ne\lambda_{2} , 那么,方程z_t =a_1 z_{t-1} + a_2 z_{t-2} 的通解为z_{t}=C_{1}\lambda_{1}^{t}+C_{2}\lambda_{2}^{t}.[证明] 由于\begin{array}{llll} a_{1}z_{t-1}+a_{2}z_{t-2}\\=a_{1}\left( C_{1}\lambda_{1}^{t-1}+C_{2}\lambda_{2}^{t-1}\right)+a_{2}\left( C_{1}\lambda_{1}^{t-2}+C_{2}\lambda_{2}^{t-2}\right)\\=C_{1}\left( a_{1}\lambda_{1}^{t-1}+a_{2}\lambda_{1}^{t-2} \right)+C_{2}\left( a_{1}\lambda_{2}^{t-1}+a_{2}\lambda_{2}^{t-2}\right)\\=C_{1}\lambda_{1}^{t}+C_{2}\lambda_{2}^{t}\\=z_{t} \end{array}所以对任意的常数 C_{1},C_{2}, 我们都有z_{t}=C_{1}\lambda_{1}^{t}+C_{2}\lambda_{2}^{t} 是方程 z_t =a_1 z_{t-1} + a_2 z_{t-2}的解.还需要验证所有的解具有这个形式. 对于给定的一组初值 z_{0},z_{1},有\begin{array}{llll}C_{1}+C_{2}=z_{0}\\C_{1}\lambda_{1}+C_{2}\lambda_{2}=z_{1}\\\end{array}这个关于 C_{1},C_{2} 的二元一次方程组的系数矩阵的行列式为\left|\begin{array}{cccc} 1 & 1 \\\lambda_{1} & \lambda_{2}\end{array}\right| \not=0所以给定初值z_{0},z_{1},就能唯一确定系数 C_{1},C_{2}. 1-2-2 相等特征根情形\bullet 如果 \lambda_{1} = \lambda_{2}= \lambda , 那么,方程 z_t =a_1 z_{t-1} + a_2 z_{t-2} 的通解为z_t =(C_1 +C_2t) \lambda^t .[证明] 由于 \lambda 是特征多项式\lambda^{2}=a_{1}x+a_{2}的二重根 ,所以它也是 \lambda^{t}=a_{1}\lambda^{t-1}+a_{2}\lambda^{t-2} 的二重根. 把\lambda^{t}=a_{1}\lambda^{t-1}+a_{2}\lambda^{t-2} 的两边对 \lambda 求导,得t\lambda^{n-1}=a_{1}\left( t-1\right)\lambda^{t-2}+a_{2}\left( t-2\right)\lambda^{t-3},因为重根求导之后仍为根,所以 \lambda 是 t\lambda^{n-1}=a_{1}\left( t-1 \right)\lambda^{t-2}+a_{2}\left( t-2 \right)\lambda^{t-3} 的根,两边乘以 \lambda 得到\lambda 也是t\lambda^{t}=a_{1}\left( t-1\right)\lambda^{t-1}+a_{2}\left( t-2\right)\lambda^{t-2} 的根,即z_{t}=t\lambda^{t} 也是特解. 容易验证z_t=(C_1 +C_2t) \lambda^t 都是方程 z_t =a_1z_{t-1} + a_2 z_{t-2} 的解.还需要验证所有的解具有这个形式. 对于给定的一组初值z_{0},z_{1},有\begin{array}{llll}C_{1}=z_{0}\\C_{1}\lambda+C_{2}\lambda=z_{1}\\\end{array}这个关于 C_{1},C_{2} 的二元一次方程组的系数矩阵的行列式为 \left|\begin{array}{cccc} 1& 0 \\ \lambda & \lambda\end{array}\right|\ne0所以给定初值z_{0},z_{1},就能唯一确定系数 C_{1},C_{2}.1-3 二阶齐次线性差分方程(容许实数解)延续上一节的记号.\bullet (i) 若特征方程有两不等实根 \lambda_1,\lambda_2 ,那么这个方程的解为z_t =C_1 \lambda_1^t+C_2 \lambda_2^t . C_1, C_2 为任意实数.\bullet (ii) 若特征方程有两相等实根 \lambda_1=\lambda_2 = \lambda ,那么这个方程的解为z_t =(C_1+C_2t) \lambda^t . C_1, C_2 为任意实数.\bullet (iii) 若特征方程有两共轭复根 \lambda_1=re^{iw}, \lambda_2=re^{-iw}, 那么两个特解为z_t=r^{t}e^{iwt} ,z'_t=r^{t}e^{-iwt},由欧拉公式有z_t=r^{t}[cos(wt)+isin(wt)],z'_t=r^{t}[cos(wt)-isin(wt)].特解含有复数部分,我们希望解是实的,可以凑出新的两个特解r^{t}cos(wt)与 r^{t}sin(wt) , 因此通解为z_t =C_1r^{t}cos(wt) +C_2 r^{t}sin(wt) .1-4 齐次线性差分方程[齐次线性差分方程] 称如下形式的方程为序列 \{z_t, \t\in \mathbb{Z} \} 的齐次线性差分方程:z_t =a_1 z_{t-1} + a_2 z_{t-2} + \cdots +a_p z_{t-p} ( )式中, p\geq 1 , a_1, a_2, \cdots a_p 为实数.[特征方程与特征根] 我们把矩阵A={ \left[ \begin{array}{cccccc} a_1 & a_2 &a_3&\cdots &a_{p-1} & a_p\\ 1 & 0 & 0&\cdots &0 & 0\\ 0 & 1 & 0&\cdots &0 & 0\\ \cdots &\cdots &\cdots&\cdots &\cdots &\cdots \\ 0 & 0 & 0&\cdots &1 & 0\end{array} \right ]} 的特征多项式\lambda^{p}=a_{1}\lambda^{p-1}+a_{2}\lambda^{p-2} +\cdots +a_p称为齐次线性差分方程 ( ) 的特征方程,而它的 p 个非零根\lambda_{1},\lambda_{2},\cdots,\lambda_{p} (可能有重根)叫做特征根.\bullet 如果 \lambda_{i} 为两两不等的实根, 那么,方程( ) 的通解为z_{t}=C_{1}\lambda_{1}^{t}+C_{2}\lambda_{2}^{t}+\cdots +C_{p}\lambda_{p}^{t}.2 线性差分方程[线性差分方程] 称如下形式的方程为序列 \{z_t, \ t\in\mathbb{Z} \} 的线性差分方程:z_t =a_1 z_{t-1} + a_2 z_{t-2} + \cdots +a_p z_{t-p}+h( t). ( )式中, p\geq 1 , a_1, a_2, \cdots a_p 为实数而 h(t) 为t 的已知函数. 并且称方程:z_t =a_1 z_{t-1} + a_2 z_{t-2} + \cdots +a_p z_{t-p} ( )为( )的导出齐次线性差分方程.\bullet 线性差分方程( )的解为导出齐次线性差分方程( )的通解和特解之和.3 例子[例1] (等差数列) 等差数列z_{t+1}=z_{t}+d 为一阶线性差分方程.它的导出齐次方程为 z_{t+1}=z_{t} , 特征根为 \lambda=1 . 于是导出齐次方程的解为 z_t=C.猜测原方程的一个特解为 z_{t} = dt , 那么全部解为 z_{t} = dt+C.[例2] z_{t}= 2 z_{t-1}+1 .它的导出齐次方程为 z_{t}=2z_{t-1} , 特征根为\lambda=2 . 于是导出齐次方程的解为 z_t=C2^t.猜测原方程的一个特解为 z_{t} = 2^t-1 , 那么全部解为z_t=C2^t-1.。
差分方程

xk = ( − a ) x0 , k = 1, 2,L
k
所以当且仅当|a|<1时 方程( 所以当且仅当|a|<1时,方程(2)的平衡点 |a|<1 从而方程( 的平衡点)才是稳定的. (从而方程(1)的平衡点)才是稳定的.
常数矩阵A构成的 常数矩阵 对于n维向量 x ( k ) 和n×n常数矩阵 构成的 对于n (3) 方程组 x ( k + 1) + Ax ( k ) = 0 其平衡点稳定的条件是A的特征根 其平衡点稳定的条件是 的特征根
g 曲线斜率 y P3 f f P4 g P4 P3 K f < Kg K f > Kg y0 y0 P0 P0 y3 P2 P2 P1 y1 P1 0 x2 x x3 x1 x 0 x0 x 0
y y2
方程模型
yk = f (xk ) x k +1 = h ( y k )
在P0点附近用直线近似曲线
yk − y0 = −α ( xk − x0 ) (α > 0) xk +1 − x0 = β ( yk − y0 ) ( β > 0)
k +1
xk +1 − x0 = −αβ ( xk − x0 ) x
− x 0 = ( −αβ ) ( x1 − x 0 )
k
αβ < 1 (α <1/ β)
xk → x0 xk → ∞
= 1,故有解 an = 2 −1
n
1.3 差分方程的平衡点及稳定性 (1) 一阶线性方程的平衡点及稳定性 一阶线性常系数分方程
x k +1 + axk = b, k = 0,1,2,L
的平衡点由 x + ax = b 当
(完整版)差分方程的常见解法

(完整版)差分方程的常见解法差分方程的常见解法差分方程是数学中的一种重要方程类型,常用于描述离散事件系统的发展规律。
在求解差分方程时,我们可以采用以下几种常见的解法。
1. 直接求解法直接求解法是最简单且常用的差分方程求解方法之一。
它的基本思想是通过观察差分方程的规律,找到解的形式,并通过代入验证得到确切的解。
举例来说,对于一阶线性差分方程$y_{n+1} = ay_n + b$,我们可以猜测解的形式为$y_n = c\lambda^n$,其中$c$和$\lambda$为待定常数。
将此解代入方程,再通过已知条件解得$c$和$\lambda$的值,从而得到原差分方程的解。
2. 特征方程法特征方程法是一种常用于求解线性齐次差分方程的方法。
对于形如$y_{n+2} = ay_{n+1} + by_n$的差分方程,我们可以通过构造特征方程来求解。
具体步骤是,我们将差分方程中的项移动到一边,得到$y_{n+2} - ay_{n+1} - by_n = 0$。
然后,假设解的形式为$y_n =\lambda^n$,将其代入方程,得到特征方程$\lambda^2 - a\lambda - b = 0$。
解这个特征方程,得到特征根$\lambda_1$和$\lambda_2$,然后通解的形式为$y_n = c_1\lambda_1^n + c_2\lambda_2^n$,其中$c_1$和$c_2$为待定常数。
3. Z 变换法Z 变换法是一种广泛应用于差分方程求解的方法,特别适用于线性时不变差分方程。
该方法的基本思想是将差分方程转化为代数方程,并利用 Z 变换的性质求解。
对于差分方程$y_{n+1} = ay_n + b$,通过取 Z 变换,我们可以得到转化后的方程$Y(z) = azY(z) + b \frac{1}{1 - z^{-1}}$,其中$Y(z)$代表$y_n$的Z 变换。
然后,将方程整理,求解得到$Y(z)$,再通过反 Z 变换将其转换为差分方程的解$y_n$。
差分方程

§1 一阶常系数线性差分方程的求解形如1()0n n y ay f n ++=≡的方程为一阶常系数线性非齐次差分方程,其中a 为非零常数,()f n 为已知函数,n 为非负整数;10n n y ay ++=为对应的齐次方程。
1. 10n n y ay ++=的通解可以由以下两种方法给出:(1) 10n n y ay ++=对应的特征方程为0a λ+=,则a λ=-为特征根,从而其通解为()n n n y C C a λ==-,于是0C y =,即通解为0()n n y a y =- 。
(2) 设0y 已知,将0,1,2,n = 依次代入1n n y ay +=-中,得10()y a y =-,220()y a y =-, ,0()n n y a y =- 。
2. 设1()0n n y ay f n ++=≡(0)a ≠有一个特解n y ,则1()0n n y ay f n ++=≡的通解为()n n n y C a y =-+其中()n n y C a =-为对应齐次差分方程10n n y ay ++=的通解。
3. 关于1()0n n y ay f n ++=≡,针对不同的()f n ,其特解的求取方法: (1) 设()f n 为关于n 的m 次已知多项式()m P n ,则特解为()k n m y n R n =其中()m R n 为n 的m 次待定多项式。
若1a ≠-,即1λ≠是特征根,则0k =;若1a =-,即1λ=是特征根,则1k = 。
(2) 设()()nm f n P n q =(1)q ≠,其中q 为已知的常数,()m R n 为n 的m 次待定多项式,则特解为()k n n m y n R n q =当q 不是特征根时,取0k =;当q 是特征根时,取1k = 。
(3) 设12()cos sin f n b n b n ωω=+,则12(cos sin )k n y n B n B n ωω=+其中12,B B 为待定系数,当cos sin i ei ωωω=+a ≠-时,取0k =;当c o s s i n i e i aωωω=+=-时,取1k = 。
差分方程模型

设特解为 an D 代入 D 0.5D 0.1 得 D 0.2 , 于是所求通解 an c(0.5) n 0.2 例3 (养老金) 解: 齐次特征方程 设特解 an D
an1 1.01an 1000
1.01 0,
* an c(1.01) n.
代入原方程得 D 100000
例 4 求非齐次差分方程
* 对应齐次方程的通解为 an c1 2n c2 n 2n
的通解
f (n) 2 中, 2 是2 重根, 设特解为
n
an A n 2 2 n
n 2 n1
代入
得 A 1 2 方法2 (化齐) :
故通解为 an c1 2 c2 n 2 n 2
Fn Fn 1 Fn 2 F1 F2 1
解:差分方程的特征方程为 x 2 x 1 0 特征根
x1
n
1 5 1 5 , x2 2 2
n
1 5 1 5 Fn c1 c2 2 2
n
2(an1 4an2 4an3 ) 2 2n1 相减得 an 6an1 12an2 8an3 0 特征方程 3 62 12 8 0 特征根 2 为三重根, 通解为:
an 4an1 4an2 2n
an c1 2n c2 n 2n c3n 2 2n
x k b1 x k 1 b2 x k 2 bk 0
称为差分方程的特征方程,其根称为特征根。 定理1(单根)若特征方程恰有k个相异的特 x1 , x2 ,, x 征根 , k 则差分方程的通解为
an c x c x ck x