n=0:50; %定義序列的長度是50
A=1; %設置信號有關的參數
a=0.4;
T=1; %采樣率
w0=2.0374;
x=A*exp(-a*n*T).*sin(w0*n*T); %pi是MATLAB定義的π，信號乘可采用“.*”
close all %清除已經繪制的x(n)圖形
subplot(3,1,1);stem(x); %繪制x(n)的圖形
title('理想信號采樣序列');

k=0:50;
Wn=exp(-i*pi/25);     %旋轉因子
X=x*(Wn.^(k'*n));       %x的fft變換－－X(N)
amgX=abs(X);            %求X的幅度
angleX=angle(X);        %求X的相位
subplot(3,1,2);stem(amgX);
title('x的fft變換的幅度譜');
subplot(3,1,3);stem(angleX);
title('x的fft變換的相位譜');