用MATLAB实现序列圆周卷积

数字信号处理实验报告实项名：

用MATLAB现列圆卷实日：实评标：1实结是正2实结分3实报是按规格

实成：A）（）（）A）（）（）A）（）（）一实目

通过本实验，掌握一些基本而且重要的离散时间信号，熟悉基本离散时间信号的MATLAB实现方法。二实器PC，MATLAB软件。三实内计两序列（n={1,2,3,4,5}（={1,2,3,4,5,4,3,2,1}圆周卷积。四实结实代码：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);1/3

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);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,N);forn=1:1:NH(n,:)=cirshift(x2,n-1,N);endy=x1*H;functiony=cirshift(x,m,N)

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(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*WNnk;function[xn]=idft(Xk,N)%计算逆离散傅里叶变换%[xn]=idft(Xk,N)n=[0:1:N-1];2/3

k=[0:1:N-1];WN=exp(-j*2*pi/N);nk=n'*k;实验结果：x1(n)

WNnk=WN.^(-nk);xn=(Xk*WNnk)/N;x2(n)

8642

8642

250200150100500

05102025

0

0510152025

0

05102025点圆周卷

点

点25020015010050

2502