主函数
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 176542
子函数
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;
结果
欢迎关注公众号:算法工程师的学习日志