實驗七卷積實驗

1、實驗七卷積實驗(2)實驗目的(1) 通過編程、上機調(diào)試程序，進一步增強使用計算機解決問題的能力。(2) 掌握線性卷與圓周卷積軟件實現(xiàn)的方法，并驗證兩者之間的關(guān)系。實驗內(nèi)容及要求編制一個計算圓周卷積的通用程序，計算兩個序列x(n)與h(n)的圓周卷積。實驗步驟1、編寫序列線性卷積函數(shù)M文件%圓周卷積函數(shù)functionyc=circonv(x1,x2,N);%6?±?6"ek?2?u?1?y%y=circonv(x1,x2,N);%y:e?3?D。aD%x1,x2e?e?DoaD%N?2?i?i?yk?3O?eiflength(x1)>Nerror('N士?D?

2、'66也卬倉6ux1?3。？倉/endiflength(x2)>Nerror('N士?D?'66也卬倉6ux2?3。？倉,end%。?e?a?o?D?d?DoaD?3O?ee?-?D?6uNx1=x1,zeros(1,N-length(x1);%i?3?D。dDx1(n)e1?3O?e?aN£?DoaDh(n)?3O?e?aN1£?DoaDx(n)?3O?e?aN2x2=x2,zeros(1,N-length(x2);%i?3?D。aDx2(n)e1?3O?e?aNn=0:1:N-1;x2=x2(mod(-n,N)+1);%eu3eDoaDx2

3、(-n)NH=zeros(N,N);forn=1:1:NH(n,:)=cirshiftd(x2,n-1,N);%?6n?KDD?ax2(k-1-n)Nendyc=x1*H'計算圓周卷積%?2?ti?i?yfunctiony=cirshiftd(x,m,N);iflength(x)>N;error('x?3O?e士?D?'66u6)RN'endx=x,zeros(1,N-length(x);n=0:1:N-1;y=x(mod(n-m,N)+1);驗證實驗結(jié)果：x1=1234;x2=1111;y=circonv(x1,x2,4)y=10101010y=circ

4、onv(x1,x2,8)y=136109740同書本結(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: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=con

5、v(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點：clearalln=0:1:4;m=0:1:3;N1=length(n);N2=length(m);xn=12345;hn=1212;y1n=conv(xn,hn);ycn=circonv(xn,hn,6);ny1=0:1:length(y1n)-1;ny2=0:1:length(ycn)-1;s

6、ubplot(2,1,1);stem(ny1,y1n);subplot(2,1,2);stem(ny2,ycn);axis(0,7,0,20);9點：clearalln=0:1:4;m=0:1:3;N1=length(n);N2=length(m);xn=12345;hn=1212;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,y1n);subplot(2,1,2);stem(ny2,ycn);axis(0,7,0,20);10點:clearalln=0:1:4;m=0:1:3;N1=length(n);N2=length(m);xn=12345;hn=1212;y1n=conv(xn,hn);ycn=circonv(xn,hn,10);ny1=0:1:length(y1n)-1;ny2=0:1:length(ycn