郑君里《信号与系统》考研考点讲义.pdf
- 文件大小: 3.63MB
- 文件类型: pdf
- 上传日期: 2025-08-25
- 下载次数: 0
概要信息:
目 录
第一章 绪论 (1)
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
第二章 连续时间系统的时域分析 (10)
!!!!!!!!!!!!!!!!!!!!!!!!!
第三章 傅里叶变换 (26)
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
第四章 拉普拉斯变换、连续时间系统的S域分析 (49)
!!!!!!!!!!!!!!!!!!!
第五章 傅里叶变换应用于通信系统-滤波、调制与抽样 (78)
!!!!!!!!!!!!!!!!
第六章 信号的空间矢量分析 (95)
!!!!!!!!!!!!!!!!!!!!!!!!!!!
第七章 离散信号与系统时域分析 (116)
!!!!!!!!!!!!!!!!!!!!!!!!!
第八章 Z变换和离散时间系统的Z域分析 (131)
!!!!!!!!!!!!!!!!!!!!!
第九章 离散傅里叶变换及其他离散正交变换 (151)
!!!!!!!!!!!!!!!!!!!!
第十一章 反馈系统 (165)
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
第十二章 系统的状态变量分析 (179)
!!!!!!!!!!!!!!!!!!!!!!!!!!
第一章 绪 论
一、信号的概念和分类
1.信号
随时间变化的物理量,消息的具体表现形式。
2.描述和分类
(1)描述:波形图和函数表达式
(2)分类:确定和随机
周期和非周期
连续和离散(注意区别模拟和数字)
一维和多维
要点:判定离散时间信号是否为数字信号
例1 判断下列图示信号的类型
(a)连续(模拟);(b)连续(量化);(c)离散(数字);(d)离散;(e)离散(数字);(f)离散(数字)
二、典型非奇异信号
1.指数信号
f(t)=Keat,τ= 1
a
—1—
郑君里《信号与系统》考点精讲
2.正弦信号
f(t)=Ksin(ωt+θ)
3.复指数信号
f(t)=Kest=Keωtcos(ωt)+jKeσtsin(ωt)
4.抽样信号
Sa(t)=sintt
5.钟形信号
f(t)=Ee-(
t
τ)
2
要点:熟记典型信号的表达式、波形特征
三、信号的运算
关于自变量:位移、反折、尺度
关于信号的幅值:加、乘(在相同时刻)、微分、积分
位移:时间轴上左右平移
f(t) →
移位
f(t)=f(t+t0)
反折:时间轴上关于纵轴镜像
f(t) →
反折
f(t)=f(-t)
尺度:时间轴上的压缩和扩展
f(t) →
尺度
f(t)=f(at)
相加
f1(t)、f2(t) →
相加
f(t)=f1(t)+f2(t)
相乘
f1(t)、f2(t) →
相乘
f(t)=f1(t)×f2(t)
微分
f(t) →
微分 d
dtf(t)
积分
f(t) →
积分∫
t
-
!
f(t)dt
掌握用上述运算进行信号波形变换的方法
(图见视频)
例2 已知f(t),求f(-4-2t)
—2—
四、典型奇异信号及其特性
斜变信号、单位阶跃、单位冲击、冲击偶信号
要点:四类奇异信号的函数定义和波形
重点:掌握单位冲击信号的筛选性,用于一般信号的单位冲击分解
1.斜变信号
(1)单位斜变
f(t)=
0 t<0
t t{ 0
(2)截平的斜变
f(t)=
k
τ
f(t) t<τ
k t
{
τ
(3)三角形脉冲
f(t)=
k
τ
f(t) t
"τ
0 t>
{
τ
—3—
郑君里《信号与系统》考点精讲
2.单位阶跃信号
(1)单位阶跃
u(t)=
0 t<0
1 t>{ 0
(2)矩形脉冲
RT(t)=u(t)-u(t-T)
(3)符号函数
sgn(t)=
1 t>0
-1 t<{ 0
用来描述符号函数:sgn(t)=2u(t)-1
描述信号的接入特性:f(t)=f1(t)u(t)
3.单位冲激信号(6种来源)
(1)矩形脉冲
δ(t)=lim
τ→0
1
τ
[u(t+τ2)-u(t-
τ
2)]
(2)三角脉冲
δ(t)=lim
τ→0
1
τ
1
τ
(1- t
τ
)[u(t+τ)-u(t-τ2{ })]
(3)双边指数脉冲
δ(t)=lim
τ→0
1
2τ
e-
t{ }τ
(4)钟形脉冲
δ(t)=lim
τ→0
1
τ
e-π(
t
τ){ }2
—4—
(5)采样信号
δ(t)=lim
k→!
k
π
Sa(kt{ })
(6)狄拉克函数
∫
!
-
!
δ(t)dt=1
δ(t)=0(当t≠0
{
)
4.冲击偶信号及其演变
δ′(t)=ddtδ(t)=
d2
dt2
u(t)
(图见视频)
5.单位冲激信号的性质
(1)抽样性
∫
!
-
!
δ(t)f(t)dt=∫
!
-
!
δ(t)f(0)dt=f(0)∫
!
-
!
δ(t)dt=f(0)
∫
!
-
!
δ(t-t0)f(t)dt=∫
!
-
!
δ(t-t0)f(t0)dt=f(t0)
(2)偶对称性
δ(t)=δ(-t)
(3)尺度变换
δ(at)= 1
aδ(t)
6.冲激偶信号的性质
(1)抽样性
∫
!
-
!
δ′(t)f(t)dt=-f′(0)
(2)乘法
f(t)δ′(t)=f(0)δ′(t)-f′(0)δ(t)
(3)积分
∫
!
-
!
δ′(t)dt=0
(4)尺度
δ′(at)= 1
a·
1
aδ′(t)
(5)卷积
—5—
郑君里《信号与系统》考点精讲
f(t)δ′(t)=ddtf(t)=f′(t)
例3 求列下表达式的值
(1)∫
!
-
!
f(t0-t)δ(t)dt
(2)∫
!
-
!
e-jωt[δ(t)-δ(t-t0)]dt
解:根据单位冲激信号的筛选性可得:
∫
!
-
!
f(t0-t)δ(t)dt=-∫
!
-
!
f(t0-t)δ(-t)d(-t)
=-∫
-
!
!
f(t0+τ)?δ(τ)dτ=∫
!
-
!
f(t0+τ)δ(τ)dτ
=∫
!
-
!
f(t)?δ(t-t0)dt=f(t0)
∫
!
-
!
e-jωt[δ(t)-δ(t-t0)]dt=∫
!
-
!
e-jωtδ(t)dt-∫
!
-
!
e-jωtδ(t-t0)dt=1-e
jωt0
五、信号的分解
直流分量与交流分量、奇分量与偶分量、脉冲分量、阶跃分量、正交分量
重点:掌握信号的单位冲击分解
1.交流分量与直流分量
f(t)→fD +fA(t)
P=1T∫
T
2
-T2
f2(t)dt
=1T∫
T
2
-T2
[fD +fA(t)]
2dt
=1T∫
T
2
-T2
[f2
D
+2fDfA(t)+f
2
A(t)]dt
=f2D +
1
T∫
T
2
-T2
f2A(t)dt
信号的功率 =直流功率 +交流功率
直流功率 =f2D
交流功率 =1T∫
T
2
-T2
f2A(t)dt
2.奇分量和偶分量
f(t)→fe(t)+fo(t)=
1
2[f(t)+f(-t)]+
1
2[f(t)-f(-t)]a
3.正交分量
用相互正交的正交函数集各个分量的线性组合来表示某些信号
—6—
注意:只用于周期信号的分解,常见的正交完备集就是三角函数和复指数函数,即周期信号的三
角级数分解和复指数分解。
4.脉冲分量
→ f(t)=lim
τ→0
f1(t)=∫
!
-
!
f(τ)δ(t-τ)dτ
六、系统模型及其分类
重点:掌握系统的分类
1.系统模型的表示方法:数学方程和框图
连续时间系统的描述:微分方程
离散时间系统的描述: }
→
差分方程
选取合适的状态变量,如电容电压、电感电流等
框图表示系统的基本运算单元:
加法
乘法
积分
例4 求解下列RLC串联电路并画出系统框图
解:根据KVL可得:
i(t)=C
duc(t)
dt
uc(t)+L
di(t)
dt +Ri(t)=e(t
{
)
整理可得该系统的微分方程为:LCd
2i(t)
dt2
+RCdi(t)dt +i(t)=C
de(t)
dt
对上式积分:i(t)=1L∫e(τ)dτ-RL∫i(τ)dτ-1LCi(τ)dτdτ
—7—
郑君里《信号与系统》考点精讲
利用框图的基本运算单元得:
2.系统的分类
1)连续时间系统、离散时间系统———输入输出信号类型
2)即时系统、动态系统———有无储能元件
3)集总参数、分布参数———状态变量手否考虑空间位置因素
4)线性系统、非线性系统———是否满足线性特性
5)时不变系统、时变系统———是否具有时不变特性
6)因果系统、非因果系统———响应与激励施加的先后关系
7)稳定系统、非稳定系统———系统响应的收敛性
8)可逆系统、不可逆系统———响应与激励是否具有
一一对应的关系
重点研究:稳定的线性时不变系统因果系统
N阶连续LTI系统
描述
→
N阶常系数线性微分方程
N阶离散LTI系统
描述
→
N阶常系数线性差分方程
七、系统模型及其分类
重点:线性时不变系统的判定
1.线性时不变系统的特点:
1)线性性质:叠加性、均匀性(齐次性)
2)时不变性:响应与激励施加的时刻无关
3)微分特性:本质即线性
2.线性性质的判定:
分解性、零状态线性、零输入线性
3.时不变特性的判定:
激励信号平移,判定响应
例5 试判定下述系统的线性特性、时变特性r(t)=x(0-)sint+te(t)
解:(1)满足分解性,易知:
rzi(t)=x(0-)sint rzs(t)=te(t)
(2)x(0-)=x1(0-)时rzi1(t)=x1(0-)sint,
x(0-)=x2(0-)时,rzi2(t)=x2(0-)sint
—8—
当x(0-)=ax1(0-)+bx2(0-)
rzi(t)=[ax1(0-)+bx2(0-)]sint
=ax1(0-)sint+bx2(0-)sint
=arzi1(t)+brzi2(t)
(3)e(t)=e1(t)时,rzs1(t)=te1(t)
e(t)=e2(t)时,rzs2(t)=te2(t)
当e(t)=ae1(t)+be2(t)时,
rzs(t)=t[ae1(t)+be2(t)]=ate1(t)+bte2(t)=arzs1(t)+brzs2(t)
满足零输入线性,综合(1)~(3),该系统为线性系统
关于时变特性的判断:当e(t)→e(t-t0)时
系统响应y(t)=x(0-)sint+te(t-t0)
而r(t-t0)=x(0-)sin(t-t0)+(t-t0)e(t-t0)≠y(t)
因此,该系统为时变系统
八、系统分析的方法
1.系统研究的内容
系统分析:给定系统在特定激励下的响应
系统综合:满足响应与激励之间关系的系统设计及实现
2.LTI系统分析的方法
1)数学描述:输入—输出描述法,适用于单变量系统
状态变量描述法,适用于多变量系统
2)模型求解:时域求解法(经典法、卷积法、数值解法)
变换域求解法(LT、FT、ZT)
→ →时域卷积 频域乘积 频域代数方程
→ →频域卷积 时域乘积 时域代数方程
—9—
郑君里《信号与系统》考点精讲
第二章 连续时间系统的时域分析
LTI系统分析的方法:
数学描述:输入—输出描述法,适用于单变量系统
状态变量描述法,适用于多变量系统
模型求解:时域求解法(经典法、卷积法、数值解法)
变换域求解法(LT、FT、ZT)
→ →时域卷积 频域乘积 频域代数方程
经典法:解微分方程,齐次通解+非齐次特解
自由响应+强迫响应
注意初始值0+、0-的区别与求法
卷积法:求卷积运算,零状态响应,物理概念明确。
一、微分方程的建立与求解
1.微分方程的建立
根据系统的物理模型、元件的约束特性、系统结构的
约束特性:KVL、KCL、欧姆定理、达朗贝尔原理。
对于机械系统的处理:机电类比法。
例1 列出下述两个系统的微分方程
对于第一个系统,选取v(t)为变量,可得元件的电流关系:
iR(t)=
1
Rv(t) iL(t)=
1
R∫
t
-
!
v(τ)dτ iC(t)=C
d
dtv(t)
根据KCL:iR(t)+iL(t)+iC(t)=iS(t)
整理得:Cd
2
dt2
v(t)+1R
d
dtv(t)+
1
Lv(t)=
d
dtiS(t)
对于第二个系统,选取v(t)为变量,根据胡克定律:
Fk(t)=k·x(t)=k∫
t
-
!
v(τ)dτ
摩擦力:Ff(t)=f·v(t) f:摩擦系数
—01—
惯性力:Fm(t)=m·
d
dtv(t) 惯性定律
力学平衡的达郎贝尔原理:Fk(t)+Ff(t)+Fm(t)=FS(t)
整理得:md
2
dt2
v(t)+fddtv(t)+kv(t)=
d
dtFS(t)
分析形式: →机械量 电气量(机电类比法)
2.微分方程的求解
线性时不变系统(LTI)以常系数线性微分方程描述:
设系统激励为e(t),响应为r(t),则系统可描述为:
C0
dn
dtn
r(t)+C1
dn-1
dtn-1
r(t)+… +Cn-1
d
dtr(t)+Cnr(t)
=E0
dm
dtm
e(t)+E1
dm-1
dtm-1
e(t)+… +Em-1
d
dte(t)+Eme(t)
根据经典法:全响应 =齐次通解 +非其次特解
C0
dn
dtn
r(t)+C1
dn-1
dtn-1
r(t)+… +Cn-1
d
dtr(t)+Cnr(t)=0
齐次解形如Aeαt的线性组合,设r(t)=Aeαt则:
C0Aα
neαt+C1Aα
n-1eαt+… +Cn-1Aαeα
t+CnAeα
t=0
得特征方程:C0α
n+C1α
n-1+… +Cn-1α+Cn =0
若αk满足上式,则Aeαk
t满足原方程
1)无重根rh(t)=∑
n
i=1
Aieαi
t
2)α1有k重根rh(t)=(∑
k
i=1
Ait
k-i
)·eα1t+∑
n-k+1
i=2
Aieαi
t
特解rp(t)形式取决于激励信号
选定形式,待定系数法
齐次解只与系统本身的结构和参数有关,而与激励信号无关,反映系统本身特征,又称为固有相
应或自由响应。而特解的形式由激励决定,称为强迫响应。
激励e(t) 响应特解rp(t)
常数E 常数B
tp B1t
p+B2t
p-1+… +Bpt+Bp+1
eαt Beαt
cosωt
sinωt
B1cosωt+B2sinωt
tp·eαt·cosωt
tp·eαt·sinωt
(B1t
p+B2t
p-1+… +Bpt+Bp+1)eα
tcosωt
+(D1t
p+D2t
p-1+… +Dpt+Dp+1)eα
tsinωt
—11—
郑君里《信号与系统》考点精讲
这样,原方程的完全解为:
r(t)=rh(t)+rp(t)=∑
n
i=1
Aie
-αit+rp(t)
为使方程有唯一解,需根据初始状态求解待定系数:
设求解区为:0+"t<+!
初始状态为:r(0+)r
、(0+)r
(n-1)(0+)
则:
r(0+)=A1+A2+…An+rp(0+)
r、(0+)=A1α1+A2α2+…Anαn+r
、
p(0+)
……
r(n-1)(0+)=A1α
n-1
1 +A2α
n-1
2 +…An+r
n-1
p (0+
)
以矩阵形式表示为:
[
r(0+)-rp(0+)
r、(0+)-r
、
p(0+)
r(n-1)(0+)-r
n-1
p (0+)
]=[
1 1 … 1
α1 α2 … αn
αn-11 αn-12 … αn-12
][
A1
A2
An
]
例2 给定微分方程:d
2
dt2
r(t)+5ddtr(t)+6r(t)=e(t)
求:(1)e(t)=2e-t,t0,r(0)=2,r′(0)=-1
(2)e(t)=e-2t,t0,r(0)=1,r′(0)=0
时系统的全响应。
解:原系统的特征方程为:λ2+5λ+6=0
求得特征根:λ1 =-2 λ2 =-3
齐次通解为:rh(t)=C1e
-2t+C2e
-3t
根据激励形式可知:
(1)e(t)=2e-t时,特解为:rp(t)=Ae
-t
将特解带入原方程得:
Ae-t+5(-Ae-t)+6Ae-t=2Ae-t
解得:A=1
可得特解为:rp(t)=e
-t
系统全响应:r(t)=rh(t)+rp(t)=C1e
-2t+C2e
-3t+e-t
现根据初始状态求待定系数C1、C2:
r(0)=C1+C2+1=2 r′(0)=-2C1-3C2-1=-1
解得:C1 =3、C2 =-2
—21—
故全响应为:r(t)=3e-2t-2e-3t+e-t
(2)e(t)=e-2t时,其指数次与特征跟之一相同:
其特解形式为:rp(t)=(A1t+A0)e
-2t
带入微分方程:A1e
-2t=e-2t
解得:A1 =1
系统全响应:
r(t)=rh(t)+rp(t)=C1e
-2t+C2e
-3t+(te-2t+A0e
-2t)
=(C1+A0)e
-2t+C2e
-3t+te-2t
代入初始条件得:
r(0)=(C1+A0)+C2 =1
r′(0)=-2(C1+A0)-3C2+1=0
解得:C1+A0 =2 C2 =-1
故全响应为:r(t)=2e-2t-e-3t+te-2t
利用初始条件无法求出C1和A0,因此分解不出自由响应和强迫响应。
二、起始点的跳变———从0-到0+状态的转换
1.关于0-状态和0+状态
状态为零输入时的初始状态,即初始值由系统的初始储能所产生。
0+状态为加入激励信号后的初始状态,即初始值不但有初始储能,还受激励信号的影响。
系统0-状态 } →
系统激励信号
0+状态
0+状态
齐次通解{ } →
非齐次特解
全响应
0-到0+ →跃变 右边包含δ(t)及其各阶导数
2.0+状态的确定———冲激函数匹配法
冲激函数匹配法:利用t=0时刻方程两边的δ(t)及各阶导数应该平衡的原理来求解(0+)的方
法。原因在于微分方程右边还有δ(t)及其各阶导数,造成(0+)和(0-)时刻的值不相等。
在经典法求全响应的积分常数时,用的是0+状态初始值。
在求系统零输入响应时,用的是0-状态初始值。
在求系统零状态响应时,用的是0+状态初始值,这时的零状态是指0-状态为零。
已知系统为:
a0y(t)+a1y′(t)+… +any
(n)(t)=b0+b1δ(t)+b2δ′(t)+… +bmδ
(m-1)(t)Q
若m
"
n则可设:
—31—
郑君里《信号与系统》考点精讲
y(n)(t)=Cmδ
(m-1)(t)+… +C1δ(t)+C0
y(n-1)(t)=Cmδ
(m-2)(t)+… +C2δ(t)+C1
y(n-m)(t)=Cm
y(n-m-1)(t)=… =y(t)=0
若m>n则可设:
y(n)(t)=Cmδ
(m-1)(t)+… +C1δ(t)+C0
y(n-1)(t)=Cmδ
(m-2)(t)+… +C2δ(t)+C1
y(t)=Cmδ
(m-n-1)(t)+… +Cn-1
将y(t)及其各阶导数带入原方程,求出C0…Cm对y(t)及各阶导数求(0-,0+)的积分.
例3 电路如图示,t<0时S处于位置1且已达到稳态,t=0时,S由1转向2。建立i(t)的微分
方程并求解其在t0+时的变化。
解:1)列写电路方程
根据KVL:R1i(t)+vC(t)=e(t)
LddtiL(t)+iL(t)R2 =vC(t)
根据KCL:CddtvC(t)+iL(t)2 =i(t)
先消去vC(t)得:
d
dti(t)=-
1
R1C
i(t)+ 1R1C
iL(t)+
1
R1
d
dte(t)
d
dtiL(t)=-
R1
Li(t)-
R2
LiL(t)+
1
Le(t)
再消去iL(t)得:
d
dtiL(t)=-
R1
Li(t)-
R2
LiL(t)+
1
Le(t)
d2
dt2
i(t)+(1R1C
+
R2
L)
d
dti(t)+(
1
LC+
R2
R1LC
)i(t)
=1R1
d2
dt2
e(t)+
R2
R1L
d
dte(t)+
1
R1LC
e(t)
—41—
带入参数得:
d2
dt2
i(t)+7ddti(t)+10i(t)=
d2
dt2
e(t)+6ddte(t)+4e(t)
2)求系统全响应
齐次解:
特征方程为:λ2+7λ+10=0
特征根:λ1 =-2 λ2 =-5
可得:ih(t)=A1e
-2t+A2e
-5t(t0+)
由题设知:t0+时,e(t)=4V为常数,
故特解形式为:ip(t)=B
带入系统方程解得:B=85
因此:t0+时,
i(t)=ih(t)+ip(t)=A1e
-2t+A2e
-5t+85
3)确定换路后的i(0+)和
d
dti(0+)
换路前:i(0-)=iL(0-)=
2
R1+R2
=45A
d
dti(0-)=0
vC(0-)=iL·R2 =
4
5×
3
2 =
6
5V
换路后,根据换路定理:
i(0+)=
1
R1
[e(0+)-vC(0+)]=
1
R1
[e(0+)-vC(0-)]=
1
1[4-
6
5]=
14
5 A
d
dti(0+)=
1
R1
[
d
dte(0+)-
d
dtvC(0+)]
=1R1
d
dte(0+)-
1
C[i(0+)-iL(0+{ })]
=1R1
d
dte(0+)-
1
C[i(0+)-iL(0-{ })]
=11[0-
1
1(
14
5-
4
5)]=-2 A/s
4)求i(t)在t0+时的全响应:
根据i(t)的全响应表达式可得:
i(0+)=A1+A2+
8
5 =
14
5
d
dti(0+)=-2A1-5A2 =-
{ 2
—51—
郑君里《信号与系统》考点精讲
解得:A1 =
4
3 A2 =-
2
15
系统全响应为:
i(t)=(43e
-2t-215e
-5t+85) A (t0+)
例4 利用冲激函数匹配法求例四中的全响应
解 前面已推得系统的微分方程为:
d2
dt2
i(t)+7ddti(t)+10i(t)=
d2
dt2
e(t)+6ddte(t)+4e(t)
考虑到激励信号在t=0时刻由2V跃变到4V,代入上式即得t=0时的微分方程如下:
d2
dt2
i(t)+7ddti(t)+10i(t)=2δ′(t)+12δ(t)+8Δu(t)
换路前:i(0-)=
4
5A
d
dti(0-)=0A /s
观察微分方程两侧可知:
d2
dt2
i(t)项中包含δ′(t)
因此可设:
(0-<t<0+)
d2
dt2
i(t)=aδ′(t)+bδ(t)+cΔu(t)
d
dti(t)=aδ(t)+bΔu(t)
i(t)=aΔu(t
)
代入前式:
[aδ′(t)+bδ(t)+cΔu(t)]+7[aδ(t)+bΔu(t)]+10aΔu(t)=2δ′(t)+12δ(t)+8Δu(t)
比较系数:
a=2
b+7a=12
c+7b+10a=
{
8
→
i(0+)-i(0-)=a=2
d
dti(0+)-
d
dti(0-)=b=-2
d2
dt2
i(0+)-
d2
dt2
i(0-)=c=
→
2
i(0+)=i(0-)+2=(
4
5+2)=
14
5A
d
dti(0+)=
d
dti(0-)-2=-2
{ A/s
—61—
三、零输入响应和零状态响应
例1 已知图示的RC电路,电容两端有起始电压vC(0-)试求时的电容t>0两端电压vC(t)。
解:可求得系统方程为:
d
dtvC(t)+
1
RCvC(t)=
1
RCe(t →)
[
d
dtvC(t)+
1
RCvC(t)]·e
t
RC = 1RCe(t)·e
t
→RC
d
dt[e
t
RC·vC(t)]=
1
RCe
t
RC·e(t →)
∫
t
0-
d
dt[e
τ
RC·vC(τ)]dτ=∫
t
0-
1
RCe
τ
RC·e(τ)d →τ
e
t
RC·vC(t)-vC(0-)=
1
RC∫
t
0-
e
τ
RC·e(τ)d →τ
vC(t)=e
-tRC·vC(0-)+
1
RC∫
t
0-
e-
1
RC(t-τ)e(τ)dτ
e-
t
RC·vC(0-)零输入响应rzi(t)
1
RC∫
t
0-
e-
1
RC(t-τ)e(τ)dτ零状态响应rzs(t)
零输入响应满足:
C0
dn
dtn
rzi(t)+C1
dn-1
dtn-1
rzi(t)+… +Cn-1
d
dtrzi(t)+Cnrzi(t)=0
齐次解的一部分:rzi(t)=∑
n
k=1
Azikeαk
t
零状态响应满足:
C0
dn
dtn
rzs(t)+C1
dn-1
dtn-1
rzs(t)+… +Cn-1
d
dtrzs(t)+Cnrzs(t)
=E0
dm
dtm
e(t)+E1
dm-1
dtm-1
e(t)+… +Em-1
d
dte(t)+Eme(t)
初始状态:r(k)(0-)=0 (k=0,1,2,…n-1)
其解:rzs(t)=∑
n
k=1
Azskeαk
t+B(t)
系统全响应:r(t)=∑
n
k=1
Akeαk
t+B(t)=∑
n
k=1
Azikeαk
t+∑
n
k=1
Azskeαk
t+B(t)
例2 将图中t<0电路看做起始状态,分别求t>0时的零输入响应和零状态响应。
—71—
郑君里《信号与系统》考点精讲
解:转换初始储能
前已求得:iL(0-)=
4
5A vC(0-)=
6
5V
d2
dt2
i(t)+7ddti(t)+10i(t)=
d2
dt2
e(t)+6ddte(t)+4e(t)
(图见视频)
1)零输入响应
→系统方程
d2
dt2
ii(t)+7
d
dtizi(t)+10ii(t)=0
→初始状态 izi(0+)
d
dtizi(0+)
(图见视频)
初始值等效电路KCL、KVL:
izi(0+)=-
vC(0+)
R1
=-
vC(0-)
R1
=-65A
iC(0+)=izi(0+)-iL(0+)=izi(0+)-iL(0-)=C
d
dt[-R1·izi(0+)]vC(0+)
d
dtizi(0+)=-
1
R1C
[izi(0+)-iL(0+)]=2 A/s
→零输入响应
d2
dt2
ii(t)+7
d
dtizi(t)+10ii(t)=0
izi(0+)=-
6
5A
d
dtizi(0+)=2
A/s
→izi(t)=(-
4
3e
-2t+215e
-5t) A (t>0)
2)零状态响应
(图见视频)
d2
dt2
izs(t)+7
d
dtizs(t)+10izs(t)=
d2
dt2
e(t)+6ddte(t)+4e(t)
已求得其全响应为:
—81—
izs(t)=Czs1e
-2t+Czs2e
-5t+85 (t>0)
将e(t)=4u(t)代入方程右端得自由项:d
2
dt2
e(t)+6ddte(t)+4e(t)=4δ′(t)+24δ(t)+16
根据冲激函数匹配法确定待定系数:
设:
d2
dt2
izs(t)=aδ′(t)+bδ(t)+cΔu(t)
d
dtizs(t)=aδ(t)+bΔu(t)
izs(t)=aΔu(t
)
(0-<t<0+)
代入微分方程,比较系数后整理可得:
izs(0+)-izs(0-)=a=4
d
dtizs(0+)-
d
dtizs(0-)=b=-
{ 4
izs(0+)=4+izs(0-)=4
d
dtizs(0+)=-4+
d
dtizs(0-)=-
{ 4
izs(t)=Czs1e
-2t+Czs2e
-5t+85 (t>0
→
)
Czs1 =
8
3 Czs1 =-
4
15
系统零状态响应:
izs(t)=(
8
3e
-2t-415e
-5t+85)A (t>0)
i(t)=43e
-2t-215e
-5
t
自由响应
暂态响应
+ 85{
强迫响应
稳态响应
=(-43e
-2t+215e
-5t
)
令输入响应
+(83e
-2t-415e
-5t+85
)
零状态响应
四、冲激响应和阶跃响应
阶跃响应:单位阶跃信号作用下的零状态响应
冲激响应:单位冲激信号作用下的零状态响应
(图见视频)
对于LTI系统,单位冲激响应满足:
C0h
(n)(t)+C1h
(n-1)(t)+… +Cn-1h′(t)+Cnh(t)
=E0δ
(m)(t)+E1δ
(m-1)(t)+… +Em-1δ′(t)+Emδ(t)
—91—
郑君里《信号与系统》考点精讲
h(k)(0-)=0
δ(t)及其各阶导数t0+时恒为0,故:
C0h
(n)(t)+C1h
(n-1)(t)+… +Cn-1h′(t)+Cnh(t)=0
n>m时:h(t)=(∑
n
k=1
Ckeαk
t)u(t)
n
"
m时:采用冲激匹配法确定系数。
例3 求例2中电流i(t)的冲激响应h(t)和单位阶跃响应g(t)。
解:1)单位冲激响应
已求得系统方程为
d2
dt2
i(t)+7ddti(t)+10i(t)=
d2
dt2
e(t)+6ddte(t)+4e(t)
故其单位冲激响应满足:
h″(t)+7h′(t)+10h(t)=δ″(t)+6δ′(t)+4δ(t)
其解为:
h(t)=C1e
-2t+C2e
-5t (t0+)
观察方程右侧,可设:
d2
dt2
h(t)=aδ″(t)+bδ′(t)+cδ(t)+dΔu(t)
d
dth(t)=aδ′(t)+bδ(t)+cΔu(t)
h(t)=aδ(t)+bΔu(t
)
(0-<t<0+)
代入方程解得:a=1,b=-1,c=1故:
h(0+)=b+h(0-)=-1
h′(0+)=c+h′(0-)=
{ 1
代入h(t)中解得:C1 =-
4
3 C2 =
1
3
由于n
"
m,故h(t)中含有aδ(t)项因此:
h(t)=δ(t)+(-43e
-2t+13e
-5t)u(t)
2)单位阶跃响应
g″(t)+7g′(t)+10g(t)=δ′(t)+6δ(t)+4u(t)
g′(0-)=g(0-)=0
全响应为:g(t)=C1e
-2t+C2e
-5t+B (t0+)
根据特解形式得:
10B=4
—02—
特解:gp(t)=
2
5
根据冲激函数匹配法设:
d2
dt2
g(t)=aδ′(t)+bδ(t)+cΔu(t)
d
dth(t)=aδ(t)+bΔu(t)
g(t)=aΔu(t
)
(0-<t<0+)
代入方程解得:a=1,b=-1,c=1故:
g(0+)=a+g(0-)=1
g′(0+)=b+g′(0-)=-
{ →
1
C1 =2/3
C2 =-1/
{ 15
全响应:g(t)=(23e
-2t-115e
-5t+25)u(t)
五、卷积
原理:将信号进行冲击分解,根据系统h(t)求任意激励信号的零状态响应
f1(t)f2(t)∫
!
-
!
f1(τ)f2(t-τ)dτ=f2(t)f1(t)
则系统的零状态响应为:r(t)=e(t)h(t)=∫
!
-
!
e1(τ)h(t-τ)dτ
计算:
1)转换自变量
2)反折
3)平移
4)相乘
5)积分
例1 计算图示信号的卷积
(图见视频)
1)-
!
<t
"
-0.5时
—12—
郑君里《信号与系统》考点精讲
e(t)h(t)=0
2)-0.5
"
t
"
1且t-2
"
-0.5时
e(t)h(t)=∫
t
-0.5
1×12(t-τ)dτ=
t2
4+
t
4+
1
16
3)t1且t-2"-0.5时,1"t"32
e(t)h(t)=∫
1
-0.5
1×12(t-τ)dτ=
3t
4+
3
16
4)-12"
t-2
"
1且t1时,32"
t
"
3
e(t)h(t)=∫
1
t-2
1×12(t-τ)dτ=-
t2
4+
t
2+
3
4
5)t-21且t1时,e(t)h(t)=0
综上可得e(t)h(t)如下图所示:
六、卷积的性质
熟练掌握卷积各种性质的数学证明
1.卷积代数
1)交换律
2)分配律
3)结合律
f1(t)f2(t)=f2(t)f1(t)
f1(t)[f2(t)+f3(t)]=f1(t)f2(t)+f1(t)f3(t)
[f1(t)f2(t)]f3(t)=f1(t)[f2(t)f3(t)]
系统串并联
(图见视频)
2)卷积的微分与积分特性
d
dt[f1(t)f2(t)]=f1(t)
d
dtf2(t)=f2(t)
d
dtf1(t)
∫
t
-
!
[f1(λ)f2(λ)]dλ=f1(t)∫
t
-
!
f2(λ)dλ=f2(t)∫
t
-
!
f1(λ)dλ
设:s(t)=f1(t)f2(t)
—22—
则:s(i)(t)=f(j)1 (t)f
(i-j)
2 (t)
3)与冲激函数或阶跃函数的卷积
f(t)δ(t)=∫
!
-
!
f(τ)δ(t-τ)dτ
f(t)δ(t)=∫
!
-
!
f(τ)δ(t-τ)dτ
筛选性 =∫
!
-
!
f(τ)δ(τ-t)dτ=f(t)
信号与单位冲激信号的卷积为本身
f(t)δ(t-t0)=∫
!
-
!
f(τ)δ(t-t0-τ)dτ=f(t-t0)
结论:与δ(t-t0)卷积的结果是延迟个时间单位。
利用卷积的微分、积分特性可得:
f(t)δ′(t)=f′(t)
f(t)u(t)=∫
t
-
!
f(τ)dτ
一般形式:
f(t)δ(k)(t)=f(k)(t)
f(t)δ(k)(t-t)=f(k)(t-t0)
例2 利用卷积特性重新计算计算例1中的卷积
解:根据微分、积分特性得
r(t)=e(t)h(t)=ddte(t)∫
t
-
!
h(τ)dτ
∫
t
-
!
h(τ)dτ=∫
t
-
!
1
2τ[u(τ)-u(τ-2)]dτ
=(∫
t
-
!
1
2τdτ)u(t)-(∫
t
-
!
1
2τdτ)u(t-2)
=14t
2[u(t)-u(t-2)]+u(t-2)=h(-1)(t)
d
dte(t)=δ(t+
1
2)-δ(t-1)
r(t)=e(t)h(t)=ddte(t)∫
t
-
!
h(τ)dτ
= 1
4(t+
1
2)
2
[u(t+12)-u(t-
3
2)]+u(t-
3
2{ })- 1
4(t-1)
2[u(t-1)-u(t-3)]+u(t-3{ })
—32—
郑君里《信号与系统》考点精讲
=
0.25(t+0.5)2 -0.5
"
t<1
0.75(t-0.25) 1
"
t<1.5
1-0.25(t-1)2 1.5
"
t<
{
3
可得:
七、微分方程的算子解法
基本算子:p=ddt
1
p=∫
t
-
!
(…)dτ
系统的微分方程可用算子法描述为:
(C0p
n+C1p
n-1+… +Cn-1p+Cn)r(t)=(E0p
m +E1p
m-1+… +Em-1p+Em)e(t)
1.基本规则
1)算子公因子不能对消
2)算子乘除顺序不能颠倒
先“乘”后“除”不抵消
p1px≠
1
ppx
先“除”后“乘”可抵消
2.算立法简历微分方程
vL(t)=L
d
dtiL(t)=Lp·iL(t) 等效电感
vC(t)=
1
C∫
t
-
!
iC(τ)dτ=
1
Cp·iC(t) 等效电容
例3 利用算子法建立系统方程。
解:用算子符号表示电路元件
—42—
(图见视频)
根据KVL:
(R1+
1
Cp)i(t)-
1
CpiL(t)=e(t)
-1Cpi(t)+
1
Cp(Lp+R2+
1
Cp)iL(t)=
{ 0
应用克莱姆法则解得:
i(t)=
(
1
R1
p+
R2
R1L
+ 1
R1LCp
)
p+(
R2
L+
1
R1C
)+(1LC+
R2
R1LC
)
1
p
允许先积分后微分,同乘p得:
i(t)=
(
1
R1
p2+
R2
R1L
p+ 1
R1LC
)
p2+(
R2
L+
1
R1C
)p+(1LC+
R2
R1LC
)
·e(t)
代入参数整理后得:
d2
dt2
i(t)+7ddti(t)+10i(t)=
d2
dt2
e(t)+6ddte(t)+4e(t)
八、小结
1.LTI连续系统的响应:
全响应=齐次解(自由响应)+特解(强迫响应)
2.关于0-和0+初始值
当系统已经用微分方程表示时,如果包含有
#
(t)及其各阶导数,说明相应的0-状态到0+状态
发生了跳变。冲激函数匹配法。
3.零输入响应和零状态响应
y(t)=yh(t)+yp(t)
自由响应+强迫响应;暂态响应+稳态响应;
零输入响应+零状态响应
4.冲激响应和阶跃响应
5.卷积积分
定义,计算的5个过程
6.卷积积分的性质
交换律、分配率、结合律
利用性质求卷积
—52—
郑君里《信号与系统》考点精讲
第三章 傅立叶变换
一、引言
通过本章的学习掌握以下内容:
1.周期信号和非周期信号频谱分析的方法
———傅立叶级数和傅立叶变换
2.周期信号和非周期信号频谱的特点与区别
3.理解非周期信号频谱密度函数的概念
4.信号时域特性与频域特性之间的关系
5.采样定理
要点:典型信号的傅立叶变换,利用傅立叶变换的
性质求信号的正、反变换,信号的频谱分析。
二、周期信号的傅立叶级数分析
1.三角函数形式的傅里叶级数
(1)狄利克雷条件
在一周期内,如果有间断点存在,则间断点的数目应该是有限个;
在一周期内,极大值和极小值的数目应该是有限个;
在一周期内,信号是绝对可积的。
周期信号可以由一系列不同角频率的三角函数的线性组合来逼近:
(2)描述
f(t)=a0+a1cos(ω1t)+b1sin(ω1t)+a2cos(2ω1t)+b2sin(2ω1t)+…+ancos(nω1t)+bnsin(nω1t)+…
=a0+∑
!
n=1
[ancos(nω1t)+bnsin(nω1t)]
—62—
a0 =
1
T1∫
t0+T1
t0
f(t)dt—直流分量幅值
an =
2
T1∫
t0+T1
t0
f(t)cos(nω1t)dt—余弦分量幅值
bn =
2
T1∫
t0+T1
t0
f(t)sin(nω1t)dt———正弦分量幅值
合并、整理可得:
f(t)=c0+∑
!
n=1
cncos(nω1t+φn)
其系数如下:
a0 =c0,cn = a2n+b
2
槡 n,φn =-arctan
bn
an
sinφn =
-bn
a2n+b
2
槡 n
,cosφn =
an
a2n+b
2
槡 n
an =cncosφn,bn =-cnsinφ
n
直流分量、基波、二次谐波、三次谐波……
2.复指数形式的傅里叶级数
欧拉公式:
cosnω0 =
1
2(e
jnω0+e-jnω0)
sinnω0 =
1
j2(e
jnω0+e-jnω0 }
)
or e±jnω0 =cosnω0±jsinnω0
代入三角函数形式的傅立叶级数可得:
f(t)=c0+∑
!
n=1
cncos(nω0t+φn)
=c0+∑
!
n=1
cn
2[e
j(nω0t+φn)+e-j(nω0t+φn)]
=c0+∑
!
n=1
cn
2e
jnω0tejφn+∑
!
n=1
cn
2e
-jnω0te-jφn
=c0+∑
!
n=1
cn
2e
jnω0tejφn+∑
-
!
n=-1
c-n
2e
-jnω0te-jφ-n
=c0+∑
!
n=1
cn
2e
jnω0tejφn+∑
-1
n=-
!
cn
2e
jnω0tejφn
整理后可得复指数形式的傅立叶级数:
f(t)=∑
!
n=-
!
cn
2e
jnω1tejφn =∑
!
n=-
!
F(nω0)e
jnω1ta=∑
!
n=-
!
Fne
jnω1t
—72—
郑君里《信号与系统》考点精讲
其中:Fn =F(nω1)=
1
T∫
t0+T
t0
f(t)e-jnω1tdt
比较上述两种形式的傅立叶级数,可得到二者系数之间的关系:
F0 =a0 =c0
Fn = Fn e
jφn =12(an+jbn)=
1
2cne
jφn
F-n =
1
2(an+jbn)=
1
2cne
-jφn
Fn =
1
2cn = F-n
φn =-arctan
bn
an
Fn+F-n =2Re[Fn]=an
j(Fn-F-n)=j2Im[Fn]=b
n
三角函数形式:cn ~ω,φn ~ω 单边频谱
复指数的形式:Fn ~ω,φn ~ω 双边频谱
二者之间关系:F(nω1) =
1
2cn(n≠0),F0 =c0 =a0
幅值谱偶对称:F(nω1) = F(-nω1)
相位谱奇对称:φ(nω1)=-φ(-nω1)
3.周期信号的频谱及其特点
幅值、相角与角频率之间的关系分别称为幅值谱、相位谱,统称频率特性或频谱。
F(nω1)= F(nω1)e
jφn
幅频特性:F(nω1) =
1
2 a2n+b
2
槡 n =
1
2cn
相频特性:φn =arctan(
-bn
an
)
离散性、谐波性、收敛性引入负频率,无物理意义
例1 试做出信号:f(t)=1+sinω1t+2cosω1t+cos(2ω1t+
π
4)的两种形式的频谱。
解:整理原式得:f(t)=1+槡5cos(ω1t-0.15π)+cos(2ω1t+
π
4)
c0 =1 c1 =槡5=2.236 c2 =1
φ0 =0 φ1 =-0.15π φ2 =0.25π
—82—
考试点(www.kaoshidian.com)名师精品课程 电话:400-6885-365
利用欧拉公示化简得:
f(t)=1+(1.12e-j0.15πejω1t+1.12ej0.15πe-jω1t)+(0.5ej0.25πej2ω1t+0.5e-j0.25πe-j2ω1t)
=∑
2
n=-2
F(nω1)e
jnω1t
4.函数的对称性与傅立叶系数的关系
当波形满足某种对称关系时,傅里叶级数中的某些特定项将不会出现,利用这种性质可以对谐波
成分迅速作出判断,以简化傅立叶系数的计算。常用的有:
周期对称:奇函数、偶函数
f(t)=-f(-t) f(t)=f(-t)
半周期对称:奇谐函数、偶谐函数
f(t)=-f(t±T2) f(t)=f(t±
T1
2)
(1)奇函数
an =0
bn =
4
T∫
T
2
0
f(t)sin(nω1t)dt
Fn =-F-n =-
1
2jbn
φn =-
π
2
不含余弦项,只含有正弦项,傅立叶系数为虚数。
(2)偶函数
an =
4
T∫
T
2
0
f(t)cos(nω1t)dt
bn =0
Fn =F-n =
1
2an =
1
2cn
φn =0
不含正弦项,只含有余弦项和直流项,傅立叶系数为实函数。
(3)奇谐函数
n为奇数时:
an =
4
T∫
T
2
0
f(t)cos(nω1t)dt
bn =
4
T∫
T
2
0
f(t)sin(nω1t)dt
—92—
郑君里《信号与系统》考点精讲
(4)偶谐函数
n为偶数时:
an =
4
T∫
T
2
0
f(t)cos(nω1t)dt
bn =
4
T∫
T
2
0
f(t)sin(nω1t)dt
5.傅里叶有限级数与最小方均误差,吉布斯现象
f(t)=a0+∑
!
n=1
[ancos(nω1t)+bnsin(nω1t)]
SN =a0+∑
N
n=1
[ancos(nω1t)+bnsin(nω1t)]
εN(t)=f(t)-SN 误差项
EN =εN
2(t)=1T1∫
t0+T1
t0
εN
2(t)dt
EN =εN
2(t)=f2(t)-[a0
2+12∑
N
n=1
(an
2+bn
2)]
有限项逼近所产生的方均误差
由上式可知:lim
N→!
SN =f(t)
对于对称方波的逼近分析发现:
(1)N越大,越接近方波
(2)快变信号,高频分量,主要影响跳变沿;
(3)慢变信号,低频分量,主要影响顶部;
(4)任一分量的幅度或相位发生相对变化时,波形将会失真。
相数越多,合成波形中出现的峰起越靠近原函数的不连续点,其值约为总跳变量的9%。称为“吉
布斯现象”。
—03—
三、典型周期信号的傅立叶级数
1.周期矩形脉冲信号
其在一个周期内的表达式为:
f(t)=[u(t+π2)-u(t-
π
2)]
偶函数,bn =0只有a0、an项。
F(nω1)=
1
T1∫
T1
2
-
T1
2
f(t)e-jnω1tdt
=1T1∫
τ
2
-τ2
Ee-jnω1tdt=EτT1
Sa(nω1
τ
2)
求得:
f(t)=EτT1+
Eτω1
π ∑
!
n=1
Sa(nω1τ2 )cos(nω1t)
f(t)=EτT1∑
!
n=-
!
Sa(nω1τ2 )e
jnω1τ
F(nω1)为实数Fn >0,相位0,Fn <0,相位 ±π
当T1→ !
时,ω1→0,
Eτ
T1
为无限小,
f(t)由周期信号→非周期信号
信号能量的主要部分称为带宽,记为:
Bω =
2π
τ
或Bf=
1
τ
—13—
郑君里《信号与系统》考点精讲
当上述矩形信号为方波时,既是偶函数,又是奇谐函数,故只有基波和奇次余弦谐波分量
f(t)=2E
π
[cos(ω1t)-13cos(3ω1t)+
1
5cos(5ω1t)-…]
谐波幅值以1/n收敛
由矩形脉冲信号的频谱特性印证了周期信号频谱的特点:离散性、谐波性、收敛性。
2.周期锯齿脉冲
f(t)=E
π∑
!
n=1
(-1)n+1 1nsin(nω1t)
只包含正弦分量,谐波幅值以1/n收敛
3.周期三角脉冲
f(t)=E2+
4E
π2∑
!
n=1
1
n2
sin2(nπ2)cos(nω1t)
包含直流分量、基波以及奇次谐波分量,谐波幅值以1/n2收敛。
4.周期半波余弦信号
f(t)=E
π
-2E
π∑
!
n=1
1
(n2-1)
cos(nπ2)cos(nω1t)
包含直流分量、基波以及各次谐波分量,
谐波幅值以1/n2收敛。
四、傅立叶变换
1.由傅立叶级数到傅立叶变换
T1→ !
lim
T→!
fT(t)=f(t)
—23—
F(nω1)=
1
T1∫
T1
2
-
T1
2
f(t)e-jnω1tdt→0
频谱间隔:ω1 =
2π
T1
→0
→离散谱 连续谱
谱系数Fn =F(nω1)无意义,但仍有相当对大小关系。
谱系数两侧同时乘以周期
T1F(nω1)=
F(nω1)
1/T =
F(nω1)
f 单位频带上的频谱
T1→ !
时,f=1T1
→0,F(nω1)→0
谱线间隔:Δ(nω1)=ω1→dω(nω1)→ω
F(ω)=lim
T1→!
T1F(nω1)=limT1→!
∫
T1
2
-
T1
2
f(t)e-jnω1tdt=∫
!
-
!
f(t)e-jnωtdt 频谱密度函数
称上式为信号的傅立叶变换
F(ω)=|F(ω)|ejφ(ω)
F(ω) ~ω:幅值谱
φ(ω)~ω:相位谱
2.傅立叶反变换
F(jω)=∫
!
-
!
f(t)e-jωtdt
fT(t)=∑
!
n=-
!
1
T[∫
T
2
-T2
fT(t)e
-jnω1tdt]ejnω1t
fT(t)=∑
!
n=-
!
ω1
2π
[∫
T
2
-T2
fT(t)e
-jnω1tdt]ejnω1t
f(t)=∫
!
-
!
1
2π
(∫
!
-
!
f(t)e-jωtdt)ejωtdω
f(t)= 12π∫
!
-
!
F(jω)ejωtdω
记:F(ω)=∫
!
-
!
f(t)e-jωtdt=F[f(t)]
f(t)= 12π∫
!
-
!
F(jω)ejωtdω=F-1[F(jω)]
分别称为傅立叶变换和傅立叶反变换
3.傅立叶变换的物理意义
f(t)= 12π∫
!
-
!
F(ω)ejωtdω
= 12π∫
!
-
!
F(ω)ejφ(ω)ejωtdω
—33—
郑君里《信号与系统》考点精讲
= 12π∫
!
-
!
F(ω)cos[ωt+φ(ω)]dω+j
1
2π∫
!
-
!
F(ω)sin[ωt+φ(ω)]dω
=1
π∫
!
0
F(ω)cos[ωt+φ(ω)]dω
=∫
!
0
F(ω)
π
dωcos[ωt+φ(ω)]
f(t)=∫
!
0
F(ω)
π
dωcos[ωt+φ(ω)]
求和 幅值 余弦信号
无穷多个振幅为(
1
π F(ω)
dω)余弦信号之和
f(t)= 12π∫
!
-
!
F(ω)ejωtdω=∫
!
-
!
F(ω)
2π
dω·ejωt
求和 幅值 复指数信号
无穷多个振幅为(
1
2π F(ω)
dω)复指数信号之和
4.傅立叶变换存在的条件
信号绝对可积∫
!
-
!
f(t)dt<!
所有能量信号均满足此条件。
充分非必要条件。借助奇异函数的概念,许多不满足此条件的信号都存在傅立叶变换。
五、典型非周期信号的傅立叶变换
1.单边指数信号
f(t)=
e-αtt>0α>0
0t<{ 0
F(ω)=F[f(t)]=∫
!
-
!
e-αtu(t)e-jωtdt=∫
!
0
e-(α+jω)tdt= 1
α+jω
F(ω) =
1
α2+ω槡
2
ω=0,F(ω) =
E
α
ω→±!
,F(ω)→
{
0
φ(ω)=-arctanωα
ω→0,φ(ω)=0
ω→+!
,φ(ω)→-π2
ω→-!
,φ(ω)→ π
2
—43—
2.双边指数信号
f(t)=e-a|t|(-
!
<t<+
!
) a>0
F(ω)=∫
!
-
!
f(t)e-jωtdt=∫
!
0
e-(a+jω)tdt
F(ω)= 1
a-jω
+ 1
a+jω
= 2a
a2+ω2
F(ω)= 2a
a2+ω2
= F(ω)
φ(ω)=
}
0
3.矩形脉冲信号
f(t)=E[u(t+τ2)-u(t-
τ
2)]
F(ω)=∫
τ
2
-τ2
Ee-jωtdt=EτSa(ωτ2)
F(ω) =EτSa(ωτ2)
φ(ω)=
0 (4nπ
τ
< ω <
2(2n+1)π
τ
)
π (2(2n+1)πτ
< ω <
4(n+1)π
τ
{
)
4.钟形脉冲信号(高斯信号)
f(t)=Ee-(
t
τ)2(-
!
<t<+
!
)
F(ω)=槡πEτe
-(ωτ2)2
5.符号函数
f(t)=sgn(t)=
+1,t>0
-1,t<{ 0
—53—
郑君里《信号与系统》考点精讲
处理方法:
衰减因子
f1(t)=sgn(t)e
-α t
求F1(ω),求极限得到F(ω)
F1(ω)=∫
0
-
!
-eαte-jωtdt+∫
!
0
e-αte-jωtdt
= -1
α-jω
+ 1
α+jω
= -j2ω
α2+ω2
F(ω)=lim
α→0
F1(ω)=lim
α→0
-j2ω
α2+ω2
=2jω
sgn(t) 2jω
=-j2
ω
= 2
ω
ej
π
2
F(ω) =( (
2
ω
)
槡
2
= 2
ω
)
F(ω)是偶函数
arctan
-2
ω
0 =
-π2,ω>0
π
2,ω<
{ 0
φ(ω)是奇函数
6.升余弦脉冲信号
f(t)=E2[1+cos(
πt
τ
)]0
" t"τ
F(ω)=∫
!
-
!
f(t)e-jωtdt
=∫
τ
-τ
E
2[1+cos(
πt
τ
)]e-jωtdt
=E2∫
τ
-τ
e-jωtdt+E4∫
τ
-τ
ej
πt
τe-jωtdt+E4∫
τ
-τ
e-j
πt
τe-jωtdt
=EτSa(ωτ)+Eτ2Sa[(ω-
π
τ
)τ]+Eτ2Sa[(ω+
π
τ
)τ]
(图见视频)
六、冲激函数和阶跃函数的傅立叶变换
1.冲击函数的傅立叶变换
F(ω)=∫
!
-
!
δ(t)e-jωtdt
—63—
2.冲激偶函数的傅立叶变换
∫
!
-
!
f(t)δ′(t)dt=-f′(0)
F[δ′(t)]=∫
!
-
!
δ′(t)e-jωtdt=-[e-jωt]′t=0 =-(-jω)=jω
3.阶跃函数的傅立叶变换
u(t)=12+
1
2sgn(t)
1
2πδ(ω)
1
2sgn(t)
1
jω
u(t)πδ(ω)+1jω
—73—
郑君里《信号与系统》考点精讲
七、傅立叶变换的基本性质
对称性、线性(叠加性)、奇偶虚实性、尺度变换特性、时移特性、频移特性、微分特性、积分特性。
重点掌握在基本的傅里叶变换对基础上利用上述八个性质对一般信号进行付傅立叶变换的方
法,了解傅里叶变换在通讯系统领域中的应用。
δ(t)1
u(t) 1jω
+πδ(ω)
u(t)e-αt 1
α+jω
ejω0t2πδ(ω-ω0)
e-at2 π
槡a
e-
ω2
4a
1.对称性
若:f(t)F(ω),则:F(t)2πf(-ω)
特殊的:f(t)为偶函数时,F(t)2πf(ω)
利用已知变换对方便的求出信号的傅里叶变换。
(图见视频)
2.线性(叠加性和齐次性)
若:f1(t)F1(ω),f2(t)F2(ω)
则:c1f1(t)+c2f2(t)c1F1(ω)+c2F2(ω)
傅里叶正、反变换公式为求积分(求和)复杂信号分解后求和。
3.奇偶虚实性
F(ω)=∫
!
-
!
f(t)e-jωtdt
=∫
!
-
!
f(t)cosωtdt-j∫
!
-
!
f(t)sinωtdt
=R(ω)+jX(ω)=F(ω)= F(ω)e-jφ(ω)
其中:
R(ω)=∫
!
-
!
f(t)cosωtdt=R(-ω)
X(ω)=-j∫
!
-
!
f(t)sinωtdt=-X(-ω)
F(ω) = R2(ω)+X2(ω槡 )
φ(ω)=arctanX(ω)R(ω)
=-φ(-ω
)
为ω的偶函数
为ω的奇函数
为ω的偶函数
为ω的奇函数
结论:
—83—
f(-t)F(-ω)
f(t)F(-ω)
f(t)为实函数时,f(-t)F(-ω)=F(ω)
4.尺度变换
若:f(t)F(ω)
则:f(at) 1
aF(
ω
a),a≠0
意义:
→时域压缩 频域扩展
→时域扩展 频域压缩
→时域对折 频域对折
信号的等效脉冲宽度与其等效带宽成反比(图见视频)
5.时移特性
若f(t)F(ω),则f(t-t0)F(ω)e
-jωt0;
设F(ω)= F(ω)ejφ(ω),则f(t-t0) F(ω)·ej[φ(ω)-ωt0]
不影响幅度频谱,只影响相位频谱
6.频移特性
若f(t)F(ω)则f(t)e±jω0tF(ωω0)
通信中调制与解调,频分复用。
时域f(t)乘以ejω0t,频谱右移ω0
时域f(t)乘以e-jω0t,频谱左移ω0
7.微分特性
若:f(t)F(ω)
则:
d(n)f(t)
dtn
(jω)nF(ω) 时域微分
d(n)F(ω)
dωn
(-jt)nf(t) {
频域微分
应用:u(t) 1jω
+πδ(ω) δ(t)1 δ′(t)F(ω)
8.积分特性
若:f(t)F(ω)
则:∫
t
-
!
f(τ)dτF(ω)jω
+πF(0)δ(ω)
例1 求图示信号的频谱
—93—
郑君里《信号与系统》考点精讲
解:已知矩形信号的频谱为采样信号,即
f0(t)F0(ω)=EτSa(
ωτ
2)
(图见视频)
而三脉冲矩形信号可表示为:
f(t)=f0(t)+f0(t+T)+f0(t-T)
根据傅里叶变换的线性特性和时移特性:
F(ω)=F0(ω)(1+e
jωT+e-jωT)=Eτ·Sa(ωτ2)[1+2cos(ωT)]
(图见视频)
与矩形信号相比:
脉冲数目增多时,幅频特性包络线不变,带宽不变。
例2 试求双Sa信号的频谱。
解:由例1和傅里叶变换的对称性知:采样信号
f0(t)=
ωc
π
Sa(ωct)的频谱为矩形。如图:
F[f0(t)]=
1(ω <ωc)
0(ω <ωc
{
)
根据傅里叶变换的时移特性知:
F[f0(t-2τ)]=
e-j2ωt(ω <ωc)
0(ω <ωc
{
)
根据傅里叶变换的线性特性得双Sa信号的频谱:
F(ω)=F[f0(t)]-F[f0(t-2τ)]=
1-e-j2ωτ(ω <ωc)
0(ω <ωc
{
)
其幅频特性为:
F(ω) =
2sin(ωτ)(ω <ωc)
0(ω <ωc
{
)
—04—
由于实际中常取τ=πωc
故:
F(ω) =
2sin(πω
ωc
)(ω <ωc)
0(ω <ωc
{
)
因此,双Sa信号及其频谱如下:
与单Sa信号相比,没有直流分量,便于传输。
例3 设矩形脉冲信号G(t)的高度为E,宽度为τ,试求矩形调幅信号f(t)=G(t)cos(ω0t),的频谱。
解:G(ω)=Eτ·Sa(ωτ2) f(t)=
1
2G(t)(e
jω0t+e-jω0t)
根据傅里叶变换的频移特性得原信号的频谱为:
(图见视频)
F(ω)=12G(ω-ω0)+
1
2G(ω+ω0)=
Eτ
2Sa[
(ω-ω0)τ
2 ]+Eτ2Sa[
(ω+ω0)τ
2 ]
例4 f(t)=
E(1-2
τ t
)(t<τ2)
0(t>τ2
{
)
,求三角脉冲信号的频谱
解:取原信号的一阶及二阶导数得:
(图见视频)
df(t)
dt =
2
τ
E(-τ2 <t<0)
-2
τ
E(0<t<τ2)
0(t>τ2
)
d2f(t)
dt2
=2E
τ
[δ(t+τ2)+δ(t-
τ
2)-2δ(t)]
二者的时域波形如下图:
—14—
郑君里《信号与系统》考点精讲
根据傅里叶变换的微分特性可得:
F[f″(t)]=∫
!
-
!
[
2E
τδ
(t+τ2)-
4E
τδ
(t)+2E
τδ
(t-τ2)]e
-jωtdt
=2E
τ
ejω
τ
2 -4E
τ
+2E
τ
e-jω
τ
2 =(jω)2F(ω)=-ω2F(ω)
整理可得:
F(ω)= 1
-ω2
[
2E
τ
ejω
τ
2 -4E
τ
+2E
τ
e-jω
τ
2]
= 1
-ω2
2E
τ
[ejω
τ
2 -2+e-jω
τ
2]
=-2E
τω2
[ejω
τ
4 -e-jω
τ
4]2 =-2E
τω2
(2jsinωτ4)
2
=8E
τω2
(sinωτ4)
2(
ωτ
4)
2
(
ωτ
4)
2 =
τE
2Sa(
ωτ
4)
2
(图见视频)
例5 求图示信号的傅里叶变换
解:(图见视频)
根据傅里叶变换的积分特性:
F2(ω)=[πδ(ω)+
1
jω
]Sa(ω2)e
-jω2 =πδ(ω)+
Sa(ω2)e
-jω2
jω
F(ω)=F[1]+F2(ω)=3πδ(ω)+
Sa(ω2)e
-jω2
jω
9.卷积特性
1)时域卷积定理
若:f1(t)F1(ω),f2(t)F2(ω)
则:f1(t)f2(t)F1(ω)·F2(ω)
F[f1(t)f2(t)]=∫
!
-
!
[∫
!
-
!
f1(τ)f2(t-τ)dτ]e
-jωtdt
=∫
!
-
!
f1(τ)[∫
!
-
!
f2(t-τ)e
-jωtdt]dτ
=∫
!
-
!
f1(τ)F2(ω)e
-jωtdt
—24—
=F1(ω)F2(ω)
2)频域卷积定理
若:f1(t)F1(ω),f2(t)F2(ω)
则:f1(t)·f2(t)
1
2π
F1(ω)F2(ω)
卷积定理揭示了时域运算与频域运算的对应关系:
时域乘积对应频域卷积,时域卷积对应频域乘积。
系统的变换域分析法的基础。
例6 求下述余弦脉冲信号的频谱。
f(t)=
Ecos(πt
τ
)(t"τ2)
0(t>τ2
{
)
解:(图见视频)
将余弦信号用窗函数截取,即可得到
余弦脉冲信号:f(t)=G(t)cos(πt
τ
)
F(ω)=F[G(t)·cos(πtτ
)]
= 12π
·[EτSa(ωτ2)][πδ(ω+
π
τ
)+πδ(ω-πτ
)]
=Eτ2Sa[(ω+
π
τ
)
τ
2]+
Eτ
2Sa[(ω-
π
τ
)
τ
2]
考点:掌握各种基本信号的傅里叶变换对,灵活运用傅里叶变换的各种特性对一般信号进行频谱
分析。
八、周期信号的傅立叶变换
满足狄利克雷条件的周期信号
fp(t)→傅里叶级数F(nω1)离散谱
满足绝对可积条件的非周期信号
f(t)→傅里叶变换F(ω)连续谱
一般的周期信号是否可进行傅里叶变换?
fp(t)→ ∑
!
n=-
!
Fne
jnω1t或c0+∑
!
n=-
!
cncos(nω1t+φn)
存在傅里叶变换
单位冲击信号的引入,使绝对可积成为非必要条件,
一般的周期信号存在FT。
—34—
郑君里《信号与系统》考点精讲
1.正、余弦信号的傅立叶变换
已知12πδ(ω)根据频移特性可得:
1·ejω1t2πδ(ω-ω1)
1·e-jω1t2πδ(ω+ω1)
应用欧拉公式,根据线性特性可得:
cosω1t
1
2[2πδ(ω-ω1)+2πδ(ω+ω1)]
=πδ(ω+ω1)+πδ(ω-ω1)
sinω1t
1
2j[2πδ(ω-ω1)-2πδ(ω+ω1)]
=jπδ(ω+ω1)-jπδ(ω-ω1)
根据上节例3,对矩形脉冲信号求极限,并根据单位冲击信号的采样信号定义:δ(t)=lim
k→!
k
π
Sa(kt)
(图见视频)
2.一般周期信号的傅立叶变换
FT(ω)=F[fT(t)]
=F[∑
!
-
!
F(nω1)e
jnω1t]
=∑
!
-
!
F(nω1)F[e
jnω1t]
=2π∑
!
-
!
F(nω1)δ(ω-nω1)
其中:F(nω1)=
1
T∫
T1
2
-
T1
2
f(t)e-jnω1tdt
从上式中可以看出:周期信号的傅里叶变换是由一系列位于信号谐波频率的冲击函数组成的,冲
击信号的幅值等于原信号傅里叶级数相应系数的2π倍。
例1 求图示周期信号的傅里叶变换
解:
Fn =
1
T∫Tf(t)·e-jω1tdt
=12∫
3
2
-12
[G1
2
(t)-G1
2
(t-1)]e-jnπtdt
—44—
=
sinnπ4
nπ
[1-(-1)n]
F(ω)=F[f(t)]
=2π∑
!
n=-
!
Fnδ(ω-nπ)
=2∑
!
n=-
!
sinnπ4
n [1-(-1)n]δ(ω-nπ)
3.单脉冲信号的频谱密度与周期信号谱系数的关系
设单个脉冲信号f0(t)周期拓展后得到:fT(t)
f0(t)=
1
2π∫
!
-
!
F0(ω)e
jωtdω
F0(ω)=∫
!
-
!
f0(t)e
-jωt{ dt
fT(t)=∑
!
n=-
!
F(nω1)e
jnω1t
F(nω1)=
1
T1∫
T
2
-T2
fT(t)e
-jnω1t{ dt
(图见视频)
比较后可得二者关系:F(nω1)=
1
T1
F0(ω)
ω=nω1
例2 求周期单位冲激序列的傅里叶变换
δT(t)=∑
!
n=-
!
δ(t-nT1)δ(t)1
可得δT(t)的傅里叶级数谱系数
F(nω1)=
1
T1
F0(ω)
ω=nω1
=1T1
δT(t)=∑
!
n=-
!
F(nω1)e
jnω1t=1T1∑
!
n=-
!
ejnω1t
F(ω)=F[δT(t)]
=1T1∑
!
n=-
!
F[ejnω1t]
=1T1∑
!
n=-
!
2πδ(ω-nω1)
=ω1∑
!
n=-
!
δ(ω-nω1)=δω(ω)
=2πT1∑
!
n=-
!
δ(ω-nω1)
(图见视频)
—54—
郑君里《信号与系统》考点精讲
例3 求周期矩形脉冲序列的傅里叶变换。
F0(ω)→F(nω1)→F(ω)
F0(ω)=EτSa(
ωτ
2)
F(nω1)=
1
T1
F0(ω)
ω=nω1
F(ω)=2π∑
!
n=-
!
F(nω1)·δ(ω-nω1)
=2π∑
!
n=-
!
Eτ
T1
Sa(
nω1τ
2 )δ(ω-nω1)
=Eτω1∑
!
n=-
!
Sa(
nω1τ
2 )δ(ω-nω1)
(图见视频)
在定义了单位冲激信号之后,周期信号的傅里叶级数和傅里叶变换得到了统一,掌握由傅里叶级
数的系数求傅里叶变换的方法。
九、抽样信号的傅立叶变换
(图见视频)
1.时域抽样
周期抽样脉冲的频谱:P(ω)=2π∑
!
n=-
!
Pnδ(ω-nωs)
其中:Pn =
1
Ts∫
Ts
2
-
Ts
2
p(t)e-jωstdt
设采样器模型为:fs(t)=f(t)p(t)
根据频域卷积定理得:
Fs(ω)=2πF(ω)P(ω)=∑
!
n=-
!
PnF(ω-nωs)
可见,采样之后所的信号的频谱是原始信号频谱的加权和,权系数为采样脉冲的傅里叶系数。
1)矩形脉冲抽样(自然抽样)
Pn =
1
Ts∫
Ts
2
-
Ts
2
p(t)e-jnωstdt=1Ts∫
Ts
2
-
Ts
2
Ee-jnωstdt=EτTs
Sa(
nωsτ
2 )
—64—
Fs(ω)=
Eτ
Ts∑
!
n=-
!
Sa(
nωsτ
2 )F(ω-nωs)
将原始信号的频谱按照采样信号的规律重复。
自然抽样的频谱关系
(图见视频)
2)冲激抽样(理想抽样)
p(t)=δT(t)=∑
!
n=-
!
δ(t-nTs)ωs∑
!
-
!
δ(ω-nωs)
fs(t)=f(t)δT(t)=∑
!
n=-
!
f(nTs)δ(t-nTs)
Fs(ω)=F[f(t)δT(t)]=
1
2π
F(ω)δT(ω)=
1
Ts∑
!
n=-
!
F(ω-nωs)
将原始信号的频谱按照等幅的规律重复。
理想抽样的频谱关系
(图见视频)
2.频域抽样
对频域信号进行理想抽样,研究二者所对应的时域信号之间的关系。
F1(ω)=F(ω)δω(ω)δω(ω)=∑
!
n=-
!
δ(ω-nω1)
F[∑
!
n=-
!
δ(t-nT1)]=ω1∑
!
n=-
!
δ(ω-nω1), (ω1=2πT1)
F-1[∑
!
n=-
!
δ(ω-nω1)]=1ω1∑
!
n=-
!
δ(t-nT1)=F-1[δω(ω)]=
δT(t)
ω1
f1(t)=f(t)F-1[δω(ω)]=
1
ω1∑
!
n=-
!
f(t-nT1)
频域理想抽样的时域对应关系(时域和频域离散的性与周期性相互对应关系)。
(图见视频)
十、抽样定理
1.时域抽样定理
—74—
郑君里《信号与系统》考点精讲
从上图可知:ωs2ωm或fsfωm时,采样所得离散信号不发生频谱混叠。鉴于傅里叶变换一一
对应的关系,离散信号中保留了连续信号的所有特征,采样有效,可通过理想低通滤波器无失真的
恢复。
2.频域抽样定理
对于时间受限信号 t"tm,对频域的采样间隔时,采样所得离散频谱所对应的时间信号不发生时
域混叠,采样有效,用窗函数从中选出单个脉冲即可无失真的恢复出原来的时域信号。
—84—
2.名词定义
状态:对于一个动态系统的装态是表示系统的一组最小变量,只要知道t=t0时这组变量和tt0
时的输入,那么就能完全确定系统在任何时间tt0的行为。
状态变量:能够表示系统状态的那些变量。
状态矢量:能够完全描述一个系统行为的k个状态变量,可以看做矢量λ(t)的各个分量的坐标。
状态空间:状态矢量λ(t)所在的空间。
状态轨迹:在状态空间中,状态矢量端点随时间变化而推出的路径。
iL(t)与uC(t)与时间t的关系
3.用状态变量分析系统的优点:
1)便于研究系统内部一些物理量在信号转换过程中的变化,这些物理量可以用状态矢量的一个
分量表现出来,从而便于研究其变化规律。
2)系统的状态分析法与系统的复杂程度没有关系,复杂系统和简单系统的数学模型形式相似,都
可表示为一些状态变量的线性组合。
3)状态方程的主要参数鲜明的表征了系统的关键性能,利用状态方程分析系统的稳定性也比较
方便。
4)状态变量分析法适用于非线性或时变系统。
5)由于状态方程都是一阶微分方程或一阶差分方程,因而便于采用数值解法,为使用计算机分析
系统提供了有效途径。
二、连续时间系统状态方程的建立
1.状态方程的一般形式和建立方法概述
一个动态连续系统的时域数学模型可利用信号的各阶导数来描述。作为连续系统的状态方程表
现为状态变量的一阶微分联立方程组,即:
状态方程
d
dtλ1(t)=f1[λ1(t),λ2(t),…λk(t);e1(t),e2(t),…em(t),t]
d
dtλ2(t)=f2[λ1(t),λ2(t),…λk(t);e1(t),e2(t),…em(t),t]
…
d
dtλk(t)=fk[λ1(t),λ2(t),…λk(t);e1(t),e2(t),…em(t),t
]
—081—
输出方程
r1(t)=h1[λ1(t),λ2(t),…λk(t);e1(t),e2(t),Λem(t),t]
r2(t)=h2[λ1(t),λ2(t),…λk(t);e1(t),e2(t),Λem(t),t]
…
rr(t)=hr[λ1(t),λ2(t),…λk(t);e1(t),e2(t),Λem(t),t
]
对于线性时不变系统(LTI),状态方程式状态变量和输入信号的线性组合,可用矢量矩阵形式表
示为:
用矢量矩阵形式可表示为:
状态方程
[
d
dtλ]k×1
=Ak×kλk×1(t)+Bk×mem×1(t)
输出方程
[r(t)]r×1 =Cr×kλk×1(t)+Dr×mem×1(t)
其中:
λ(t)=
λ1(t)
λ2(t)
…
λk(t
)
A=
a11 a12 … a1k
a21 a22 … a2k
… … …
ak1 ak2 … a
kk
C=
c11 c12 … c1k
c21 c22 … c2k
cr1 cr2 … c
rk
d
dtλ(t[ ]) =
d
dtλ1(t)
d
dtλ2(t)
…
d
dtλk(t
)
B=
b11 b12 … b1m
b21 b22 … b2m
bk1 bk2 … b
km
D=
d11 d12 … d1m
d21 d22 … d2m
dr1 dr2 … d
rm
—181—
郑君里《信号与系统》考点精讲
r(t)=
r1(t)
r2(t)
…
rr(t
)
e(t)=
e(t
e(t
…
e(
t
2.由电路图直接建立状态方程
状态变量(均为独立变量):
电容两端电压与流经电感电流或电容电荷与电感磁链。
系统阶数(k):状态变量个数。
(图见视频)
注意选取独立变量
在选定状态变量之后,即可利用KCL和KVL列写电路方程,经化简消去一些不需要的变量,只留
下状态变量和输入信号经整理给出状态方程。
3.由系统的输入-输出方程或流图建立状态方程
假定某一物理系统可用如下微分方程表示
dk
dtk
r(t)+a1
dk-1
dtk-1
r(t)+… +ak-1
d
dtr(t)+ak
dk
dtk
r(t)
=b0
dk
dtk
e(t)+b1
dk-1
dtk-1
e(t)+… +bk-1
d
dte(t)+bke(t)
表示为算子形式为
(pk+a1p
k-1+Λ+ak-1p+ak)r(t)=(b0p
k+b1p
k-1+Λ+bk-1p+bk)e(t)
为便于选择状态变量,将其传输算子表示成
H(p)=
b0+b1/p+… +bk-1/p
k-1+bk/p
k
1+a1/p+… +ak-1/p
k-1+ak/p
k
这样当用积分器来实现该系统时,有流图如下:
(图见视频)
将系统的状态方程和输出方程表示成矩阵简化形式:
λ
·
(t)=Aλ(t)+Be(t)
r(t)=Cλ(t)+De(t{
)
对应的A、B、C、D矩阵分别是:
—281—
A=
0 1 0 …
0 0 1 …
0 0 0 … 1
-ak -ak-1 -ak-2 … -a
1
B=
0
0
0
1
C=[(bk-akb0),(bk-1-ak-1b0),…(b2-a2b0),(b1-a1b0)]
D=b0
三、连续时间系统状态方程的求解
1.用拉普拉斯变换法求解状态方程
若给定方程:
d
dtλ(t)=Aλ(t)+Be(t)
r(t)=Cλ(t)+De(t
{
)
两边取拉氏变换得
sΛ(s)-λ(0-)=AΛ(s)+BE(s)
R(s)=CΛ(s)+DE(s)
λ(0-)=
λ1(0-)
λ2(0-)
λk(0-
)
式中:λ(0-)为初始条件
整理得
Λ(s)=(sI-A)-1λ(0-)+(sI-A)
-1BE(s)
R(s)=C(sI-A)-1λ(0-)+[C(sI-A)
-1B+D]E(s{
)
因而时域表示式为:
λ(t)=L-1[(sI-A)-1λ(0-)]+L
-1[(sI-A)-1B]L-1E(s)
r(t)=CL-1[(sI-A)-1λ(0-)]+{CL
-1[(sI-A)-1B]+Dδ(t)}L-1E(s
{
)
2.用时域法求解状态方程(矢量微分方程求解)
对给定的状态方程,若已知:
d
dtλ(t)=Aλ(t)+Be(t)
并给定起始状态矢量λ(0-)=
λ1(0-)
λ2(0-)
λk(0-
)
—381—
郑君里《信号与系统》考点精讲
可解得方程的一般解为
λ=eAtλ(0-)+e
AtBe(t)
将此结果带入输出方程得到
r(t)=CeAtλ(0-)+[Ce
AtB+Dδ(t)]e(t)
无论状态方程的解或输出方程的解都由零输入解(由λ(0-)引起)和零状态解(由e(t)引起)相
加组成。两部分的变化规律都与矩阵eAt有关,因此eAt反映了系统状态变化的本质。
eAt称为“状态转移矩阵”,它的拉氏变化(sI-A)-1称为“特征矩阵”。
其它直接求eAt的方法:
1)计算机求解:将给定的A带入eAt=∑
!
k=0
1
k!A
ktk
利用计算机求无穷级数之和。
2)化对角阵计算法:若将A矩阵化为对角阵 A^,则eA^t也是对角阵,当 A^的对角各元素为ai时,e
A^t对
应元素为eait,即:
A^=diag[a1,a2,…,ak]
eA^t=diag[ea1t,ea2t,…,eakt]
3)化eAt为有限项之和求解,此处需要利用凯莱 -哈密顿定理。
3.由状态方程求系统函数H(s)
当给定系统的状态方程时,可利用已知的A,B,C,D矩阵表示系统转移函数H(s)。
设给定状态方程
λ(t)=Aλ(t)+Be(t)
r(t)=Cλ(t)+De(t{
)
(A,B,C,D为常数阵)
可得H(s)=R(s)E(s)=C(sI-A)
-1B+D
对其取逆变换的冲击响应:h(t)=CeAtB+Dδ(t)
此外转移函数的一般表达式为:H(s)=
b0s
k+b1s
k-1+… +bk-1s+bk
sk+a1s
k-1+… +ak-1s+ak
比较上面转移函数表达式可以看出|(sI-A)|即为H(s)分的特征多项式,所以称(sI-A)-1为
系统的特征矩阵。
四、离散时间系统状态方程的建立
1.状态方程的一般形式和建立方法概述
对于一个动态的时域离散系统,它的数学模型是用各阶差分方程形式描述的。作为离散系统的
状态方程表现为一阶差分联立方程组的形式。
即:
—481—
状态方程
λ1(n+1)=f1[λ1(n),λ2(n),…λk(n);x1(n),x2(n),…xm(n),m]
λ2(n+1)=f2[λ1(n),λ2(n),…λk(n);x1(n),x2(n),…xm(n),m]
…
λk(n+1)=fk[λ1(n),λ2(n),…λk(n);x1(n),x2(n),…xm(n),m
]
输出方程
y1(n)=h1[λ1(n),λ2(n),…λk(n);x1(n),x2(n),…xm(n),n]
y2(n)=h2[λ1(n),λ2(n),…λk(n);x1(n),x2(n),…xm(n),n]
…
yr(n)=hr[λ1(n),λ2(n),…λk(n);x1(n),x2(n),…xm(n),n
]
如果是线性时不变系统,则状态方程和输出方程是状态变量和输入信号的线性组合,即:
状态方程
λ1(n+1)=a11λ1(n)+a12λ2(n)+… +a1kλk(n)+b11x1(n)+b12x2(n)+… +b1mxm(n)
λ2(n+1)=a21λ1(n)+a22λ2(n)+… +a2kλk(n)+b21x1(n)+b22x2(n)+… +b2mxm(n)
…
λk(n+1)=ak1λ1(n)+ak2λ2(n)+… +akkλk(n)
+bk1x1(n)+bk2x2(n)+… +bkmxm(n
)
输出方程
y1(n)=c11λ1(n)+c12λ2(n)+… +c1kλk(n)+d11x1(n)+d12x2(n)+… +d1mxm(n)
y2(n)=c21λ1(n)+c22λ2(n)+… +c2kλk(n)
+d21x1(n)+d22x2(n)+… +d2mxm(n)
…
yr(n)=cr1λ1(n)+cr2λ2(n)+… +crkλk(n)+dr1x1(n)+dr2x2(n)+… +drmxm(n
)
2.由系统的输入 -输出差分方程建立状态方程
对于离散系统通常用下列k阶差分方程描述
y(n)+a1y(n-1)+a2y(n-2)+… +ak-1y[n-(k-1)]+aky(n-k)
=b0x(n)+b1x(n-1)+b2x(n-2)+… +bk-1x[n-(k-1)]+bkx(n-k)
表示成算子形式为
(Ek+a1E
k-1+a2E
k-2+… +ak-1E+ak)y(n)
=(b0E
k+b1E
k-1+b2E
k-2+…bk-1E+bk)x(n)
传输算子为
H(E)=
b0+b1/E+b2/E
2+… +bk-1/E
k-1+bk/E
k
1+a1/E+a2/E
2+… +ak-1E
k-1+ak/E
k
—581—
郑君里《信号与系统》考点精讲
(图见视频)
如图,选延时单元作为状态变量,则有:
λ1(n+1)=λ2(n)
λ2(n+1)=λ3(n)
λk-1(n+1)=λk(n)
λk(n+1)=-akλ1(n)-ak-1λ2(n)-…
-a2λk-1(n)-a1λk(n)+x(n
)
y(n)=(bk-akb0)λ1(n)+(bk-1-ak-1b0)λ2(n)+… +(b2-a2b0)λk-1(n)+(b1-a1b0)λk(n)
+b0x(n)
表示成矢量方程的形式为
λ(n+1)=Aλ(n)+Be(n)
y(t)=Cλ(n)+De(n{
)
其中
A=
0 1 0 …
0 0 1 …
0 0 0 … 1
-ak -ak-1 -ak-2 … -a
1
B=
0
0
0
1
C=[(bk-akb0),(bk-1-ak-1b0),…(b2-a2b0),(b1-a1b0)]
D=b0
可见根据离散系统的传输算子来列写系统的状态方程其步骤和结果与连续系统是一样的。
3.由给定系统的方框图或流图建立状态方程
给定系统的方框图或流图,只要取延时单元的输出作为状态变量便可建立系统的状态方程。
作为状态方程与输出方程,在形式上是相同的,就是由输入量、输出量、状态变量以及联系它们之
间关系A、B、C、D矩阵构成。
对连续系统状态方程和输出方程形式为
λ
·
(t)=Aλ(t)+Be(t)
r(t)=Cλ(t)+De(t{
)
对离散系统状态方程和输出方程形式为
λ(n+1)=Aλ(n)+Be(n)
y(t)=Cλ(n)+De(n{
)
—681—
五、离散时间系统状态方程的求解
离散系统状态方程的求解和连续系统的求解方法相似,包括时域和变化域。
1.矢量差分方程的时域求解
离散系统的状态方程为在n=n0有λ(n0),则有
λ(n0+1)=Aλ(n0)+Be(n0)
当n>n0利用迭代法可归结为:
λ(n)=An-n0λ(n0)+∑
n-1
i=n0
An-1-iBx(i)
若n0 =0,并将上述对n值得限制以阶跃信号的形式写入表达式
λ(n)=Anλ(0)u(n)+[∑
n-1
i=n0
An-1-iBx(i)]u(n-1)
解得输出为:
y(n)=Cλ(n)+Dx(n)=CAnλ(0)u(n)+[∑
n-1
i=n0
CAn-1-iBx(i)]u(n-1)+Dx(n)u(n)
若令x(n)=δ(n),则系统的单位脉冲响应为
h(n)=CAn-1Bu(n-1)+Dδ(n)
2、An的计算
与连续系统的情况类似,离散系统状态方程时域求解的关键步骤是计算状态转移矩阵,即An
利用凯莱 -哈密顿定理,有
Aj=c0I+c1A+c2A
2+… +ck-1A
k-1(jk)
分别用A的特征值代入上式,解联立方程即可求出系数
c0,c1,Λ,ck-1
3.离散系统状态方程的z变换解
由离散系统的状态方程和输出方程
λ(n+1)=Aλ(n)+Be(n)
y(t)=Cλ(n)+De(n{
)
两边取z变换并整理得
Λ(z)=(zI-A)-1zλ(0)+(zI-A)-1BX(z)
Y(z)=C(zI-A)-1zλ(0)+C(zI-A)-1BX(z)+DX(z{
)
取其逆变换即得时域表示式为
λ(n)=Z-1[(zI-A)-1z]λ(0)+Z-1[(zI-A)-1B]Z-1[X(z)]
y(n)=Z-1[C(zI-A)-1]zλ(0)+Z-1[C(zI-A)-1B+D]Z-1[X(z{
)]
状态转移矩阵为
An-1u(n-1)=Z-1[C(zI-A)-1]
—781—
郑君里《信号与系统》考点精讲
两边取z变换并整理得
Λ(z)=(zI-A)-1zλ(0)+(zI-A)-1BX(z)
Y(z)=C(zI-A)-1zλ(0)+C(zI-A)-1BX(z)+DX(z{
)
取其逆变换即得时域表示式为
λ(n)=Z-1[(zI-A)-1z]λ(0)+Z-1[(zI-A)-1B]Z-1[X(z)]
y(n)=Z-1[C(zI-A)-1]zλ(0)+Z-1[C(zI-A)-1B+D]Z-1[X(z{
)]
状态转移矩阵为
An-1u(n-1)=Z-1[C(zI-A)-1]
六、状态矢量的线性变换
1.在线性变换下状态方程的特性
设一组状态变量λ与另一组状态变量γ之间有
γ1 =p11λ1+p12λ2+… +p1kλk
γ2 =p21λ1+p22λ2+… +p2kλk
γk =pk1λ1+pk2λ2+… +pkkλ
k
表示为矢量式为γ
→
=P
→
λ
→
其中λ和γ为列矢量
γ=
γ1
γ2
γ
k
λ=
λ1
λ2
λ
k
P=
p11 p12 … p1k
p21 p22 … p2k
…
pk1 pk2 … p1
k
λ经线性变换变成γ,则原状态方程也做相应的改变。设源基底下状态方程表示为
d
dtλ(t)=Aλ+Be(t)
经γ
→
=P
→
λ
→
变换后
P-1 ddtγ(t)=AP
-1γ(t)+Be(t)
因而在新的状态变量下,状态方程与输出方程中的系数矩阵与原方程的系数矩阵有如下关系
A^=PAP-1
B^=PB
C^=CP-1
D^
=D
—881—
2.系统转移函数阵在线性变换下是不变的
状态方程和系统转移函数是描述系统的两种方法,当状态矢量用不同基底时,并不影响系统的物
理本质,因此对同一系统不同状态变量的选择,系统函数应是不变的。
证明:
H(s)=C^(sI-A)-1B^+D^=CP-1(sI-PAP-1)-1PB+D
=C[(sI-PAP-1)P]-1PB+D
=C[P-1(sI-PAP-1)P]-1B+D
=C(sI-A)-1B+D=H(s)
3.由状态方程判断系统的稳定性
特征值决定了系统的自由运动情况,因此可以根据A矩阵的特征值来判断系统的稳定情况。
1)连续系统
判断连续系统是否稳定,只要定出A的特征值Re[ai]是否小于零,即对稳定系统Re[ai]<0。
判断上述情况需要求解方程|aI-A|=0
方程|sI-A|=0的根在s平面上的位置决定了系统的稳定情况,因为其与求解方程|aI-A|=
0的特征根来判定系统稳定性是一致的。即只要知道根落在s平面的左半平面,便可确定系统的稳定情
况,即判定Re[ai]<0
2)离散系统
对于离散系统,根据A矩阵的特征值|ai|<1判定。即系统的特征根位于单位圆内系统稳定。和连
续系统相似,A矩阵的特征值和离散系统转移函数特征多项式的根位置相同,判定准则也相同。
七、系统的可控制性与可观测性
1)控制作用是否必然可以使系统在有限时间内,从起始状态指引到所要求状态,这就是可控制性
问题。
2)是否可以做到,通过观测有限时间内输出量,而识别系统的起始状态。
1.系统可控制性定义、可控阵满秩判别法
1)可控性:
当系统用状态方程描述时,给定系统的任意初始态,可以找到容许的输入量(即控制矢量),在有
限时间内把系统的所有状态引向状态空间的原点(即零状态),如果能做到这一点,则称系统是完全可
控制的。如果只对部分状态变量适用,则称系统不完全可控。
2)可控满秩判别法
设连续系统有k个状态变量,则可用如下状态方程
λ
·
k×1(t)=Ak×kλk×1(t)+Bk×ke(t)
它的解为
—981—
郑君里《信号与系统》考点精讲
λ(t)=eAtλ(0-)+∫
t
0
eA(t-τ)Be(τ)dτ
只要方程
eAtλ(0-)+∫
t
0
eA(t-τ)Be(τ)dτ=0有解即可满足要求。
将上式改写为:λ(0-)=-∫
t
0
e-AτBe(τ)dτ
考虑到凯莱 -哈密顿定理,则有:
e-Aτ =∑
k-1
i=0
ci(τ)A
i
代入上式有
λ(0-)=-∑
k-1
i=0
AiB∫
t1
0-
ci(τ)e(τ)dτ
令ωi(t1) =∫
t1
0-
ci(τ)e(τ)dτ
则上式可写成
λ(0-)=-∑
k-1
i=0
AiBωi(t1)
=-(B AB … Ak-1B)
ω0(t1)
ω1(t1)
…
ωk-1(t1
)
A为非奇异矩阵,若系统状态完全可控,也就是给定一组起始状态 λ(0-),满足之前方程,则必
须有:
M =[B AB A2B … Ak-1B]
满秩。这就是连续系统完全可控的充要条件。
M称为系统的可控性判别矩阵,简称可控阵,在给定系统状态方程时,只要M阵满秩,系统即为完
全可控系统
2.系统的可观性定义,可观阵满秩判别法
1)系统的可观测性:
如果系统用状态方程描述,在给定控制后,能在有限时间间隔内根据系统输出唯一地确定系统的
所有起始状态,则称系统完全可观;若只能确定部分起始状态,则称系统不完全可观。
2)可观阵满秩判别法
设连续时间系统有k个状态变量,则系统的输出可写作:
r(t)=CeAtλ(0-)+[Ce
AtB+Dδ(t)]e(t)
令e(t)=0,这样有:
—091—
r(t)=CeAtλ(0-)
或按凯莱 -哈密顿定理可写成
r(t)=C(c0I+c1A+Λ+ck-1A
k-1)λ(0-)
上式表明输出量是系统所有起始状态的线性组合。
r(t)=[c0 c1 Λ ck-1]
C
CA
M
CAk-
1
λ(0-)
因而要在0<t<t0时间间隔内,根据输出量唯一确定起
始状态必须使矩阵
N=
C
CA
…
CAk-
1
有K个线性无关列向量,亦即上阵满秩。这是连续系统可观的充要条件。N称为系统的可观阵,在
给定系统时,只要N阵满秩,则系统为完全可观系统。
—191—
郑君里《信号与系统》考点精讲
当前页面二维码