dfs_sample_sub.m

很多matlab的源代码

xn=[ones(1,L),zeros(1,N-L)];
Xk=dfs(xn,N);

if rem(N,2)==0
   magXk=abs([Xk(N/2+1:N) Xk(1:N/2+1)]);
   k=[-N/2:N/2];
else
   magXk=abs([Xk((N+3)/2:N) Xk(1:(N+1)/2)]);
   k=[-(N-1)/2:(N-1)/2];
end   

x_g=[xn,xn];
n_g=[0:1:length(x_g)-1];

subplot(2,2,1);
h1=stem(n_g,x_g);
set(h1,'MarkerSize',2);
axis([min(n_g),max(n_g),0,1.1])
xlabel('N');
ylabel('x(n)');

subplot(2,2,3);
h2=stem(k,magXk);
xlabel('k');
ylabel('X(k)');
axis([min(k),max(k),0,max(magXk)*1.1])
set(h2,'MarkerSize',2);