资源描述:
《通信与数字信号处理实验指导书(11电信)》由会员上传分享,免费在线阅读,更多相关内容在工程资料-天天文库。
1、福建狡林丈席东方修院DONGFANGCOLLEGE,FUJIANAGRICULTUREANDFORESTRYUNIVERSITY《通信与數芳信号处理》实验指專不热创卿偏写福建农林丈修京方修院针尊机科曇系(部丿2013弄®月实验一常用序列的MATLAB表示及其基本运算1实验二利用MATLAB进行时域离散系统的时域分析6实验三利用MATLAB进行时域离散信号和系统的频域分析8实验四IIR数字滤波器的设计10实验五FIR数字滤波器的设计13实验一.常用序列的MATLAB表示及其基本运算一、实验目的:熟悉MATLAB软件的使用,掌握MA
2、TLAB的程序设计方法;掌握数字信号处理小儿种常用序列的表达式和波形、序列的儿种基本运算;学会用MATLAB画出序列的波形图和实现序列基本运算的方法。二、实验仪器实验的硬件环境是:微机一•台实验的软件环境是:MATLAB软件三、实验原理MATLAB软件屮捉供了很多的内部函数,可川于实现通信和数字信号处理的相关功能。本实验屮用到的内部函数主要有:(1)zeros(m,n):用于产生一个m行n列的0矩阵。(2)ones(m,n):用于产生一个m彳亍n列的1矩阵。(3)stem(n,xn):用于绘制离散序列xn随n变化的波形图。(4)
3、subplot(m,n,p):用于画子图,m表示子图排成几行,n表示子图排成儿列,p表示是从左到右第几个位置。(5)xlabeKlexf):用于给出X轴的坐标名text。(6)ylabelCtext*):用于给出Y轴的坐标名text。(7)title('texf):用于给所绘制的波形取一个波形名text。(8)linspace(xl,x2,N):用于在xl和x2之间进行N等分。(9)conv(a,b):用于对序列a和b进行卷积和运算。(10)length(x):用于给出矢量x的长度。卩q、实验内容1、利用MATLAB软件分别画出
4、单位釆样序列》(〃)、单位阶跃序列u(n)、矩形序列RN(n).实指数序HJT列anu(n)和正弦序列sin(——n)的波形。画出各序列波形的MATLAB参考程序分别如下。b(1)单位采样序列<5(77)n=-50:50;stem(n,xn)%该函数用于绘制离散序列的波形xlabelCn*)ylabelCdelta(n)^)(2)单位阶跃序列u(n)n=-50:50;xn=[zeros(1,50),ones(l,51)];%函数ones(1,n)用于产生n个1stem(n,xn)xlabel('n')ylabel('u(n)')
5、(3)矩形序列RN(n)N=4n=-20:20;xn=[zeros(l,20),ones(l,N),zeros(1,21-N)];stem(n,xn)xlabel('n')ylabel(4RN(n)')(4)实指数序列anu(n)a=0.6;n=0:10;xn=a.An;stem(n,xn)xlabel('n')(5)正弦序列sin(—n)ba=2;b=5;n=linspace(0,2*pi,60);%60个点的n坐标,将(0,2pi)分成60等分xn=sin(a*pi*n/b);stem(n,xn)xlabelCn')2、利用
6、MATLAB软件实现序列的相加、相乘和卷积和运算,要求分別価出运算前后各序列的波形。实现这三种运算的MATLAB参考程序分别如下。(1)相加n=linspace(O,2*pi,60);xl=sin(pi*n/3);x2=sin(2*pi*n/3);xn=xl+x2;subplot(3,l,l)stem(n,xl)xlabel('n')ylabelCxl(nnsubplot(3,l,2)stem(n,x2)xlabel('n')ylabel('x2(n)')subplot(3,l,3)stem(n,xn)xlabel('n')yl
7、abel('x(n)')(2)相乘n=linspace(0,2*pi,60);xl=sin(pi*n/3);x2=sin(2*pi*n/3);xn=xl.*x2;subplot(3,lJ)stem(n,xl)xlabel('n1)ylabel(,xl(n)^)subplot(3,l,2)stem(n,x2)xlabel('n')ylabel(,x2(n)^)subplot(3,l,3)stem(n,xn)xlabel('n')ylabel('x(n)')(3)卷积和n=0:4;xl二[1,2,3,4,5];subplot(3,l
8、,l)stem(n,xl)xlabel('n')ylabel('x1(n)')n=0:3;x2=[1,1,1,1];subplot(3,l,2)stem(n,x2)xlabel('n')ylabcl('x2(n)‘)y=conv(xl,x2);N=length(