%《數(shù)字信號(hào)處理教程——MATLAB釋義與實(shí)現(xiàn)》第四章例5.3.1程序hc531
% 用逆z變換求離散系統(tǒng)的輸出
% 電子工業(yè)出版社出版  陳懷琛編著 2004年9月
%
x=[2,4,3,5];		% 輸入序列及初始時(shí)間
nfx=length(x)-1;		% 計(jì)算序列終止時(shí)間
b=-3;a=[2,-2.2,0.5]; 		% 系統(tǒng)函數(shù)的分子分母系數(shù),
B=conv(-3,x);A=a;			% 輸入與分子z變換的多項(xiàng)式乘積,分母不變
[r,p,k]=residuez(B,A)		% 求留數(shù)r,極點(diǎn)p及直接項(xiàng)k		
nf=input('終點(diǎn)時(shí)間nf= ');	% 要求鍵入終點(diǎn)時(shí)間
n = 0:nf;		            % 生成總時(shí)間數(shù)組
% 求無限序列yi和直接序列yd
yi=(r(1)*p(1).^n+r(2)*p(2).^n).*stepseq(0,0,nf);
yd =k(1)*impseq(0,0,nf)+k(2)*impseq(1,0,nf);
y=yi+yd;										 	% 合成輸出
subplot(2,1,1),
stem(n,[x,zeros(1,nf-nfx)],'.')             % 畫輸入序列
xlabel('n'),ylabel('x(n)')
subplot(2,1,2),stem(n,y,'.')                % 畫輸出序列
xlabel('n'),ylabel('y(n)')
set(gcf,'color','w')                        % 置圖形背景色為白