%x(n)的[0~6*pi]DTFT
clc;clear;
xn=[1,1,1,1];
N=4;
M=1000;
w=6*pi/M*[0:M-1];
Xejw=dtft(xn,w)
plot(w/pi,abs(Xejw));
title('R4(n)的DTFT波形图');
xlabel('w/pi');
ylabel('幅度');
grid on
%x(n)X(k)与X(ejw)比较
clc;clear;
N1=4;
N2=32;
xn=[1,1,1,1];
k1=(2*pi/N1)*[0:N1-1];
k2=(2*pi/N2)*[0:N2-1];
xn=[xn,zeros(1,N1-length(xn))];
Xk=dft(xn,N1);
%x(n)的4点DFT
subplot(221);
stem(k1/pi,abs(Xk));
title('X(k) n=4点DFT');
xlabel('w/pi');
ylabel('幅度');
grid on
M=1000;
w=2*pi/M*[0:M-1];
Xejw=dtft(xn,w);
%x(n)的4点DTFT
subplot(222);
plot(w/pi,abs(Xejw));
title('X(ejw) n