资源描述：
《matlab序列的运算》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。
1、实验三序列的运算一、实验内容1、序列的乘积：仿真并完成例题1.222、序列的移位：编写seqshift()函数，实现序列的移位，并作图验证。3、序列的卷积：编写convwthn()函数，实现具体序列的卷积，并且进行理论计算，比较结果。4、序列的相关：仿真并完成例题1.24二、MATLAB仿真1、序列的乘积解：function[y,n]=seqadd(x1,n1,x2,n2)n=min((n1),min(n2)):max((n1):max(n2));y1=zeros(1,length(n));y2=y1;y1(find((n>=min(n1))&(n<=max(n1))=
2、=1))=x1;y2(find((n>=min(n2))&(n<=max(n2))==1))=x2;y=y1+y2;function[y,n]=seqmult(x1,n1,x2,n2);n=min((n1),min(n2)):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;代码：clc;clearall;x1=[1,3,5,7,6,4,2,1]
3、;ns1=-3;x2=[4,0,2,1,-1,3];ns2=1;nf1=ns1+length(x1)-1;nf2=ns2+length(x2)-1;n1=ns1:nf1;n2=ns2:nf2;n=min(ns1,ns2):max(nf1:nf2);y1=zeros(1,length(n));y2=y1;y1(find((n>=ns1)&(n<=nf1)==1))=x1;y2(find((n>=ns2)&(n<=nf2)==1))=x2;ya=y1+y2;ym=y1.*y2;7subplot(221);stem(n1,x1,'*','r');ylabel('x1(n)')
4、;grid;subplot(223);stem(n2,x2,'*','r');xlabel('n');ylabel('x2(n)');grid;subplot(222);stem(n,ya,'*','r');ylabel('y1(n)+y2(n)');grid;subplot(224);stem(n,ym,'*','r');xlabel('n');ylabel('y1(n)*y2(n)');grid;如图1所示。clc;clearall;x1=[1,3,5,7,6,4,2,1];ns1=-3;x2=[4,0,2,1,-1,3];ns2=1;nf1=ns1+length(
5、x1)-1;nf2=ns2+length(x2)-1;n1=ns1:nf1;n2=ns2:nf2;[ya,n]=seqadd(x1,n1,x2,n2);[ym,n]=seqmult(x1,n1,x2,n2);subplot(221);stem(n1,x1,'*','r');ylabel('x1(n)');grid;subplot(223);stem(n2,x2,'*','r');xlabel('n');ylabel('x2(n)');grid;subplot(222);stem(n,ya,'*','r');ylabel('y1(n)+y2(n)');7grid;subp
6、lot(224);stem(n,ym,'*','r');xlabel('n');ylabel('y1(n)*y2(n)');grid;图172、序列的移位解：function[y,ny]=seqshift(x,nx,m)ny=nx+m;y=x;代码：clc;clearallx=[1,2,3,4,5];nx=-2:2;m=2;%向右移2个单位[y,ny]=seqshift(x,nx,m);subplot(211);stem(nx,x,'*','r');title('原序列x(n)');xlabel('n');ylabel('x(n)');subplot(212);stem
7、(ny,y,'*','r');title('移位后序列y(n)');xlabel('n');ylabel('y(n)');grid;如下图2所示：图273、序列的卷积解：function[y,ny]=convwthn(x,nx,h,nh)ny1=nx(1)+nh(1);ny2=nx(end)+nh(end);y=conv(x,h);ny=(ny1:ny2);代码：clc;clearall;x=[1,2,3,-1,-2];nx=-1:3;h=[2,2,1,-1,4,-2];nh=-3:2;[y,ny]=convwthn(x,nx,h,nh)
