1、實驗七 卷積實驗（ 2）實驗?zāi)康?(1) 通過編程、上機調(diào)試程序，進(jìn)一步增強使用計算機解決問題的能力。(2) 掌握線性卷與圓周卷積軟件實現(xiàn)的方法，并驗證兩者之間的關(guān)系。 實驗內(nèi)容及要求 編制一個計算圓周卷積的通用程序，計算兩個序列x(n)與h(n)的圓周卷積。實驗步驟 1、 編寫序列線性卷積函數(shù)M 文件%圓周卷積函數(shù)function yc=circonv(x1,x2,N);%0 ?± ?6飛卩？?2? u ? i ?y%y=circonv(x1,x2,N);%y: e ?3?DodD%x1,x2 e ? e ?DodD%N?2U ? ?y ii ?3O ?eif length(x1)

2、>Nerror( 'N ± ?D?，6ouieou x1 i ?3。? e'endif length(x2)>Nerror( 'N ± ?D?，66uie6u x2 i ?3。? e'end%o ?e ?d ? 6 ?D?d ?Dod Di ?3O ?ee ? ?D?6uNx1=x1,zeros(1,N-length(x1);%i ?3?Dod Dx1( n) e 1?3 O ?e ?aN£ ?Dod Dh( n) i ?3O ?e ?aNl£ ?Dod Dx(n) i ?3O ?e ?aN2 x2=x2,ze

3、ros(1,N-length(x2);%i ?3?Dod Dx2( n) e 1?3 O ?e ?aNn=0:1:N-1;x2=x2(mod(-n,N)+1);%eu 3 e Dod Dx2(-n)NH=zeros(N,N);for n=1:1:NH(n,:)=cirshiftd(x2,n-1,N);%?6i ?KDD?ax2(k-1-n)Nendyc=x1*H'計算圓周卷積%?2?u ? i ?yfunction y=cirshiftd(x,m,N);if length(x)>N;error( 'x i ?3O ?e± ?D?，6 6uo)u N'en

4、dx=x,zeros(1,N-length(x);n=0:1:N-1;y=x(mod(n-m,N)+1);驗證實驗結(jié)果： x1=1 2 3 4; x2=1 1 1 1;y=circonv(x1,x2,4) y =10 10 10 10y=circonv(x1,x2,8)1 3 6 10 9 7 4 0 同書本結(jié)果一致。研究圓周卷積與線性卷積兩者之間的關(guān)系 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:leng

5、th(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 al

6、l n=0:1:4;m=0:1:3;N1=le ngth( n);N2=le ngth(m);xn=1 2 3 4 5;hn=1 2 1 2;y1n=con v(x n,hn);ycn=circ on v(x n,hn, 6); ny1= 0:1:le ngth(y1 n)-1; ny2=0:1:le ngth(yc n)-1; subplot(2,1,1);stem( ny1,y1 n); subplot(2,1,2);stem (n y2,yc n); axis(0,7,0,20);20LLLiL15-10-一5L-0Cr12345679點：clear all n=0:1:4;m=0:1:

7、3;N1=le ngth( n);N2=le ngth(m);xn=1 2 3 4 5;hn=1 2 1 2;y1n=con v(x n,hn);ycn=circ on v(x n,hn, 9); ny1= 0:1:le ngth(y1 n)-1; ny2=0:1:le ngth(yc n)-1; subplot(2,1,1);stem( ny1,y1 n); subplot(2,1,2);stem (n y2,yc n); axis(0,7,0,20);10點：clear alln=0:1:4;m=0:1:3;N1=le ngth( n);N2=le ngth(m);xn=1 2 3 4 5;hn=1 2 1 2;y1n=con v(x n,hn);ycn=circ on v(x n,hn ,10); ny1= 0:1:le ngth(y1 n)-1; ny2=0:1:le ngth(