eg10_4.m

来自「详细的介绍MATLAB 7 的功能」· M 代码 · 共 13 行

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))