《matlab序列地运算.docx》由会员分享，可在线阅读，更多相关《matlab序列地运算.docx（8页珍藏版）》请在三一文库上搜索。
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=min(n2)&(n=min(n1)&(n=min
2、(n2)&(n=ns1)&(n=ns2)&(n=nf2)=1)=x2;ya=y1+y2;ym=y1.*y2;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);grid；subplot(224);stem(n,ym,*,r);xlabel(n);ylabel(y1(n)*y2(n);grid；如图1所示。clc;clear all;x1=1,3,5
3、,7,6,4,2,1;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;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
4、(n)+y2(n);grid；subplot(224);stem(n,ym,*,r);xlabel(n);ylabel(y1(n)*y2(n); grid；文档大全(豆+cHA实用标准2、序列的移位解：function y,ny = seqshift(x,nx,m) ny = nx+m;y = x;代码:clc;clear allx=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)
5、;stem(ny,y,*,r);title(移位后序列y(n);xlabel(n);ylabel(y(n);grid;图5文档大全如下图2所示:原.序列Mn)1.5*0 n 鼬后序列y(n)22.533.54n3、序列的卷积解：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;clear all;x = 1,2,3,-1,-2;nx = -1:3;h = 2,2,1,-1,4,-2;nh = -3:2;y,ny = con
6、vwthn(x,nx,h,nh);stem(ny,y,*,r);xlabel(n);ylabel(y(n);grid;x(n)与h(n)相卷积得到的图形如图3所示:4、序列的相关解：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;clear allx = 2,1,3,2,1,5,1;nx = 1:7;y = 2,1,3,4;ny = 1:4;N1 = length(x)-1;N2 = length(y)-1;rxy = convwthn(fliplr(y),-ny,x,nx);k = -(N2):N1;stem(k,rxy,*,r);xlabel(n);ylabel(rxy(n);title( x(n)与y(n)之间的互相关);axis(-N2,N1,0,35);grid;实用标准三、实验小结通过这次实验使我知道序列的卷积可以通过matlab实现。Matlab里本身有丰富的库函数，画图，求卷积各种运算只要调用相应的库函数 即可。本次实9的难点就在conv()函数和各序列长度的计算，如果不 注意其中的计算关系很容易出错。而 conv()函数只要写错，后面的图 形就不能准确的绘制出来。文档大全