《数字信号处理》上机实验指导.docVIP

《数字信号处理》上机实验指导 实验一、Z变换及离散时间系统分析 （一）、实验目的 1、通过本实验熟悉Z变换在离散时间系统分析中的地位和作用。 2、掌握并熟练使用有关离散系统分析的MATLAB调用函数及格式，以深入理解离散时间系统的频率特性。 （二）、实验内容及步骤 对于一个给定的LSI系统，其转移函数H(z)习惯被定义为H(z)=B(z)/A(z),即： 公式中和分别是H(Z)分子与分母多项式的阶次，在有关MATLAB 的系统分析的文件中，分子和分母的系数被定义为向量，即 并要求=1,如果≠1，则程序将自动的将其归一化为1。 1、系统的阶跃响应 调用格式为： y=filter(b,a,x),其中x,y,a,b都是向量。 例1 令 求该系统的阶跃响应（y（n））。实现该任务的程序如下： clear; x=ones(100);% x(n)=1,n=1~100； t=1:100;% t用于后面的绘图； b=[.001836，.007344，.011016，.007374，.001836]; % 形成向量b； a=[1，-3.0544，3.8291，-2.2925，.55075]; % 形成向量a； y=filter(b,a,x); % 求所给系统的输出，本例实际上是求所给系统的阶跃响应； plot(t,x,r.,t,y,k-);grid on;% 将x(n)（绿色）y(n)（