eg10_4.m

方便初学者学习

N=12;%序列长度
n=0:N-1;%时域取样
xn=cos(pi*n/6);%产生序列
k=0:N-1;%频域取样
wn=exp(-j*2*pi/N);
nk=n'*k;
wnnk=wn.^nk;
xk=xn*wnnk %计算DFT
figure(1)%画图
stem(n,xn)
figure(2)
stem(k,abs(xk))