数字信号处理实验离散序列的基本运算

实验卡和实验报告信息科学与工程学院课程名称数字信号处理实验课时4学时实验项目名称和编号离散序列的基本运算同组者姓名实验目的学会运用MATLAB产生常用离散信号，并掌握程序的编写方法，运用find和fliplr子函数，验证型实验，综合型实验实验环境MATLAB实验内容和原理离散序列的时域运算包括信号的相加、相乘，信号的时域变换包括信号的移位、反折、尺度变换等。在MATLAB中，离散序列的相加、相乘等运算是两个向量之间的运算是两个向量之间的运算，因此参加运算的两个序列向量必须具有相同的维数，否则应进行相应的处理。实验卡和实验报告信息科学与工程学院实验步骤方法关键代码一、序列的移位x(n)=x(n-m)1.n1=-10;n2=10;k0=0;k1=-6;k2=4;n=n1:n2;x0=[n=k0];x1=[(n-k1)=0];x2=[(n-k2)=0];subplot(3,1,1),stem(n,x0,'filled','k');axis([n1,n2,1.1*min(x0),1.1*max(x0)]);ylabel('u(n)');subplot(3,1,2),stem(n,x1,'filled','k');axis([n1,n2,1.1*min(x1),1.1*max(x1)]);ylabel('u(n+6)');subplot(3,1,3),stem(n,x2,'filled','k');axis([n1,n2,1.1*min(x2),1.1*max(x2)]);ylabel('u(n-4)');2.n=-2:10;n0=2;n1=-2;x=2*sin(2*pi*n/10);x1=2*sin(2*pi*(n-n0)/10);x2=2*sin(2*pi*(n-n1)/10);subplot(3,1,1),stem(n,x,'filled','k');ylabel('x(n)');subplot(3,1,2),stem(n,x1,'filled','k');ylabel('x(n-2)');subplot(3,1,3),stem(n,x2,'filled','k');ylabel('x(n+2)');二、序列的相加1.n1=0;n2=10;n01=2;n02=4;n=n1:n2;x1=[(n-n01)==0];x2=[(n-n02)==0];x3=x1+x2;subplot(3,1,1);stem(n,x1,'filled','k');axis([n1,n2,1.1*min(x1),1.1*max(x1)]);ylabel('δ(n-2)');subplot(3,1,2);stem(n,x2,'filled','k');实验卡和实验报告信息科学与工程学院axis([n1,n2,1.1*min(x2),1.1*max(x2)]);ylabel('δ(n-4)');subplot(3,1,3);stem(n,x3,'filled','k');axis([n1,n2,1.1*min(x3),1.1*max(x3)]);ylabel('δ(n-2)+δ(n-4)');2.n1=-4:6;n01=-2;x1=[(n1-n01)=0];n2=-5:8;n02=4;x2=[(n2-n02)=0];n=min([n1,n2]):max([n1,n2]);N=length(n);y1=zeros(1,N);y2=zeros(1,N);y1(find((n=min(n1))&(n=max(n1))))=x1;y2(find((n=min(n2))&(n=max(n2))))=x2;x=y1+y2;subplot(3,1,1);stem(n1,x1,'filled','k');ylabel('x1(n)');axis([min(n),max(n),1.1*min(x),1.1*max(x)]);subplot(3,1,2);stem(n2,x2,'filled','k');ylabel('x2(n)');axis([min(n),max(n),1.1*min(x),1.1*max(x)]);subplot(3,1,3);stem(n,x,'filled','k');ylabel('x(n)');axis([min(n),max(n),1.1*min(x),1.1*max(x)]);实验卡和实验报告信息科学与工程学院(接上页)实验步骤方法关键代码三、序列相乘n1=-4:10;x1=3*exp(-0.25*n);n2=-2:6;n02=-1;x2=[(n2-n02)=0];n=min([n1,n2]):max([n1,n2]);N=length(n);y1=zeros(1,N);y2=zeros(1,N);y1(find((n=min(n1))&(n=max(n1))))=x1;y2(find((n=min(n2))&(n=max(n2))))=x2;x=y1.*y2;subplot(3,1,1),stem(n1,x1,'filled','k');ylabel('x1(n)');axis([min(n),max(n),1.1*min(x1),1.1*max(x1)]);subplot(3,1,2),stem(n2,x2,'filled','k');ylabel('x2(n)');axis([min(n),max(n),1.1*min(x2),1.1*max(x2)]);subplot(3,1,3),stem(n,x,'filled','k');ylabel('x(n)');axis([min(n),max(n),1.1*min(x),1.1*max(x)]);四、序列反折n=-4:4;x=exp(-0.3*n);x1=fliplr(x);n1=-fliplr(n);subplot(1,2,1),stem(n,x,'filled','k');title('x(n)');subplot(1,2,2),stem(n1,x1,'filled','k');title('x(-n)');五、序列的尺度变换n=0:40;tn=n./20;x=sin(2*pi*tn);x1=sin(2*pi*tn*2);x2=sin(2*pi*tn/2);实验卡和实验报告信息科学与工程学院subplot(3,1,1),stem(tn,x,'filled','k');ylabel('x(n)');axis([0,2,1.1*min(x),1.1*max(x)]);subplot(3,1,2),stem(tn,x1,'filled','k');ylabel('x(2n)');axis([0,2,1.1*min(x),1.1*max(x)]);subplot(3,1,3),stem(tn,x2,'filled','k');ylabel('x(n/2)');axis([0,2,1.1*min(x),1.1*max(x)]);三、流程图测试记录分析结论一、序列的移位x(n)=x(n-m)实验卡和实验报告信息科学与工程学院1.实验卡和实验报告信息科学与工程学院2.二、序列的相加实验卡和实验报告信息科学与工程学院1.实验卡和实验报告信息科学与工程学院2.三、序列相乘实验卡和实验报告信息科学与工程学院四、序列反折实验卡和实验报告信息科学与工程学院五、序列的尺度变换实验卡和实验报告信息科学与工程学院小结经本次实验我了解了序列的移位、序列的相加、序列相乘、序列反折及序列的尺度变换的实现方法，并通过绘出它们的图形掌握了各个函数的特点，此次实验编写让我对序列的基本应用有所认识以下由实验教师填写记事评议实验卡和实验报告信息科学与工程学院成绩评定平时成绩_______实验报告成绩________综合成绩_________指导教师签名：