资源描述：
《数字信号处理实验二》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、WORD格式可编辑实验二离散时间系统的时域分析实验室名称：实验时间：姓名：学号：专业：指导教师：成绩教师签名：年月日一、实验目的1.利用MATLAB仿真简单的离散时间系统，研究其时域特性；2.对线性时不变系统进行重点分析研究，掌握其特性。二、实验内容2.1对M=2，运行上述程序，生成输入x[n]=s1[n]+s2[n]的输出信号。输入x[n]的哪个分量被该离散时间系统抑制？2.2若线性时不变系统由y[n]=0.5(x[n]+x[n-1])变成y[n]=0.5(x[n]-x[n-1])，对输入x[n]=s1[
2、n]+s2[n]的影响是什么？2.4修改程序P2.1，用一个长度为101、最低频率为0、最高频率为0.5的扫频正弦信号作为输入信号（见程序P1.7），计算其输出信号。你能用该系统对扫频信号的响应来解释习题Q2.1和习题Q2.2的结果吗？2.7运行程序P2.3，对由加权输入得到的y[n]与在相同权系数下输出y1[n]和y2[n]相加得到的yt[n]进行比较，这两个序列是否相等？该系统是线性系统吗？2.12运行程序P2.4并比较输出序列y[n]和yd[n-10]。这两个系列之间有什么关系？该系统是时不变系统吗？
3、2.19运行程序P2.5，生成式（2.15）所给离散时间系统的冲激响应。2.20修改程序P2.5，产生如下因果线性时不变系统的冲激响应的前45个样本：2.23运行程序P2.6，计算输出序列y[n]和y2[n]以及差值信号d[n]。Y[n]和y2[n]相等吗？2.28运行程序P2.7，对序列h[n]和x[n]求卷积，生成y[n]，并用滤波器h[n]对输入x[n]滤波，求得y1[n]。y[n]和y1[n]有差别吗？为什么要使用对x[n]补零后得到的x1[n]作为输入来产生y1[n]？三、实验器材及软件1.微型计
4、算机1台2.MATLAB7.0软件专业知识整理分享WORD格式可编辑四、实验原理1.三点平滑滤波器是一个线性时不变的有限冲激响应系统，将输出延时一个抽样周期，可得到三点平滑滤波器的因果表达式，生成的滤波器表示为归纳上式可得此式表示了一个因果M点平滑FIR滤波器。2.对线性离散时间系统，若y1[n]和y2[n]分别是输入序列x1[n]和x2[n]的响应，则输入的输出响应为则系统称为线性系统。3.对于离散时不变系统，若y1[n]是x1[n]的响应，则输入x[n]=x1[n-n0]的输出响应为y[n]=y1[n-
5、n0]则称系统为时不变系统。五、实验步骤2.1首先利用MATLAB产生一个高频正弦信号和一个低频正弦信号，利用两个信号生成一个输入信号，接着利用filter函数生成输出信号，最后利用plot函数画出4个信号，再对输出信号进行分析。2.2在2.1的基础上编写num=ones[1-1]，运行程序得出结论。2.4分别用扫频信号通过2.1、2.2的系统，进行比较分析。2.7分别计算出y1[n]和y2[n]，得到yt[n]；再利用filter函数求得y[n]，计算差值输出，比较y[n]和yt[n]。2.20根据imp
6、z函数的调用方式，得到num=[0.9-0.450.350.002]den=[10.71-0.46-0.62]，再调用impz函数，画出图像。2.23首先产生序列x[n]，把它作为四阶系统的输入，生成y[n]。然后将同样的输入x[n]应用到第一级得到y1[n]。接着用相同的方法得到y2[n]。最后求得两者的差，并画出图像。2.28分别用conv函数和filter函数求得输出，进行图像比较。六、实验记录（数据、图表、波形、程序等）专业知识整理分享WORD格式可编辑2.1对M=2，运行上述程序，生成输入x[n]
7、=s1[n]+s2[n]的输出信号。输入x[n]的哪个分量被该离散时间系统抑制？%ProgramP2_1clf;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');s
8、ubplot(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('时间序号n');ylabel('振幅');title('输入信号');subplot(2,2,4);plot(n,y);axis([0,1
