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

2、ose allclcx1=1,2,3,4,5,6,7,8;x2=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=x1 zeros(1,N-length(x1);x2=x2 zeros(1,N-length(x2);Xf1=dft(x1,N);Xf2=dft(x2,N);Xf=Xf1.*Xf2;x=idft(Xf,N);x=real(x);subpl

3、ot(2,3,1)stem(n,x1);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點(diǎn)圓周卷積);subplot(2,3,5);stem(n1,y2);title(N-1點(diǎn)圓周卷積);subplot(2,3,6);stem(n2,y3);title(N-2點(diǎn)圓周卷積);function y=circonvt(x1,x2,N)if length(x1)N error(N 必須 = x1的長(zhǎng)

4、度)endif length(x2)N error(N 必須 = x2的長(zhǎng)度)endx1=x1 zeros(1,N-length(x1);x2=x2 zeros(1,N-length(x2);m=0:1:N-1;x2=x2(mod(-m,N)+1);H=zeros(N,N);for n=1:1:N H(n,:)=cirshift(x2,n-1,N);endy=x1*H;function y=cirshift(x,m,N)if length(x)N error(N 必須 = x的長(zhǎng)度)endx=x zeros(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.nk;Xk= xn * WNnk;function xn=idft(Xk,N)%計(jì)算逆離散傅里葉變換%xn=idft(Xk,N)n=0:1:N-1;k=0:1:N-