资源描述：
《离散时间信号的分析》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、武汉工程大学实验报告实验课程数字信号处理姓名张鹏学号0704140228专业及班级07通信02班实验地点实验组号实验日期篡•月」9••旦实验项目•…―离散时间信号的分析实验目的认识常用的各种信号，理解其数字表达式和波形表示掌握在计算机中生成及绘制数字信号波形的方法—掌握序列的简单运算及计算机实现与作用理解J®散时JW里吐变軽…Z…变换及宣们的性质.和值号的.频域分析实验原理，实验步骤，实验仪器设备（名称，型号，功能，量程，在本次试验中的用途）二、实验设备计算机，MATLAB语言环境。三、实验基础理论1•序列的相关概念2•常见序列（哪几个？）3•序列的基本运算（哪些?4•离散傅
2、里叶变换的相关概念5.Z变换的相关概念（4、5的关系？）U1实验内容与步骤1•离散时间信号（序列）的产生利用MATLAB语言编程产生和绘制单位样值信号、单位阶跃序列、指数序列、正弦序列及随机离散信号的波形表示。2•序列的运算（1）利用MATLAB语言编程实现信号平滑运算。（2）利用MATLAB语言编程实现信号的调制。利用MATLAB语言编程实现信号卷积运算。（4）利用MATLAB语言编程实现信号离散傅立叶的正反变换。利用MATLAB语言编程实现信号的圆周移位、圆周卷积,验证DFT的圆周时移、圆周卷积性质和圆周卷积与线性卷积的关系。验证一个周期实序列奇偶部分的DFT与此序列本身
3、的DFT之间的关系。利用MATLAB语言编程实现信号的Z变换及其反变换、Z变换的零、极点分布。1.单位值信号程序x=2;Y=l;stem(x,y);titleC单位值J2.单位阶跃序列程序nO=O;nl=-10;n2=10;n=[nl:n2];x=[(n-nO)>=O];stem(nzx);指数序列n=[0:10];x=(0.5)•^n;stem(n,x)xlabel(1n1);ylabel(1x(n)1);title(1实指数序列「实指数序列正弦序列n=[-50:50];x=2*sin(0・05*pi*n);stem(nfx);xlabel(1n1);ylabel(1x(n
4、)1);title(1正弓玄序歹U1);n随机序列程序n=[l:10];x=rand(1,10);subplot(221)stem(n,x);xlabel(1n1);ylabel(1x(n)!)/title(1随机序歹U1)隨机序列平滑运算R=51;d=0.8*(rand(R,1)-0.5);m二0:RT;s=2*m.*(0.9.F);x二s+d';subplot(2,1,1);plot(m,d',，i•-'，m,s,，g一，，m,x,，b~.，);xlabelTimeindexn');ylabel('Amp1itude');legendCd[nj，，‘s[n]','x[nj
5、，);xl=[00x];x2=[0x0];x3二[x00];y二(xl+x2+x3)/3;subplot(2,1,2);plot(m,y(2:R+1),',m,s,'g一-');legend('y[n]'，‘s[n]');xlabelCTimeindexn);ylabel(*Amplitude，)d[n]s[n]刈n]505■Q)pn=-dluv°051015202S3035404550Timeindexn8642apnu-dluv两个序列相加减函数function[yzn]=sigadd(xl,nlzx2,n2)n=min(min(nl),min(n2)):max(max(
6、nl),max(n2));yl=zeros(1,length(n));y2=yl;yl(find((n>=min(nl))&(n<=m&x(nl))==1))=xl;y2(find((n>=min(n2))&(n<=max(n2))==1))=x2;y=yl+y2;两个序列相乘函数function[y,n]=sigmult(xl,nl,x2,n2)n=min(min(nl),min(n2)):max(max(nl),max(n2));yl=zeros(1,length(n));y2=yl;yl(find((n>=min(nl))&(n<=mmx(nl))==1))=xl;y2
7、(find((n>=min(n2))&(n2))==1))=x2;y=yi.*y2；两序列相乘和相加程序clc;clear;xl=[0zl,2z3z4z3z2rl,0];nl=-2:6;x2=[2f2f0f0f0f-2f-2];n2=2:8;[yl,n]=sigmult(xl,nl?x2‘n2);[y2,n]=sigadd(xl,nl,x2,n2);subplot(2,2,1);stem(nl,xl)/title(fDdABxl1)xlabel(1n1);ylabel(1xl(n)!);subplot(
