第三章 LTI离散系统的响应
合集下载
信号与线性系统分析-(第四版)第三章

(2) 特解 yp(k) p(2)k,k 0
p(2)k 4 p(2)k1 4 p(2)k2 2k
p 4 p(2)1 4 p(2)2 1
p
1 4
特解
yp
(k)
1 4
(2)k
(3) 全解
y(k
)
(C1k
C2
)(2)k
1 4
(2)k,k
0
根据初始条件
1 y(0) C2 4 0
1 y(1) 2C1 2C2 4 2 1
y(k) 4 y(k 1) 4 y(k 2) f (k) 已知初始条件y(0)=0,有y(1)= - 1,激励 f (k) 2k , k 0。
求方程的全解。
解: (1) 齐次解 特征方程
齐次解
2 4 4 0 特征根 1 2 2
yh(k) (C1k C2 )(2)k 代入差分方程
10cos(0.5 k)
P Q 1
yp (k) cos(0.5 k) sin(0.5 k)
2 cos(0.5 k )
4
y(k) yh (k) yp (k)
C1
1 2
k
C2
1 3
k
2 cos(0.5 k )
4
y(0) C1 C2
2 cos( ) 0
4
y(1) C1 C2 2 cos(0.5 ) 1
y(2) 3 y(1) 2 y(0) f (2) 2
y(3) 3y(2) 2y(1) f (3) 10
y(4) 3 y(3) 2 y(2) f (4) 10
便于计算机求解
二、差分方程的经典解
LTI系统的数学模型:n阶常系数线性差分方程
y(k) an1 y(k 1) a0 y(k n) bm f (k) bm1 f (k 1) b0 f (k m)
信号与系统第三章

y (4) 3 y (3) 2 y (2) f (4) 10 ...
特点:便于用计算机求解
2、差分方程的经典解
• 若单输入-单输出的LTI系统的激励为 f(k),全响应为y(k),则描述系统激 励与响应之间关系的数学模型是n阶 常系数线性差分方程,一般可写为:
a y (k i ) b
例3.1-1
• 解:将差分方程中除y(k)以外的各项都移到等 号右端,得
y(k ) 3 y(k 1) 2 y(k 2) f (k )
对k=2,将已知初始值y(0)=0,y(1)=2代入上式,得
y(2) 3 y(1) 2 y(0) f (2) 2
依次迭代可得 y (3) 3 y (2) 2 y (1) f (3) 10
位移单位序列:
运算:
• 加: (k) 2 (k) =3(k)
乘:(k) (k) (k)
延时:
0
取样性质:f (k)(k) f (0)(k)
2. 单位阶跃序列: (k)
(1)定义: (2)运算:
3) δ(k)与ε(k)的关系:
δ(k)=△ε(k)= ε(k)-ε(k-1) 差分表示,对应 的微分δ(t)=dε(t)/dt ε(k)=
第三章 离散系统的时域分析
连续系统与离散系统的比较
时域连续系统
f (t ) y(t )
常系数线性微分方程 卷积积分
时域离散系统
f (k ) y (k )
常系数线性差分方程 卷积和
y(t ) yzi (t ) yzs (t )
yzs (t ) f (t ) h(t )
y(k ) yzi (k ) yzs (k )
信号与线性系统分析第三章

系统描述 分析方法
连续系统 微分方程 卷积积分 变换域(傅氏、s) 系统函数
离散系统 差分方程 卷积和 变换域(离散傅氏、z) 系统函数
第 2页
§2.1 LTI离散系统的响应
• 差分与差分方程 —前向差分、后向差分以及差分方程
• 差分方程解 —数值解、经典解,以及不同特征根对应的齐 次解和不同激励对应的特解
yzi (-2) = y(-2)
-----------
yzi (n) = ?
----------------yzi (-n) = y(-n)
第 13 页
零输入举例
例1:系统方程为 y(k) + 3y(k –1) + 2y(k –2) = f(k) 已知激励f(k)=2k , k≥0;初始状态 y(–1)=0, y(–2)=1/2 求系统的零输入响应
解:yzi(k)零输入响应满足:
yzi(k) + 3yzi(k –1)+ 2yzi(k –2)= 0
yzi(–1)= y(–1)= 0 yzi(–2) = y(–2) = 1/2 递推求 yzi(0)、 yzi(1) yzi(k)= – 3yzi(k –1) –2yzi(k –2)
yzi(0)= –3yzi(–1) –2yzi(–2)= –1
yzs(0)、yzs(1)、---yzs(n)=? 借助微分方程
n
若其特征根均为单根: yzk (k ) Czsjkj y p (k ) j 1
第 16 页
零状态举例
例1:系统方程为 y(k) + 3y(k –1) + 2y(k –2) = f(k) 已知激励f(k)=2k , k≥0;求系统的零状态响应 解:零状态响应yzs(k) 满足
3_10离散时间LTI系统的单位脉冲响应

等效初始条件由差分方程和h[-1] = h[-2] = = h[-n] = 0 递推求出。
此方法称为等效初始条件法
2. 单位脉冲响应的求解
[例] 某离散因果LTI系统的差分方程为 y[k] 3y[k 1] 2y[k 2] x[k] 求系统的单位脉冲响应h[k]。
解:h[k]满足方程 h[k] 3h[k 1] 2h[k 2] d[k] (1) 确定h[k]的形式 特征方程为 r 2 3r 2 0 特征根为 r1 1, r2 2
h[k ] C1 (1) k C2 (2) k , k 0
2. 单位脉冲响应的求解
[例] 某离散因果LTI系统的差分方程为 y[k] 3y[k 1] 2y[k 2] x[k] 求系统的单位脉冲响应h[k]。
✓选择初始条件基本原则是必须将d[k]的作用体现在初始条件中 解:h[k]满足方程 h[k] 3h[k 1] 2h[k 2] d[k]
2. 单位脉冲响应的求解
[例] 某离散因果LTI系统的差分方程为 y[k] 3y[k 1] 2y[k 2] x[k] 求系统的单位脉冲响应h[k]。
分析: 如何确定系统的初始条件?
由于d [k]信号在k>0后函数值都为零。对于因果系统, 将d [k]对系统的瞬时作用转化为系统的等效初始条件。
北京交通大学 信号处理课程组
离散时间LTI系统的单位脉冲响应
单位脉冲响应的定义 单位脉冲响应的求解
1.பைடு நூலகம்单位脉冲响应的定义
在系统初始状态为零的条件下,以单位脉冲序列δ [k]激 励系统所产生的输出响应,称为系统的单位脉冲响应,以符 号h[k]表示。
d [k]
离散时间
此方法称为等效初始条件法
2. 单位脉冲响应的求解
[例] 某离散因果LTI系统的差分方程为 y[k] 3y[k 1] 2y[k 2] x[k] 求系统的单位脉冲响应h[k]。
解:h[k]满足方程 h[k] 3h[k 1] 2h[k 2] d[k] (1) 确定h[k]的形式 特征方程为 r 2 3r 2 0 特征根为 r1 1, r2 2
h[k ] C1 (1) k C2 (2) k , k 0
2. 单位脉冲响应的求解
[例] 某离散因果LTI系统的差分方程为 y[k] 3y[k 1] 2y[k 2] x[k] 求系统的单位脉冲响应h[k]。
✓选择初始条件基本原则是必须将d[k]的作用体现在初始条件中 解:h[k]满足方程 h[k] 3h[k 1] 2h[k 2] d[k]
2. 单位脉冲响应的求解
[例] 某离散因果LTI系统的差分方程为 y[k] 3y[k 1] 2y[k 2] x[k] 求系统的单位脉冲响应h[k]。
分析: 如何确定系统的初始条件?
由于d [k]信号在k>0后函数值都为零。对于因果系统, 将d [k]对系统的瞬时作用转化为系统的等效初始条件。
北京交通大学 信号处理课程组
离散时间LTI系统的单位脉冲响应
单位脉冲响应的定义 单位脉冲响应的求解
1.பைடு நூலகம்单位脉冲响应的定义
在系统初始状态为零的条件下,以单位脉冲序列δ [k]激 励系统所产生的输出响应,称为系统的单位脉冲响应,以符 号h[k]表示。
d [k]
离散时间
LTI离散系统的时域分析

二、差分方程的解
1、用迭代法求差分方程的数值解 差分方程是具有递推关系的代数方程,当已知 初始条件和激励时可以利用迭代法求得差分方 程的数值解 当差分方程阶次较低时可以使用此法
例3.1-1若描述某离散系统的差分方程为
y(k ) 3 y(k 1) 2 y(k 2) f (k )
3.有一对共轭复根λ1 、2=a+jb =ρe±jβ Yh()=ρk[Ccos(βk)+Dsin(βk)]
(或Aρkcos(βk-θ),其中Aejθ=C+jD )
y(k) + an-1y(k-1) +…+ a0y(k-n) = bmf(k)+…+ b0f(k-m)
特解yp(k):
表3-2典型激励对应的特解
已知初始条件y(0)=0,y(1)=2,激励f(k)=2k(k),求y(k) • 解:将差分方程中除y(k)以外的各项都移到等号 右端,得
y(k ) 3 y(k 1) 2 y(k 2) f (k )
对k=2,将已知初始值y(0)=0,y(1)=2代入上式,得
y(2) 3 y(1) 2 y(0) f (2) 2
k r (Pmk m Pm1k m1 P 1的特征根) 1k P 0 )(有r重为
cos k 或 sin k P1 cos k P2 sin k (特征根不等于 e j )
选定特解后代入原差分方程,求出待定系数就得 出方程的特解。 3)全解
j n bj ( 1 ) ( ,2....n j),j 0,1
2. 差分方程
包含未知序列y(k)及其各阶差分的方程式称为差 分方程。 将差分展开为移位序列,得一般形式 y(k) + an-1y(k-1) +…+ a0y(k-n) = bmf(k)+…+ b0f(k-m) 差分方程的阶数:未知序列最高与最低序数的差 上式各移位序列的系数均为常系数,即常系数差分 方程,用来描述LTI离散系统; 若系数是变量K的函数则为变系数差分方程
02-1 LTI离散系统响应的求解课件

特解
(1)齐次解:齐次差分方程的解。
y(k ) an1 y(k 1) L a0 y(k n) 0
(a)通一信阶方原 程理
y(k) a0 y(k 1) 0 y(k)
y(k 1) a0 y(k) C(a0 )k
特征方程: a0 0 特征根: a0 差分方程齐次解形式: Ck
解:(1)求齐次解。 2 4 4 0
可解得特征根 1 2 2 为二重根,齐次解形式为:
y (k ) [C k C ](2yp (k )
P
2 k
,
k
0
将 yp (k ), yp (k 1), yp (k 2) 代入微分方程中得
P 2k 4P 2k 1 4P 2k 2 2k 4P 1
y(k):响应信号
初始状态: y(-1) ,y(-2),…, y(-n)
二、L通T信I离原散理系统响应的求解 1、迭代法:差分方程是递推的代数方程,若已知初
始条件和激励,利用迭代法可求得其数值解。
例:若描述某系统的差分方程为: y(k ) 3 y(k 1) 2 y(k 2) f (k)
已知初始条件 y(0)=0, y(1)=2, 激励 f (k) 2k(k) , 求 y(k)。 解: y(k ) 3 y(k 1) 2 y(k 2) f (k)
4
4
通信原理
通信原理
LTI离散系统响应的求解
主讲人:曹红梅 通信与信息工程学院
1
1
一、L通T信I离原散理系统的描述
LTI离散系统描述方法:n阶线性常系数差分方程
y(k ) an1 y(k 1) L a0 y(k n) bm f (k ) bm1 f (k 1) L b0 f (k m)
f(k):激励信号
离散时间LTI系统零状态响应

x[k]* h[n] ( x[n])* h[k] y[n]
n
n
n
3. 卷积和的性质
[例]
计算
x[k]
{1,
0,
2,
4}与
h[k]
{1,
4,
5,
3}的卷积和。
解: x[k] [k 2] 2[k] 4[k 1]
利用分配律和位移特性
x[k]*h[k] {[k 2] 2[k] 4[k 1]}*h[k]
x[1] h[1]
x[1] h[2]
x[1] h[3]
对角斜线上各数值就 是 x[n]h[k-n]的值。
x[ 2 ] x[2] h[0]
x[ 3 ] x[3] h[0]
x[2] h[1] x[3] h[1]
x[2] h[2] x[3] h[2]
x[2] h[3] x[3] h[3]
对角斜线上各数值的 和就是y[k]各项的值。
h[k] h [ n ] 1
解:
n
n
0
3
k
02
k
y[k]
综上可知: 3
2 1
k 0 1 2 3 4 56
h[-n] 1
-2 0
n
两个信号的卷积和,卷积和结果仍为一个信号。该信号的起点 等于那两个信号起点之和,终点等于那两个信号的终点之和。
2. 卷积和的计算
列表法: 设x[k]和h[k]都是因果序列,则有
2. 卷积和的计算
[例]
计算x [ k ]
{1,
2,
0,
3,
2}
与h [ k ]
{1,
4,
2,
3 } 的卷积和。
解:
h [ -1 ] 1
离散时间LTI系统的单位脉冲响应

系统分析和设计
通过单位脉冲响应可以分析系统 的稳定性、频率响应和因果性等 特性,用于系统的设计和优化。
信号处理
单位脉冲响应可以用于信号的滤 波、预测和合成等处理,提高信 号的质量和性能。
控制工程
单位脉冲响应可以用于控制系统 的分析和设计,优化控制性能和 稳定性。
BIG DATA EMPOWERS TO CREATE A NEW ERA
IIR系统
系统的输出不仅与当前的输入有关, 还与过去的输入有关,因此其单位脉 冲响应在时间上是无限的。
系统的表示方法
差分方程
离散时间LTI系统的动态行为通常由差分方 程描述,如 $y(n) = f(n) + g(n)u(n)$。
传递函数
通过将差分方程转换为传递函数的形式,可以更方 便地分析系统的频率响应和稳定性。
仿真分析的步骤与过程
建立数学模型
根据系统定义,建立离散时间LTI系统的数学模型,包括差分方程或传递函数。
生成单位脉冲信号
在仿真中,生成一个单位脉冲信号,用于输入到离散时间LTI系统中。
计算单位脉冲响应
将单位脉冲信号输入到系统中,并记录系统的输出,即单位脉冲响应。
分析单位脉冲响应
对单位脉冲响应进行分析,包括幅度和相位特性,以及稳定性等。
性质
单位脉冲响应是线性时不变系统的内 部动态特性,具有稳定性、因果性和 可预测性。
单位脉冲响应的求解方法
直接法
根据系统函数或差分方程,直接计算单位脉冲响 应的数值解。
迭代法
根据系统函数或差分方程,通过迭代计算单位脉 冲响应的数值解。
逆系统法
通过求解系统的逆系统,得到单位脉冲响应的数 值解。
单位脉冲响应的应用
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
yzi(k)= – 3yzi(k –1) –2yzi(k –2)
yzi(0)= –3yzi(–1) –2yzi(–2)= –1
yzi(1)= –3yzi(0) –2yzi(–1)=3 方程的特征根为λ1= –1 ,λ2= – 2 , 其解为 yzi(k)=Czi1(– 1)k+Czi2(–2)k 将初始值代入 并解得
3.1 LTI离散系统的响应
例3:若描述某系统的差分方程为 y(k)+ 4y(k – 1) + 4y(k – 2) = f(k) 已知初始条件y(0)=0,y(1)= – 1;激励f(k)=2k,k≥0。求全解。
解: 特征方程为 λ2 + 4λ+ 4=0 可解得特征根λ1=λ2= – 2,其齐次解为 yh(k)=(C1k +C2) (– 2)k 特解为 yp(k)=P (2)k , k≥0 代入差分方程得 P(2)k+4P(2)k –1+4P(2)k–2= f(k) = 2k ,P=1/4 特解 yp(k)=2k–2 , k≥0
3.1 LTI离散系统的响应 二、差分方程的经典解 y(k) + an-1y(k-1) +…+ a0y(k-n) = bmf(k)+…+ b0f(k-m) 与微分方程经典解类似,y(k) = yh(k) + yp(k) 1. 齐次解yh(k)
例2:一阶齐次方程的解 y(k 1) a0 y(k ) 0
yh(t)(自由响应) yp(t) (强迫响应)
当 1, 则自由响应是衰减变化 ,系统稳定。
1, 则自由响应是增长变化 ,系统不稳定。
3.1 LTI离散系统的响应
三、零输入响应和零状态响应 y(k) = yzi (k) + yzs(k) 设激励f(k)在k=0时接入系统, 初始状态:y(–1), y(–2) , …,y(–n)
y(k )
非移变系统 If
f (k )
then
f (k i ) y (k i )
3.1 LTI离散系统的响应 一、差分与差分方程 设有序列f(k),则…,f(k+2),f(k+1),…,f(k-1),f(k-2)…等 称为f(k)的移位序列。 1. 差分运算
d f (t ) f (t ) f (t t ) f (t ) f (t ) f (t t ) lim lim lim t 0 t 0 t 0 dt t t t
3.1 LTI离散系统的响应 求解常系数线性差分方程的方法一般有以下几种
1、迭代法 2、时域经典法
逐次代入求解, 概念清楚, 比较简便, 适用于计算机, 缺点是不能得出通式解答。
全响应=齐次通解 自由响应
+
特解
强迫响应
求解过程比较麻烦, 不宜采用。 3、全响应=零输入响应+零状态响应 零输入响应求解与齐次通解方法相同。 零状态响应求解利用卷积和法求解,十分重要。 4、变换域法(Z变换法)
y (k 1) a0 y (k )
y (k )
是个公比为 a0的级数
y(k ) c(a0 )k
y(0) c(a0 )0 c
c是待定常数,有初始条件决定
3.1 LTI离散系统的响应
齐次方程 y(k) + an-1y(k-1) + … + a0y(k-n) = 0 其特征方程为 λn + an-1λn– 1 + … + a0 = 0 其根λi( i = 1,2,…,n)称为差分方程的特征根。
i 1 自由响应 强迫响应 i 1 零输入响应 i 1 零状态响应
Ci i Czi i Czsi i
k k i 1 i 1 i 1
n
n
n
k
Czii ——仅由系统的初始状态决定。
Ci ——由初始状态和激励决定。
3.2 单位序列响应和阶跃响应
不同特征根所对应的齐次解 特征根 单实根 齐次解
y (k )
j
ck
cr 1k r 1k cr 1k r 2 k c1kk c0 k
p k [c cos k D sin k ] 或Apk cos(k )
r 重实根
1,2 a jb pe
k
f (k ) (k i) f (i)
3.2 单位序列响应和阶跃响应
( 2)单 位 阶跃 序 列 1 k 0 (k ) 0 k 0 (k )
移位单位阶跃序列 (k i ) 1 k i 0 k i
(k 2)
1
1
0
1 2 3
k
k
0
1 2 3 4 5
k
显然: (k ) (k ) (k ) (k 1)
(k )
i
(i) (k j)
j 0
3.2 单位序列响应和阶跃响应
3.2 单位序列响应和阶跃响应
二、单位序列响应 由单位序列δ(k)所引起的零状态响应称为单位序列响应或 单位样值响应或单位取样响应,或简称单位响应,记为h(k)。 方法一:若方程右端只有f(k),而无移位项---经典法。 例1: 求如图所示
y(k ) yzi (k ) yzs (k )
零输入响应
n k
零状态响应
n k n k
y(k ) yzi (k ) yzs (k ) Ci i y p (k ) Czii i Czsi i y p (k )
一、常用序列
(1)单位脉冲序列 定义: (k )
( k )
1 k 0 0 k0
移位单位脉冲序列 定义: (k i )
( k i )
1 k i 0 k i
1
1 1 2
0
k
0
1 2 i
k i
k
(k )的性质:f (k) (k i) f (i) (k i) f (i)
yzi(j)和yzs(j) ( j = 0, 1, 2 , … ,n – 1)
3.1 LTI离散系统的响应 例4:若描述某离散系统的差分方程为 y(k) + 3y(k –1) + 2y(k –2) = f(k) 已知激励f(k)=2k , k≥0,初始状态y(–1)=0, y(–2)=1/2, 求系统的 零输入响应、零状态响应和全响应。 解:(1)yzi(k)满足方程 yzi(k) + 3yzi(k –1)+ 2yzi(k –2)= 0 其初始状态 yzi(–1)= y(–1)= 0, yzi(–2) = y(–2) = 1/2 首先递推求出初始值yzi(0), yzi(1), yzi(k)= – 3yzi(k –1) –2yzi(k –2)
Czi1=1 , Czi2= – 2
所以 yzi(k)=(– 1)k – 2(– 2)k , k≥0
3.1 LTI离散系统的响应 (2)零状态响应yzs(k) 满足 yzs(k) + 3yzs(k –1) + 2yzs(k –2) = f(k) yzs(–1)= yzs(–2) = 0 递推求初始值 yzs(0), yzs(1), yzs(k) = – 3yzs(k –1) – 2yzs(k –2) + 2k , k≥0 yzs(0) = – 3yzs(–1) – 2yzs(–2) + 1 = 1 yzs(1) = – 3yzs(0) – 2yzs(–1) + 2 = – 1 分别求出齐次解和特解,得 yzs(k) = Czs1(–1)k + Czs2(–2)k + yp(k) = Czs1(– 1)k + Czs2(– 2)k + (1/3)2k 代入初始值求得 Czs1= – 1/3 , Czs2=1 所以 yzs(k)= – (– 1)k/3+ (– 2)k + (1/3)2k , k≥0
1
r重等于1的特征根
a 特征根 a 特征单根
pak
a
k
p1ka k p0 a k
Pcos(βk)+Qsin(βk)
pr k r a k pr 1k r 1a k p1ka k p0 a k a r 重特征根
cos(βk)或
sin(βk)
所有特征根均不等于e±jβ
第三章:LTI离散系统的时域分析
Chapter3
Discrete systems
本章要点
F F
LTI离散时间系统的响应 单位序列响应和阶跃响应 卷积和
F
引言
f (k )
激励是离散
时间信号 离散系统
响应是离散
时间信号
y (k )
连续时间系统 微分方程 卷积积分 拉氏变换 连续傅立叶变换 卷积定理 连续系统模拟
y (k) f (k) ∑ D 1 D 2
3.1 LTI离散系统的响应 2. 差分方程
包含未知序列y(k)及其各阶差分的方程式称为差分方程。 将差分展开为移位序列,得一般形式 y(k) + an-1y(k-1) +…+ a0y(k-n) = bmf(k)+…+ b0f(k-m)
例1:若描述某系统的差分方程为 y(k) + 3y(k – 1) + 2y(k – 2) = f(k) 已知初始条件y(0)=0, y(1)=2, 激励f(k)=2kε(k), 求y(k)。 解: y(k) = – 3y(k – 1) – 2y(k – 2) + f(k) y(2)= – 3y(1) – 2y(0) + f(2) = – 2 y(3)= – 3y(2) – 2y(1) + f(3) = 10 …… 一般不易得到解析形式的(闭合)解。