%離散傅里葉變換及x(n),FT[x(n)]的繪圖函數
 function [c,l]=DFT(x,N,str)
 n=0:N-1;
 k=-200:200;
 w=(pi/100)*k;
 l=w;
 c=x*(exp(-j*pi/100)).^(n'*k);  %計算DFT[x(n)]
 magX=abs(c);                    %幅值
 angX=angle(c);                 %輔角
 
 subpolt(2,2,1);
 t=max(x);
 n=0:N-1;
 stem(n,x,'.');
 xlabel('n')
 ylabel('xa(n)');
 title('信號的原型');
 text((0.3*N),(0.8*t),str);
 hold on
 n=0:N-1;
 m=zeros(N);
 plot(n,m);
 
 subpolt(2,2,3);
 plot(w/pi,magX);
 xlabel('w/pi');
 ylabel('|Y(jw)|');
 title('上圖信號的傅氏變換');
 
 %產生信號xa(n)
 function c=FF(A,a,w,fs)
 n=0:50-1;
 c=A*exp((-a)*n/fs).*sin(w*n/fs).*stepseq(0,0,49);
 
 %產生脈沖信號
 function [x,n]=impseq(n0,n1,n2)
 n=[n1:n2];
 x=[(n-n0)==0)];
 
 %產生矩形框信號
 function [x,n]=stepseq(n0,n1,n2)
 n=[n1:n2];
 x=[(n-n0)>=0];