002第2章ch_2_2离散时间系统分析end_(修复的)_(自动保存的)

合集下载
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

2.2离散系统的时域分析
离散系统分析与连续系统分析在于多方面是相互平行的,有许多类似之处。

连续系统可用微分方程描述,离散系统可用差分方程描述。

差分方程和微分方程的许多求解方法在许多方面是相互对应的。

在连续系统分析中,卷积积分具有重要的意义;在离散系统中,卷积和也具有同等重要的地位。

连续系统分析与离散系统分析的相似性为读者学习本章节提供了有力条件,不过,读者应该十分注意他们之间存在着的重要差异。

在离散系统分析中,激励(输入)用()f n 表示,响应(输出)用()y n 表示,其中n 为整数;初始状态用{}0()x n 表示,其中0n 为正常数,通常取00n =。

下面,从离散系统的差分方程(或系统框图)及其求解开始,研究LTI 离散系统的时域分析。

2.2.1 LTI 离散系统的响应 差分与差分方程
与连续时间信号的微分与积分运算相对应,离散时间信号有差分及序列求和运算。

设有序列()f n ,则称
(2)f n +,(1)f n +,
(1)f n -,(2)f n -等为()f n 的位移序
列。

序列的差分可分为前向差分和后向差分。

一阶前向差分定义为
()
(1)()def
f n f n f n ∆+-
(2.2-1)
一阶后向差分定义为
()
()(1)def
f n f n f n ∇--
(2.2-2)
式中∆和∇称为差分算子。

由式(2.2-1)和(2.2-2)可见,前向差分与后向差分的关系为
()(1)f n f n ∇=∆-
(2.2-3)
二者仅位移不同,没有原则的差别,因而它们的性质也相同。

本书主要采用后向差分,并简称其为差分。

由差分的定义,若有序列1()f n 、2()f n 和常数1a 、2a ,则
112211221122[()()][()()][(1)(2)]a f n a f n a f n a f n a f n a f n ∇+=+--+- 111222[()(1)][()(1)]a f n f n a f n f n =--+-- 1122()()
a f n a f n =∇+∇
(2.1-4)
这表明差分运算具有线性性质。

二阶差分可定义为
2(n)[()][()(1)]()(1)def
f f n f n f n f n f n ∇∇∇=∇--=∇-∇-
()2(1)(2)f n f n f n =--+-
(2.2-5)
类似的,可定义三阶、四阶、五阶···差分。

一般的,k 阶差分
k
1
0()[()](1)()
k
def
k j j k f n f n f n j j -=⎛⎫
∇∇∇
=-- ⎪⎝⎭∑
(2.2-6)
式中
!
,0,1,2,...()!!
k k j k
j k j j ⎛⎫== ⎪
-⎝⎭
(2.2-7)
为二项式系数。

序列()f n 的求和运算为
()n
i f i =-∞

(2.2-8)
差分方程是包含关于变量n 的未知序列()f n 及各阶差分方程的方程式,它的一般形式可写为
[,(),(),
(),(),
,()]0r N F n f n f n f n y n y n ∇∇∇∇=
(2.2-9a )
式中差分的最高阶为n 阶,称为n 阶差分方程。

由式(2.2-6)可知,各阶差分均可写成
(n)y 及其各位移序列的线性组合,故上式常写为
[,(),(1),
(),(),(1),
,()]0F n f n f n f n r y n y n y n N ----=(2.2-9b )
通常所说的差分方程是指式(2.2-9b )形式的方程。

差分方程是具有递推关系的代数方程,若已知初始条件和激励,利用迭代法可求得差分方程的数值解。

例2.2-1 若描述某离散系统的差分方程为
()3(1)2(2)()y n y n y n f n +-+-=
已知初始条件(0)0,(1)2,y y ==激励()2()n
f n n ε=,求()y n 。

解:将差分方程中除()y n 以外的各项都移到等号右端,得
()3(1)2(2)()y n y n y n f n =----+
对于2,n =将已知初始值(0)0,(1)2y y ==代入上式,得
(2)3(1)2(0)(2)2y y y f =--+=-
类似地,一次迭代可得
(3)3(2)2(1)(3)10y y y f =--+=
(4)3(3)2(2)(4)10
y y y f =--+=-
由上例可见,用迭代法求解差分方程思路清楚,便于编写计算机程序,能得到方程的数值解,但它常常不宜得出解析形式(或称闭式)的解。

2.2.2差分方程的经典解
一般而言,如果单输入-单输出的LTI 系统的激励为()f n ,其全响应为()y n ,那么,描述该系统激励()f n 与响应()y n 之间的关系的数学模型是n 阶常系数线性差分方程,它
可以写为
10()(1)()N y n a y n a y n N -+-++-
10()(1)()M M b f n b f n b f n M -=+-+
+-
(2.2-10a )
式中(0,1,2,
,1)i a i N =-,(0,1,
,)j b j M =都是常数。

上式可缩写为
N
M
N M 00
()()i j i j a y n i b f n j --==-=-∑∑ (式中N 1a =) (2.2-10b )
与微分方程的经典解相类似,上述差分方程的解由齐次节和特解两部分组成。

齐次解用()h y n 表示,特解用()p y n 表示。


(n)()()
h p y y n y n =+
(2.2-11)
齐次解
当式(2.2-10)中的(n)f 及其各位移项均为零时,齐次方程
10()(1)()0
n y n a y n a y n N -+-++-=
(2.2-12)
(2.2-12)的解称为齐次解。

首先分析最简单的一阶差分方程,若一阶差分方程的齐次方程为
()(1)0y n ay n +-=
(2.2-13)
它可改写为
()
(1)
y n a y n =--
()y n 与(1)y n -之比等于()a -表明,序列(n)y 是一个公比为()a -的等比序列,因此()y n 应有如下形式
n ()()y n C a =-
(2.2-14)
式中C 为常数,由初始条件确定。

对于n 阶齐次差分方程,它的齐次解由形式为n C λ的序列组合而成,将n
C λ代入到时
(2.2-12),得
111100n n n N n N n C a C a C a C λλλλ-----++++=
由于C 0≠,消去C ;且0λ≠,以n N λ-除上式,得
11100N N n a a a λλλ--++++=
(2.2-15)
上式称为差分方程(2.1-10)和(2.1-12)的特征方程,它有n 个根(1,2,
,)i i n λ=,称
为差分方程的特征根。

显然,形式为C n
i i λ的序列都满足式(2.2-12),因而它们是式(2.1-10)方程的解。

依特征根的不同取值,差分方程齐次解的形式如表2.2-1,其中C i 、
D i 、A i 、i θ等为待定常数。

表2.2-1 不同特征根所对应的齐次解
特征根λ 齐次解
单实根
n C λ
r 重实根
121210C C C C r n r n n n r r k k n λλλλ----++
++
一对共轭复根
1,2j a jb e βλρ±=+=
[cos()sin()]n C n D n ρββ+或cos()n A n ρβθ-,其中
j Ae C jD θ=+
特解
特解的函数形式与激励的函数形式有关,表2.2-2列出了几种典型的激励()f n 所对应的特解的(n)p y 。

选定特解后代入原差分方程,求出其待定系数i P (或A 、θ)等,就得出方程的特解。

表2.2-2 不同激励所对应的齐次解
激励(n)f
齐次解
n m
111n m m m m P P n Pn P --++++所有特征根均不等于1时;
1110[n ]r m m m m n P n P n P P --++
++当有r 重等于1的特征根
时。

n a
n Pa ,当a 不等于特征根时;
1
10n m n P a
P a -+当a 是特征根时; 111
0r n r n n
n r r P n a P n a Pna P a --++
++,
当a 是r 重特征根
时。

cos()n β或sin()n β Pcos()sin()n Q n ββ+当所有的特征根均不等于j e θ±
或Acos()n βθ-,其中=P+j j Ae Q θ
全解
式(2.2-10)的线性差分方程的完全解是其齐次解与特解之和。

如果方程的特征根均为单根,则差分方程的全解为
1
()()()()N
n h p i i p i y n y n y n C y n λ==+=+∑ (2.2-16)
如果特征根1λ为r 重根,而其余()n r -个特征根为单根时,差分方程的全解为
1
1
1
()()r
N
r i
n
n i j
j p i j r y n C n C y n λλ
-==+=+
+∑∑(2.2-17)
式中系数i C 、j C 由初始条件确定。

如果激励信号是在0n =时接入的,差分方程的解适用于0n ≥。

对于N 阶差分方程,用给定的n 初始条件(0),(1),
,(1)y y y N -就可确定全部待定系数i C 和j C 。

如果差分方
程的特征根都是单根,则方程的全解为式(2.2-16),将给定的初始条件(0),(1),...,(
y y y N -分别代入到式(2.2-16),可得
121122111
1122(0)(0)
(1)(1)(1)(1)n p n n p N N N N n p y C C C y y C C C y y N C C C y N λλλλλλ---=++++⎫⎪=++++⎪⎬⎪
⎪-=++++-

(2.2-18)
由以上方程可求得全部待定系数(1,2,,)i C i n =。

例2.2-2 若描述某系统的差分方程为
()4(1)4(2)(n)y n y n y n f +-+-=
(2.2-19)
已知初始条件(0)0y =,(1)1y =-;激励()2,0n
f n n =≥。

求方程的全解。

解:首先求齐次解。

上述差分方程的特征方程为
2440λλ++=
可解得特征根122λλ==-,为二重根,由表2.2-1可知,其齐次解
12()(2)(2)n n h y n C n C =-+-
其次求特解。

由表2.2-2,根据()f n 的形式可知特解
()2,0n p y n P n =≥
将()p y n ,(1)p y n -和(2)p y n -代入到式(2.2-19),得
1224242()2n n n n P P P f n --++==
上式中消去2n
,可解得1
4
P =
,于是得特解 1
()(2),04
n p y n n =≥
差分方程的全解
()()()()()()121222,04
n n
n
h p y n y n y n C n C n =+=-+-+
≥ 将已知的初始条件代入上式,有
21
(0)04
y C =+
= 121
(1)22.214
y C C =--+=-
由上式可求得11C =,21
4
C =-
,最后得方程的全解为 11
()(2)(2)(2)n 044
n n n y n n =-+-+≥自由响应
强迫响应

(2.2-20)
差分方程的齐次解也称为系统的自由响应,特解也称为强迫响应。

本例中由于1λ>,故
其自由响应随着n 的增大而增大。

例 2.2-3 若描述某离散系统的差分方程为
6()5(1)(2)()y n y n y n f n --+-= (2.2-21)
已知初始条件(0)0y =,(1)1y =,激励为有始的周期序列()10cos 02
n f n π
⎛⎫
=≥ ⎪⎝⎭
,n ,求其全解。

解:首先求齐次解。

差分方程的特征方程为
26510λλ-+=
其特征根为112λ=
,21
3
λ=,方程的齐次解 ()121123n
n
h y n C C ⎛⎫⎛⎫
⎪ ⎪⎝⎭⎝⎭
=+
其次求特解。

由表2.2-2可知,特解
()cos sin 2
2
p n n y n P Q π
π⎛⎫⎛⎫
=+
⎪ ⎪⎝⎭⎝⎭
其移位序列
()()()111cos[
]sin[]cos sin 2
2
22p n n n n y n P Q P Q ππππ--⎛⎫⎛⎫-=+=- ⎪ ⎪⎝⎭
⎝⎭
(2)(2)(2)cos[
]sin[]cos sin 2222
p n n n n y n P Q P Q πππ
π
--⎛⎫⎛⎫
-=+=-- ⎪ ⎪⎝⎭⎝⎭
将()p y n ,(n 1)p y -,(n 2)p y -代入到式(2.2-21)并稍加整理,得
n (65)cos (65)sin ()10cos 2
22
n n P Q P Q P Q f n π
ππ
⎛⎫⎛⎫⎛⎫
+-+--==
⎪ ⎪ ⎪⎝⎭⎝⎭⎝⎭
由于上式对于任何0n ≥成立,因为等号两端的正、余弦序列的系数应相等,于是有
6510P Q P +-= 650Q P Q --=
由上式可解得1P Q ==,于是特解
()cos sin 2cos ,n 02224p n n n y n ππ
ππ⎛⎫⎛⎫⎛⎫
⎪ ⎪ ⎪⎝⎭⎝
⎭⎝⎭=+=-≥
(2.2-22)
方程的全解
()()()1211cos sin ,02322n n
p h n n y n y n y n C C n ππ⎛⎫⎛⎫⎛⎫⎛⎫
⎪ ⎪ ⎪ ⎪⎝⎭⎝⎭⎝⎭⎝⎭
=+=+++≥
将已知的初始条件代入上式,有
()12010y C C =++= ()1211
11123
y C C =++=
由上式可解得12C =,23C =-,最后解得全解
()()
()
1123cos sin 232211232cos 02324n n
n
n
n n y n n n ππππ⎛⎫⎛⎫⎛⎫⎛⎫
⎪ ⎪
⎪ ⎪⎝⎭⎝⎭
⎝⎭⎝⎭
⎛⎫⎛⎫⎛⎫
⎪ ⎪ ⎪⎝⎭⎝⎭⎝⎭
=-++=-+-≥强迫响应自由响应
稳态响应
瞬态响应
, (2.2-23)
由上式可见,由于本例中特征根1,21λ<,因而其自由响应是衰减的。

一般而言,如果差分方程所有特征根均满足
()1,2,,1i i n λ=<那么自由响应将随着n 的增大而逐渐衰减趋近
于零。

这样系统称为稳定的系统,这时自由响应也称为瞬态响应。

稳定系统在阶跃序列或有始周期序列作用下,其强迫响应也称为稳态响应。

2.2.3零输入响应和零状态响应
LTI 系统的全响应()n y 也可分为零输入响应和零状态响应。

零输入响应是激励为零
时仅由初始状态所引起的响应,用()x y n 表示;零状态响应是系统的初始状态为零时,仅
由输入信号
()f n 所引起的响应,用()f y n 表示。

这样LTI 系统的全响应将是零输入响
应与零状态响应之和,即
()()()
x f y n y n y n =+
(2.2-24)
在零输入条件下,式(2.2-10)等号右端为零,化为齐次方程。

若其特征根均为单根,则其零状态响应
()1
N
n
x xi i i y k C λ==∑
(2.2-25)
式中x i C 为待定常数。

若系统的初始状态为零,这时方程(2.2-10)仍是非齐次方程,若其特征根均为单根,则其零状态响应
()()
1
N
n p i f fi i y n C y n λ==+∑
(2.2-26)
式中f i C 为待定常数。

系统的全响应可分为自由响应和强迫响应,也可分为零输入响应和零状态响应,它们的关系是
()()()1
1
1
n N N
n n n i i
p x i i
i p f i i i i y n C y n C C y n λλλ====+=++∑∑∑自由响应
零输入响应
零状态响应
强迫响应
(2.2-27)
式中
1
1
1
=
N
N
N
n n n i i xi i i f i i i i C C C λλλ===+
∑∑∑ (2.2-28)
可见,两种分解方式有明显的区别。

虽然自由响应与零输入响应都是齐次解的形式,
但它们的系数并不相同,xi
C 仅由系统的初始状态所决定,而i C 是由初始状态和激励共同
决定。

在用经典法分别求解系统的零输入响应和零状态响应时,需要已知各响应的初始值,用以分别确定常数x i C 和f i C 。

由式(2.2-25)可知,在时刻n j =
()()()0,1,2,
,1x f y j y j y j j N =+=-
(2.2-29)
这时各初始值()y j 中不仅包含有零输入响应的初始值()x y j ,也包含零状态响应的初始

()f y j ,这常常不便区分。

为此,如果激励是在0n =时接入的,通常以
()()()1,2,
,y y y N ---描述系统的初始状态。

因为在0n <时,激励尚未接入,显然
零状态响应在这些时刻的值为零,即
()()()120f f f y y y N -=-==-=
(2.2-30)
由式(3.1-25)可知,这时
()()()()()()11,22,,x x x y y y y y N y N -=--=--=-(2.2-31)
它们给出了该系统已往历史的全部信息,称其为系统的初始状态。

这样,在求解时可分别根据()()1,
x
x y y N -和()()1,,f f y y N --,利用该系统所满足的差分方程,用迭代
法分别求得零输入响应和零状态响应的初始值()x y j 和()()0,1,2,
,1x y j j N =-,并
进一步确定系数x i C 和f i C 。

例2.2-4 若描述某离散系统的差分方程为
()()()()3122n y n y n y n f +-+-=
(2.2-32)
已知激励
()2,0n f n n =≥,初始状态()()1
10,22
y y -=-=
,求系统的零输入响应、零状态响应和全响应。

解:
(1) 零输入响应
()()()31220x x x y n y n y n +-+-=
(2.2-33)
由式(2.2-32),其初始状态()()()()1110,222
x
x y y y y -=-=-=-=。

首先求出初始值()()0,1x
x y y ,式(2.2-33)可写为
()()()3122x x x y n y n y n =----
令0,1n =,并将()()12x
x y y --,
代入,得 ()()()031221x x x y y n y n =----=- ()()()130213x x x y y y =---=
式(2.2-34)的特征值为12
1,2λλ=-=-,由表2.2-1,其齐次解
()()()1212n n
x x x y n C C =-+-
(2.2-34)
将初始值代入得
()1201x x x y C C +=-= ()12123x x x y C C =--=
可解得1
21,2x
x C C ==-,于是得该系统的零输入响应
()()()122,0n
n
x y n n =---≥(2.2-35)
实际上,式(2.2-34)满足齐次方程式(2.2-33),而初始值()0x y ,()1x y 也是有该方
程递推出的,因而直接用()()12x
x y y --,
确定待定常数12x x C C ,将更加简便。

即在式(2.2-34)中令1,2n =--,有
()()12121
10
211
242
x x x x x x y C C y C C -=--=-=+=
可解得121,2x x C C ==-,与前述结果相同。

(2) 零状态响应
()()()()3122f f f y n y n y n f n +-+-=(2.2-36)
和初始状态()()120f f y y -=-=。

首先求出初始值()()0,1f f y y ,将式(2.2-37)写为
()()()()3122f f f y n y n y n f n =----+
令0,1n =,并代入()()120f f y y -=-=和()()0,1f f ,得
()()()()()()()()03122011302111f f f f f f y y y f y y y f ⎫=----+=⎪
⎬=---+=-⎪⎭
(2.2-37)
系统的零状态响应是非齐次差分方程式(2.2-36)的全解,分别求出方程的齐次解和特解,得
()()()()()()()n
12121121223
n n n n
f f f p f f y n C C y n C C =-+-+=-+-+
将式(2.2-37)的初始值代入上式,有
()()1212101;3
2
121
3f f f f f f y C C y C C =++
==--+=-
可解得121,13
f f C C =-=,于是得零状态响应
()()()()11122,0
33n n
n f y n n =-
-+-+≥
(2.2-38) 系统的全响应是零输入响应与零状态响应之和,即
()()()()()()()()1112212233
n
n
n
n
n
x f n y n y n y =-----+-+=+自由响应强迫响应
自由响应
零输入响应
零状态响应
()()()211220
33n n
n n =
---+≥,
(2.2-39)
2.2.4单位序列和单位序列响应
单位序列响应
当LTI 离散系统的激励为单位序列()n δ
时,系统的零状态响应称为单位序列响应(或
单位样值响应、单位取样响应、单位函数响应),用()h n 表示,它的作用与连续系统的
冲激响应()h
t 相类似。

例2.2-5 ()()()()122y n y n y n f n ----=,求离散系统的单位序列响应()n h 。

解:
()()()()122n y n y n y n f ----=
(2.2-40)
根据单位序列响应()h
n 的定义,它应满足方程
()()()()
122h n h n h n n δ----=
(2.2-41)
且初始状态()()120h
h -=-=。

将上式移项有
()()()()122h n h n h n n δ=-+-+
令0,1n =,并考虑到()()01,10δδ==,可求得单位序列响应()h n 的初始值
()()()()()()()()012201102111h h h h h h δδ=-+-+=⎫⎪
⎬=+-+=⎪⎭
(2.2-42)
(2)求()h
n
对于0n >,由式(2.2-41)知()h
n 满足齐次方程
()()()1220h n h n h n ----=
其特征方程为
()()22120λλλλ--=+-=
其特征值121,2λλ=-=,得方程的齐次解
()()()1212,0n n
h n C C n =-+>
将初始值(式2.2-42)代入,有
()1201h C C =+= ()12121h C C =-+=
请注意,这时已将()0,(1)h h 代入,因为方程的解也满足0n =。

由上式可解得
1212
,33C C ==。

于是得系统的单位序列响应
()()()1212,033
n
n h n n =
-+≥ 或写为
()()()()n 1
212n 33n h n ε⎡⎤=-+⎢⎥⎣⎦
阶跃响应
当LTI 离散系统的激励为单位阶跃序列()n ε时,系统的零状态响应称为单位阶跃响应
或阶跃响应,用()g n 表示。

若已知系统的差分方程,那么利用经典法可以求得系统的单
位阶跃响应()g
n 。

例2.2-6求()()()()122y n y n y n f n ----=系统的单位阶跃响应。

解:(1)经典法
系统的差分方程为
()()()()122y n y n y n f n ----=
根据阶跃响应的定义,()g
n 满足方程
()()()()
122g n g n g n n ε----=
(2.2-43)
和初始状态()()120g g -=-=。

上式可写为
()()()()122g n g n g n n ε=-+-+
将0,1n =和()()011ε
ε==代入上式,得初始值
()()()()012201g g g ε=-+-+=
()()()()102112g g g ε=+-+=
式(2.2-4)的特征根121,2λλ=-=,容易求得它的特解()1
,02
p g n n =-
≥。

于是得 ()()()12112,
02
n n
g n C C n =-+-≥
将初始值代入上式,可求得1214
,63
C C =
=,最后得该系统的阶跃响应 ()()()()1
41126
32n n g n n ε⎡⎤=-+-⎢⎥⎣⎦
2.2.5卷积和 2.2.5.1卷积和
在LTI 连续时间系统中,把激励信号分解为一系列冲激函数,求出各冲激函数单独作用于系统时的冲激响应,然后将这些响应相加就得到系统对于该激励信号的零状态响应。

这个相加的过程表现为求卷积积分。

在LTI 离散系统中,可用与上述大致相同的方法进行分析。

由于离散信号本身是一个序列,因此,激励信号分解为单位序列的工作很容易完成。

如果系统的单位序列响应为已知,那么,也不难求得每个单位序列单独作用于系统的响应。

把这些序列相加就得到系统对于该激励信号的零状态响应,这个相加的过程表现为求卷积和。

任意离散时间序列()(
),2,1,0,1,2,
f n --可表示为
()()()()()()()22110f n f n f n f n δδδ=+-++-++
()()()()11f n f i n i δδ+-+
+-+
()()i f i n i δ∞
=-∞
=
-∑
(2.2-44)
如果LTI 系统的单位序列响应为()h
n ,那么,由线性系统的齐次性和时不变系统的
移位不变性可知,系统对
()()f i h n i -。

根据系统的零状态线性性质,式(2.2-44)的序
列()f n 作用于系统所引起的零状态响应()n f y 应为
()()()()()()()22110f y n f h n f h n f h n =+-++-++
()()()()11f h n f i h n i +-+
+-+
()()i f i h n i ∞
=-∞
=
-∑
(2.2-45)
式(2.1-45)称为序列()f n 和()h n 的卷积和。

式(2.2-45)表明,LTI 系统对于任何
激励
()f n 的零状态响应是激励与系统单位序列响应()h n 的卷积和。

一般而言,若有两个序列()1
f n 和()2f n ,和式
()()()12i f n f i f n i ∞
=-∞
=
-∑
(2.2-46)
称为()1
f n 与()2f n 的卷积和,也简称为卷积。

卷积常用符号“*”表示,即
()()()
()()1212def
i f n f n f n f i f n i ∞
=-∞
=*-∑ (2.2-47)
如果序列()1
f n 是因果序列,即有0n <,()10f n =,则式(2.2-47)中求和下限
可改写为零,即若()1
0f n =,0n <,则
()()()()
12120
i f n f n f i f n i ∞
=*=-∑ (2.2-48)
如果()1
n f 不受限制,而()2f n 为因果序列,那么式(2.2-47)中,当0n i -<,即
i n >时,()0f n i =-,因而求和的上限可改写为n ,即:
若()2
0f n =,0n <,则
()()()()1212n
i f n f n f i f n i =-∞
*=
-∑ (2.2-49)
如果()1f n ,()2f n 均为因果序列,即:
若()()1
20,0f n f n n ==<,则
()()()()12120
n
i f n f n f i f n i =*=-∑
(2.2-50)
例2.2-7 如()()112n
f n n ε⎛⎫

⎝⎭
=,()2n 1f =,n -∞<<∞,()()3f n n ε=,求:
(1)()()12f n f n *;(2)()()13f n f n *
解:
(1) 由卷积和的定义式(2.2-46),考虑到()2
1f k i -=得
()(
)()12112i
i f k f k i ε∞
=-∞⎛⎫

⎝⎭
*=⨯∑
()()()()1312i
i f n f n i n i εε∞
=-∞⎛⎫
⎪⎝⎭
*=-∑
上式中,当0i <时()0i ε
=,故求和下限可改写为
0;当0k i -<,即i k >时
()0n i ε-=,因而从1n +到∞的和为零,故求和上限可改写为k ;而在0i n ≤≤区间
()()1i n i εε=-=,于是上式可写为
()()1
1130112112112212
n i n n
i f n f n ++=⎛⎫ ⎪⎡⎤
⎛⎫
⎛⎫⎝⎭
⎢⎥ ⎪ ⎪⎢⎥⎝⎭⎝⎭⎣⎦
-*===--

显然,上式中0n ≥,故应写为
()()()()()113112122n
n f n f n n n n εεε+⎡⎤
⎛⎫
⎛⎫⎢⎥ ⎪
⎪⎢⎥⎝⎭
⎝⎭⎣⎦
*=*=-
2.2.5.2卷积和的性质
离散信号卷积和的运算也服从某些代数运算规则,
()()12f n f n *=
()()21f n f n *
(2.2-51)
即离散序列的卷积和也服从交换律,类似地,二序列的卷积和也服从分配率和结合律,即有
()()()()()()()1231213f n f n f n f n f n f n f n ⎡⎤⎢⎥


*+=*+* (2.2-52)
()()()()()()123123f n f n f n f n f n f n ⎡⎤⎡⎤⎢⎥⎢⎥




**=** (2.2-53)
如果两序列之一是单位序列,由于()n δ仅当0n =时等于1,0n ≠时全为零,因而

()()()()()()()i f n n n f n i f n i f n δδδ∞
=-∞
*=*=
-=∑ (2.2-54) 即序列
()f n 与单位序列()n δ的卷积和就是序列()f n 本身。

将式(2.2-54)推广,
()f n 与移位序列()1n n δ-的卷积和
()()()()11()i f n i f n n n f i n i n δδ∞
=-∞
=-*-=
--∑
考虑到交换律,有
()()()()11f n n n n n f n δδ=*--*
(2.2-55)
此外还有
()()()()()122112f n n n n f n n n n f n n n δδ==--*--*--
(2.2-56)

()()()12f n f n f n =*

()()()()()1212122112f n n f n n f n n f n n f n n n ==--*--*--(2.2-57)
以上各式中1n 、2n 均为整常数,各式的证明和图示与连续系统相似,不多赘述。

例 2.2-8复合系统由两个子系统级联组成,已知子系统的单位序列响应分别为
()()()()12,n n h n a n h n b n εε==,(,a b 为常数)求复合系统的单位序列响应()h n 。

解: 根据单位序列响应的定义,复合系统的单位序列响应()h n 是激励()()f n n δ=时系
统的零状态响应,即()()f
y n h n =。


()()f n n δ=,则子系统1的零状态响应
()()()()()()111f x n h n h h n f n n n δ===**
当子系统2的输入为()f
x k 时,子系统2的零状态响应亦即复合系统的零状态响应
()()()()()()
122f f x n h n h n h n y n h n ===**
即复合系统的单位序列响应
()()()()()12i i n i h n h n a i b n i h n εε∞
=-∞
-==∙-*∑
考虑到当0i <时()0i ε
=,i n >时()0n i ε-=以及在0i n ≤≤区间
()()i n i εε=-=1,以上各式可写为:
当a b ≠时
()()()100i
n
n
n n i n n i i a h n a n b n a b b b εε-==⎛⎫
⎪⎝⎭
=*==∑

1
11
11n n n n a b
b a b a b a
b
+++⎛⎫

⎝⎭
--==
--
当a b =时
()()0
11n
n
n
i h n b
n b ===+∑ 显然上二式在0n ≥成立,故得
()()()()()()11
,1,n n n n n b a n a b b a h n a n b n n b n a b εεεε++⎧⎪
⎨⎪⎩
-≠-=*=+=
上式中,若1,1a b ≠=,则有
()()()111n n
a a n n n a εεε+-*=- 若1a
b ==,有
()()()()1n n n n εεε*=+
最后举例说明时域分析求解LTI 离散系统全响应的有关问题。

例 2.2-9 如图 2.2-9所示离散系统,已知初始状态
()10y -=,()126y -=,激励()()()()()cos 1n f n n n n πεε==-,求系统的全响应()y n 。

图2.2-9离散系统
解: 按图2.2-9,不难写出描述系统的差分方程为
()()()()122y n y n y n f n ----=
求零输入响应
根据零输入响应的定义,它满足方程
()()()1220x x x y n y n y n ----=
和初始状态()()110x y y -=-=,()()1
226x y y -=-=。

可推得其初始条件
()()()1
01223
x x x y y y =-+-= ()()()110213
x x x y y y =+-= 差分方程的特征根为121,2λλ=-=-,故有
()()()
1212n n x x x y n C C =-+
将初始条件代入,有 ()12103
x x x y C C =+= ()121123x x x y C C =-+=
解得1212,99
x x C C ==,得零输入响应 ()()()1212,099
n n x y n n =-+≥ 求单位序列响应和零状态响应
根据单位序列响应的定义,系统的单位序列响应()h n 满足方程
()()()()122h n h n h n n δ----=
和初始状态()()120h h -=-=。

()()()()121233n n h n n ε⎡⎤⎢⎥⎣⎦
=-+ 系统的零状态响应等于激励()f n 与单位序列响应()h n 的卷积和,即
()()()()()()()()1
212133n n n f y n h n f n n n εε⎡⎤⎢⎥⎣⎦=*=-+
*- ()()()()()()()()12112133
n n n n n n n n εεεε=-*-+*- (2.2-58) 由于
()()()()()()()1111n
n n n n n n εεε-*-=+- ()()()()()
()()11111212121123n n n n n n n n εε++++⎡⎤⎢⎥⎣⎦
--*-==---- ()()()212133n n n ε⎡⎤+-⎢⎥⎣⎦
= 将它们代入到式(2.2-58),得
()()()()()()()()()()()122111213333154112399n n n f n n n y n n n n n n εεε⎡⎤⎢⎥⎣⎦
⎡⎤⎢⎥⎣⎦=+-++-=-+-+
最后得到系统的全响应
()()()()()()()()()()()()1215412112993991221233n n n n n x f n n y n y n y n n n n ε⎡⎤⎢⎥⎣⎦
=+=-++-+-+=+-+
习题二
2.1 已知描述系统的微分方程和初始状态如下,求系统的零输入响应。

(1)
(2) ,
2.2已知描述系统的微分方程和初始状态如下,求系统的全响应。

(1) ,
(2) ,
2.3已知描述系统的微分方程和初始状态如下,求系统的零输入响应,零状态响应与系统的
全响应。

(1) ,
(2) ,
2.4已知描述系统的微分方程如下,求系统的冲激响应和h(t)与阶跃响应g(t)。

(1)
(2)
(3)
(4)
2.5计算下列函数卷积积分
(1)
(2)
(3)
2.6 已知
(1)
(2)
求, 。

2.7 已知系统: , 求
2.8 某LTI系统,当激励时,零状态响应为:
求冲激响应。

2.9 计算
2.10 求系统的零状态响应
2.11 已知描述某离散时间LTI系统的差分方程
, 求系统在作用下的固有响应,强迫响应,完全响应。

2.12求系统的零输入响应
(1)a=1/2; (2)a=-1; (3)a=5/4
2.13求系统的零状态响应
(1)
(2)
(3)
2.14 求下列方程描述的离散时间系统的单位脉冲响应
(1)
(2)
(3)
2.15求系统的差分方程
(1)
(2)
(3)
2.16计算卷积
(1){1,2,1}{1,2,1}
(2){1,0,2,0,1}{1,0,2,0,1}
(3){1,2,1} , , , ,
2.17计算卷积
(1)
(2)
(3)
2.18已知离散时间系统的差分方程为求全响应。

相关文档
最新文档