文档介绍：. .
优选
实验一离散时间系统的时域分析
一、实验目的
１. 运用MATLAB仿真一些简单的离散时间系统，并研究它们的时域特性。
２. 运用MATLAB中的卷积运算计算系统的输出序列，加深x2(n)时的输出y2(n)
yn=filter(num,den,x,ic); %计算输入为x (n)时的输出y(n)%画出输出信号
subplot(2,2,1)
stem(n,y1);
ylabel('振幅');
title('y1输出');
subplot(2,2,2)
stem(n,y2);
ylabel('振幅');
title('y2输出');
subplot(2,2,3)
stem(n,yn);
ylabel('振幅');
title('yn输出');
实验2：
clf;
n=0:40;n1=0:50;
a=2;b=-3;
x1= cos(2*pi**n);
x2=cos(2*pi**n);
x3=a*x1+b*x2;
x4=[zeros(1,10), x3];
den=[1, -,];
num=[ ];
ic=[0 0]; %设置零初始条件
y3=filter(num,den,x3,ic);
y4=filter(num,den,x4,ic);%计算输入为x (n)时的输出y(n)
%画出输出信号
subplot(2,1,1)
stem(n,y3);
ylabel('振幅');
title('yn输出');
subplot(2,1,2)
stem(n1,y4);
ylabel('振幅');
title('y1输出');
. .
优选
实验3：
clf;
x=[1 3 2];%冲激
u=[1 1 1]; %输入序列
y=conv(u,x);
n=0:4;
stem(n,y);
xlabel('时间序号n');ylabel('振幅');
title('用卷积得到的输出');grid;
实验二〔1〕离散时间信号的DTFT
一、实验目的
１. 运用MATLAB理解Z变换及其绘制H(z)的零极点图。
２. 运用MATLAB计算逆Z变换。
二、实验原理
〔一〕、MATLAB在ZT中的应用。
线性时不变离散时间系统的冲激响应h(n)的z变换是其系统函数H(z)，在MATLAB中可以利用性质求解Z变换，例如可以利用线性卷积求的Z变换。假设H(z)的收敛域包含单位圆，即系统为稳定系统，即系统在单位圆上处计算的是系统的频率响应。
〔二〕、逆Z变换
Z变换对于分析和表示离散线性时不变系统具有重要作用。但是在MATLAB中不能直接计算Z变换，但是对于一些序列可以进展逆Z变换。
序列的Z变换及其收敛域，求序列称为逆Z变换。序列的Z变换及共逆Z变换表示如下：
通常，直接计算逆Z变换的方法有三种：围线积分法、长除法和局局部式展开法。在实际中，直接计算围线积分比较困难，往往不直接计算围线积分。由于序列的Z变换常为有理函数，因此采用局局部式展开法比较切合实际，它是将留数定律和常用序列的Z变换相结合的一种方法。
设x(n)的Z变换X(z)是有理函数，分母多项式是N阶，分子多项式是M阶，将X(z)展成一些简单的常用的局局部式之和，通过常用序列的Z变换求得各局部的逆变换，再相加即得到原序列x(n)。在MATLAB中提供了函数residuez来实现上述过程，调用格式如下：
[R，P，K]= residuez〔B，A〕
其中B、A分别是有理函数分子多项式的系数和分母多项式的系数，输出R是留数列向量，P是极点列向量。如果分子多项式的阶数大于分母多项式的阶数，那么K返回为常数项的系数。
三、实验容与步骤
选做一个实验：
1、.运行下面程序并显示它，验证离散时间傅立叶变换DTFT的时移性。
两个线性时不变的因果系统，系统函数分别为
，
分别令N=8，a=，计算并图示这两个系统的零、极点图及幅频特性。
. .
优选
程序：
2、运行下面程序并显示它，验证离散时间傅立叶变换DTFT的频移性。
四、实验仪器设备
计算机，MATLAB软件
五、实验本卷须知
课前预先阅读并理解实验程序；
六、实验结果
clear
num1=[1 0 0 0 0 0 0 0 -1];%分子系数高阶到低阶
den1=[1 0 0 0 0 0 0 0 0];
subplot(2,2,1)
zplane(num1, den1)
grid