exa2_4.m

来自「离散信号及其Matlab实现,包括单位冲激序列 ,单位阶跃序列,斜坡序列等」· M 代码 · 共 17 行

% exa2-4_dft.m, for example 2-4
% to test dft.m

N=32; %length of sequence
n=0:N-1; %time sample
xn=cos(pi*n/6); %generate sequence
k=0:N-1; %frequency sample
WN=exp(-j*2*pi/N);
nk=n'*k;
WNnk=WN.^nk;
Xk=xn*WNnk; %compute DFT
%plot sequence and its DFT
figure(1)
stem(n,xn)
figure(2)
stem(k,abs(Xk))