用MATLAB實現(xiàn)序列的圓周卷積

數(shù)字信號處理實驗報告實驗項目名稱：用MATLAB實現(xiàn)序列的圓周卷積實驗日期：2012-11-28 實驗成績： 實驗評定標準：1）實驗結(jié)果是否正確A（ ）B（ ）C（）2）實驗結(jié)果分析A（ ）B（ ）C（）3）實驗報告是否按照規(guī)定格式A（ ）B（ ）C（）一、 實驗?zāi)康?gt;通過本實驗，掌握一些基本而且重要的離散時間信號，熟悉基本離散時間信號的MATLAB實現(xiàn)方法。二、 實驗器材PC機，MATLAB軟件。三、 實驗內(nèi)容計算兩序列x1（n）={123,4,5}，x2（n）={1,2,3,4,5,4,3,2,1}的圓周卷積。四、 實驗結(jié)果>實驗代碼：clearallclccloseallx1=[1,2,3,4,5,6,7,8];

closeallx2=[1,2,3,4,5,6,7,8,7,6,5,4,3,2,1];N=length(x1)+length(x2);n=0:N-1n1=0:N-2;n2=0:N-3;y1=circonvt(x1,x2,N);y2=circonvt(x1,x2,N-1);y3=circonvt(x1,x2,N-2);x1=[x1zeros(1,N-length(x1))];x2=[x2zeros(1,N-length(x2))];Xf1=dft(x1,N);Xf2=dft(x2,N);Xf=Xf1.*Xf2;x=idft(Xf,N);x=real(x);subplot(2,3,1)stem(n,x1);title('x1(n)');subplot(2,3,2)stem(n,x2);title('x2(n)')title('x1(n)');subplot(2,3,2)stem(n,x2);title('x2(n)')subplot(2,3,3);stem(n,x);title('x(n)=IDFT(X(k))');subplot(2,3,4);stem(n,y1);title('N點圓周卷積');subplot(2,3,5);stem(n1,y2);title('N-1點圓周卷積');subplot(2,3,6);stem(n2,y3);title('N-2點圓周卷積');iflength(x1)>Nerror('N必須>=x1的長度')endiflength(x2)>Nerror('N必須>=x2的長度')endx1=[x1zeros(1,N-length(x1))];x2=[x2zeros(1,N-length(x2))];m=[0:1:N-1];x2=x2(mod(-m,N)+1);H=zeros(N,N);forn=1:1:NH(n,:)=cirshift(x2,n-1,N);endy=x1*H;functiony=cirshift(x,m,N)iflength(x)>Nerror('N必須>=x的長度')x=[xzeros(1,N-length(x))];n=[0:1:N-1];x=[xzeros(1,N-length(x))];n=[0:1:N-1];n=mod(n-m,N);y=x(n+1);function[Xk]=dft(xn,N)n=[0:1:N-1];k=[0:1:N-1];WN=exp(-j*2*pi/N);nk=n'*k;WNnk=WN.Ank;Xk=xn*WNnk;function[xn]=idft(Xk,N)%計算逆離散傅里葉變換%[xn]=idft(Xk,N)n=[0:1:N-1];k=[0