离散时间系统的时域分析
实验五 离散时间系统的时域分析

实验五 离散时间系统的时域分析一、实验目的:(1)理解离散时间信号的系统及其特性。
(2)对简朴的离散时间系统进行分析,研究其时域特性。
(3)运用MATL AB对离散时间系统进行仿真,观测成果,理解其时域特性。
二、具体实验:1、离散时间系统的仿真——滑动平均系统s1s2xFigur e 5-1 T he wave form of s1,s2,x由图5-1所示及其运算可知,s1=cos(2*pi*0.05*n),s 2=cos(2*pi*0.47*n ),s1周期T1=1/0.05=20,s 2周期T2=1/0.47=100/47。
x=s 1+s 2,x 的周期为T1、T2的最小公倍数,因此x的周期为100。
Time Serial n A m p l i t u d eSignal #1Time Serial n A m p l i t u d eSignal #2Time Serial nA m p l i t u d eInput SignalTime Serial nA m p l i t u d eOutput SignalTime Serial n A m p l i t u d eSignal #1Time Serial n A m p l i t u d eSignal #2Time Serial nA m p l i t u d eInput SignalTime Serial nA m p l i t u d eOutput SignalF ig ure 5-2 Fi gure 5-3(1)如图5-2,当M=2时,第一种图显示的是一种低频信号,第二个是高频信号,第三个图是信号一和信号二的合成的输入,第四个是通过函数Y 的得出的输出。
成果是低频信号,前后对比得出是高频信号被克制了。
本系统是滑动平均滤波器,为低通滤波系统,功能就是从信号中滤除高频分量,因此输入的高频分量s2[n]被该系统克制了。
离散时间信号的时域分析实验报告

离散时间信号的时域分析实验报告实验报告:离散时间信号的时域分析一、实验目的本实验旨在通过MATLAB软件,对离散时间信号进行时域分析,包括信号的显示、基本运算(如加法、减法、乘法、反转等)、以及频域变换(如傅里叶变换)等,以加深对离散时间信号处理的基本概念和原理的理解。
二、实验原理离散时间信号是在时间轴上离散分布的信号,其数学表示为离散时间函数。
与连续时间信号不同,离散时间信号只能在特定的时间点取值。
离散时间信号的时域分析是研究信号的基本属性,包括幅度、时间、频率等。
通过时域分析,我们可以对信号进行各种基本运算和变换,以提取有用的信息。
三、实验步骤1.信号生成:首先,我们使用MATLAB生成两组简单的离散时间信号,一组为正弦波,另一组为方波。
我们将这些信号存储在数组中,以便后续分析和显示。
2.信号显示:利用MATLAB的绘图功能,将生成的信号在时域中显示出来。
这样,我们可以直观地观察信号的基本属性,包括幅度和时间关系。
3.基本运算:对生成的信号进行基本运算,包括加法、减法、乘法、反转等。
将这些运算的结果存储在新的数组中,并绘制出运算后的信号波形。
4.傅里叶变换:使用MATLAB的FFT(快速傅里叶变换)函数,将信号从时域变换到频域。
我们可以得到信号的频谱,进而分析信号的频率属性。
5.结果分析:对上述步骤得到的结果进行分析,包括比较基本运算前后的信号波形变化,以及傅里叶变换前后的频谱差异等。
四、实验结果1.信号显示:通过绘制图形,我们观察到正弦波和方波在时域中的波形特点。
正弦波呈现周期性的波形,方波则呈现明显的阶跃特性。
2.基本运算:通过对比基本运算前后的信号波形图,我们可以观察到信号经过加法、减法、乘法、反转等运算后,其波形发生相应的变化。
例如,两个信号相加后,其幅度和时间与原信号不同。
反转信号则使得波形在时间轴上反向。
3.傅里叶变换:通过FFT变换,我们将时域中的正弦波和方波转换到频域。
正弦波的频谱显示其频率为单一的直流分量,方波的频谱则显示其主要频率分量是直流分量和若干奇数倍的谐波分量。
信号与系统 第8章 离散时间系统的时域与变换域分析

4
8.1.1 线性时不变离散时间系统
例8.1-1 设某离散系统激励x[n]与响应y[n]之间的关系为
y[n] = nx[n],判断该系统是否为线性时不变系统。
1 M2
1
M2
k M1
(x1[n
k]
x1[n
k ])
1
M2
1
M2
M1 M 2 1 kM1 x1[n k] M1 M 2 1 kM1 x2[n k]
y1[n] y2[n]
该系统满足叠加性,所以该系统是线性系统。
(3)假设输入信号为x[n]= x1[n-m],则输出信号为
y[n]
y[n] = x[n] + ay[n-1] = a n
此范围仅限于n ≥ 0,
故
y[n] = anu[n] 12
8.2 常系数线性差分方程的求解
N
M
ak y[n k] br x[n r]
k 0
r0
(8.2-2)
8.2.1 线性常系数差分方程的时域经典法求解
一般地,常系数线性差分方程的解由齐次解和特解组成。
的完全解。
其中激励信号为x[n] n2,且边界条件为 y[1] 1
解:(1)齐次解为 yh[n] C(2)n
(2)将 x[n] n代2 入差分方程的右端,得自由项为
2n 1
从而特解为 yp[n] D1n D2
其中,D1和D2为待定系数,代入原方程得
3D1n 3D2 2D1 2n 1
离散时间系统的时域分析

称为混叠。 常称作折叠频率。 2
信号频率
fa nfs fm
fa fs / 2
假频
Fδ(jω)
抽样频率
ω Ω-ωm ωm Ω
例如:当抽样率为5kHz对3kHz的余弦信号 抽样,然后用截止频率为2.5kHz的低通滤波 器进行滤波,输出的频谱只包含2kHz的频率, 这是原信号中所没有的。
对一个低通滤波器的冲激响应进行抽样,抽 样后低频通带将在整个频率轴上周期的重复出现, 这种现象称为“伪门”。在设计数字滤波器时要 适当选择抽样率,使得伪门在干扰频率之外。
H(jω)
ω 0 数字滤波器的伪门
例1:对于频率为150Hz的正弦时间序列,分别以4ms 和8ms采样结果会如何?
100HZ 25HZ
在实际工作中应用抽样定理时,还应考虑下 面两个实际问题:
1、在理论上讲,按照奈奎斯特抽样率抽样, 通过理想低通滤波器以后,就可以恢复原信 号。但理想低通滤波器在物理上是不可实现 的,实际滤波器都存在一个过渡带,为了保 证在滤波器过渡带的频率范围内信号的频谱 为零,必须选择高于2fm的抽样率。
u (n) 0, n 0
...
n -1 0 1 2 3
(n) u(n) u(n) u(n 1)
u(n) (n m) (n) (n 1) (n 2) m0
3.矩形序列 R N (n )
1, R N (n) 0,
0 n N 1 其他n
RN (n) u(n) u(n N )
第五章 离散时间系统 的时域分析
§5.1 离散信号与抽样定理
一、离散信号及其表示
1、离散时间信号是指只在一系列离散的时刻 tk (k = 0,1,2,…)时,信号才有确定值,在其它时 刻,未定义; 2、离散时间信号是离散时间变量 tk 的函数; 3、抽样间隔可以是均匀的,也可以非均匀。
3.2.3离散时间LTI系统的时域分析 - 离散时间LTI系统的时域分析(精品文档)

ci可由初始状态 yzi (1),yzi (2), ,yzi (k) 确定
10
信号处理与系统
DLTI系统零输入响应通解
y(n) yzi (n) yzs (n)
故有: yzi (1) y(1), , yzi (k) y(k)
n0
yzi (1) y(1), ,yzi (k) y(k), n 0
y(1) y(2)
c1 c1
(3)1 c2 (3)2 c2
0 1/
2
cc12
3/4 9 / 4
yzi
(n)
3 4
9 4
(3)n ,
n0
12
信号处理与系统
DLTI系统零输入响应分析
DLTI系统零输入响应通解为:
yzi (n) c1(1 )n ckr (kr )n ckr1nr1(0 )n ck1n(0 )n ck (0 )n
其中 1 2 kr ,即k-r个单根,0为r个重根
(i )n ,i 1, , k r
例2. 一信号处理过程是:每当收到一个数据,就将此 数据与前一步的处理结果平均。求这一信号处理 过程的输入输出关系。
解:
y(n) 1 y(n 1) 1 x(n)
2
2
一阶后向差分方程系统模拟框图
4
信号处理与系统
一、DLTI系统方程的建立
离散时间线性时不变 (discrete-time ,linear, time-invariant, 记作DLTI) 系统:用常系数差分方程来描述
-
试从微分方程推导其差分方程。
解: d y(t) 1 y(t) 1 x(t)
6.离散时间信号与系统的时域分析

0, n 1 1 z ( n) x ( n) y ( n) , n 1 2 1 n 1 ( 2 )( n 1)( 2 ) , n 0
6 线性时不变离散系统的时域分析
5. 累加 设某一序列为x(n),则x(n)的累加序列 y(n)定义为
y ( n)
k
x(k ) x(n) * u(n)
n
根据上述性质可以推得以下结论:
f (n n1 ) * (n n2 ) f (n n1 n2 )
6 线性时不变离散系统的时域分析
例 已知 x1 (n) (n) 3 (n 1) 2 (n 2) x2 (n) u(n) u(n 3) 试求信号 x (n) ,它满足 x(n) x1 (n) x2 (n) 解:可利用上面讲述的性质求解。
1 1/ 2 1/4 -2 -1 0 1 1/8 ... 2
n
x(-n) 1 1/2 1/8 1/4 ... -2 -1 0
1
2
n
6 线性时不变离散系统的时域分析
3.序列的加减 两序列的加、减是指同序号(n)的序列值逐项对 应相加得一新序列。
6 线性时不变离散系统的时域分析
例:
x(n) 1 1/2 1/4 -2 -1 0 y(n) 2 1 1/4 1/2 1 2 n …
6 线性时不变离散系统的时域分析
2.单位阶跃序列
u(n)
1, u ( n) 0,
n0 n0
u(n)
...
-1 0 1 2 3 n
(n) u (n) u (n) u (n 1)
m 0
u (n) (n m) (n) (n 1) (n 2)
离散时间系统的时域分析

第七章离散时间系统的时域分析§7-1 概述一、离散时间信号与离散时间系统离散时间信号:只在某些离散的时间点上有值的信号。
离散时间系统:处理离散时间信号的系统。
混合时间系统:既处理离散时间信号,又处理连续时间信号的系统。
二、连续信号与离散信号连续信号可以转换成离散信号,从而可以用离散时间系统(或数字信号处理系统)进行处理:三、离散信号的表示方法:1、 时间函数:f(k)<——f(kT),其中k 为序号,相当于时间。
例如:)1.0sin()(k k f =2、 (有序)数列:将离散信号的数值按顺序排列起来。
例如:f(k)={1,0.5,0.25,0.125,……,}时间函数可以表达任意长(可能是无限长)的离散信号,可以表达单边或双边信号,但是在很多情况下难于得到;数列的方法表示比较简单,直观,但是只能表示有始、有限长度的信号。
四、典型的离散时间信号1、 单位样值函数:⎩⎨⎧==其它001)(k k δ 下图表示了)(n k −δ的波形。
这个函数与连续时间信号中的冲激函数)(t δ相似,也有着与其相似的性质。
例如:)()0()()(k f k k f δδ=,)()()()(000k k k f k k k f −=−δδ。
2、 单位阶跃函数:⎩⎨⎧≥=其它001)(k k ε这个函数与连续时间信号中的阶跃函数)(t ε相似。
用它可以产生(或表示)单边信号(这里称为单边序列)。
3、 单边指数序列:)(k a k ε比较:单边连续指数信号:)()()(t e t e t a at εε=,其底一定大于零,不会出现负数。
(a) 0.9a = (d) 0.9a =−(b) 1a = (e) 1a =−(c) 1.1a = (f) 1.1a =−4、 单边正弦序列:)()cos(0k k A εφω+双边正弦序列:)cos(0φω+k A五、离散信号的运算1、 加法:)()()(21k f k f k f +=<—相同的k 对应的数相加。
离散时间系统的时域分析

x(2n)
6
4
2
O 123456 n
已知x(n)波形,请画出 x(2n), x n 波形。
2
x n 6 2 5 4 3 2 1 O 1 2 3 4 5 6 7 8 9 10 12 n
三.常用离散信号
•单位样值信号 •单位阶跃序列 •矩形序列 •斜变序列 •单边指数序列 •正弦序列 •复指数序列
z(n) x(n + m) 左移位
x(n)
x(0) x( 1) x(1) x(3)
2 1 o 1 3 n
x(n 1) x(0)
x( 1) x(1) x(3)
3 1 o 1 2 4 n
x(2)
x(2)
5.倒置: z(n) x(n) 6.差分:前向差分:x(n) x(n + 1) x(n)
1.迭代法
2.时域经典法:齐次解+特解
3.零输入响应+零状态响应 利用卷积求系统的零状态响应
4. z变换法反变换y(n)
解差分方程的基础方法
一.迭代法 差分方程本身是一种递推关系,
(n)在n 0取有限值(不是面积)。
利用单位样值信号表示任意序列
x(n) x(m) (n m) m f (n)
1.5 2
1 o 1
34 n
3
f (n) 1,1.5,0,3,0,0, (n + 1) + 1.5 (n) 3 (n 2)
a nu(n)
1 a 0
1 1 O
1 2 3 4n
1 1 O 1
23
4n
x(n) sin(nω0 )
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
(f) a = −1.1
比较:单边连续指数信号: eatε (t) = (ea )tε (t) ,其 底一定大于零,不会出现负数。
4、 单边正弦序列: Acos(ω0k + φ )ε (k)
双边正弦序列: Acos(ω0k + φ )
五、离散信号的运算 1、 加法: f (k) = f1(k) + f2 (k) <—相同的 k 对应
似,只不过将其中的积分器变成延时(移序)器。
离散时间系统的初始状态可以包含在延时
(移序)器中。
e(k) y(k +1) y(k)
∑
D
−a
一阶离散时间系统的模拟框图
n 阶离散时间系统模拟框图
§7-4 离散时间系统的零输入响应
离散差分方程的解法: 1) 时域经典法
与微分方程一样,将解分为通解(齐次解) 和特解两部分。首先确定形式解,再代入初始 条件(或边界条件),确定其中的待定系数。 优点:物理概念清晰,可以一次得到全部解; 缺点:特解有时很难求,不实用。 2) 近代时域法:
(b) 原信号的频谱 F ( jω)
(c)单位冲激序列 δT (t)
(d)单位冲激序列的频谱ωsδωs (ω )
(ωs
=
2π T
)
(e)
fδ (t)
=
1 τ
fs (t)
=
f (t)δT (t)
(f) fδ (t) 的频谱
如果原来信号最大频率分量为的谱 ωm ,抽样 频率 ωs > 2ωm ,则周期化后的各个频谱不会相互 重叠。将抽样信号通过一个截止频率为 ωs / 2 、增 益为 T 的 ILPF,可以不失真地还原原来的信号。 此低通滤波器的冲激响应:
将解分为零输入响应 rzi (k ) 和零状态响应 rzs (k ) 两部分。对零输入响应 rzi (k ) 仍然用时域 经典法;零状态响应 rzs (k ) 用卷积和求解。
这种方法是求解差分方程的主要方法; 3) 变换域解法:Z 变换( Z.T.),相当于连续时
间系统中的 L.T.变换法。在第八章中介绍。 4) 数值解法:利用前向预测形式的差分方程,
例 1:人口(或虫口)问题:
z 假设人口的年出身率为 a,则 k 年人口 y(k)和
下一年的人口 y(k+1)之间的关系为:
y(k + 1) = (1+ a) y(k) <—前向(预测)方程;
或:
y(k)
=
1 (1 +
a)
y(k
+
1)
<—后向(滤波)方程;
或: y(k +1) − (1+ a) y(k ) = 0 <—一般差分方程。
ωs > 2ωm 时,用信号的一些离散的时间点上的数 值来代替这个信号可以不损失任何信息。
能够完全不失真地还原信号所需要的最小的
抽 样 频 率 ωs = 2ωm 称 为 Nyquist 抽 样 频 率 ,或 Shannon 抽样频率。
z 在实际工程中的做法与取样中的过程正好相 反:首先测量得到 f(kT),然后再构成抽样信号。 工程上的采样就是指测量到 kT 时刻 f(t)的值。
z 在构成抽样信号时,不可能产生冲激信号,这 时候可以用任意的周期性脉冲信号代替,其结 果不变。
z 恢复信号时,ILPF 是不可能实现的,只能用其 它的 LPF,所以抽样频率必须进一步增加,一 般取 ωm 的 3~5 倍。
抽样信号经过非理想低通滤波器
z 如果原来的信号是一个带限信号,则 Nyquist 抽样定理还可以做适当修改。
四、典型的离散时间信号
1、
单位样值函数:
δ
(k
)
=
⎧1 ⎩⎨0
k =0 其它
下图表示了δ (k − n) 的波形。
这个函数与连续时间信号中的冲激函数 δ (t) 相似,也有着与其相似的性质。例如:
f (k)δ (k ) = f (0)δ (k ) , f (k )δ (k − k0 ) = f (k0 )δ (k − k0 ) 。
3、 线性移不变离散时间系统 同时满足线性和移不变性的系统。
七、离散时间系统的描述方法:见§7-3。
§7-2 抽样信号与抽样定理
离散信号可以通过对连续信号抽样得到;连 续信号可以通过抽样转化为离散信号,从而可以 用离散时间系统进行处理。但是,这牵涉到两个 问题:
1)怎样进行抽样? 2)如 何 抽 样 才 能 不 损 失 原 来 信 号 中 的 信
本章重点介绍近代时域法。 首先,在本节中介绍近代时域法中零输入
月必然都长成大兔子 所以,第 k+2 月兔子的总对数为: y(k+2)=y(k)+y(k+1) 或者:y(k+2)- y(k+1)-y(k)=0
差分方程的一般形式:
r(k + n) + an−1r(k + n −1) + ... + a1r(k + 1) + a0r(k ) = bme(k + m) + bm−1e(k + m −1) + ... + b1e(k + 1) + b0e(k ) z 差分方程在形式上与微分方程相似,只不过微
τ →0
δ
k =−∞
(t
−
kT
)
=
limτ
τ →0
⋅δT
(t)
可见,开关函数近似成为一个幅度为无穷小的
周期性冲激序列。这个“无穷小”会给我们分析
带来不便,所以一般直接用幅度为 1 的周期性冲 激序列代替它,即:
+∞
∑ s(t) = δ (t − kT ) = δT (t) k =−∞
这样,抽样以后的信号为:
第七章 离散时间系统的时域分析 §7-1 概述
一、离散时间信号与离散时间系统 离散时间信号:只在某些离散的时间点上有值的
信号。 离散时间系统:处理离散时间信号的系统。 混合时间系统:既处理离散时间信号,又处理连
续时间信号的系统。
二、连续信号与离散信号 连续信号可以转换成离散信号,从而可以用 离散时间系统(或数字信号处理系统)进行 处理:
二、 抽样定理 显然,利用原来的信号在某些离散的时间点
上的值构成的信号,是否会损失信息?或者,在
何条件下,可以用抽样后的信号,不失真地还原 出原来的信号?
1、 抽样信号的频谱:
+∞
∑ fs (t) = f (t) δ (t − kT ) k =−∞
∑ Fs ( jω)
=
1 2π
F
(
jω
)
*
⎢⎣⎡ω
s
k
+∞
δ
=−∞
(ω
−
kω
s
⎤ )⎥⎦
∑ =
ωs 2π
+∞
F(
k =−∞
jω) *δ
(ω
− kωs )
∑ =
1 T
+∞
F(
k =−∞
jω
)
*δ
(ω
−
kω s
)
其中 ωs
=
2π T
,称为抽样(角)频率;T
称为抽样
(取样)周期。
可见,抽样后信号的谱是抽样以前的谱按抽
样(角)频率周期化的结果。
(a) 原信号 f (t)
的数相加。 2、 乘法: f (k ) = f1(k ) ⋅ f2 (k ) 3、 标量乘法: f (k ) = a ⋅ f1(k ) 4、 移序: f (k ) = f1(k − n) 当 n>0 时,信号向右移(后移)——>称为减序; 当 n<0 时,信号向左移(前移)——>称为增序。
离散信号的移序计算相当于连续时间信号 的时间平移计算。
通过迭代计算的方法,得到数值解。这种方 法用计算机求解比较方便,但是无法得到通 式。 例如:对于 Fibonacci 问题,有差分方程:
y(k+2)- y(k+1)-y(k)=0; Î y(k+2)= y(k+1)+y(k); Î y(k)= y(k-1)+y(k-2); 现在已知:y(0)=0,y(1)=1,则可以得到: y(2)=1,y(3)=2,y(4)=3,y(5)=5,y(6)=8,……
z 抽样也是一个线性处理过程,它满足齐次性和 叠加性。这是我们通过它达到用离散时间系统 处理连续信号的基础。
z 通过抽样可以将连续信号转化为离散数字信 号,从而可以用数字信号处理系统进行处理, 达到模拟信号处理无法达到的效果。
采 e(t) 样
A/D DSP D/A LPr(tF) 转换 处理 转换 滤波
2、
单位阶跃函数:
ε
(k
)
=
⎧1 ⎩⎨0
k ≥0 其它
这个函数与连续时间信号中的阶跃函数 ε (t) 相似。用它可以产生(或表示)单边信号 (这里称为单边序列)。 3、 单边指数序列: akε (k)
(a) a = 0.9
(d) a = −0.9
(b) a = 1
(e) a = −1
(c) a = 1.1
分计算变成了移序计算; z 差分方程也有阶,差分方程的阶定义为其中最
大移序与最小移序之差; z 求解差分方程也必须有初始条件,初始条件的
个数必须等于差分方程的阶数; z 与连续时间系统中的结论相似,线性移不变系
统可以用一个常系数差分方程描述。 z 因为差分方程可以很方便地用计算机求其数
值解,所以很多微分方程可以近似为差分方程 求近似数值解。
六、线性移不变离散时间系统 1、 线性离散时间系统
系统的激励和响应之间满足齐次性和叠加性