%《MATLAB及其在理工課程中的應用指南》(十一五規劃版)第九章例9－2－1程序
% 離散信號的MATLAB表述方法
% 西安電子科技大學出版社出版  陳懷琛編著 2007年7月
%
clear,n0=0;nf=10;ns=3;
n1=n0:nf;x1=[zeros(1,ns-n0),1,zeros(1,nf-ns)];  %單位脈沖序列的產生
% n1 = n0:nf; x1=[(n1-ns)==0]; 	% 顯然,用邏輯式是比較高明的方法
n2=n0:nf;x2=[zeros(1,ns-n0),ones(1,nf-ns+1)];  %單位階躍序列的產生
% 也有類似的用邏輯比較語句的方法,留給讀者思考
n3 = n0:nf; x3=(0.9).^n3;			% 實數指數序列
n4 = n0:nf; x4=exp((-0.2+0.5j)*n3);	% 復數指數序列
subplot(2,2,1),stem(n1,x1);title('單位脈沖序列')
subplot(2,2,2),stem(n2,x2);title('單位階躍序列')
subplot(2,2,3),stem(n3,x3);title('實指數序列')
% 注意sunplot的輸入變元
subplot(4,2,6),stem(n4,real(x4));line([0,10],[0,0])
title('復指數序列'),ylabel('實部')
subplot(4,2,8),stem(n4,imag(x4));line([0,10],[0,0]),			%畫橫坐標
ylabel('虛部')