matlab离散系统差分求全响应,实验四 离散LTI系统的时域和z域分析

实验四离散LTI 系统的时域和z 域分析

【实验目的】

1. 掌握利用MATLAB 计算离散系统响应的数值方法，包括冲激响应、全响应等。

2. 掌握离散信号z 变换和逆z 变换的MATLAB 实现方法；

3. 掌握离散系统的系统函数零极点分布与系统频率特性分析的MATLAB 实现方法。【实验原理】

1..单位序列δ(k)

单位序列的定义：

下面为绘制δ(k-k0)波形图的子程序：

function impseq(k1,k2,k0)

k=k1:k2; %k1,k2 为时间序列的起始及终止时间序号

fk=[(k-k0)==0]; %k0 为单位序列在时间轴上的位移量

stem(k,fk)

axis([k1,k2,0,1.1])

xlabel('k')

title('单位序列')

输入命令impseq(-1,5,3)，则可获得单位序列δ(k-3)的波形图，如图1 所示。

2..单位阶跃序列ε(k)

单位阶跃序列的定义：

下面为绘制ε(k-k0)波形图的MATLAB 子程序。

function stepseq(k1,k2,k0)

k=k1:k2; %k1,k2 为时间序列的起始及终止时间序号

fk=[(k-k0)>=0]; %k0 为阶跃序列在时间轴上的位移量

stem(k,fk)

axis([k1,k2,0,1.1])

xlabel('k')

title('单位阶跃序列')