离散傅立叶变换dft实验报告.doc

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