N=64;
n=0:63;
t=0.01*n;
q=n*2*pi/N;
x=cos(3*pi/20*t);
                                      %生成信號序列
y=fft(x,N);
subplot(3,1,1);
plot(t,x);
title('source signal');
                                      %繪制源曲線信號
subplot(3,1,2);
plot(q,abs(y));            %幅度譜
title('magnitude');
subplot(3,1,3);
plot(q,angle(y));          %相位譜
title('phase');