山东大学 信息科学工程学院5-3 常系数线性差分方程的时域求解法
常系数线性差分方程的解
常系数线性差分方程的解 方程)(...110n b x a x a x a n k k n kn =+++-++（1）其中k a a a ,...,,10为常数，称方程（1）为常系数线性方程。
又称方程0...110=+++-++n k k n kn x a x a x a（2）为方程（1）对应的齐次方程。
如果（2）有形如nnx λ=的解，带入方程中可得：0 (11)10=++++--k k k k a a a a λλλ（3）称方程（3）为方程（1）、（2）的特征方程。
显然，如果能求出（3）的根，则可以得到（2）的解。
基本结果如下：（1） 若（3）有k 个不同的实根，则（2）有通解：nkk nnn c c c x λλλ+++=...2211，（2） 若（3）有m 重根λ，则通解中有构成项：nm m nc n c c λ)...(121----+++（3）若（3）有一对单复根βαλi ±=，令：ϕρλi e±=，αβϕβαρarctan,22=+=，则（9）的通解中有构成项：nc n c nnϕρϕρsin cos 21--+（4） 若有m 重复根：βαλi ±=，φρλi e±=，则（2）的通项中有构成项：n nc n c c n nc n c c nm m m m nm m ϕρϕρs i n )...(c o s )...(1221121---++---+++++++综上所述，由于方程（10）恰有k 个根，从而构成方程（2）的通解中必有k 个独立的任意常数。
通解可记为：-n x 如果能得到方程（1）的一个特解：*n x ，则（1）必有通解：=n x -nx +*n x （4）特解可通过待定系数法来确定。
差分方程的法
差分方程常用解法1、 常系数线性差分方程的解方程)(...110n b x a x a x a n k k n k n =+++-++ （1）其中k a a a ,...,,10为常数，称方程（1）为常系数线性方程。
又称方程0...110=+++-++n k k n k n x a x a x a （2）为方程（1）对应的齐次方程。
如果（2）有形如n n x λ=的解，代入方程中可得： 0...1110=++++--k k k k a a a a λλλ （3） 称方程（3）为方程（1）、（2）的特征方程。
显然，如果能求出方程（3）的根，则可以得到方程（2）的解。
基本结果如下：（1） 若（3）有k 个不同的实根，则（2）有通解：nk k n n n c c c x λλλ+++=...2211，（2） 若（3）有m 重根λ(即m 个根均为λ)，则通解中有构成项：n m m n c n c c λ)...(121----+++（3）若（3）有一对单复根 βαλi ±=，令：ϕρλi e ±=，αβϕβαρarctan,22=+=，则（2）的通解中有构成项： n c n c n nϕρϕρsin cos 21--+ （4） 若有m 重复根：βαλi ±=，φρλi e ±=，则（2）的通项中有构成项:n n c n c c n n c n c c n m m m m n m m ϕρϕρsin )...(cos )...(1221121---++---+++++++综上所述，由于方程（3）恰有k 个根，从而构成方程（2）的通解中必有k 个独立的任意常数。
通解可记为：-n x如果能得到方程（1）的一个特解：*n x ，则（1）必有通解： =n x -n x +*n x （4） 方程（4） 的特解可通过待定系数法来确定。
例如：如果)(),()(n p n p b n b m m n =为n 的m 次多项式，则当b 不是特征根时，可设成形如)(n q b m n 形式的特解，其中)(n q m 为n 的m 次多项式；如果b 是r 重特征根时，可设特解：r n n b )(n q m ，将其代入（1）中确定出系数即可。
差分方程的解法分析及MATLAB实现
差分方程的解法分析及MATLAB实现差分方程是描述离散时序系统行为的数学工具。
在离散时间点上，系统的行为由差分方程给出，这是一个递归方程，其中当前时间点的状态取决于之前的状态和其他外部因素。
解差分方程的方法可以分为两类：直接解法和转化为代数方程的解法。
直接解法通过求解差分方程的递归形式来得到解析或数值解。
转化为代数方程的解法则将差分方程转化为代数方程进行求解。
一、直接解法的步骤如下：1.将差分方程表示为递归形式，即将当前时间点的状态表示为之前时间点的状态和其他外部因素的函数。
2.根据初始条件，确定初始时间点的状态。
3.根据递归形式，计算出后续时间点的状态。
以下是一个简单的差分方程的例子：y(n)=2y(n-1)+1，其中n为时间点。
按照上述步骤求解该差分方程：1.将差分方程表示为递归形式：y(n)=2y(n-1)+12.根据初始条件，假设y(0)=1，确定初始时间点的状态。
3.根据递归形式，计算出后续时间点的状态：y(1)=2y(0)+1=2*1+1=3y(2)=2y(1)+1=2*3+1=7y(3)=2y(2)+1=2*7+1=15...依此类推计算出所有时间点的状态。
二、转化为代数方程的解法的步骤如下：1.假设差分方程的解具有指数形式，即y=r^n，其中r为待定参数。
2.将差分方程代入上述假设中，得到r的方程。
3.解得r的值后，再根据初始条件求解出常数值。
4.得到差分方程的解析解。
以下是一个复杂一些的差分方程的例子：y(n)=2y(n-1)+3y(n-2)，其中y(0)=1，y(1)=2按照上述步骤求解该差分方程：1.假设差分方程的解具有指数形式：y=r^n。
2.代入差分方程得到：r^n=2r^(n-1)+3r^(n-2)。
3.整理得到：r^2-2r-3=0。
4.解得r的值为：r1=-1，r2=35.根据初始条件求解出常数值：y(0)=c1+c2=1，y(1)=c1-c2=2、解得c1=1.5，c2=-0.56.得到差分方程的解析解：y(n)=1.5*(-1)^n+-0.5*3^n。
1.3 常系数线性差分方程
同一个差分方程，边界条件不同，所求的h(n)表达 式不同。即：
同一个差分方程，边界条件不同，其对应的系统 是不同的。
二、常系数线性差分方程的求解
解得：此系统不是线性系统，也不是移不变系统。 结论：常系数线性差分方程，其所对应的
系统并不一定是线性移不变的。
一些关于差分方程的结论
一个差分方程不能唯一确定一个系统 常系数线性差分方程描述的系统不一定
利用查找表实现4bit x 4bit
4bit x 4bit 乘法器
用D触发器实现延时器
1. 己知差分方程，作出系统运算结构
2. 己知系统运算结构，求差分方程表达式
例：已知某系统结构如下所示，求此系统所对应的 差分方程。
四、系统运算结构的实现
当输入x(n)=nR10(n)，求输出y(n)。
(输入和输出信号均为8件编程
课后自训
某线性移不变离散时间系统的单位抽样响应序 列h(n)=R3(n)，
1.求此系统对应的差分方程； 2.作出此系统的运算结构； 3. 分别用硬件电路和软件编程实现此系统，
当输入信号x(n)＝nR10(n)时，求出输出 信号y(n)。（输入输出信号均用8bit表示）
1.3 常系数线性差分方程
(3) y(n) x(n) h(n) x(m)h(n m) m
(使用3之前要证明此系统是线性移不变系统)
一、常系数线性差分方程的定义
二、常系数线性差分方程的求解
二、常系数线性差分方程的求解
二、常系数线性差分方程的求解
解得：
二、常系数线性差分方程的求解
解得：
是线性移不变的 不一定是因果的 不一定是稳定的
在今后的讨论中，通常假设常系数 线性差分方程就代表线性移不变系统， 且多数代表可实现的因果系统。
常系数线性差分方程的求解28页PPT
1、战鼓一响，法律无声。——英国 2、任何法律的根本；不，不成文法本 身就是பைடு நூலகம்讲道理 ……法 律，也 ----即 明示道 理。— —爱·科 克
3、法律是最保险的头盔。——爱·科 克 4、一个国家如果纲纪不正，其国风一 定颓败 。—— 塞内加 5、法律不能使人人平等，但是在法律 面前人 人是平 等的。 ——波 洛克
31、只有永远躺在泥坑里的人，才不会再掉进坑里。——黑格尔 32、希望的灯一旦熄灭，生活刹那间变成了一片黑暗。——普列姆昌德 33、希望是人生的乳母。——科策布 34、形成天才的决定因素应该是勤奋。——郭沫若 35、学到很多东西的诀窍，就是一下子不要学很多。——洛克
3常系数线性差分方程
该“线性”与线性系统的“线性”含义不 同
2、常系数差分方程的求解：
① 经典解法：类似于模拟系统求解微分方程的方法，要求 齐次解、特解，并由边界条件求待定系数。 由于计算复杂，较少使用。 ② 递推(迭代)法：简单、适于用计算机进行求解。但只能 得到一系列数值解，不易得到封闭式(公 式)解答。 ③ 变换域法：将差分方程变换到z域求解。 ④ 卷积法：由差分方程求出系统的h(n)，再与已知的x(n) 进行卷积，得到y(n)。
前面已经证明： 当 x1(n)=(n) 时，y1(n) = anu(n) 当 x2(n)=(n-1) 时， y2(n) = anu(n)+ an-1u(n-1) 令：x3(n)=(n)+(n-1)， y3(0)=1 y3(1) = ay3(0)+x3(1) = a+1 y3(2) = ay3(1)+x3(2) = a2+a … y3(n) = ay3(n-1)+x3(n) = an+an-1 ∴ y3(n) = anu(n)+ an-1u(n-1) ∵ 当x3(n)=x1(n)+x2(n)时，y3(n)≠y1(n)+y2(n)， 所以，该系统也不是线性系统。
B、当边界条件为y(0)=0时，为线性、移变系统 C、当边界条件为y(-1)=0时，为线性、移不变系统 证：(这里只证明A，B和C留给大家课后思考证明。)
令：x1(n)=(n)， y1(0)=1 y1(1) = ay1(0)+x1(1) = a y1(2) = ay1(1)+x1(2) = a2 … y1(n) = ay1(n-1)+x1(n) = an ∴ y1(n) = anu(n)
§1.3 常系数线性差分方程
§7.4 常系数线性差分方程的求解
(
) u(n)
2、若把初值y(0)=1,看作激励加入后系统的初始样值y+(0)， 若把初值y(0)=1,看作激励加入后系统的初始样值 (0)， 始样值y (0)=1应满足方程 应满足方程： 则y+(0)=1应满足方程： y(n)-3y(n-1)= u(n) <0时 迭代法得： 当n<0时，由迭代法得： y+(n)=0 当n ≥ 时，则有： 0 则有： y+(0)= 1 +3y y+(1)= u(1) +3y+(0)=1+3*1=4
y − (− 1) = 1 1 y − (0 ) = 3 3 2 1 1 y − (− 2 ) = y − (− 1) = 3 3
…...
1 1 y − (n ) = y − (n + 1) = 3 3
−n
假设系统是因果系统， 假设系统是因果系统， 由于激励u n=0 由于激励u(n)在n=0接 那么,此解就是n 入，那么,此解就是n<0 时系统的零输入响应。 时系统的零输入响应。
如果系统起始样值 如果系统起始样值y-(n) ≠ 0，则系统差分方程的完全 起始样值y 0，则系统差分方程的完全 解将不满足线性时不变的特性。 解将不满足线性时不变的特性。 今后我们规定，所有初值如无下标 值如无下标， 今后我们规定，所有初值如无下标，则一律按初始 样值处理。 样值处理。
返回
种方法） 二、差分方程的解法（前3种方法） 差分方程的解法（
y+(2)= u(2) +3y+(1)=1+3+32=13 +3y …... 1 2+……+3n = (3 n +1 − 1) y+(n)= u(n) +3y+(n-1)=1+3+3 +3y 2 1 n +1 则方程的解为： 则方程的解为： y(n)= (3 − 1) u(n)
常系数线性差分方程
同步骤1），由
y3 (n
1)
1 a
[
y3 (n)
x3 (n )]
得y3(n) an1，n 1
y3(n) (1 a) (n) (1 a a2 )an1u(n 1)
an1u(n 1)
4）结论：
当输入x1(n) (n)时，输出
y1(n) (1 a)anu(n) an1u(n 1) 当输入x2 (n) (n 1)时，输出
y2 (n) an1(a2 1)，n 1
同步骤1），由
y2 (n
1)
1 a
[
y2
(n)
x2 (n)]
得y2 (n) an1，n 1
y2 (n) a (n) (1 a2 )an1u(n 1) an1u(n 1)
3）令输入x3(n) x1(n) x2 (n) (n) (n 1)，
例1：已知常系数线性差分方程
y(n) ay(n 1) x(n)
若边界条件
y(1) 0
求其单位抽样响应。
解：令输入x(n) (n)，则输出y(n) h(n)，
又已知y(1) 0
由y(n) ay(n 1) x(n)，得
y(0) ay(1) x(0) 1 y(1) ay(0) x(1) a y(2) ay(1) x(2) a2
y1(n) an (a 1)，n 0
由y1 (n
1)
1 a
[
y1(n)
x1(n)]，得
y1 ( 2)
1 a
[
y1(1)
x1(1)]
a 1
y1 ( 3)
1 a
[
y1(2)
x1(2)]
a 2
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
y zi ( n ) y h ( n )
由初始状态
后向差分方程： 前向差分方程： y ( 1 ), y ( 2 ), , y ( N ) y ( 0 ), y ( 1 ), , y ( N 1 )
形同
可确定齐次解中的待定系数
)=0,y(1)=2,激励x(n)=2nu(n),求y(n) 解： y(n) = – 3y(n – 1) – 2y(n – 2) + x(n)
y(2)= – 3y(1) – 2y(0) + x(2) = – 2
y(3)= – 3y(2) – 2y(1) + x(3) = 8 …
二、经典解法
§５.3
常系数线性差分方程 的时域解
迭代法、时域经典法、时域零输入零状态法、变换域法
一、迭代法
特点：概念清晰 方法简单 易于编程 但一般不易得到差分方程解析形式的(闭合)解
差分方程本质上是递推的代数方程，若已知边 界条件和激励，利用迭代法可求得其数值解。 例：若描述某系统的差分方程为
y(n) + 3y(n – 1) + 2y(n – 2) = x(n)
k
a 不等于特征根 a是一特征单根 a是一个r 阶特征根
k
所有特征根均不等于1 k r ( B 0 B 1 n B k n ) n 有r 阶等于1的特征根
B 0 sin n B 1 cos n
例题：见讲义p35 例5.3-1
或 cos n
三、零输入响应
由定义推知：
见讲义 p36 5.5.3
与微分方程经典解类似， y(n) = yh(n) + yp(n) 齐次解 齐次解
（见讲义 p33）
特解
特解： 线性时不变系统输入与输出有相同的形式
输入
常数
Ba
n
特解
B
( B 0 B 1 n )a
n
r n
a n
sin n
n
( B 0 B 1 n B r n )a
B 0 B1 n B k n