《数字信号处理》实验报告汇总.docx
HuaiyinNormalIlniversity物理与电子电气工程学院试验报告课程名称：数字信号处理院系：物电学院专业：电子信息科学与技术班级：学号：姓名：试验报告(D试验名称常见高散信号产生与实现试验日期2016年9月13日指导老师普凤莲一、试验目的1 .驾驭MAT1.AB语言的班本操作，学习基本的编程功能。2 .驾驭几种基本的离散时间信号(包括单位抽样序列,单位阶跃序列,单频正弦序列.单版复指数序列，实指数序列等).3 .能弱娴热利用MAT1.AB产生这些基本的离故时间信号。二、Sl习要求1,熟识各种常用序列，驾驭序列matlab实现。4 .侦习网T1.AB中序列产生的词用函数及绘制图形函数stc三()三、管验内容利用YATIAB编程产生和绘制下列有限长序列：(I)单位脓冲序列3()(2)单位阶跃序列”()(3)矩形序列45)(4)正弦型序列M)=Asin(qn+q)(5)随意序列x(n)=6()+2J(-1)+第(-2)+4J(?-3)+55(-4)/:()=8()+25(-1)+6(-2)+2Z>(n-3)四、祓3报告1.实现上述各序列(I)单位冲激序列functiony=impDT(n)%定义inpDT函数y=(n=0):n=-3:3;x=inpDT(n);stcm(n.x：fiir):%离散时间信号的波形绘制在MAT1.AB中一殷用S(Cm函数xlabel(,n'):gridon;Ele(即位冲激序列姓名学号')axis(33-0.11.1)(2)单位阶跃序列functiony=uDT(n>y=(n>=0);n=-3:5；x=uDT(n);slcm(n.'fil);xlabel(*n,);gridonEle('单位阶跃序列姓名学号')axis(35-0.1l.l)(3)矩形序列&()functiony=uDT(n)y=(n>=O);n=-5:IO:x=uDT(n>-uDT(n-8):%此处调用了uDT.m函数stcm(n,xfi);xlabeC11,);gridonaxis(-51()-0.11.1J)IiUC('矩形序列姓名学号')(4)正弦里序列=ASin(三+工)53n=20r20;x=sin(pi5*n÷pW);s<em(n.'fil);Xlabclfnt):gridonlhle('正弦序列姓名学号)axis(-2020-1.51.5)(三)随意序列1：x(n)=(n)+2(n-1)+3(n-2)+4(n-3)+5J(-4)functionx.n=impscq(n.nl.n2)if(n<nI)(n>n2)(nI>n2)errrf参数必需涓痣nl<=n<=n2,)endn=(nkn2J;x=(n-nO)=O;%以上,定义inpcq.m解数n=IOJO;x=impseq(0.-l0.10)+2*inpseq(1.-10.10H3巧mpstq(2.-10.IO)+4*irnpscq(3.-10.10)+5*impseq(4.-10.10);s(e11M11.fillr);xlabel(an>>gridonCitIe(,随意序列1姓名学号)随意序列2：h(r)=b(")+2(w-l)+5(-2)+25(-3)n=-IOJO:h-impscq(O,-O,I()÷2*impscq(l,-OJO)÷impscq(2,-l()JO)÷2*impscq(3t-IOJO);%同样的调用5.1中的inpseq.m函数stcm(n.h,fil);xlabcl(n,);ylabcCh,);gridonthle('随怠序列2姓名学号)axis(-10IO-I3|)2.画出各序列的图形，并对结果进行分析。单位冲激序列李杰171413275111QU.o11CU.OAU-*O011-5r)3单位阶跃序列李杰171413275正弦序列李杰171413275任意序列2李杰1714132752.521.5U10.500.5F-i1Q，.i-1-103-6-4-20246810n分析总结，MatIab的函数定义要放在M文件中，弁以文件名与函数名一样才能运行.此次试脸大致巩固了MIalab基本的语法与表达式。同时，我了解了单位阶跃序列、单位冲激序列、矩形序列、正弦序列的实现方法,并且通过他们的绘图广解这函数.从代码到图像.洋细形象地理解/这些信号的特征.试验报告(2)忒验名称离散时间系统的时域分析试验日期2016年9月20日指导老师曹凤莲一、试验目的1 .学会运用时T1.RB求解及散时间系统的零状态响应.2 .学会运用见V1.B求解禹散时间系统的单位取样响应.3 .学会运用MAT1.AB求好康散时间系统的卷积和.二、Sl习要求1 .预习试验中基础学问,熟识MAT1.AB指令及y=fHtcNbax),ImPZ(b.a.N)、y=conv(x.h>函数.2 .结合试验内容，提前编制相应的程序。3 .思索变更差分方程的形式，单位抽样响应将如何变更。三、管验内容I.试用MAT1.AB吩咐求解以下圉散时间系统的单位取样响应. 1) I)3>(J)+4如-1)+y(n-2)=x(n)+<11-I) 2) y(11)+61y(/:-1)+10y(11-2)=x(r)2.已知某系统的单位取样响应为叫)=()”“()-“(-10),试用MAT1.AB求当激励信号为8.v()=w(11)-w(n-5)Bt.系统的零状态响应.四、忒3报告1.依据内容要求，写出调试好的MT1.H诺言程序及对应的耿T1.AB运算结果.1.试用MNr1.AB吩咐求解以下离散时间系统的单位取抑响应.%在m文件中定义两个函数以便调用funcliony=impDT(n);y=(n=O);functiony=DT(n);y=(n>=O(1)3.v(")+4y(j-l)+y(n-2)=x(n)+x(n-1)a=34l;b=(lU；n=O:IO;x=inDT(n);h=fiher(b.3Xs(em(n.h.'fi).gridon;xlabcl(*n,):Iit1«,系统单位取样响应姓名学号)系统单位取样响应李杰171413275U.oO0.3U.Q0.2fU.IuC40.1U.Uz)AiIU-u.Ut>114u.I。13)3I5678910n(2)Jy(n)+6y(i-l)+lOy(j-2)=X(M)a=(2.56101；Mlkn=0:50;x=impDT(n);hfiltcrtb,a,x);stem(nJ)Jfil).gridon:xlabel('n'):IitICC系统单位取样响应姓名学号')x10系统单位取择响应李杰1714132752.已知某系统的单位取样响应为M)=(N)"M")i(一IO)1.试用MAT1.AB求当激励信号为x(n)=()(-5)时,系统的零状态响应，nh=-IO:2O:nx=-5:5;x=uDT(nx)-uDT(nx-5);h=(778).Anh.*(uDT(nh)uDT(nh-l();y=conv(x.h);nyl=nx(l)+nh(l);ny=nyI+():(lcngth(nx)+lcngtMnh)-2);suhplo<(3II);s(em(nx.x.'fil);gridon:xlabcl('n,):mIC(Wn)姓名学号)axis(-511-0.1l.l);subpkx(312):stcm(nh.h.,fil):gridon;xlabel(,n,);tide(,h(n)姓名学号)axis(-511-0.1I.1J);SUbPkX(313);stem(ny,y.,fil);gridon;xlabcCn'):tit1cCy(n)=x(n)*h(n)姓名学号');axis(-518-15):Mn)李杰171413275-50510分析总结：函数的调用。编辑好M文件后调用单位冲激序列和堆位阶跃序列，并利用这些荷洁序列的调用描述其他的序列，从而简化工作。离散时间信号的波形绘制在MT1.B中一般用Slem函数，由于MAT1.AB中矩阵元素的个数有限，所以MAT1.AB只能表示肯定时间范围内有限长度的序列；而对于无限序列，也只能在肯定时间范围内表示出来。试验报告(3)试验名称离散时间1.Tl系统的z域分析试验日期2016年9月27日指导老师普凤莲一、试验目的1 .学会运用时T1.RB分析离散时间系统的系统函数的零极点。2 .学会运用WVr1.AB分析系统函数的零极点分布与其时域特性的关系.3 .学会运用MAT1.AB进行疡散时间系统的领率特性分析.二、5(习央求1 .预习试验中基础学问,熟识MAT5B指令及freqz()»abs().angleO函数.2 .结合试验内容，提前编制相应的程序，3,熟识系统侵定的充要条件及学过的梗定判据。三、管验内容1 .试用MAT1.AB柯出下列因果系统的系统函数零极点分布图，并推断系统的稳定性。(1) H(Z)=-r-2T1609zj-2.52+1.96-0.48(2) H(z)=-r,z4-0.9za-0.65z2+0,873z2 .试用MAT1.AB绘制系统H(Z)=r的频率响应曲线。>3IZ2十一
