资源描述:
《数字信号处理实验二().doc》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、实验二离散时间系统的时域分析实验室名称:实验时间:姓名:学号:专业:指导教师:成绩教师签名:年月日一、实验目的1.通过MATLAB仿真一些简单的离散时间系统,并研究它们的时域特性。2.通过MATLAB仿真滤波器,研究它们的性质。3.用MATLAB仿真研究离散时间系统的性质。二、实验内容Q2.1M=2运行内容一的程序代码,生成输入x[n]=s1[n]+s2[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,用一个长度为
2、101、最低频率为0、最高频率为0.5的扫频正弦信号作为输入信号(见程序P1.7),计算其输出信号。Q2.7运行程序P2.3,对由加权输入得到的y[n]与在相同权系数下输出y1[n]和y2[n]相加得到的yt[n]进行比较,这两个序列是否相等?该系统是线性系统吗?Q2.12运行程序P2.4并比较输出序列y[n]和yd[n-10]。Q2.19运行程序P2.5,生成式(2.15)所给离散时间系统的冲击响应。Q2.20修改程序P2.5产生如下因果线性时不变系统的冲激响应的前45个样本:y[n]+0.71y[n-1]-0.46y[n-2]-0.62y[n-3]=0.9x[n]-0
3、.45x[n-1]+0.35x[n-2]+0.002x[n-3]。Q2.23运行程序P2.6,计算输出序列y[n]和y2[n]以及差值信号d[n]。y[n]和y2[n]相等吗?Q2.28运行程序P2.7,对序列h[n]和x[n]求卷积,生成y[n],并用FIR滤波器h[n]对输入x[n]滤波,求得y1[n]。三、实验器材及软件1.微型计算机1台2.MATLAB7.0软件四、实验原理1.离散时间系统中信号运算方法和基本性质;2.离散时间系统的MATLAB表示和基本运算;3.MATLAB处理离散时间系统的方法。五、实验步骤1.根据实验题目要求进行分析;2.运用所学知识用MAT
4、LAB编程实现题目要求;3.对结果进行分析总结。六、实验记录(数据、图表、波形、程序等)Q2.1程序代码:%产生输入信号n=0:100;s1=cos(2*pi*0.05*n);%一个低频正弦s2=cos(2*pi*0.47*n);%一个高频正弦x=s1+s2;%滑动平均滤波器的实现M=input('Desiredlengthofthefilter=');num=ones(1,M);y=filter(num,1,x)/M;clf;%显示输入和输出信号subplot(2,2,1);plot(n,s1);axis([0,100,-2,2]);xlabel('Timeindexn
5、');ylabel('Amplitude');title('Signal#1');subplot(2,2,2);plot(n,s2);axis([0,100,-2,2]);xlabel('Timeindexn');ylabel('Amplitude');title('Signal#2');subplot(2,2,3);plot(n,x);axis([0,100,-2,2]);xlabel('Timeindexn');ylabel('Amplitude');title('InputSignal');subplot(2,2,4);plot(n,y);axis([0,100,-
6、2,2]);xlabel('Timeindexn');ylabel('Amplitude');title('OutputSignal');axis;显示结果:Q2.2代码(将输出信号的代码做出修改,其余代码与Q2.1相同):n=0:100;s1=cos(2*pi*0.05*n);s2=cos(2*pi*0.47*n);x=s1+s2;M=input('Desiredlengthofthefilter=');num=(-1).^[0:M-1];y=filter(num,1,x)/M;clf;%显示输入和输出信号subplot(2,2,1);plot(n,s1);axis([
7、0,100,-2,2]);xlabel('Timeindexn');ylabel('Amplitude');title('Signal#1');subplot(2,2,2);plot(n,s2);axis([0,100,-2,2]);xlabel('Timeindexn');ylabel('Amplitude');title('Signal#2');subplot(2,2,3);plot(n,x);axis([0,100,-2,2]);xlabel('Timeindexn');ylabel('Amplitude');title