方波 离散傅里叶级数 MATLAB

%方波 离散时间傅里叶变换
L = 5; N = 10;
k = [-N/2:1:N/2];    %占空比 基本周期 离散时间的参数
xn = [ones(1,L),zeros(1,N-L)]; %生成方波序列
XK = dfs(xn,N);
magXK = abs([XK(N/2+1:N),XK(1:N/2+1)]);
subplot(2,2,3);
stem(k,magXK);
axis([-N/2,N/2,-0.5,5.5]);
xlabel('k');
ylabel('X(k)');
title('DFS of SQ.wave:L = 5,N = 10');
% 零点为占空比的倒数 即 N/L