实验七 卷积实验.docVIP

实验七 卷积实验（2） [实验目的] 通过编程、上机调试程序，进一步增强使用计算机解决问题的能力。 掌握线性卷与圆周卷积软件实现的方法，并验证两者之间的关系。 [实验内容及要求] 编制一个计算圆周卷积的通用程序，计算两个序列x(n)与h(n)的圆周卷积。 [实验步骤] 编写序列线性卷积函数M文件 %圆周卷积函数 function yc=circonv(x1,x2,N); %ó??±?ó·¨êμ???2?ü?í?y %y=circonv(x1,x2,N); %y:ê?3?DòáD %x1,x2 ê?è?DòáD %N?2?ü?í?yμ?3¤?è if length(x1)N error(N±?D?′óóúμèóúx1μ?3¤?è); end if length(x2)N error(N±?D?′óóúμèóúx2μ?3¤?è); end %ò?é?á???ó????D??á?DòáDμ?3¤?èê?·?D?óúN x1=[x1,zeros(1,N-length(x1))]; %ì?3?DòáDx1(n)ê1??3¤?è?aN￡?DòáDh(n)μ?3¤?è?aN1￡?DòáDx(n)μ?3¤?è?aN2 x2=[x2,zeros(1,N-length(x2))]; %ì?3?DòáDx2(n)ê1??3¤?è?aN n=[0:1:N-1]; x2=x2(mod(-n,N)+1); %éú3éDòáDx2((-n))N H=zeros(N,N); for n=1:1:N H(n,:)=cirshiftd(x2,n-1,N); %?????óμ?KDD?ax2((k-1-n))N end yc=x1*H; 计算圆周卷积 %?????2?ü?í?y function y=cirshiftd(x,m,N); if length(x)N; error(xμ?3¤?è±?D?′óóúóúN); end x=[x,zeros(1,N-length(x))]; n=[0:1:N-1]; y=x(mod(n-m,N)+1); 验证实验结果： x1=[1 2 3 4]; x2=[1 1 1 1]; y=circonv(x1,x2,4) y = 10 10 10 10 y=circonv(x1,x2,8) y = 1 3 6 10 9 7 4 0 同书本结果一致。 研究圆周卷积与线性卷积两者之间的关系 n=[0:1:11]; m=[0:1:5]; N1=length(n); N2=length(m); xn=0.8.^n; hn=ones(1,N2); y1n=conv(xn,hn); ycn=circonv(xn,hn,N1); ny1=[0:1:length(y1n)-1]; ny2=[0:1:length(ycn)-1]; subplot(2,1,1); stem(ny1,y1n); subplot(2,1,2); stem(ny2,ycn); N2=length(m); xn=[1,2,3,4,5]; hn=[1,2,1,2]; yln=conv(xn,hn); ycn=circonv(xn,hn,N1); ny1=[0:1:length(yln)-1]; ny2=[0:1:length(ycn)-1]; subplot(2,1,1); stem(ny1,yln); subplot(2,1,2); stem(ny2,ycn); axis([0,7,0,20]); 6点： clear all n=[0:1:4]; m=[0:1:3]; N1=length(n); N2=length(m); xn=[1 2 3 4 5]; hn=[1 2 1 2]; y1n=conv(xn,hn); ycn=circonv(xn,hn,6); ny1=[0:1:length(y1n)-1]; ny2=[0:1:length(ycn)-1]; subplot(2,1,1); stem(ny1,y1n); subplot(2,1,2); stem(ny2,ycn); axis([0,7,0,20]); 9点： clear all n=[0:1:4]; m=[0:1:3]; N1=length(n); N2=length(m); xn=[1 2 3 4 5]; hn=[1 2 1 2]; y1n=conv(xn,hn); ycn=circonv(xn,hn,9); ny1=[0:1:length(y1n)-1]; ny2=[0:1:length(ycn)-1]; subplot(2,1,1); stem(ny1,y