二、离散信号运算MATLAB实现.docxVIP

实验题目离散信号运算的MATLAB实现　实验时间　实验台号　实验目的序列的和、积、移位运算；序列的卷积和运算；　实验仪器PC台式机，软件MATLAB实验预习1、序列移位；2、序列之和；　3、序列之积；4、时间尺度变换；5、序列翻褶；实验步骤应用函数conv();flipl();1、直接定义法计算序列之和、序列之积。n1=[-5:4];%序列x1的时间起始和终止位置n1s=-5;n1f=4;x1=[2,3,1,-1,3,4,2,1,-5,-3];%序列x1不同时间幅度n2=[0:9];%序列x2的时间起始和终止位置n2s=0;n2f=9;x2=[1,1,1,1,1,1,1,1,1,1];%序列x2不同时间幅度ns=min(n1s,n2s);nf=max(n1f,n2f);%求取新信号的时间起始以及终止位置n=ns:nf;y1=zeros(1,length(n));%延拓序列初始化y2=zeros(1,length(n));y1(find((n=n1s)(n=n1f)==1))=x1;%给延拓序列y1赋值x1y2(find((n=n2s)(n=n2f)==1))=x2;%给延拓序列y2赋值x2ya=y1+y2;yp=y1.*y2;subplot(4,1,1),stem(n,y1,.);ylabel(y1);line([n(1),n(end)],[0,0]);subplot(4,1,2),stem(n,y2,.);ylabel(y2);line([n(1),n(end)],[0,0]);subplot(4,1,3),stem(n,ya,.);ylabel(ya);line([n(1),n(end)],[0,0]);subplot(4,1,4),stem(n,yp,.);ylabel(yp);line([n(1),n(end)],[0,0]);%结果2-12、题目2-1：N=8,a=0.8,求:1、X1(n)=a^n[u(n)-u(n-N)];2、x2(n)=x1(n+3);3、x3(n)=x1(n-2);4、x4(n）=x1(-n)画出x1,x2,x3,x4对应图像。a=0.8;N=8;n=-12:12;Rn=[zeros(1,12),ones(1,N),zeros(1,5)];x=a.^n.*Rn;n1=n;n2=n1-3;n3=n1+2;n4=-n1;subplot(4,1,1);stem(n1,x,fill);grid on;ylabel(x1(n));axis([-15,15,0,1]);subplot(4,1,2);stem(n2,x,fill);grid on;ylabel(x2(n));axis([-15,15,0,1]);subplot(4,1,3);stem(n3,x,fill);grid on;ylabel(x3(n));axis([-15,15,0,1]);subplot(4,1,4);stem(n4,x,fill);grid on;ylabel(x1(n));xlabel(n);axis([-15,15,0,1]);%图像2-23、定义函数法进行序列运算：定义序列相加： function[y,n]=sigadd(x1,n1,x2,n2);n=min(min(n1),min(n2)):max(max(n1),max(n2));y1=zeros(1,length(n));y2=y1;y1(find((n=min(n1))(n=max(n1))==1))=x1;y2(find((n=min(n2))(n=max(n2))==1))=x2;y=y1+y2定义序列相乘：function[y,n]=sigmult(x1,n1,x2,n2);n=min(min(n1),min(n2)):max(max(n1),max(n2));y1=zeros(1,length(n));y2=y1;y1(find((n=min(n1))(n=max(n1))==1))=x1;y2(find((n=min(n2))(n=max(n2))==1))=x2;y=y1.*y2定义序列移位：function[y,ny]=sigshift(x,nx,m);ny=nx+m;y=x;定义序列翻褶：function[y,ny]=sigfold(x,nx);ny=fliplr(x);ny=-fliplr(x);题目2-2：X1=[1,2,3,4],n1=[0:3],x2=[3,4,2,1,3,2],n2=[-2:3]分别计算x1和x2相加、相乘；计算x2移位、翻褶。x1=[1,2,3,4];n1=[0:3];x2=[3,4,2,1,3,2];n2=[-2:3];[y1,n]=sigadd(x1,n1,x2,