第5章LTI系统的频域分析
实验2连续LTI系统的频域分析

1.按所通过信号的频段分为低通、高通、带通和带阻滤波器四种。
低通滤波器:它允许信号中的低频或直流分量通过,抑制高频分量或干扰和噪声。
高通滤波器:它允许信号中的高频分量通过,抑制低频或直流分量。
带通滤波器:它允许一定频段的信号通过,抑制低于或高于该频段的信号、干扰和噪声。
带阻滤波器:它抑制一定频段内的信号,允许该频段以外的信号通过。
实验前,一定要针对下面的实验项目做好相应的实验准备工作,包括事先编写好相应的实验程序等事项。
给定三个连续时间LTI系统,它们的微分方程分别为
系统1: Eq.2.1
系统2: Eq.2.2
系统3:
Eq.2.3
1.修改程序Program2_1,并存盘,使之能够能够接受键盘方式输入的微分方程系数向量。并利用该程序计算并绘制由微分方程Eq.2.1、Eq.2.2和Eq.2.3描述的系统的幅度响应特性、相位响应特性、频率响应的实部和频率响应的虚部曲线图。
xlabel('Frequency in rad/sec')
input the coefficient vector of the right side of the differential equationb =
1 -1
the coefficient vector of the left side of the differential equationa =
262
the coefficient vector of the left side of the differential equationa =
phai = angle(H); % Compute the phase response phai
信号频域分析LTI系统频域分析

信号的频域分析LTI系统频域分析学院:班级:姓名:学号:指导教师:目录一.实验目的 (1)二.实验原理与方法 (1)三.实验内容 (4)四.实验结果 (5)五.实验总结 (8)一、 实验目的:1、 掌握信号的MATLAB 表示及其可视化方法.2、 掌握信号基本时域运算的MATLAB 实现方法.3、利用MATLAB 分析常用信号,加深对信号时域特性的理解二、实验原理与方法1、连续周期信号的频谱分析实验原理:如果信号满足狄里赫力条件,就可以展开为傅立叶级数形式,即其中T 0 表示基波周期,w 0 为基波频率。
上式定义为周期信号复指数形式的傅里叶级数,系数c k 称为x(t)的傅里叶系数。
同时,傅里叶级数还可以用三角函数的线性组合来表示,即展开为三角函数形式的傅里叶级数。
可见,任何满足狄里赫力条件的周期信号都可以表示成一组谐波关系的复指数函数或三角函数的叠加。
对于第一题中的周期矩形脉冲信号,其傅里叶级数展开式为:该题目中,当A=1,T=1,=0.5时X(t)=1/2 + 2/(n π)sin(n π/2)cos(2πnt);因此可以将信号的频谱化为各次谐波的幅度和相位来表示,根据该信号傅立叶展开式编写MATLAB 命令如下:(1) 该周期矩形信号的傅立叶级数表示为 :(2)利用MATLAB 绘出前N 次谐波波形,观察N 变化合成信号的变化规律 >> t=-1.5:0.01:1.5;>> N=7; //绘出前7次谐波合成的信号波形 ⎪⎪⎩⎪⎪⎨⎧==⎰∑-∞-∞=T dt e t x T c e c t x t jk k k t jk k 000)(1)(0ωω>> x=zeros(size(t));>> for n= 1:1:Nx=x+2/(n*pi)*sin(n*pi/2)*cos(2*pi*n*t); //计算x的傅里叶级数前N次谐波end>> plot(t,x+0.5); //绘出图形>> title(['N=7']); //设定图形名称>> xlabel('Time(sec)'); //添加坐标轴标注时间单位为秒(3)利用MATLAB汇出周期矩形脉冲信号的频谱,观察参数T和τ变化时对频谱波形的影响。
第五章1-连续LTI系统频域分析

连续时间LTI系统的频域分析 离散时间LTI系统的频域分析 信号的幅度调制和解调
时域分析的要点是,以冲激函数为基本信号,
任意输入信号可分解为一系列冲激函数;而系统零 状态响应yzs(t) = x(t)*h(t)。 由单位冲激函数δ (t)所引起的零状态响应称为单位 冲激响应,简称冲激响应,记为h(t)。
解: 利用H(j)与h(t)的关系
H ( j) F[h(t)] 1 1 j 1 j 2
1
( j)2 3( j) 2
只有当连续系统是稳定的LTI系统时,才存在H(j), 且可以由h(t)计算出H(j)。
电路系统的频率响应:
分析电路系统的频率响应,主要有两种方法。
H ( j) Yzs ( j)
( j) 3
X ( j) ( j)2 3( j) 2
在实际应用中, 只有当连续系统是稳定的LTI系统时,
才存在H(j),且频响函数才有意义。
例 已知某LTI系统的冲激响应为
h(t) = (e-t-e-2t) u(t),求系统的频率响应H(j)。
vR (t) RiR (t)
VR ( jw) R IR ( jw)
ZR
VR ( IR(
jw) jw)
R
vL
(t)
L
diL (t) dt
VL ( jw) jwLIL ( jw)
ZL
VL ( jw) IL ( jw)
jwL
iC
(t)
C
d
vC (t) dt
IC ( jw) jwCVC ( jw)
例 已知某LTI系统的动态方程为 y"(t) + 3y'(t) + 2y(t) = x(t),
理学LTI系统的频域分析PPT学习教案

本节小结
第12页/共13页
例1 例2
第2页/共13页
例3:当激励是ej 0t时,用频域分析法求系统 响应:
Y( j) H( j)2( 0)
H( j0)2( 0)
0
y (t)= H(j 0) ej
t
第3页/共13页
上例的两个推广
☆ 周期信号通过系统的响应
对于fT周(t期) 信号有F:n e jn t n
y(t) FnH ( jn) e jn t n
θ (ω)
H ( j) g2C () e j td
第7页/共13页
• 理想低通filter的冲激响应h(t)
d
c h(t)= ℱ-1[g 2
( ) e -j t
]=
t 1
h(t )
c π
sin c (t td (t td )
)
c
Sa[c
(t
td
)]
t
可见,它实际上是不可实现的非因果 系统。
td
t
π
c
第8页/共13页
• 理想低通filter的阶跃响应g(t)
g(t)=h(t)*(t)
g t
1
1
2
ππ
c c
O
t0Байду номын сангаас
t
tr
上升时间:
tr
第9页/共13页
理想低通滤波器在物理上是不可实现 的,近 似理想 低通滤 波器的 实例
• 一种可实现的低通滤波器
L
v1(t) C
H j
1
R v2(t)
系统无失真传输条件
第5页/共13页
无失真传输系统的频响图: H(j)=Ke-jωtd
LTI系统的频域分析

*
LTI hHale Waihona Puke t)傅里叶变换法=
变换 y(t) ③傅氏 Y (j ω ) 反变换
②傅氏 H(jω )
×
=
频率响应H(j)可定义为系统零状态响应的傅里叶变 换Y(j)与激励f(t)的傅里叶变换F(j)之比,即
Y ( j ) j[ y ( ) f ( )] Y ( j ) j ( ) H ( j ) e e H ( j ) F ( j ) F ( j )
|H(jω)| 1 ω - ωC 0 ωC θ (ω)
j t d , C e j t d H ( j ) g 2C ( ) e C 0,
(1)冲激响应
h(t ) F g 2c ( )e
1 jtd
c Sa[ ( t t )] c d
f (t )
FT FT
F ( )
1 2
相乘
1 Fs ( ) F ( ns ) Ts n
s
相卷积
T (t )
n
(t nT )
y(t ) K f (t td ) Y ( j ) Ke
其频谱关系为
jtd
F ( j)
(2)无失真传输条件:
系统要实现无失真传输,对系统h(t),H(j)的要求是: (a)对h(t)的要求: |H(jω)| K h(t)=K(t–td) (b)对H(j)的要求: ω H(j)=Y(j)/F(j)=Ke-jtd 0 即 θ (ω)
例:某LTI系统的H(j)和θ()如图,
若f(t)= 2 + 4cos(5t) + 4cos(10t),求系统的 响应。
θ(ω) π 1 -10 0 -π |H(jω)|
SignalsSystems_Chapter5

信号与系统天津大学电子信息工程学院第五章连续系统的复频域分析一、拉普拉斯变换(LT)(一)从傅里叶变换到拉普拉斯变换z1、从FT到双边LT信号f(t)的傅里叶变换(FT)为z许多函数不满足绝对可积条件,其F( jω)中一般都含有冲激函数。
用衰减因子e-σt乘以f(t),适当选择σ的值,使f(t)·e-σt绝对可积,从而可求得其FT:如果令s=σ+jω——称为f(t)的双边LT3z根据FT-1反变换式,可得:——F(s)的双边拉普拉斯反变换z F(s)称为f(t)的象函数,f(t)称为F(s)的原函数。
z记作:F(s)=_{f(t) },f(t)=_-1{F(s) },或者简52、收敛域(ROC)使双边LT 的象函数F b (s )存在的s 平面的区域称为双边LT 的收敛域z (1)因果信号z (2)反因果信号z(3)双边函数73、单边拉普拉斯变换单边拉普拉斯变换单边拉氏逆变换4、单边LT的收敛域——F(s)存在的充分条件对于双边LT,必须认真研究收敛域问题,须由F(s)和收敛域共同确定原函数f(t)9LT的收敛域分为以下三种情况:z①收敛域是整个s平面根据收敛条件:推广:凡时宽有限且幅度有限的信号(满足绝对可11②F (s )在s 平面的部分区域收敛z 一般而言,单边LT 的收敛域是在s 平面上σ>σ0的区域。
z 收敛域的横坐标σ0(=α)称为收敛坐标,直线σ=σ0称为收敛轴。
③在整个平面上,F (s )都不收敛,即F (s )不存在z 如:、t t 等函数,其随t 上升而增加的速度超过指数阶函数,F (s )不存在。
2t e 如因果信号f (t )满足:(1)在有限区间a <t <b ()内可积,(2)对于某个有,则对于,拉普拉斯积分式绝对且一致收敛。
(教材P214定理)0a b ≤<<∞0σ0lim |()|0,t t f t e σσσ−→∞=>0Re()s σσ=>(二)常用函数的单边LT变换z1、复指数函数13可推出一些函数的LT:15z2、f(t)=t n·ε(t),n为正整数17 3、冲激函数δ(t)与冲激偶δ’(t)二、Laplace变换的性质z1、线性性质19 2、尺度变换(比例性)注:a < 0不适用于单边LT213、时移(延时)特性说明:①注意f (t -t 0)·ε(t -t 0) 与f (t -t 0)·ε(t )的区别z②注意延时性与比例性综合应用的情况例123有始周期函数的拉氏变换等于其第一周期的拉氏变换-Ts25z例2(教材P219例5.2-3)试求在t =0-时接入的周期性冲激序列的象函数。
LTI系统的频域分析

y(t ) h(t )* fT (t ) Fn [h(t )*e jnt ] Fn H ( jn) e jnt n n 若
则可推导出
A0 y(t ) H (0) An | H ( jn) | cos[nt n (n)] 2 n 1
h( ) e j d
y(t ) H ( j) e
j t
H ( j )反映了响应y(t)的幅度和相位。
二、一般信号f(t)作用于LTI系统的响应
e
1 2
j t
H(j ) ej t
1 2
齐次性
1 j t F ( j ) e d 2
F(j )H(j ) ej t d
FT [TS (t )] S
n
( n
S
)
如果f(t)是带限信号[即f(t)的频谱在(- m,m) 为有限值,而其余区间为0]。
设f(t)←→F(j),取样信号fS(t)的频谱函数
1 FS ( j ) F ( j )* S ( nS ) 2 n 1 F[ j ( n S )] TS n
f (t )
FT FT
F ( )
1 2
相乘
1 Fs ( ) F ( ns ) Ts n
s
相卷积
T (t )
n
(t nT )
FT
FT
p( ) s
n Βιβλιοθήκη ( ns
)
在画取样信号fS(t)的频谱时,设定ωS ≥2ωm ,这时 其频谱不发生混叠,因此能设法(如利用低通滤波器), 从FS(j)中取出F(j),即从fS(t)中恢复原信号f(t)。 否则将发生混叠,而无法恢复原信号。
自动控制原理第5章

jY (ω )
ω =∞
X (ω )
ω
积分环节的Nyquist图 积分环节的Bode图
幅频特性与角频率ω成反比,相频特性恒为-90° 成反比, 90° 对数幅频特性为一条斜率为 - 20dB/dec的直线,此 线通过L(ω)=0,ω=1的点
三、微分环节 微分环节的频率特性为
G ( jω ) = jω = ωe
奈奎斯特(N.Nyquist)在1932年基于极坐标图阐述 奈奎斯特(N.Nyquist)在1932年基于极坐标图阐述 了反馈系统稳定性。 极坐标图(Polar 极坐标图(Polar plot) =幅相频率特性曲线=幅相曲线 幅相频率特性曲线=
G ( jω )
可用幅值 G( jω ) 和相角ϕ (ω ) 的向量表示。
当输入信号的频率 ω → 0 ~ ∞ 变化时,向量 G ( jω ) 的幅值和相位也随之作相应的变化,其端点在复平面 上移动的轨迹称为极坐标图。
jY (ω )
ω →∞
ϕ (ω ) A(ω )
ω = 0 X (ω )
ω
RC网络对数频率特性 RC网络频率特性
5.2 典型环节的频率特性
用频域分析法研究控制系统的稳定性和动态 响应时,是根据系统的开环频率特性进行的, 响应时,是根据系统的开环频率特性进行的, 而控制系统的开环频率特性通常是由若干典 型环节的频率特性组成的。 型环节的频率特性组成的。 本节介绍八种常用的典型环节。 本节介绍八种常用的典型环节。
频率响应: 正弦输入信号作用下, 系统输出的稳态分量。 频率响应 : 正弦输入信号作用下,系统输出的稳态分量。 (控制系统中的信号可以表示为不同频率正弦信号的合成) 控制系统中的信号可以表示为不同频率正弦信号的合成) 频率特性: 系统频率响应和正弦输入信号之间的关系, 频率特性 : 系统频率响应和正弦输入信号之间的关系,它 和传递函数一样表示了系统或环节的动态特性。 和传递函数一样表示了系统或环节的动态特性。 数学基础:控制系统的频率特性反映正弦输入下系统响应 数学基础:控制系统的频率特性反映正弦输入下系统响应 的性能。研究其的数学基础是Fourier变换。 的性能。研究其的数学基础是Fourier变换。 频域分析法:应用频率特性研究线性系统的经典方法。 频域分析法:应用频率特性研究线性系统的经典方法。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
A
e jn
1 aej
1 ae j
- 当n趋向无穷大时含有 的前亮相趋向于零 - 前两项为暂态响应,后一项为稳态响应 - 在实际应用中,系统的暂稳态响应是不重要的,
因此,暂稳态响应被忽略 13
5.1.3 周期输入信号的稳态响应
➢基本周期为N的周期信号x(n)输入到稳定的线性
时不变系统时,系统在任意时刻n的响应是稳态响 应
其中, Vk () e j zk
,
k () e j zk
Uk () e j pk
k () e j pk
23
M
e j zk
M
Vk ()e jk ()
H () b0e j(NM )
k 1 N
b0e j ( N M )
k 1 N
e j pk
U k ()e jk ()
➢根据频域特性,滤波器分为低通、高通、带通、
带阻、全通滤波器
30
31
z平面极点、零点位置和数字滤波器设计
➢理想滤波器的特性:
- 常数增益(通常是单位增益)的带通特性,而 在带阻部分的增益为零
- 线性相位特性
➢理想滤波器在物理上是不可实现!
➢ 例:
hlp
(n)
sin cn n
,
n
不是因果,也不是绝对可和 → 不稳定 32
-输入信号的叠加性:xn
1 2
x1 n
x2
n
A
cosn
10
-叠加信号的输出响应:
yn
1 2
y1n
y2 n
A Hcosn
➢ 如果输入信号是
xn
1 j2
x1
n
x
2
n
A sin n
- 输出响应:
yn
1 j2
y1
n
y
2
n
yn A H sinn
11
➢ 正弦输入信号的输出响应
L- 输入信号:xn Ai osin i i1➢线
低通滤波器的零-极点模型
高通滤波器的零-极点模型
34
单极点低通滤波器
➢系统函数
1 a H1(z) 1 az1
1 a 1 z 1 H 2 (z) 2 1 az1
幅度响应
相位响应
35
单极点低通滤波器
➢通过将低通滤波器的
零极点位置在z平面关于 虚轴进行反转(折叠), 就可以得到简单的 高通滤波器
z平面极点、零点位置和数字滤波器设计
➢放置极、零点的基本原则:
- 在单位圆上对应于需要加强频率的点附近放 置极点
- 在需要拉低的频率点处附近放置零点
➢约束:
1:所有的极点必须放置在单位圆内,零点可以 放在任何位置
2:所有复值的极、零点必须以共轭的形式出现 33
5.4.2 低通、高通和带通滤波器
➢ 如果知道了 H 和在 0 上的值,
就可得到这两个函数在 0 上的值。
9
➢ 正弦函数的复共轭指数函数
-输入信号:x1n A e jn
x 2 n A e jn
-输出响应 y1n A H e je jn
y2 n A H e je jn A H e je jn
➢ 线性系统系统可能会改变周期输入信号的形状
(如:收缩、放大、相位移动) 但,不能改变输入信号的周期
14
5.1.4 非周期输入信号的响应
➢ 根据卷积定理来具体计算 ➢ 如xn表示输入序列,yn 表示输出序列,
hn 表示系统的单位采样响应
• Y() ,X(),H() 分别是 yn ,xn,hn 对应
输入-输出关系的卷积公式:
yn hn xn hkxn k k
- 单位采样响应: hn, n
➢ 激励系统的复指数信号表示:
xn Aejn , n
A是幅度, 是限制在频率区间 , 上的任意频率
3
➢ 代入之后可以得到输出响应:
yn h k A e jnk k
A
k
n
- 输出响应:
L
xn Ai Hi cosin i i , i1 12
5.1.2 正弦输入信号的稳态和暂稳态响应
➢线性时不变系统对作用于 n 的指数和
正弦输入信号的响应时,输出响应是稳态响应, 而没有暂态响应。为什么?
➢例:输入是复指数信号的一阶差分方程
y n a n1y 1 Aa n1e jn1 e jn
45
➢系统函数的描述: Hz
b0
1 a1z 1 a 2 z 2
➢差分方程:yn a1yn 1 yn 2 b0n
➢定义:对所有频率具有常数幅度响应的系统
H 1, 0
➢应用范围:延迟系统
➢特征:原始信号不改变,只延迟k单位
Hz z k
➢全通滤波器系统函数的一般形式:
N
A(z) ak zk , a0 1 k 0
H (z) zN A(z1)
A( z )
43
➢全通滤波器的极点和零点互为倒数 ➢系统的相位响应达不要求时此滤波器来补偿
的傅立叶变换
• Y () H ()X () • 输出信号的频谱等于输入信号的频谱乘以系统
的频率响应 15
➢ 如果 Y() , X (),H() 表示为点斜式 ➢ 输出信号的幅度:Y() H() X () ➢ 输出信号的相位:Y () H () X ()
16
➢ 输出信号的能量密度谱: ➢ Y () H ()X () 等式两边取平方
M
H L (z) h(k)zkL k 0
➢梳状滤波器的频率响应:
M
重复性! H L () h(k)e jkL H (L) k 0 40
➢HL(ω)的频率特性仅仅是在区间[0, 2π]内的L
阶重得
➢例:L=5
41
➢例:M=3, L=3的FIR梳状滤波器的结构图
42
5.4.6 全通滤波器
H
3
(
z
)
1
2
a
1 z 1 1 az1
幅度响应
相位响应
36
5.4.3 数字谐振器
➢可利用两极点带通滤波
器的两个极点以复共轭对 的形式来表示。如:语音信号
幅度响应
相位响应
37
5.4.4 槽口滤波器
➢包含n个深槽口的滤波器,在特定点的频率响应
为零(如:0 ,1)
➢单位圆的角 0 处引入复共轭的零点方式来产生 ➢应用范围:语音信号处理领域(共振)
7
➢ 其中实部和虚部定义为:
HR hkcosk k
HI hksin k k
➢ H 的幅度: H
H
2 R
H
2 I
➢ H
的相位:
arctan
H I H R
8
➢ 收敛范围
因为HR 是 的偶函数,HR HR 因为HI是 的奇函数, HI HI
所以,H 是 的偶函数,是 的奇函数
槽口滤波器的频率响应特性
38
➢在零点附近引入极点
的方式来调整槽口宽带
➢极点产生实际的共振 ➢多次调整的方式来
幅度响应
产生共振
➢优点:结构简单 ➢缺点:调试难
相位响应
39
5.4.5 梳状滤波器
➢FIR滤波器的系统函数:
M
H (z) h(k)zk k 0
➢用zL来代替z,就可以构造梳状滤波器
➢极点 pk 和零点 zk 位于z平面上的A和B处,计算
L点的傅里叶变换
CL CA AL CL CB BL
(a)
25
➢然而 CL e j , CA pk , CB zk ➢因此 AC e j pk
BL e j zk
AL e j pk Uk e jk BL e j zk Vk e jk
k 1 N
1 pke j
k 1
M
e j zk
或等介于
H () b0e j(N M )
k 1 N
e j pk
k 1
22
M
e j zk
H () b0e j(N M )
k 1 N
e j pk
k 1
➢ 将上式中的各复值因子写成点斜式
e j zk Vk ()e jk () e j pk Uk ()e jk ()
1 2
H
2
Sxx
d
28
5.4 作为频率选择滤波器的线性时不变系统
➢滤波器:用来描述一个设备,根据作用于输入
端的对象的某些属性进行分辨过滤,以让某部分 通过它
➢频率选择滤波器:在某些频段的频率分量的信
号可以通过,而包含在其他频段的频率分量中的 信号将被衰减
➢线性时不变系统和滤波器是同义的
29
5.4.1 理想滤波特性
➢若H(z)是有理函数形式,且 H(z)=B(z)/A(z)
M
M
H ()
B() A()
bk e jk
k 0
N
1 ak e jk
b0
1 zk e j
k 1
N
1 pk e j
k 1
k 1
18
M 1 zk e j
H () b0
k 1 N
1 pk e j
k 1
M 1 zke j
H () b0
k 1 N
1 pke j
k 1
➢ H*(ω)是H* (1/z*)在单位圆上的值
M
1
z
k
z
H (1/ z) b0
k 1 N
1 pk z
k 1
19
➢当{h(n)}是实数,或等价的系数{ak}和{bk}是
实数时 H (1/ z) H (z1), H () H ()