离散线性时不变系统分析
线性时不变系统的稳定性分析
线性时不变系统的稳定性分析稳定性是控制系统理论中的重要概念，对于线性时不变系统来说尤其重要。
稳定性分析可以帮助我们确定系统的输出是否会在输入变化或干扰的情况下产生不受控制的波动或发散。
本文将探讨线性时不变系统的稳定性分析方法。
一、线性时不变系统的定义线性时不变系统（Linear Time-Invariant System，LTI系统）是指满足叠加性和时移不变性两个性质的系统。
叠加性指系统对输入的响应是可加的，时移不变性指系统对延时输入的响应是不变的。
线性时不变系统可以用微分方程或差分方程来描述。
二、稳定性的定义在系统稳定性分析中，我们关注的是系统的零输入响应或者零状态响应。
稳定性可以分为BIBO稳定性和渐近稳定性两种类型。
1. BIBO稳定性BIBO稳定性（Bounded-Input Bounded-Output Stability）是指当输入有界时，系统的输出也是有界的。
如果对于任意有界的输入信号，系统的输出都有界，则系统是BIBO稳定的。
2. 渐近稳定性渐近稳定性是指当输入信号趋于稳定时，系统的输出也趋于稳定。
如果对于任意渐近稳定的输入信号，系统的输出也渐近稳定，则系统是渐近稳定的。
三、稳定性分析方法稳定性分析的常用方法包括传输函数法、状态空间法和频域法。
下面将分别介绍这三种方法。
1. 传输函数法传输函数法是用传输函数来描述系统的稳定性。
传输函数是输入和输出的关系，它是Laplace变换或Z变换的比值。
对于连续时间系统，传输函数可以表示为H(s)；对于离散时间系统，传输函数可以表示为H(z)。
通过分析传输函数的极点（Pole）可以判断系统的稳定性。
对于连续时间系统，如果传输函数的极点都位于左半平面，则系统是BIBO稳定的；如果传输函数的极点有位于右半平面的，则系统是不稳定的。
对于离散时间系统，如果传输函数的极点都位于单位圆内部，则系统是BIBO稳定的；如果传输函数的极点有位于单位圆外部的，则系统是不稳定的。
离散线性时不变系统分析
实验六 离散线性时不变系统分析一、 实验目的1. 掌握离散LSI 系统的单位序列响应、单位阶跃响应和任意激励下响应的MATLAB 求解方法。
2. 掌握离散LSI 系统的频域分析方法；3. 掌握离散LSI 系统的复频域分析方法；4. 掌握离散LSI 系统的零极点分布与系统特性的关系。
二、实验原理及方法1.离散LSI 系统的时域分析描述一个N 阶线性时不变离散时间系统的数学模型是线性常系统差分方程，N 阶LSI 离散系统的差分方程一般形式为)()(0i n x b k n y a Mi i N k k -=-∑∑== （6.1） 也可用系统函数来表示12001212120()()()()()1MiM ii M NNkN k k b zb b z b z b z Y z b z H z X z a z a z a z a z a z----=----=++++====++++∑∑ （6.2）系统函数()H z 反映了系统响应和激励间的关系。
一旦上式中k a ，i b 的数据确定了，系统的性质也就确定了。
特别注意0a 必须进行归一化处理，即01a =。
对于复杂信号激励下的线性系统，可以将激励信号在时域中分解为单位序列或单位阶跃序列的线性叠加，把这些单元激励信号分别加于系统求其响应，然后把这些响应叠加，即可得到复杂信号作用于系统的零状态响应。
因此，求解系统的单位序列响应和单位阶跃响应尤为重要。
由图6-1可以看出一个离散LSI 系统响应与激励的关系。
()h n ()H z ()x n ()X z ()()()Y z X z H z =()()*()y n x n h n =图6-1 离散LSI 系统响应与激励的关系(1) 单位序列响应（单位响应）单位响应()h n 是指离散LSI 系统在单位序列()n δ激励下的零状态响应，因此()h n 满足线性常系数差分方程(6.1)及零初始状态，即()()N Mkik i a h n k b n i δ==-=-∑∑， (1)(2)0h h -=-== (6.3)按照定义，它也可表示为()()()h n h n n δ=* (6.4) 对于离散LSI 系统，若其输入信号为()x n ，单位响应为()h n ，则其零状态响应()zsy n为()()*()zs y n x n h n = (6.5)可见，()h n 能够刻画和表征系统的固有特性，与何种激励无关。
系统的时域分析 线性时不变系统的描述及特点 连续时间LTI系统的响应
y x (t ) K1e 2t K 2 e 3t
y(0)=yx(0)=K1+K2=1 y' (0)= y'x(0)= 2K13K2 =3
解得 K1= 6，K2= 5
y x (t ) 6e 2t 5e 3t , t 0
18
[例] 已知某线性时不变系统的动态方程式为: y" (t)+4y ' (t) +4y (t) = 2f ' (t )+3f(t), t>0 系统的初始状态为y(0) = 2，y'(0) = 1， 求系统的零输入响应yx(t)。 解: 系统的特征方程为 系统的特征根为
2t
Be
4t
1 y (0) A B 1 3 解得 A=5/2，B= 11/6 1 y ' (0) 2 A 4 B 2 3
5 2t 11 4t 1 t y(t ) e e e , t 0 2 6 3
12
1 t e 3
系统的几个概念：
9
[例] 已知某二阶线性时不变连续时间系统的动态方程
y" (t ) 6 y' (t ) 8 y(t ) f (t ), t 0
初始条件y(0)=1, y '(0)=2, 输入信号f (t)=et u(t)， 求系统的完全响应y(t)。
解:
(1) 求齐次方程y''(t)+6y'(t)+8y(t) = 0的齐次解yh(t)
11
[例] 已知某二阶线性时不变连续时间系统的动态方程
y" (t ) 6 y' (t ) 8 y(t ) f (t ), t 0
离散时间系统的时域特性分析实验报告
信号、系统与信号处理实验报告实验一、离散时间系统的时域特性分析姓名:学号:班级：专业：一．实验目的线性时不变（LTI）离散时间系统在时域中可以通过常系数线性差分方程来描述，冲激响应列可以刻画时域特性。
本次实验通过使用MATLAB函数研究离散时间系统的时域特性，以加深对离散时间系统的差分方程、冲激响应和系统的线性和时不变性的理解。
二．基本原理一个离散时间系统是将输入序列变换成输出序列的一种运算。
离散时间系统中最重要、最常用的是“线性时不变系统”。
1.线性系统满足叠加原理的系统称为线性系统，即若某一输入是由N个信号的加权和组成的，则输出就是系统对这几个信号中每一个输入的响应的加权和。
即那么当且仅当系统同时满足和时，系统是线性的。
在证明一个系统是线性系统时，必须证明此系统同时满足可加性和比例性，而且信号以及任何比例系数都可以是复数。
2.时不变系统系统的运算关系在整个运算过程中不随时间（也即序列的先后）而变化，这种系统称为时不变系统（或称移不变系统）。
若输入的输出为，则将输入序列移动任意位后，其输出序列除了跟着位移外，数值应该保持不变，即则满足以上关系的系统称为时不变系统。
3.常系数线性差分方程线性时不变离散系统的输入、输出关系可用以下常系数线性差分方程描述：当输入为单位冲激序列时，输出即为系统的单位冲激响应。
当时，是有限长度的，称系统为有限长单位冲激响应（FIR）系统；反之，则称系统为无限长单位冲激响应（IIR）系统。
三．实验内容及实验结果1.实验内容考虑如下差分方程描述的两个离散时间系统：系统1：系统2：输入：（1）编程求上述两个系统的输出，并画出系统的输入与输出波形。
（2）编程求上述两个系统的冲激响应序列，并画出波形。
（3）若系统的初始状态为零，判断系统2是否为时不变的？是否为线性的？2.实验结果（1）编程求上述两个系统的输出和冲激响应序列，并画出系统的输入、输出与冲激响应波形。
clf;n=0:300;x=cos((20*pi*n)/256)+cos((200*pi*n)/256);num1=[0.5 0.27 0.77];den1=[1];num2=[0.45 0.5 0.45];den2=[1 -0.53 0.46];y1=filter(num1,den1,x);y2=filter(num2,den2,x);subplot(3,1,1);stem(n,x);xlabel('时间信号');ylabel('信号幅度');title('输入信号');subplot(3,1,2);stem(y1);xlabel('时间信号n');ylabel('信号幅度');title('输出信号');subplot(3,1,3);stem(y2);xlabel('时间序号n ');ylabel('信号幅度');title('冲激响应序列');（2）N=40;num1=[0.5 0.27 0.77];den1=[1];num2=[0.45 0.5 0.45];den2=[1 -0.53 0.46];y1=impz(num1,den1,N);y2=impz(num2,den2,N);subplot(2,1,1);stem(y1);xlabel('时间信号n ');ylabel('信号幅度');title('³冲激响应');subplot(2,1,2);stem(y2);xlabel('时间信号n ');ylabel('信号幅度');title('³冲激响应');1.应用叠加原理验证系统2是否为线性系统：clear allclcn = 0 : 1 : 299;x1 = cos(20 * pi * n / 256);x2 = cos(200 * pi * n / 256);x = x1 + x2;num = [0.45 0.5 0.45];den = [1 -0.53 0.46];y1 = filter(num, den, x1);y2 = filter(num, den, x2);y= filter(num, den, x);yt = y1 + y2;figuresubplot(2, 1, 1);stem(n, y, 'g');xlabel('时间信号n');ylabel('信号幅度');axis([0 100 -2 2]);grid;subplot(2, 1, 2);stem(n, yt, 'r');xlabel('时间信号n');ylabel('信号幅度');axis([0 100 -2 2]);grid;2.应用时延差值来判断系统2是否为时不变系统。
随机序列通过离散线性系统分析_-两种分析方法
问题：给定输入和离散线性系统的特性；求解：输出的统计特征。
✓冲激响应法✓频谱法]离散线性时不变系统的描述方法：✓单位样值响应✓系统传递函数⏹冲激响应法⏹频谱法⏹常用时间序列模型][][][][][]k Y n h k X n k h n X n +∞=−∞=−=∗∑系统输出：均值:[]([])([()])[][]Y X m n E Y n E L X n h n m n ===∗X m []Y n[,][][][,]121212Y X R n n h n h n R n n =∗∗[,]([][])[][,]1212212XY X R n n E X n Y n h n R n n ==∗相关函数:][][][][][]k Y n h k X n k h n X n +∞=−∞=−=∗∑系统输出：[][][]XY X R m h m R m =−∗[][][][]Y X R m h m h m R m =−∗∗[]()0Y X X k m m h k m H +∞=−∞=∑若X [n ]平稳：则输出平稳，而且输入输出联合平稳2. 频谱法若X [n ]平稳：2()()()|()|()Y XY X G H G H G ωωωωω==*()()()XY X G H G ωωω=1()()()XY X G z H z G z −=1()()()()()()Y XY X G z H z G z H z H z G z −==若用z 变换表示：则输出平稳，而且输入输出联合平稳只适用于平稳随机序列的分析例1：一个平稳随机序列X[n]的自相关函数为，线性系统的单位样值响应是, 求输出Y(n)的自相关函数及功率谱密度。
2[]m σδ[],nh n r =0,||1n r ≥<[][][][]Y X R m h m h m R m =−∗∗2[][]h m h m =−∗σ222201mm k k r r r r ∞==−∑σσ[解]20k m kk r r ∞+=⋅∑σ当m ≧0时，例1：一个平稳随机序列X[n]的自相关函数为，线性系统的单位样值响应是, 求输出Y(n)的自相关函数及功率谱密度。
信号分析与处理第5章 离散时间线性时不变系统的分析
yzs (n) 3yzs (n 1) 2 yzs (n 2) x(n)
对于因果系统，有 yzs (i) 0 i 0 ，因此，
y(i) yzi (i)
i0
y(i) yzi (i) yzs (i) i 0
求 N 阶差分方程，需用到 N 个初始条件。求零输
入响应的 N 个初始条件为 yzi (i) i 1,2, , N ；求零状态
响应的 N 个初始条件为 yzs (i) i 0,1, , N -1 ；求全响应的 N
5.3 离散时间线性时不变系统的 Z 域分析
5.3.1 Z 变换解差分方程 5.3.2 系统函数 5.3.3 由系统函数的零极点分布确定时域特性 5.3.4 系统的因果性和稳定性
5.4 离散时间线性时不变系统的频域分析
5.4.1 系统的频率特性 5.4.2 离时间信号通过系统的频域分析 5.4.3 理想低通数字滤波器
yzi (n) Czi1 (1)n Czi2 (2)n
n 0
代入初始条件
yzi
(1)
y(1)
0,
yzi
(2)
y(2)
1 2
，有
yzi (1)
Czi1
1 2 Czi2
0
yzi (2)
Czi1
1 4
Czi
2
1 2
解得 Czi1 1， Czi2 2。故系统的零输入响应为
个初始条件为 y(i) i 0,1, , N -1。
2、离散时间线性时不变系统响应的迭代求解
差分方程可用迭代法求解。差分方程可改写为
第二章 线性时不变系统
9
例5 y[n] 6,5,24,13,22,10,n 0,1,2,3,4,5 h[n] 3,1,4,2 n 0,1,2,3
y[n] x[n]h[n] 求 x[n]
2 t 5t2 x(t)
x[n] x[k] [n k] 离散的信号分解成脉冲
k
信号的 线性组合的形式
把任意一个序列表示成一串移位的单位脉冲序列 [n k]
的线性组合，而这个线性组合式中的权因子就是 x[k]
4
二. 离散时间线性时不变系统卷积和表示
[n] h[n]
[n k] h[n k]
时不变
x[k] [n k] x[k]h[n k] 齐次性
11
二. 连续时间线性时不变系统的卷积积分表示
(t) h (t)
(t k)
x(k) (t k)
x(k) (t k)
k
h (t k)
时不变
x(k
)h
(t
k
)
齐次性
x(k)h (t k) 可加性
k
xˆ(t)
yˆ (t )
y(t) x( )h(t )d x(t) h(t)
12
卷积的计算
（1）由定义计算卷积积分
例：设某一线性时不变系统的输入为x(t)，其单位冲
激响应为h(t) x(t) eatu(t) , a 0 h(t) u(t)
试求 x(t) h(t)
x(t) h(t) ea u( )u(t )d
t ea d ,
0
t0
0,
t0
1 1 eat u(t) a
1
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)
离散时间系统稳定的充要条件
离散时间系统稳定的充要条件离散时间系统是指系统的输入和输出在时间上是离散的情况下进行的系统分析和设计。
而离散时间系统的稳定性是一个重要的性质，它决定了系统是否能够在一定范围内保持稳定的输出。
本文将介绍离散时间系统稳定性的充要条件。
一、离散时间系统的稳定性概念稳定性是指系统在有限时间内是否能够保持有限的幅值，而不会出现无限增长或发散的情况。
对于离散时间系统而言，其稳定性可以分为两类：绝对稳定和相对稳定。
绝对稳定是指系统的输出在有限时间内始终保持有限的幅值，不会发散或无限增长。
相对稳定是指系统的输出在有限时间内保持有限的幅值，但可能会在无穷时间后发散或无限增长。
二、离散时间系统的稳定性充要条件1. 线性时不变系统对于线性时不变系统而言，其稳定性充要条件是系统的传递函数的极点都位于单位圆内。
也就是说，系统的所有极点的模长都小于1。
2. 有限冲激响应系统对于有限冲激响应系统而言，其稳定性充要条件是系统的冲激响应是绝对可和的。
也就是说，系统的冲激响应的绝对和是有限的。
3. 时变系统对于时变系统而言，其稳定性充要条件是系统的输入和输出序列都是绝对可和的，并且系统的输入和输出序列的绝对和都是有界的。
4. 有限差分方程系统对于有限差分方程系统而言，其稳定性充要条件是系统的差分方程的根都位于单位圆内。
也就是说，系统的所有根的模长都小于1。
5. 正态系统对于正态系统而言，其稳定性充要条件是系统的所有特征值的实部都小于等于零。
6. 离散时间系统的Lyapunov稳定性对于离散时间系统而言，其稳定性充要条件是系统的状态方程存在一个正定矩阵，使得系统的状态的Lyapunov函数是递减的。
三、离散时间系统的稳定性判定方法除了以上充要条件外，还可以通过以下方法判断离散时间系统的稳定性：1. 构造系统的Lyapunov函数。
通过构造系统的Lyapunov函数来判断系统的稳定性。
如果系统的状态的Lyapunov函数是递减的，则系统是稳定的。
实验三 线性时不变(LTI)连续系统的时域分析
执行结果
实验任务 1：LTI系统的微分方程y''(t)+2y'(t)+y(t)=f'(t)+2f(t)，激励f (t)=e-2tε(t)，
（1） 利用 impulse 函数获得冲激响应； （2） 利用 lsim 函数求取零状态响应； （3） 用卷积分析法计算其零状态响应； 要求：在一个图形窗口里以 3 个子图形式绘制冲激响应和两种方法得到的零状 态响应的波形。 （4）改变系统的 a 系数矩阵，观察冲激响应和零状态响应时域波形的变化情况。 建议 a 系数向量分别如下取值讨论。
用。
一、实验目的
1. 掌握系统时域分析常用函数的使用方法； 2. 理解系统特征根对系统时域特性的影响；
二、实验原理及内容 2.1 连续系统的时域分析 2.1.1 连续系统时域分析的几个常用函数
设 LTI 连续时间系统的微分方程为
a 2 y''(t)+a 1 y'(t)+a 0 y(t)=b 2 f''(t)+b 1 f'(t)+b 0 f(t)
将积分变量离散化即将用n替代d用替代只要时域取样间隔足够小上式可近似为再把观察响应时刻离散化即将t用k替换只要足够小通常将fn简记为fnhkn简记为hknyk简记为yk这样上式便可表示为因此两个连续信号ft和ht的卷积yt可用ft和ht的取样信号f均取整数
实验三 线性时不变（LTI）连续系统的时域分析
2.1.2 连续信号卷积的近似计算
连续系统的零状态响应 y(t)可通过输入信号 f(t)与系统冲激响应 h(t)的卷积求 得；但是计算机只能处理数字信号，不能直接处理模拟信号，因此，卷积积分不 能直接用计算机计算。为了解决这个问题，可以将连续信号用取样信号来近似表 示，利用卷积和近似求得卷积积分。下面就连续信号卷积积分的近似计算进行简 单推导。
信号处理与系统分析 第2章线性时不变系统
从波形的角度来观察离散时间信号，它可以 看成是由许多加权了的单位冲激信号组合 而成的
x[n] x[1] [n 1] x[0] [n] x[2] [n 2]
对于任意的离散时间信号：
累加序号 自变量
加权值 移位的冲激信号
x[n]
k
x[k ] [n k ]
n
卷积公式是无穷多项求和，而我们实际遇到的常 常是有限长度序列，特别是在计算机离线处理的场 合，因为计算机不可能处理无穷多的信息。 在进行有限长度的序列的卷积时候，长度为N和M 的2个序列作卷积时，反转序列从左到右进入重叠 直至移出重叠，只有存在重叠项时，卷积和才可能 非零。 卷积序列的长度为M+N-1。
求解系统响应的卷积方法是系统分析的重要工具。
单位冲激响应h[n]完全描述了线性时不变系统的变换 规律。不同的系统输入，都在h[n]的作用下产生相应的 响应，因此，给定了一个LTI系统的单位冲激响应h[n]就 等于给定了该系统。
从计算某一个特定点的角度来看
yy [n [n 0]
k k
第2章 线性时不变系统
线性时不变（简称LTI，Linear， Time-invariant）系统
为什么引入LTI ？
如果不对系统的性质加以限制，那么分析 一个系统将是十分困难的。 给系统加上线性和时不变性的限制，那么 系统的分析将变得十分简便。 LTI系统的分析还为非线性系统的分析方法 提供了思路。例如，线性时不变系统可以 用冲激响应来表达，非线性系统可以用 Volterra级数来表达。
上式应该理解为许多以为n自变量的函数的相 加，而不是数值相加。
许多移了位的冲激信号的加权和，构成了x[n] 。
特别地，我们有
数字信号处理讲义线性时不变系统的变换分析
第５章线性时不变系统的变换分析[教学目的]1．了解LTI 系统频率响应的概念；2．掌握线性常系数差分方程所表征系数的系统函数的方法；3．掌握有理系统频率响应分析方法4．理解线性相位系统、广义线性相位系统与因果广义线性相位系统的概念，几类线性相位系统。
[教学重点与难点]重点：1．线性常系数差分方程所表征系数的系统函数的方法；2．有理系统频率响应分析方法；3．几类线性相位系统。
难点：1． 有理系统频率响应分析方法几类线性相位系统5.1 LTI 系统的频率响应前面已经讨论过，在时域中，一个线性时不变系统完全可以由它的单位脉冲响应h (n )来表示。
对于一个给定的输入x (n )，其输出y (n )为对等式两端取Z 变换，得则 (5-1)两边做离散傅立叶变换有:|Y (e j ω)|=|H (e j ω)|·|X (e j ω)| (5-2)|Y (e j ω)|=|H (e j ω)|·|X (e j ω)|arg ［Y (e j ω)］=arg ［H (e j ω)］+arg ［X (e j ω)］|H (e j ω)| 幅度响应 : 增益/幅频特性调整输入信号各频率分量的相对强度(幅度)关系Arg[H (e j ω)] 频率响应的相位响应 : 相移/相频特性调整输入信号各频率分量的相对位置(相位)关系H (e j ω) 调整输入信号各频率分量的相对大小(幅度)及位置(相∑∞-∞=-=*=m m n h m x n h n x n y )()()()()()()()(z X z H z Y =)()()(z X z Y z H =位)关系5.1.1理想低通滤波器的选择性5.1.2相位失真与延时线性相位 : 不会改变信号的相对位置,时延相同线性相位的效应 : 时延 非线性相位：改变信号的相对位置时延不相同⎩⎨⎧≤<≤=πωωωωω||,0||,1)(c c j H n n n h c F πωsin ][=−→←()()|()|j H j H j H j eωωω=0 : ()near Phase H j t ωω=-0()H j t ωω≠-5.2 用线性常系数差分方程所表征系统的系统函数一个线性时不变系统也可以用常系数线性差分方程来表示，其N 阶常系数线性差分方程的一般形式为若系统起始状态为零，这样就可以直接对上式两端取Z 变换，利用Z 变换的线性特性和移位特性可得这样就得到系统函数为（5-3）由此看出系统函数分子、分母多项式的系数分别就是差分方程的系数。
信号与系统-第二章线性时不变系统
n
1
k
f1 (k )
f2 (0
k)
3,
k
f1 (k )
f2 (1 k)
3,
n0 n 1
k
f1 (k )
f2(2 k)
1,
0,
n2 n14 3
三. 卷积和的计算：(3)列表法
分析卷积和的过程，可以发现有如下特点：
① x(n与) 的h(所n)有各点都要遍乘一次；
② 在遍乘后，各点相加时，根据 x(k)h(n k)， k
x (t) x(t)
20
x(t) x (t)
x(k)
t
0
k (k 1)
引用 (t，) 即：
(t)
1
/ 0
0t otherwise
则有：
(t
)
1 0
0t otherwise
21
第 个k 矩形可表示为： x(k) (t k)
这些矩形叠加起来就成为阶梯形信号 x，(t)
即： x (t) x(k) (t k) k 当 时0 ， k d
un 4 ak
an3
1un 4
k 0
a 1
9
例4： x(n) nu(n) 0 1 h(n) u(n)
x(k) ku(k)
1
0
k ...
h(n k) u(n k)
1
k
0
n
y(n) x(n) h(n)
x(k)h(n k) ku(k)u(n k)
k
k
u(n) n k 1 n1 u(n)
例2 :
1 x(t) 0
h( )
2T
0t T otherwise
信号与系统 第二章 线性时不变系统的时域分析
外加信号 常数A
特解 常数B
r 1i k t i r 1 i 1
tr
sin t或cos t
eλt
k1 cost k2 sin t keλt， λ不是方程的特征根 kteλt， λ是方程的特征根
k t
i 1 i
r 1
r 1i t
e ， λ是方程的r阶特征重根
一、微差分方程的建立以及经典解法
'' 1
di1 (t ) 1 t L i2 ( )d R2i2 (t ) f (t ) dt C
一、微差分方程的建立以及经典解法
1 (2) Li (t ) i2 (t ) R2i2 ' (t ) f ' (t ) C 1 ( R2i2 i2 ( )d ) 1 U C i2 (t ) y (t ) (3) i1 i2 i2 (4) R2 R1 R1
(1)
t
i ( )d
1 (2) Li (t ) i2 (t ) R2i2 ' (t ) f ' (t ) C 1 ( R2i2 i2 ( )d ) 1 U C i2 (t ) y (t ) (3) i1 i2 i2 (4) R2 R1 R1
例题，已知线性时不变系统方程如下： y˝(t)+6y΄(t)+8y(t)= f(t), t>0. 初始条件y(0)=1， y΄(0)=2，输入信号f(t)=e-tu(t) , Q求系统的完全响应y(t)。
解：1）求方程的齐次解 特征方程为：m2+6m+8=0 显然特征根为：m1=-2，m2=-4
故原方程的齐次解为：yn(t)= Ae-2t+Be-4t
信号与系统matlab实验线性时不变系统的时域分析(最新整理)
答案
1. x n hn u n u n 4 ;
nx=0:9;x=ones(1,length(nx)); nh=0:4;h=ones(1,length(nh)); y=conv(x,h); % 下限=下限1+下限2 ny_min=min(nx)+min(nh); % 上限=上限1+上限2 ny_max=max(nx)+max(nh); ny=ny_min:ny_max; subplot(3,1,1);stem(nx,x); xlabel('n');ylabel('x(n)');axis([ny_min ny_max 0 max(x)]); subplot(3,1,2);stem(nh,h); xlabel('n');ylabel('h(n)');axis([ny_min ny_max 0 max(h)]); subplot(3,1,3);stem(ny,y); xlabel('n');ylabel('x(n)*h(n)');axis([ny_min ny_max 0 max(y)]);
到连续卷积的数值近似，具体算法如下：
y=conv(x,h)*dt
% dt 为近似矩形脉冲的宽度即抽样间隔
例 2-2：采用不同的抽样间隔 值，用分段常数函数近似 x t u t u t 1 与
h t sin t u t u t π 的 卷 积 ， 并 与 卷 积 的 解 析 表 达 式
x(t)
h(t)
1 0.5
0 0 0.5 1 1.5 2 2.5 3 3.5 t
1 0.5
0 0 0.5 1 1.5 2 2.5 3 3.5 t
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
实验六 离散线性时不变系统分析一、 实验目的1. 掌握离散LSI 系统的单位序列响应、单位阶跃响应和任意激励下响应的MATLAB 求解方法。
2. 掌握离散LSI 系统的频域分析方法；3. 掌握离散LSI 系统的复频域分析方法；4. 掌握离散LSI 系统的零极点分布与系统特性的关系。
二、实验原理及方法1.离散LSI 系统的时域分析描述一个N 阶线性时不变离散时间系统的数学模型是线性常系统差分方程，N 阶LSI 离散系统的差分方程一般形式为)()(0i n x b k n y a Mi i N k k -=-∑∑== （6.1） 也可用系统函数来表示12001212120()()()()()1MiM ii M NNkN k k b zb b z b z b z Y z b z H z X z a z a z a z a z a z----=----=++++====++++∑∑ （6.2）系统函数()H z 反映了系统响应和激励间的关系。
一旦上式中k a ，i b 的数据确定了，系统的性质也就确定了。
特别注意0a 必须进行归一化处理，即01a =。
对于复杂信号激励下的线性系统，可以将激励信号在时域中分解为单位序列或单位阶跃序列的线性叠加，把这些单元激励信号分别加于系统求其响应，然后把这些响应叠加，即可得到复杂信号作用于系统的零状态响应。
因此，求解系统的单位序列响应和单位阶跃响应尤为重要。
由图6-1可以看出一个离散LSI 系统响应与激励的关系。
()()()z X z H z =()()*()n x n h n图6-1 离散LSI 系统响应与激励的关系(1) 单位序列响应（单位响应）单位响应()h n 是指离散LSI 系统在单位序列()n δ激励下的零状态响应，因此()h n 满足线性常系数差分方程(6.1)及零初始状态，即()()N Mkik i a h n k b n i δ==-=-∑∑， (1)(2)0h h -=-== (6.3)按照定义，它也可表示为()()()h n h n n δ=* (6.4) 对于离散LSI 系统，若其输入信号为()x n ，单位响应为()h n ，则其零状态响应()zsy n为()()*()zs y n x n h n = (6.5)可见，()h n 能够刻画和表征系统的固有特性，与何种激励无关。
一旦知道了系统的单位响应()h n ，就可求得系统对任何输入信号()x n 所产生的零状态响应()zsy n 。
MATLAB 提供了专门用于求连续系统冲激响应的函数impz()，其调用格式有 [h,n]=impz(b,a)求解离散系统的单位响应，其中012[,,,,]M b b b b b =，12[1,,,,]Na a a a =，[0,1,2,]n '=；[h,n]=impz(b,a,N)求解离散系统的单位响应，采样点数由N 确定，[0,1,2,,N-1]n '=；impz(b,a) ：在当前窗口，用stem(n,h)绘出图形。
（2）单位阶跃响应单位阶跃响应()s n 是指离散LTI 系统在单位阶跃序列()u n 激励下的零状态响应，它可以表示为()()()()nm s n u n h n h m =-∞=*=∑ (6.6)上式表明，离散LSI 系统的单位阶跃响应是单位响应的累加和，系统的单位阶跃响应和系统的单位响应之间有着确定的关系，因此，单位阶跃响应也能完全刻画和表征一个LSI 系统。
MATLAB 提供了专门用于求离散系统单位阶跃响应的函数stepz( )，其调用格式有 [s,n]=stepz(b,a) ：求解离散系统的单位阶跃响应，其中012[,,,,]Mb b b b b =，12[1,,,,]N a a a a =，[0,1,2,]n '=；[s,n]=stepz(b,a,N) ：求解离散系统的单位阶跃响应，采样点数由N 确定，[0,1,2,,N-1]n '=；stepz(b,a) ：在当前窗口，用stem(n,s)绘出图形。
(3) 任意激励下的零状态响应已经知道，离散LSI 系统可用常系数线性差分方程(6.1)式来描述，Matlab 提供的函数dlsim( )能对上述差分方程描述的离散LSI 系统的响应进行仿真，该函数不仅能绘制指定时间范围内的系统响应波形图，而且还能求出系统响应的数值解。
其调用格式有dlsim(b,a, x) ：求解输入序列为x 的零状态响应需要特别强调的是，Matlab 总是把由分子和分母多项式表示的任何系统都当作是因果系统。
所以，利用impz (b,a)，stepz(b,a)，d lsim(b,a,x)函数求得的响应总是因果信号。
同时，卷积和也是LSI 系统求解零状态响应的重要工具之一。
假设系统的输入信号为()x n ，单位响应为()h n ，则系统的零状态响应()zs y n 可由（6.5）式求解。
Matlab 提供了专门用于求离散系统卷积和的函数conv( )，其调用格式有y=conv(x,h) ：求解序列x ，h 的卷积和，若序列x 的长度为n1，序列h 的长度为n2，卷积和y 的长度为n1+n2-1。
这一点需要特别注意，否则，作图时容易造成横纵坐标长度不匹配。
（4）带初始状态的任意激励下的全响应任意激励下的离散LSI 系统的全响应为零输入响应和零状态响应之和，表示为()()()zi zs y n y n y n =+ （6.7）在理论学习的过程中，同学们对低阶差分方程的求解已颇为头痛，高阶差分方程直接求解几乎不可能。
Matlab 提供了用于求离散系统全响应的函数filter( )，其调用格式有y=filter( b,a,x) ：求解零状态响应；y=filter( b,a,x,zi) ：求解初始条件为zi 的系统的全响应，zi 向量的长度为max(length(a),length(b))-1，返回值为系统的全响应。
z = filtic(b,a,y,x)：将初始状态转换为初始条件，其中[(1),(2),(3),,()]x x x x x m =----，[(1),(2),(3),,()]y y y y y n =----；z = filtic(b,a,y)：将初始状态转换为初始条件，其中0x =，[(1),(2),(3),,()]y y y y y n =----。
2 离散LSI 系统的复频域（Z 域）分析 （1）利用Z 变换解差分方程在前面图6-1中表示了离散系统的响应与激励的关系，由图可知，系统的响应既可以用时域的方法求解，也可以用Z 域的方法求解。
当已知系统输入序列的Z 变换()X z ，系统函数()H z 时，系统响应序列的Z 变换可由()()()Y z X z H z =求出。
Matlab 提供了用于求序列Z 变换和Z 反变换的函数，其调用格式有X=ztrans(x)：求无限长序列x 的Z 变换，返回Z 变换的表达式，注意这里x ，X 都是符号表达式；x=iztrans(X)：求X （z ）的Z 反变换x(n)，返回Z 反变换的表达式，注意这里x ，X 都是符号表达式；[r,p,c]=residuez(b,a)：把b(z)/a(z)展开成部分分式；[b,a]=residuez(r,p,c):根据部分分式的r 、p 、c 数组，返回有理多项式。
（2）系统的零极点分布与系统因果性和稳定性的关系因果系统的单位响应)(n h 一定满足当0<n 时，)(n h =0，那么其系统函数)(z H 的收敛域一定包含∞点，即∞点不是极点，极点分布在某个圆的圆内，收敛域是圆外区域。
系统稳定要求 ∞<∑∞-∞=n n h |)(| ，对照z 变换定义，系统稳定要求收敛域包含单位圆。
如果系统因果且稳定，收敛域包含∞点和单位圆，那么收敛域可表示为：10,||<<∞≤<r z r (6.8)MATLAB 提供了用于求系统零极点的函数，其调用格式有 roots()：利用多项式求根函数来确定系统函数的零极点位置；roots(a)：求极点位置，a 为系统函数)(z H 分母多项式所构成的系数向量； roots(b)：求零点位置，b 为系统函数)(z H 分子多项式所构成的系数向量；zplane(b,a)：绘制由行向量b 和a 构成的系统函数的零极点分布图；zplane(z,p)：绘制由列向量z 确定的零点、列向量p 确定的极点构成的零极点分布图。
（3）系统的零极点分布与系统频率响应的关系 将式(6.2)因式分解，得到∏∏=-=---=Nk k Mmm z d z c A z H 1111)1()1()( (6.9) 式中，00a bA =，m c 是)(z H 的零点，k d 是其极点。
A 参数影响频率响应的幅度大小，影响系统特性的是零点m c 和极点k d 的分布。
下面采用几何方法研究系统零极点分布对系统频率特性的影响。
将式(6.9)的分子、分母同乘以M N z +，得到：∏∏∏∏==-=-=---=--=N k kMm mMN Nk kMm mdz cz Az z dz cA z H 111111)()()1()1()( (6.10)假设系统稳定，将ωj e z =代入上式，得到频率响应)](arg[11)()()()()(ωωωωωωj eH j j N k kjMm m j M N j j e e H de c e Ae e H =--=∏∏==- (6.11)设M N =，由式(6.11)得到∏∏==--=Nk kjNm m j j de c eA e H 11)()()(ωωω (6.12)在z 平面上，mj c e -ω用一根由零点mc 指向单位圆（ωj e ）上任一点B 的向量B c m 表示，同样kj d e -ω用一根由极点kd 指向单位圆（ωje ）上任一点B 的向量B d k 来表示，B c m 和Bd k 分别称为零点矢量和极点矢量，用极坐标表示为：m j m m e c B c α=，kj k k e d B d β=。
将B c m 和Bd k 表示式代入式(6.12)，得到)(1111)()()()(ωϕωωωωj j N k kNm mNk kjNm mj j e e H BdBcA de ce A e H ==--=∏∏∏∏====∏∏===Nk kNm mj dcA e H 11|)(|ω (6.13)∑∑==-=Nk kNm m 11)(βαωϕ (6.14)系统或者信号的频率特性由式(6.13)和式(6.14)确定。
按照式(6.13)，知道零极点的分布后，可以很容易地确定零极点位置对系统特性的影响。