MATLAB做离散傅里叶变换DFT

 

主函数

N=16; n=0:N-1; x1n=exp(j*pi*n/8); X1k=dft(x1n,N); x2n=cos(pi*n/8); X2k=dft(x2n,N); x3n=sin(pi*n/8); X3k=dft(x3n,N); subplot(2,3,1);stem(n,x1n,'.'); title('序列x1(n)'); xlabel('k');ylabel('x1(n)'); subplot(2,3,2);stem(n,x2n,'.'); title('序列x2(n)'); xlabel('k');ylabel('x2(n)'); subplot(2,3,3);stem(n,x3n,'.'); title('序列x3(n)'); xlabel('k');ylabel('x3(n)'); subplot(2,3,4);stem(n,abs(X1k),'.'); title('16点DFT[x1(n)]'); xlabel('k');ylabel('|X1(k)|'); subplot(2,3,5);stem(n,abs(X2k),'.'); title('16点DFT[x2(n)]'); xlabel('k');ylabel('|X2(k)|'); subplot(2,3,6);stem(n,abs(X3k),'.'); title('15点DFT[x3(n)]'); xlabel('k');ylabel('|X3(k)|'); % % 定制开发需求：qq 1762016542

子函数

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;

结果

欢迎关注公众号：算法工程师的学习日志

 

总结

以上是生活随笔为你收集整理的MATLAB做离散傅里叶变换DFT的全部内容，希望文章能够帮你解决所遇到的问题。

如果觉得生活随笔网站内容还不错，欢迎将生活随笔推荐给好友。