clear all
clc
close all
%實(shí)現(xiàn)R4(n)的DFT
x1n=[1 1 1 1];
M=[3,4,5];
N=2.^M;
for i=1:3
x1k=fft(x1n,N(i));
subplot(3,2,i);
stem(0:N(i)-1,abs(x1k),'.');
grid on
xlabel('k');
ylabel('幅度');
title(['X1(n)  N=' num2str(N(i))]);

subplot(3,2,2+i)
stem(0:N(i)-1,angle(x1k),'.')
grid on
xlabel('k');
ylabel('相位');
title(['X1(n)  N=' num2str(N(i))]);

end