实验十一 离散信号时域分析的MATLAB 实现
一、实验目的
1. 熟悉MATLAB 编程方法、常用语句和可视化绘图技术;
2. 掌握序列时域运算的MATLAB 编程方法。
3. 掌握离散信号时域分析的MATLAB 方法。
二、实验原理
1.离散时间系统
LTI 离散系统中,其输入和输出的关系由差分方程描述:
00()()n m
i j
i j a y k i b
f k j ==+=+∑∑ (前向差分方程)
()()n
m
i j
i j a y k i b
f k n j ==-=-+∑∑ (后向差分方程)
当系统的输入为单位序列δ(k )时产生的零状态响应称为系统的单位函数响应,用h (k )表示。当输入为 ε(k )时产生的零状态响应称为系统的单位阶跃应,记为:g (k ),如下图所示。
如果系统输入为e (k ),冲激响应为h (k ),系统的零状态响应为y(k ),则有:
()()()y k h k f k =*。与连续系统的单位冲激响应h (t )相类似,离散系统的单位函数响应h (k )
也包含了系统的固有特性,与输入序列无关。我们只要知道了系统的单位函数响应,即可求得系统在不同激励信号作用下产生的响应。因此,求解系统的单位函数响应h (k )对我们进行离散系统的分析也同样具有非常重要的意义。
MATLAB 中为用户提供了专门用于求解离散系统单位函数响应, 并绘制其时域波形的函数impz( )。同样也提供了求离散系统响应的专用函数filter( ),该函数能求出由差分方程所描述的离散系统在指定时间范围内的输入序列作用时,产生的响应序列的数值解。当系统初值不为零时,可以使用dlsim( )函数求出离散系统的全响应,其调用方法与前面连续系统的lsim( )函数相似。另外,求解离散系统阶跃响应可以通过如下两种方法实现:一种是直接调用专用函数dstep( ),其调用方法与求解连续系统阶跃响应的专用函数step( )的调用方法相似;另一种方法是利用求解离散系统零状态响应的专用函数filter( ),只要将其中的激励信号看成是单位阶跃信号ε(k )即可。 函数的调用格式分别如下: ⑴ impz( )函数
impz(b,a) 以默认方式绘出由向量a 和b 所定义的离散系统单位函数响应的时域波形。