数字信号处理实验2解说.doc

实验二 离散时间系统的时域分析 实验室名称：计算机基础实验室（信息学院2202） 实验时间： 2015 年 9 月 24 日 姓 名： 王凤琼 学号：20131060114 专业： 电子信息工程 指导教师：柏正尧 成绩 教师签名： 年 月 日 一、实验目的 1. 通过MATLAB仿真一些简单的离散时间系统，并研究它们的时域特性。 2. 了解滑动平均系统的特性，仿真滑动平均系统，与理论课相结合，分析系统特性。 3. 运用MATLAB程序得出线性时不变离散时间系统的冲激响应。 4. 熟悉相应的MATLAB命令，如filter，impz。 5. 实验与理论相结合，能达到对理论的更好理解。 二、实验内容 滑动平均系统的仿真； 线性和非线性系统的仿真； 线性时不变离散时间系统的仿真； Q2.1对 M=2，运行上述程序，生成输入 x[n]=s1[n]+s2[n]的输出信 号。输入 x[n]的哪个分量被该离散时间系统抑制？ Q2.2若 线 性 时 不 变 系 统 由 y[n]=0.5(x[n]+x[n-1]) 变 成 y[n]=0.5(x[n]-x[n-1])，对输入 x[n]= s1[n]+s2[n]有什么影响？ Q2.4修改程序P2.1，用一个长度为101，最低频率为0，最高频率为0.5的扫频正弦信号作为输入信号，计算其输出信号。 Q2.7运行程序P2.3，对由加权输入得到的y[n]与在相同权系数下输出y1[n]和 y2[n]相加得到的yt[n]进行比较，这两个序列是否相等？该系统是线性系统吗？ Q2.12采用三个不同的延时变量D的值运行程序P2.4并比较输出序列y[n]和yd[n-10]。这两个序列之间有什么关系？该系统是时不变系统吗？ Q2.19运行程序 P2.5，生成2.15所给离散时间系统的冲激响应。 y[n]-0.4y[n-1]+0.75y[n-2]=2.2403x[n]+2.4908x[n-1]+2.2403x[n-2] Q2.20修改程序 P2.5,产生如下因果线性时不变系统的冲激响应的前 45 个样本： y[n]+0.71y[n-1]-0.46y[n-2]-0.62y[n-3]=0.9x[n]-0.45x[n-1]+0.35x[n2]+0.002x[n-3] Q2.23运行程序p2.6，计算输出序列y[n]和y2[n]以及 差值信号d[n],y[n]和y2[n]相等 Q2.28运行程序p2.7，对序列h[n]和y[n]求卷积，生成y[n],并用FIR 滤波器h[n]对输入x[n]滤波，求得y1[n]。y[n]和y1[n]有差别吗？为什么要使用对x[n]补零后得到的x1[n]作为输入来产生y1[n]？ 计算线性时不变离散时间系统的冲激响应，运用MATLAB命令y=impz(num，den，N)，得出冲激响应的前N个样本。 三、实验器材及软件 1. 微型计算机1台 2. MATLAB 7.0软件 四、实验原理 1.若y1[n]和y2[n]分别是因果离散时间系统输入信号u1[n]和u2[n]的响应，则当 u1[n]=u2[n] nN时，有y1[n]=y2[n] nN 2.对于离散时不变系统，若y1[n]是x1[n]的响应，则输入x[n]=x1[n=n0]的输出响应为y[n]=y1[n-n0]。 3.描述因果线性时不变离散时间系统，可用命令 filter进行仿真。 计算冲激响应可用函数impz表示计算冲击响应钱N个样本可用y=impz(num,den,N)来表示。 五、实验步骤 按照课本上的例题书写程序。 六 实验记录（数据、图表、波形、程序等） Q2.1 clf; n=0:100; s1=cos(2*pi*0.05*n); s2=cos(2*pi*0.47*n); x=s1+s2; M=input(滤波器所需的长度=); num=ones(1,M); y=filter(num,1,x)/M; subplot(2,2,1); plot(n,s1); axis([0,100,-2,2]); xlabel(时间序号n);ylabel(振幅); title(信号#1); subplot(2,2,2); plot(n,s2); axis([0,100,-2,2]); xlabel(时间序号n);ylabel(振幅); title(信号#2); subplot(2,2,3); plot(n,x); axis([0,100,-2,2]); xlabel