%《MATLAB在電子信息課程中的應用》第七章例7.4程序q704
% 求時域離散系統對常用序列的響應
% 電子工業出版社出版  陳懷琛 吳大正 高西全合著 2001年10月

%N: 輸入數據長度 
%B: 差分方程系數
%A: 差分方程系數
clear;close all;
N=64;n=0:N-1;m=10;
% 設定系統參數A,B
B=0.0003738* conv([1,1],conv([1,1],conv([1,1],conv([1,1],conv([1,1],[1,1])))))
A=conv([1,-1.2686,0.7051],conv([1,-1.0106,0.3583],[1,-0.9044,0.2155]))
x1=[n==0]; 			% 產生輸入信號x1(n)
y1=filter(B,A,x1); 	% 對x1(n)的響應
x2=[(n-m)==0];		% 產生輸入信號x2(n)
y2=filter(B,A,x2); 	% 對x2(n)的響應
x3=[n>=0];				% 產生輸入信號x3(n)
y3=filter(B,A,x3); 	% 對x3(n)的響應
x4=[(n>=0)&(n<32)];	% 產生輸入信號x4(n)
y4=filter(B,A,x4); 	% 對x4(n)的響應
x5=exp(j*pi*n/8);		% 產生輸入信號x5(n)
y5=filter(B,A,x5); 	% 對x5(n)的響應
subplot(3,2,1);stem(n,y1,'.');line([0,N],[0,0])
axis([0,N,min(y1),max(y1)]);ylabel('y1(n)')
subplot(3,2,2);stem(n,y2,'.');line([0,N],[0,0])
axis([0,N,min(y2),max(y2)]);ylabel('y2(n)')
subplot(3,2,3);stem(n,y3,'.');line([0,N],[0,0])
axis([0,N,min(y3),max(y3)]);ylabel('y3(n)')
subplot(3,2,4);stem(n,y4,'.');line([0,N],[0,0])
axis([0,N,min(y4),max(y4)]);ylabel('y4(n)')
subplot(3,2,5);stem(n,real(y5),'.');line([0,N],[0,0])
axis([0,N,-1,1]);ylabel('Re[y5(n)]')
subplot(3,2,6);stem(n,imag(y5),'.');line([0,N],[0,0])
axis([0,N,-1,1]);ylabel('Im[y5(n)]')
set(gcf,'color','w')