matlab中idfs,DFS IDFS 离散傅里叶级数

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