离散时间系统的时域分析

课程名称：数字信号处理实验

任课教师： 周浩

实验名称：离散时间系统的时域分析

年级、专业：07级通信工程

学  号：20071060110

姓  名： 刘灵芝

日期：2009年10月16日

                       云南大学 信息学院

一、实验目的

1. 学会使用MATLAB仿真简单的离散时间系统，并研究这些系统的时域特性。

2. 学会对实验结果进行分析比较。

二、实验内容

1. 离散时间系统的仿真。

2. 线性时不变离散时间系统的冲击响应的计算和卷积。

习题：Q2.1   Q2.2  Q2.3      Q2.7  Q2.8  Q2.9   Q2.19~Q2.22  Q2.28  Q2.29

三、主要算法与程序

Q2.1:对ｍ＝２，运行程序。输入ｘ［ｎ］的哪各分量被抑制？

x[n]的s2部分被抑制，即高频部分被滤除。

Ｑ２.２:

Q.2.3 :对滤波器长度m和正弦信号s1和s2的频率取其他值。M=2，s1=0.01zH，s2=0.47zH

Q2.4:S输入扫频信号

Q.2.7

这两个序列是相等的，该系统是线性的

Q2.8

a=3,b=4 x1=100zh  x2=250zh

Q2.9

当初始条件非零时，这两个序列是相等的，该系统是线性的。

Q2.19运行程序P2.5，生成离散时间系统的冲击响应。

Q2.20

Q2.21

Q2.22

Q2.28

clf;

h = [3 2 1 -2 1 0 -4 0 3];    

x = [1 -2 3 -4 3 2 1];        

y = conv(h,x);

n = 0:14;

subplot(2,1,1);

stem(n,y);

xlabel('Time index n'); ylabel('Amplitude');

title('Output Obtained by Convolution'); grid;

x1 = [x zeros(1,8)];

y1 = filter(h,1,x1);

subplot(2,1,2);

stem(n,y1);

xlabel('Time index n'); ylabel('Amplitude');

title('Output Generated by Filtering'); grid;  

Q2.29 修改程序P2.7，重做问题Q2.28

h = [ 2.2403    3.3869    1.9148   -1.7743   -2.1458    0.4724 1.7983 0.3651   -1.2027   -0.7549    0.6001    0.8062 -0.1276   -0.6557 0.1666 ] 

x = [1 -2 3 -4 3 2 1 4 5 6 ];    

y = conv(h,x);

n = 0:24;

subplot(2,1,1);

stem(y);

xlabel('Time index n'); ylabel('Amplitude');

title('Output Obtained by Convolution'); grid;

x1=[x zeros(1,14)];

y1=filter(h,1,x1);

subplot(2,1,2);

stem(y1);

xlabel('Time index n'); ylabel('Amplitude');

title('Output Obtained by Convolution');

grid;

四、实验小结

在三点滑动平均系统中当y[n]=0.5(x[n]+x[n-1])变为y[n]=0.5(x[n]-x[n-1])时是把信号翻转后再叠加。Fiter命令和impz命令都能求出系统的冲击响应，用fiter命令时把输入函数改为u（t）。卷积后序列的长度为h=a+b-1。