矩形序列离散时间傅里叶变换matlab实现

n=0:6;
x=R(n,7);
w=[0:1:500]*pi/500;
[X]=dtft(x,n,w);
mX=abs(X);aX=angle(X);
re=real(X);im=imag(X);
subplot(2,2,1);
plot(w/pi,mX);
xlabel('以pi为单位的频率');
title('幅度部分');
ylabel('幅度');
subplot(2,2,3);
plot(w/pi,aX);
xlabel('以pi为单位的频率');
title('相角部分');
ylabel('弧度');
subplot(2,2,2);
plot(w/pi,re);
xlabel('以pi为单位的频率');
title('实部');
ylabel('实部');
subplot(2,2,4);
plot(w/pi,im);
xlabel('以pi为单位的频率');
title('虚部');