WN=exp(-j*2*pi/N); nk=n'*k; WNnk=WN.nk; Xk=xn*WNnk; subplot(2,1,1) stem(n,xn); subplot(2,1,2) stem(k,abs(Xk)) 2.矩形序列x(n)=R(n),求N分别取8,32时的DFT,最后绘出结果图形。 参考程序如下： function[Xk]=dft (xn,N) n=[0:1:N-1]: %n的行向量 k=[0:1:N-1]; k的行向量 WN=exp (-j*2*pi/N) %旋转因子 nk=n'*k; %产生一个含k值的N乘N维矩阵 WNnk=WN.nk; DFT矩阵 Xk=xn*WNnk %DFT系数的行向量 调用上面函数解题。 N=8;x=[ones(1,5),zeros(1,N-5)]: n=0:N-1: X=dft (x,N) N=8点离散傅立叶变换 magX=abs (X):phaX=angle(X)*180/pi; k=(0:length (magX)'-1)*N/length (magX); subplot(2,2,1)stem(n,x);ylabel('x(n)'); subplot(2,2,2)stem(k,magX);axis([0,10,0,5]);ylabel('X(k)|'); N=32:x=[ones(1,5),zeros(1,N-5)]: n=0:-1: X=dft (x,N) N=32点离散傅立叶变换 magX=abs(X)phaX=angle(X)*180/pi; k=(0:length (magX)'-1)*N/length (magX) 1111 WN=exp(-j*2*pi/N); nk=n'*k; WNnk=WN.^nk; Xk=xn*WNnk; subplot(2,1,1) stem(n,xn); subplot(2,1,2) stem(k,abs(Xk)); 2.矩形序列 x(n)=R5(n),求 N 分别取 8，32 时的 DFT，最后绘出结果图形。 参考程序如下： function[Xk]=dft(xn,N) n=[0:1:N-1]; %n 的行向量 k=[0:1:N-1]; %k 的行向量 WN=exp(-j*2*pi/N); %旋转因子 nk=n'*k; %产生一个含 nk 值的 N 乘 N 维矩阵 WNnk=WN.^nk; %DFT 矩阵 Xk=xn*WNnk; %DFT 系数的行向量 调用上面函数解题。 N=8;x=[ones(1,5),zeros(1,N-5)]; n=0:N-1; X=dft(x,N); %N＝8 点离散傅立叶变换 magX=abs(X);phaX=angle(X)*180/pi; k=(0:length(magX)'-1)*N/length(magX); subplot(2,2,1);stem(n,x);ylabel('x(n)'); subplot(2,2,2);stem(k,magX);axis([0,10,0,5]);ylabel('|X(k)|'); N=32;x=[ones(1,5),zeros(1,N-5)]; n=0:N-1; X=dft(x,N); %N＝32 点离散傅立叶变换 magX=abs(X);phaX=angle(X)*180/pi; k=(0:length(magX)'-1)*N/length(magX);