数字信号处理A实验指导书全解.doc

数字信号处理（A） 实验指导书 柏逢明 吉林城建学院 2012年10月20日 实验一 离散时间信号分析 实验目的：加深对常用离散信号的理解，掌握序列运算和序列构造方法。 实验原理：连续信号采样获得离散时间序列，对序列进行运算和构造典型序列过程。 三、实验内容： 信号是信息的载体，它可以传递不同的信息。离散时间信号是时间变量离散化的信号。它是对连续时间信号经过抽样而获得的。经过采样得到的多个离散点的有序排列的信号可称为序列x(n)。 （一）、 序列的运算 序列相加与相乘： 已知序列 ，分别求两序列相加及相乘 % sequence1.m w0=pi/20;a=1.05; n1=[-20:20];n2=[-20:20]; subplot(2,2,1),stem(n1,sin(w0.*n1),.k);axis([-20,20,-1,3]); text(-10,2,x1(n)); xlabel(n) subplot(2,2,2),stem(n2,a.^n2,.k);axis([-20,20,-1,3]); text(-10,2,x2(n)); xlabel(n) n=[min(min(n1),min(n2)):max(max(n1),max(n2))]; x1=zeros(1,length(n)); x2=zeros(1,length(n)); x1([find((n=min(n1))(n=max(n1)))])=sin(w0.*n1); x2([find((n=min(n2))(n=max(n2)))])=a.^n2; y1=x1+x2; y2=x1.*x2; subplot(2,2,3);stem(n,y1,.k);axis([-20,20,-1,3]); text(-10,2,y1(n)); xlabel(n) subplot(2,2,4);stem(n,y2,.k);axis([-20,20,-1,3]); text(-10,2,y2(n)); xlabel(n) 2.序列移位（延迟）： 设x(n)＝｛1,2,3,4,5,5,4,3,2,1｝且n＝｛0,1,2,3,4,5,6,7,8,9｝，求出序列x(n)及其移位x(n-2)、x(n+2)的图形。 % sequence2.m n1=[0:9]; x=[1,2,3,4,5,5,4,3,2,1]; subplot(3,1,1);stem(n1,x,.b),axis([-3,12,0,5]) text(0,3,x(n)) n2=n1-2; subplot(3,1,2);stem(n2,x,.b),axis([-3,12,0,5]) text(9,3,x(n+2)) n3=n1+2; subplot(3,1,3);stem(n3,x,.b),axis([-3,12,0,5]) text(0,3,x(n-2)) 3.序列反褶（翻转）： 已知序列，求x(n)的反褶序列。 解：x(n)的反褶序列为，设，a＝1.1，－10≤n≤20。序列x(n)的反褶序列x(-n)的Matlab程序实现： % sequence3.m a=1.1; n=[-10:20];x=a.^n; subplot(2,1,1);stem(n,x,.k),axis([-20,20,0,8]) text(0,5,x(n)) n=fliplr(-n);x=fliplr(x); subplot(2,1,2);stem(n,x,.k),axis([-20,20,0,8]) text(0,5,x(-n)) （二）、典型序列 单位冲激序列（单位抽样序列）： 采用Matlab编程实现： function[x,nx]=delta(n) x=(n==0); nx=n; 采用Matlab的 zeros函数实现： x=zeros(1,N); x(1)=1 单位阶跃序列： 采用Matlab编程实现： function u=u(n) u=[n=0]; nu=n; 采用Matlab的 ones函数实现： x=ones(1,N); 矩形序列： function[x,nx]=RNn(N,n); % 设N=5，此时R5(n) x=[(n=0)(n=N-1)]; nx=n; 实指数序列： % sequence4.m a1=1.1;a2=0.9;a3=-1.1;a4=-0.9; n=[-5:15]; x1=(a1.^n);x2=(a2.^n);x3=(a3.^n);x4=(a4.^n); subplot(2,2,1);stem(n,x1,.k);axis([-5,15,-0.5,5]); xlabel(n); text(3.5,3.5,a1); subplot(2,2,2)