[foxstar]

下面是程序:
N=16;
n=0:N-1;
x=cos(n*pi/8+pi/3);
X=1/N*fft(x,N);
omega=2*pi/N*(n-N/2);%这个为什么要这样写w哦?请高手指导一下
subplot(211);
stem(omega,abs(fftshift(X)));%还有这一条
axis([-pi,pi,0,1]);
ylabel('magnitude');xlabel('frequency(rad)');
subplot(212);
stem(omega,angle(fftshift(X)));
axis([-pi,pi,-4,4]);
ylabel('phase');xlabel('frequency(rad)');

请高手们说一下,先谢谢了!