资源描述：
《用MATLAB实现序列的圆周卷积》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、数字信号处理实验报告实验项目名称：用MATLAB实现序列的圆周卷积实验日期：2012-11-28实验成绩：实验评定标准：1）实验结果是否正确A（）B（）C（）2）实验结果分析A（）B（）C（）3）实验报告是否按照规定格式A（）B（）C（）一、实验目的Ø通过本实验，掌握一些基本而且重要的离散时间信号，熟悉基本离散时间信号的MATLAB实现方法。二、实验器材ØPC机，MATLAB软件。三、实验内容Ø计算两序列x1（n）={1,2,3,4,5}，x2（n）={1,2,3,4,5,4,3,2,1}的圆周卷积。四、实验结果
2、Ø实验代码：clearallcloseallclcx1=[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=[x1zeros(1,N-length(x1))];x2=[x2zeros(1,N-length(x2))];Xf1=
3、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)')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);ti
4、tle('N-1点圆周卷积');subplot(2,3,6);stem(n2,y3);title('N-2点圆周卷积');functiony=circonvt(x1,x2,N)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,
5、N);forn=1:1:NH(n,:)=cirshift(x2,n-1,N);endy=x1*H;functiony=cirshift(x,m,N)iflength(x)>Nerror('N必须>=x的长度')endx=[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.^nk;Xk=xn*W
6、Nnk;function[xn]=idft(Xk,N)%计算逆离散傅里叶变换%[xn]=idft(Xk,N)n=[0:1:N-1];k=[0:1:N-1];WN=exp(-j*2*pi/N);nk=n'*k;WNnk=WN.^(-nk);xn=(Xk*WNnk)/N;Ø实验结果：一、实验结果分析
