%《MATLAB及在電子信息課程中的應用》第七章例7.2程序q702
% 離散信號序列的移位、相加和相乘
% 電子工業出版社出版  陳懷琛 吳大正 高西全合著 2001年10月初版，2003年7月第二版
%
clear;close all
n0=input('輸入序列起點:n0=');
N=input('輸入序列長度:N=');
n1=input('輸入位移:n1=');
n=n0:n1+N+5; 		                    % 生成自變量數組
u=[(n-n1)>=0];                          % 產生單位階躍序列（u(n-n1)）
x1=[(n-n1)>=0]-[(n-n1-N)>=0]	        % 用階躍序列之差產生矩形序列
x2=[(n>=n1)&(n<(N+n1))];	            % 用邏輯式產生矩形序列
x3=exp(j*n*pi/8).*x2;		            % 對復正弦序列加矩形窗（元素群乘）
subplot(2,2,1);stem(n,x1,'.');	
set(gcf,'color','w')                    % 設置圖形背景色為白色
xlabel('n');ylabel('x1(n)');			% 標注
axis([n0,max(n),0,1]);			        % 定坐標范圍
subplot(2,2,3);stem(n,x2,'.');
xlabel('n');ylabel('x2(n)');			% 標注
axis([n0,max(n),0,1]);			        % 定坐標范圍
subplot(2,2,2);stem(n,real(x3),'.');
xlabel('n');ylabel('x3(n)的實部');		 % 標注
line([n0,max(n)],[0,0]);		        % 畫橫軸
axis([n0,max(n),-1,1]);			        % 定坐標范圍
subplot(2,2,4);stem(n,imag(x3),'.');
xlabel('n');ylabel('x3(n)的虛部');		 % 標注
line([n0,max(n)],[0,0]);		        % 畫橫軸
axis([n0,max(n),-1,1]);			        % 定坐標范圍