怀疑你给的公式有错,结果与你给的图像不一致。
经分析,我认为指数项的系数应该是π而不是2π。把取π和2π的图都画出来作为对比如下:
参考代码:
syms n
N = 16;
L = 32;
k = 0 : L-1;
x = sin(2*pi*0.135*n)+cos(2*pi*(0.135+0.06)*n);
subplot 211
a = @(k)symsum(x.*exp(i*(-2*pi*n*k/L)),n,0,N-1);
p = abs(a(k)).^2 / N;
stem(k,p,'fill')
xlim([0 31])
title('指数项为 e^{-j2\pink/L}')
subplot 212
a = @(k)symsum(x.*exp(i*(-pi*n*k/L)),n,0,N-1);
p = abs(a(k)).^2 / N;
stem(k,p,'fill')
xlim([0 31])
title('指数项为 e^{-j\pink/L}')