离散时间系统的时域分析
实验五 离散时间系统的时域分析
实验五 离散时间系统的时域分析一、实验目的:(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 )
new第三章离散时间系统的时域分析
3. 举例 • 例1 已知 x(n)=(n),y(-1)=0, 用迭代法解方程:
y(n) ay(n 1) x(n)
• 解:y(0)=ay(-1)+1=1 • y(1)=ay(0)+0=a • y(2)=ay(1)+0=a2 • • y(n)=ay(n-1)+0=an • y(n)=ay(n-1)+0=anu(n)
n y(n) 0.45(0.9) u(n) 0.5u(n) 自由响应 强迫响应
• 零输入响应和零状态响应
用边界条件求系数
C1
5
1
, C2
n
5
1
最终解
1 1 5 1 1 5 y ( n) 5 2 5 2
n
例3 求 y(n)+6y(n-1)+12y(n-2)+8y(n-3)=x(n) 的齐次解 • 解(有重根)
差分方程特解的形式 • • • • • • • • • 激励 x(n) 特解 yp(n)的形式 A(常数) C(常数) An C1n+C2 nk C1 nk+ C2 nk-1++ Ck+1 nkan an(C1 nk+ C2 nk-1++ Ck+1 ) sin(bn)或 C1sin(bn)+C2cos(bn) con(bn) an [sin(bn)或 an[C1sin(bn)+C2cos(bn)] cos(bn)]
– 常系数线性差分方程(递归关系式) – 后向(或右移) 差分方程;前向(或左移) 差分方程
例2 已知离散时间系统如图示,写出 系统的差分方程。
数字信号处理实验离散时间 LTI 系统的时域分析与 Z 域分析
实验一离散时间LTI系统的时域分析与Z域分析一、实验目的1、掌握用MATLAB求解离散时间系统的零状态响应、单位脉冲响应和单位阶跃响应;2、掌握离散时间系统系统函数零极点的计算方法和零极点图的绘制方法,并能根据零极点图分析系统的稳定性。
二、实验原理1、离散时间系统的时域分析(1)离散时间系统的零状态响应离散时间LTI系统可用线性常系数差分方程来描述,即MATLAB中函数filter可对式(1-1)的差分方程在指定时间范围内的输入序列所产生的响应进行求解。
函数filter的语句格式为:y=filter(b,a,x)其中,x为输入的离散序列;y为输出的离散序列;y的长度与x的长度一样;b与a分别为差分方程右端与左端的系数向量。
(2)离散时间系统的单位脉冲响应系统的单位脉冲响应定义为系统在 (n)激励下系统的零状态响应,用h(n)表示。
MATLAB求解单位脉冲响有两种方法:一种是利用函数filter;另一种是利用函数impz。
impz函数的常用语句格式为impz(b,a,n),其中b和a的定义见filter,n表示脉冲响应输出的序列个数。
(3)离散时间系统的单位阶跃响应系统的单位阶跃响应定义为系统在ε(n)激励下系统的零状态响应。
MATLAB求解单位脉冲响应有两种方法:一种是利用函数filter,另一种是利用函数stepz。
stepz函数的常用语句格式为stepz(b,a,N)其中,b和a的定义见filter,N表示脉冲响应输出的序列个数。
2、离散时间系统的Z域分析(1)系统函数的零极点分析离散时间系统的系统函数定义为系统零状态响应的z变换与激励的z变换之比,即如果系统函数H(z)的有理函数表示式为那么,在MATLAB中系统函数的零极点就可通过函数roots得到,也可借助函数tf2zp得到。
roots的语法格式为:Z=roots(b)%计算零点b=[b1b2…bmbm+1]P=roots(a)%计算极点a=[a1a2…anan+1]tf2zp的语句格式为[Z,P,K]=tf2zp(b,a)其中,b与a分别表示H(z)的分子与分母多项式的系数向量。
离散时间系统的时域分析
§7.1 引言
离散时间信号通过将连续时间信号进行取样得到
f t 4.2
3.1
采样(sampling)过程就是对模拟信号的 时间取离散的量化值过程——得到离 散信号。
1.5 0.9 2T 3T
o
3
f q t
T
4
t
幅值量化——幅值只能分级变化。
2 1
o
T
2T
3T
t
§7.1 引言
• 经过量化的离散时间信号称 为数字信号(digital signal)
经典法:齐次解 特解 时域分析 零输入响应 零状态响应 变换域分析: 拉氏变换法
离散时间系统——差分方程描述 差分方程的解法与微分方程类似
经典法:齐次解 特解 时域分析 零输入响应 零状态响应 变换域分析: z变换法
§7.2 取样信号与取样定理
• 取样定理(抽样定理)
• 通常将这种方程形式称为前向预测差分方程 (forward difference equation)
§7.3 离散时间系统的描述和模拟
• 差分方程与微分方程相比 在取样间隔Ts足够小时
dy( t ) y[( k 1)Ts ] y( kTs ) 微分方程 dt Ts 也可写做 dy( t ) y( kTs ) y[( k 1)Ts ] dt Ts
x n
3 4 5
1 2
9 10 11 6 7 8
22
n
一个周期
§7.1 引言
信号xn sin0.4n是否为周期信号?
0 0.4
2π
0
5π是无理数 所以为非周期的序列
§7.1 引言
• 离散信号 sin n0与连续信号 sin 0 t 的关系 2 对连续信号 f t sin2πf 0 t sinΩ0 t Ω0 T 离散点(时刻)nT’上的正弦值
实验三 离散时间系统的时域分析(附思考题程序)
实验三 离散时间系统的时域分析1.实验目的(1)理解离散时间信号的系统及其特性。
(2)对简单的离散时间系统进行分析,研究其时域特性。
(3)利用MATLAB 对离散时间系统进行仿真,观察结果,理解其时域特性。
2.实验原理离散时间系统,主要是用于处理离散时间信号的系统,即是将输入信号映射成的输出的某种运算,系统的框图如图所示:][n x ][n y Discrete-timesystme(1)线性系统线性系统就是满足叠加原理的系统。
如果对于一个离散系统输入信号为12(),()x n x n 时,输出信号分别为12(),()y n y n ,即:1122()[()]()[()]y n T x n y n T x n ==。
而且当该系统的输入信号为12()()ax n bx n +时,其中a,b 为任意常数,输出为121212[()()][()][()]()()T ax n bx n aT x n bT x n ay n by n +=+=+,则该系统就是一个线性离散时间系统。
(2)时不变系统如果系统的响应与激励加于系统的时刻无关,则该系统是时不变系统。
对于一个离散时间系统,若输入()x n ,产生输出为()y n ,则输入为()x n k -,产生输出为()y n k -,即:若()[()]y n T x n =,则[()]()T x n k y n k -=-。
通常我们研究的是线性时不变离散系统。
3.实验内容及其步骤(1)复习离散时间系统的主要性质,掌握其原理和意义。
(2)一个简单的非线性离散时间系统的仿真 参考:% Generate a sinusoidal input signalclf; n = 0:200; x = cos(2*pi*0.05*n); % Compute the output signal x1 = [x 0 0]; % x1[n] = x[n+1] x2 = [0 x 0]; % x2[n] = x[n] x3 = [0 0 x];% x3[n] = x[n-1]y = x2.*x2-x1.*x3; y = y(2:202); % Plot the input and output signalssubplot(2,1,1) plot(n, x)xlabel('Time index n'); ylabel('Amplitude'); title('Input Signal')subplot(2,1,2) plot(n,y)xlabel('Time index n'); ylabel('Amplitude');title('Output signal');(3)线性与非线性系统的仿真参考:% Generate the input sequencesclf; n = 0:40; a = 2; b = -3;x1 = cos(2*pi*0.1*n); x2 = cos(2*pi*0.4*n);x = a*x1 + b*x2;num = [2.2403 2.4908 2.2403];den = [1 -0.4 0.75];ic = [0 0]; % Set zero initial conditionsy1 = filter(num,den,x1,ic); % Compute the output y1[n]y2 = filter(num,den,x2,ic); % Compute the output y2[n]y = filter(num,den,x,ic); % Compute the output y[n]yt = a*y1 + b*y2; d = y - yt; % Compute the difference output d[n] % Plot the outputs and the difference signalsubplot(3,1,1) stem(n,y); ylabel('Amplitude');title('Output Due to Weighted Input: a \cdot x_{1}[n] + b \cdot x_{2}[n]');subplot(3,1,2) stem(n,yt); ylabel('Amplitude');title('Weighted Output: a \cdot y_{1}[n] + b \cdot y_{2}[n]');subplot(3,1,3) stem(n,d); xlabel('Time index n'); ylabel('Amplitude');title('Difference Signal');(4)时不变与时变系统的仿真参考:% Generate the input sequencesclf; n = 0:40; D = 10; a = 3.0; b = -2;x = a*cos(2*pi*0.1*n) + b*cos(2*pi*0.4*n);xd = [zeros(1,D) x]; num = [2.2403 2.4908 2.2403]; den = [1 -0.4 0.75];ic = [0 0]; % Set initial conditions% Compute the output y[n]y = filter(num,den,x,ic);% Compute the output yd[n]yd = filter(num,den,xd,ic);% Compute the difference output d[n]d = y - yd(1+D:41+D);% Plot the outputssubplot(3,1,1) stem(n,y); ylabel('Amplitude'); title('Output y[n]'); grid;subplot(3,1,2) stem(n,yd(1:41)); ylabel('Amplitude');title(['Output due to Delayed Input x[n - ', num2str(D),']']); grid;subplot(3,1,3) stem(n,d); xlabel('Time index n'); ylabel('Amplitude');title('Difference Signal'); grid;4.实验用MATLAB函数介绍在实验过程中,MATLAB函数命令plot, figure, stem, subplot, axis, grid on, xlabel, ylabel, title, clc等在不同的情况下具体表述也有所不同,应该在实验中仔细体会其不同的含义。
离散时间系统的时域分析
离散时间系统的时域分析离散时间系统是指系统输入和输出信号都是在离散的时间点上进行采样的系统。
时域分析是分析系统在时域上的性质和特征。
在离散时间系统的时域分析中,常用的方法包括冲击响应法、单位样值法和差分方程法等。
冲击响应法是通过对系统施加单个冲击信号,观察系统在输出上的响应来分析系统的时域特征。
冲击响应法的基本思想是将系统的输出表示为输入信号与系统的冲击响应之间的卷积运算。
冲击响应法适用于线性时不变系统,在实际应用中可以使用软件工具进行计算。
单位样值法是通过将系统输入信号取为单位样值序列,观察系统在输出上的响应来分析系统的时域特征。
单位样值法的基本思想是将系统的输出表示为输入信号与系统的单位样值响应之间的卷积运算。
单位样值法适用于线性时不变系统,可以用来计算系统的单位样值响应和单位样值响应序列。
差分方程法是通过建立系统输入和输出之间的差分方程来分析系统的时域特征。
差分方程法的基本思想是根据系统的差分方程,利用系统的初始条件和输入序列,递推计算系统的输出序列。
差分方程法适用于线性时不变系统,可以用来计算系统的单位样值响应和任意输入信号下的输出序列。
以上所述的方法是离散时间系统时域分析中常用的方法,通过这些方法可以获得系统的冲击响应、单位样值响应和任意输入信号下的输出序列,进而分析系统的时域特征和性质。
在实际应用中,根据系统的具体情况和需求,选择合适的方法进行时域分析,能够更好地理解离散时间系统的动态行为和响应特性。
离散时间系统的时域分析是研究系统在离散时间上的动态行为和响应特性的关键方法。
通过分析系统的时域特征,可以深入了解系统的稳定性、响应速度、频率选择性和滤波特性等方面的性能。
冲击响应法是离散时间系统常用的时域分析方法之一。
它通过施加一个单个的冲击信号,即输入信号序列中只有一个非零元素,然后观察系统在输出上的响应。
这样可以得到系统的冲击响应序列,它描述了系统对单位幕函数输入信号的响应情况。
冲击响应法的核心思想是将系统的输出表示为输入信号序列与系统的冲击响应序列之间的卷积运算。
信号与系统(第三版)第五章离散时间系统的时域分析
连续时间系统的信号在任意时刻都有取值,而离散时间系统的信 号只在离散时刻上取值。
离散时间系统的数学描述
02
差分方程
定义
差分方程是描述离散时间信号变化的数学方程,通常表示为y[n] = f(n) + g(n),其中y[n]是离散时间信号,f(n)和g(n)是已知的 离散时间信号。
类型
差分方程可以分为线性和非线性两种类型。线性差分方程是指方程中未知数的系数为常数且方程中未知数次数不超过1的差分方 程。
稳定性判据
通过判断系统的极点位置,确定系统的稳定性。
稳定性分析的意义
对于实际应用中的系统,稳定性是非常重要的性能指标。
系统的动态性能分析
动态性能的定义
描述系统在输入信号激励下,输出信号随时间变 化的特性。
动态性能的参数
包括超调和调节时间、上升时间和峰值时间等。
动态性能的分析方法
通过系统函数的Leabharlann 点和零点位置,以及时间常数等参数进行分析。
04 离散时间系统的时域响应 单击添加文本具体内容
离散时间系统 的定义与特点
离散时间系统的定义
离散时间系统
在时间上离散取样,信号在离散时刻上变化的系统。
离散时间信号
只在离散时刻上取值的信号。
离散时间系统分析
通过数学模型对离散时间信号和系统进行描述和分析 的方法。
离散时间系统的特点
时域离散
01
离散时间系统的状态变量和信号只在离散时刻上取值,时
定义
分类
稳定性判据
劳斯判据 通过求解劳斯表,判断系统的极点和稳定性。
赫尔维茨判据 通过判断系统的特征方程的根的性质,判断系统的 稳定性。
波波夫判据 通过求解波波夫矩阵,判断系统的稳定性。
离散系统时域分析_OK
例:设 y(k)+3y(k-1)+2y(k-2)=2k (k),y(0)=0, y(1)=2,求y(k)。
f(k)=ak(k)
|a| >
1
f(k)=ak(k)
|a| <
11
1
-2 -1 0 1 2 3
k
-2 -1 0 1 2 3
k
3
发散
收敛
5.正弦序列
f (k) Acos(kω0 )
0序列依次重复出现的频率。
2
ω 0
为有理数,正弦序列为周期序列。
f (k N ) A cosω[ 0(k N ) ] A cosω[ 0k ω0 N ]
any(k)+an-1y(k-1)+…+a1y(k-n+1)+a0y(k-n)=0(后向)
any(k+n)+an-1y(k+n-1)+…+a1y(k+1)+a0y(k)=0(前向)
对应的特征方程为:ann+an-1n-1+ + …+a1 + a0=0
1.特征根均为单根: 则齐次通解为:
1≠2≠…≠n
10
§5–2 离散时间系统的数学模型
一、线性时不变离散时间系统
1.离散系统:激励和响应都是离散信号的系统
f(k)
y(k)
离散时间系统
2.分类:亦可分为线性与非线性;时不变与时变;因果与非 因果等。
时不变: f(k) → y(k) f(k-m) → y(k-m)
因果系统:响应总是出现在激励之后。即: 当k < k0 ,f(k)
(2) 初始条件y(0), y(1),…, y(n-1)(与外施激励有关)代入完全解,可确 定待定常数Ci 。
离散时间系统的时域特性分析
离散时间系统的时域特性分析离散时间系统是指输入和输出均为离散时间信号的系统,如数字滤波器、数字控制系统等。
时域分析是研究系统在时间上的响应特性,包括系统的稳定性、响应速度、能否达到稳态等。
在时域分析中,我们通常关注系统的单位采样响应、阶跃响应和脉冲响应。
1. 单位采样响应单位采样响应是指当输入信号为单位脉冲序列时,系统的输出响应。
在时间域上,单位脉冲序列可以表示为:$$ u[n] = \begin{cases}1 & n=0\\ 0 & n \neq 0\end{cases} $$系统的单位采样响应可以表示为:$$ h[n] = T\{ \delta[n]\} $$其中,$T\{\}$表示系统的传输函数,$\delta[n]$表示单位脉冲序列。
通常情况下,我们可以通过借助系统的差分方程求得系统的单位采样响应。
对于一种具有一阶差分方程的系统,其单位采样响应可以表示为:2. 阶跃响应其中,$\alpha$为系统的传递常数。
3. 脉冲响应脉冲响应是指当输入信号为任意离散时间信号时,系统的输出响应。
其主要思路是通过将任意输入信号拆解成单位脉冲序列的线性组合,进而求得系统的输出响应。
设输入信号为$x[n]$,系统的脉冲响应为$h[n]$,则系统的输出信号$y[n]$可以表示为:$$ y[n] = \sum_{k=-\infty}^{\infty} x[k]h[n-k] $$在实际计算中,通常采用卷积算法实现脉冲响应的计算,即将输入信号和脉冲响应进行卷积运算。
总之,时域特性分析是对离散时间系统进行分析和设计时的基础。
对于实际工程应用中的系统,需要综合考虑其时域和频域特性,进而选择合适的滤波器结构、控制算法等来实现系统的优化设计。
《信号、系统与数字信号处理》第四章 离散时间系统的时域分析
当系统的阶数不高,并且激励不复杂时,可以用迭代 (递推)法求解差分方程。
例4-6 已知 yn ayn 1 xn ,且 yn 0, n 0 ,xn n ,求 yn 。
P
F[
p(t)]
F
Pn
e
jns
t
n
Pn F e jnst 2
Pn ns
n
n
由
fs t
f
tpt
Fs
1
2
F P
得:
Fs
1
2
F 2
n
Pn ns
Pn F ns
n
当开关函数 pt 是周期冲激序列时也称理想抽样,
此时
p t T t t nTs
折叠后位移:
yn x n m
m0 m0
x n 逐项右移(时延)m位; x n 逐项左移(时延)m位。
xn
3
2
1
3
0 12
n
-2 -1
x n
3
2
1
0
n
x n 1
3
2
-2
1
-1 0
n
(5)尺度变换
yn xmn
是 xn 序列每m点取一点形成的,即时间轴n压缩了m倍。
例m=2时
xn
3 2
1 01
-1
频谱组成的。
f t
F
0
t
pt
m 0 m
P
…
……
…
0 T 2T
t
fs t
s
0 s
Fs
1/T
…
…
Hale Waihona Puke …0 T 2Tt
- 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、 线性离散时间系统
系统的激励和响应之间满足齐次性和叠加性