关于MATLAB的FFT转换问题

N=16;
n=0:1:N-1; %n是一个,一行N列的向量
xn=sin(n*pi/8)+sin(n*p1/4); %由n产生的一个序列
k=0:1:N-1; %k是一个,一行N列的向量
WN=exp(-j*2*pi/N);
nk=n'*k; %n'是转置的意思.也就是N行一列*一行N列,得到的是N*N的矩阵
WNnk=WN.^nk %这个涉及点运算,就是对矩阵的每一个元素进行nk次幂运算.不加点,变成矩阵的nk次幂是不一样的.
Xk=xn*WNnk;
subplot(2,1,1)
stem(n,xn); %画火柴杆那种图像.就是函数的未端是一个大点.
stem(k,abs(XK));