clf ;
L = 5 ;
N = 20 ;
xn =[ones(1,L) zeros(1,N-L)];
k = [-N/2:N/2];
Xk = dfs(xn, N);
magXk = abs([Xk(N/2 + 1:N ), Xk(1:N/2+1)]);
subplot(221);
stem(k, magXk);
xlabel('k');
ylabel('Xk');
title('L = 5 , N = 20 wave ');
L = 5 ;
N = 40 ;
xn =[ones(1,L) zeros(1,N-L)];
k = [-N/2:N/2];
Xk = dfs(xn, N);
magXk = abs([Xk(N/2 + 1:N ), Xk(1:N/2+1)]); % 这个的目的就有点类似于fftshift
subplot(222);
stem(k, magXk);
xlabel('k');
ylabel('Xk');
title('L = 5 , N = 40 wave ');
L = 6 ;
N = 40 ;
xn =[ones(1,L) zeros(1,N-L)];
k = [-N/2:N/2];
Xk = dfs(xn, N);
magXk = abs([Xk(N/2 + 1:N ), Xk(1:N/2+1)]);
subplot(223);
stem(k, magXk);
xlabel('k');
ylabel('Xk');
title('L = 6 , N = 40 wave ');
L = 7;
N = 20 ;
xn =[ones(1,L) zeros(1,N-L)];
k = [-N/2:N/2];
Xk = dfs(xn, N);
magXk = abs([Xk(N/2 + 1:N ), Xk(1:N/2+1)]);
subplot(224);
stem(k, magXk);
xlabel('k');
ylabel('Xk');
title('L = 7 , N = 20 wave ');