实验三_线性卷积与圆周卷积的计算

通信班课程考试报告学院：基础信息工程系统：电子信息工程课程名称：数字信号处理名字：学位：日期：实验3实验名称：线性卷积和圆周卷积的计算一、实验目的(1)不通过编程、调试程序，进一步加强使用计算机解决问题的能力。(2)掌握线性卷积和圆周卷积软件实现的方法，并对两者的关系进行实验。二、主要函数的介绍对于无限长序列，MATLAB无法直接计算线性卷积。 在MATLAB内部，只能使用一个conv函数计算两个有限长度序列的线性卷积。 由于圆周卷积MATLAB内部不提供现成的函数，所以可以基于表达式直接计算。计算两个有限长序列的线性卷积函数y(n)=conv(x，h )输入参数： x，h是已知的两个有限长序列向量。输出参数： y是线性卷积的序列向量。三、实验内容和要求已知两个有限长序列x (n )=(n ) 2(n-1 ) 3(n-2 ) 4(n-3 ) 5(n-4 )。h(n)=(n) 2(n-1) (n-2) 2(n-3 )。(1)实验前，事先计算这两个序列的线性卷积和以下情况下的圆周卷积x (n )h (n )x (n )h (n )x (n )h (n )x (n )h (n )(2)创建一个计算两个序列的线性卷积的通用程序，并计算x(n)*h(n )。(3)创建用于计算圆周卷积的通用程序，并且计算上述四种情况下的两个序列x(n )和h(n )的圆周卷积。(4)在桌子上调试实验结果打印或记录。注意：上述两个卷积的计算可以通过一个程序以菜单形式实现。(5)将实验结果与预先计算的结果进行比较，验证其真实性。首先，定义圆周卷积的函数，并保存为. m文件。函数YC=circo NV (x1，x2，n )一个if长度(x1) nerror(N必须大于或等于x1 )结束2个if长度(n )error(N必须大于或等于x2 )结束x1=x1，zeros(1，n长度(x1 ) ) ；x2=x2，zeros(1，n长度(x2 ) ) ；n=0:N-1；x2=x2 (模型(-n，N) 1)H=zeros(N，n )for n=1:1:NH(n，)=cirshiftd(x2，n-1，n )；结束yc=x1*H；重新定义圆周位移函数，并保存为. m文件。函数y=cirshiftd (x，m，n )if长度(x ) nerror(x的长度必须小于n )结束x=x，zeros(1，N-length(x)；n=0:1:N-1；y=x (模数(n-m，N) 1)在此基础上执行MATLAB函数(1) x(n)y(n )。清除全部；N1=5；N2=4；xn=1 2 3 4 5；hn=1 2 1 2；yln=conv(xn，hn )ycn=circonv(xn，hn，5 )ny1= 0:603360长度(yln )-1 ；ny2= 0333 601:长度(ycn )-1 ；subplot (2，1，1 )stem(ny1、yln )ylabel (线性卷积)subplot (2，1，2 )stem(ny2，ycn )ylabel (圆周卷积) x(n)y(n )。清除全部；N1=5；N2=4；xn=1 2 3 4 5；hn=1 2 1 2；yln=conv(xn，hn )ycn=circonv(xn，hn，6 )ny1= 0:603360长度(yln )-1 ；ny2= 0333 601:长度(ycn )-1 ；subplot (2，1，1 )stem(ny1、yln )ylabel (线性卷积)subplot (2，1，2 )stem(ny2，ycn )ylabel (圆周卷积) x(n)y(n )。清除全部；N1=5；N2=4；xn=1 2 3 4 5；hn=1 2 1 2；yln=conv(xn，hn )PS=PS K (xn，hn，9 )ny1= 0:603360长度(yln )-1 ；ny2= 0333 601:长度(ycn )-1 ；subplot (2，1，1 )stem(ny1、yln )ylabel (线性卷积)subplot (2，1，2 )stem(ny2，ycn )ylabel (圆周卷积) x(n)y(n )清除全部；N1=5；N2=4；xn=1 2 3 4 5；hn=1 2 1 2；yln=conv(xn，hn )ycn=circonv(xn，hn，10 )ny1= 0:603360长度(yln )-1 ；ny2= 0333 601:长度(ycn )-1 ；subplot (2，1，1 )stem(ny1、yln )ylabel (线性卷积)subplot (2，1，2 )stem(ny2，ycn )ylabel (圆周卷积)(2)(3)。清除全部PSn=五六九十； %圆周卷积的长度向量xn=1 2 3 4 5；hn=1 2 1 2；y1n=conv(xn，hn)%计算线性卷积ny1=0:长度(y1n )-1；%计算x(n )和h(n )五点、六点、九点和十点的圆周卷积yc1=circonv(xn，hn，N(1) )yc2=circonv(xn，hn，N(2) )yc3=circonv(xn，hn，N(3) )yc4=circonv(xn，hn，N(4) )%分别取与线性卷积不同的分数，进行了圆周卷积的图像比较subplot (1，2，1 )stem(ny1，y1n )xlabel (时间编号n )ylabel (信号振幅)title (线性卷积)subplot (1，2，2 )stem(0:N(1)-1，yc1 )xlabel (时间编号n )ylabel (信号振幅)title(5点圆周卷积)花样滑冰。subplot (1，2，1 )stem(ny1，y1n )xlabel (时间编号n )ylabel (信号振幅)title (线性卷积)subplot (1，2，2 )stem(0:N(2)-1，yc2 )xlabel (时间编号n )ylabel (信号振幅)title(6点圆周卷积)花样滑冰。subplot (1，2，1 )stem(ny1，y1n )xlabel (时间编号n )ylabel (信号振幅)title (线性卷积)subplot (1，2，2 )stem(0:N(3)-1，yc3 )xlabel (时间编号n )ylabel (信号振幅)title(9点圆周卷积)花样滑冰。subplot (1，2，1 )stem(ny1，y1n