资源描述：
《实验4线性卷积与圆周卷积的计算》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、题目：已知两个有限长序列x(n)=8(n)+28(n-1)+3S(n-2)+4§(n-3)+56(n-4)h(n)=5(n)+26(n~l)+5(n~2)+25(n~3)计算以下两个序列的线性卷积和圆周卷积(1)x(n)(§)y(n)(2)x(n)@y(n)(3)x(n)(9)y(n)(4)x(n)⑩y(n)•调用函数circonvfunctionyc=circonv(xlzx2,N)令用直接法实现圆周卷积%y=circonv(xlzx2,N)%y:输岀序列%xl,x2:输入序列%N:圆周卷积的长度iflength(xl)>Nerror;endiflength
2、(x2)>Nerror;end马以上语句判断两个序列的长度是否小于Nxl=[xl,zeros(lzN-length(xl))];辛填充序列xl(n)使其长度为N,序列h(n)的长度为Nl,序列x(n)的长度为N2x2=[x2/zeros(1,N-length(x2))];令填充序列x2(n)使其长度为Nn=[0:1:N-1];x2=x2(mod(-n,N)+1);%生成序列x2((-n))N,镜像，可实现对x(n)以N为周期的周期延拓，力111是因为MATLAB向量下标只能从1开始。H=zeros(N,N);%生成N行N列的零矩阵forn=l:1:NH(n,:
3、)=cirshiftd(x2,n-1fN);%该矩阵的k彳了为x2((k-l-n))Nendyc=xl*H*;%计算圆周卷积•调用函数cirshiftdfunctiony=cirshiftd(xfmAN)饗直接实现序列x的圆周移位%y=cirshiftd(xzm,N)%x:输入序列，冃它的长度小于N令m:移位位数%N：圆周卷积的长度%y:输出的移位序列iflength(x)>Nerror(*x的氏度必须小丁"N*);endx=[x,zeros(1,N-length(x))];n=[0:1:N-l];y=x(mod(n-m,N)+1);•函数(1)x(n^y(n
4、)clearall;Nl=5;N2=4;xn=[12345];%4<成x(n)hn=[l212];%牛成h(n)yln=conv(xnzhn);%宜接用函数convil"算线性卷积ycn=circonv(xn,hn,5);%用函数circonv计算N1点圆周卷积ny1=[0:1:length(yin)-1];ny2=[0:1:length(yen)-1];subplot(2,1r1);%lDll图stem(nylzyin);ylabel「线性卷积T;subplot(2,1,2);stem(ny2,yen);ylabel(，圆周卷积，);更報匝画•函数(2)x(
5、n)(§)y(n)clearall;Nl=5;N2=4;xn=[l2345];%d成x(n)hn=[l212];%生成h(n)yln=conv(xn,hn);%直接用函数conv计算线性卷积ycn=circonv(xn,hn,6);%用函数circonv计算N1点圆周卷积nyl=[O:l:length(yln)-l];ny2=[0:1:length(yen)-1];subplot(2,1,1);stem(nyl,yin);ylabel(1线性卷积T;subplot(2,1/2);stem(ny2,yen);ylabel(*圆周卷积201510502015105
6、0.51.522.5J4<)G5•函数(3)x(n2y(n)clearall;Nl=5;N2=4;xn=[12345];%4<成x(n)hn=[l212];%牛.成h(n)vln=conv(xnzhn);%Tl接用函数conv计算线性卷积ycn=circonv(xn,hn,9);%用函数circonv计算N1点圆周卷积ny1=[0:1:length(yin)-1];ny2=[0:1:length(yen)-1];subplot(2,1,1);stem(nyl,yin);ylabel(*线性卷枳T;subplot(2Z2);stem(ny2,yen);ylabe
7、K*圆周卷积T;1510502015105
8、012345678•函数(4)x(n)(®y(n)clearall;Nl=5;N2=4;xn=[l2345];%d成x(n)hn=[l212];%生成h(n)yln=conv(xn,hn);%直接用函数conv计算线性卷积ycn=circonv(xn,hn,10);%用函数circonv计算N1点圆周卷积ny1=[0:1:length(yin)-1];ny2=[0:1:length(yen)-1];subplot(2Z1,1);stem(nyl,yin);ylabel(*线性卷积T;subplot(2,1,2);st
9、em(ny2,yen);ylabel(
