%《MATLAB在電子信息課程中的應用》第七章例7.8程序q708
% 求z的多項式分式的逆z變換
% 電子工業出版社出版  陳懷琛 吳大正 高西全合著 2001年10月

clear, close all
x=[2,4,3,5];nsx=-2;		% 輸入序列及初始時間
nfx=nsx+length(x)-1;		% 計算序列終止時間
Bw=-3; nsbw=-1;				% 系統函數的分子系數,及z的最高次數
Aw=[2,-2.2,0.5]; nsaw=0;	% 系統函數的分母系數,及z的最高次數
B=conv(-3,x);				% 輸入與分子z變換的多項式乘積
A=Aw;						% 分母不變
nsy=nsaw-(nsbw-nsx)			% 分子比分母z的次數高nsy
[r,p,k]=residuez(B,A)		% 求留數r,極點p及直接項k		
nf=input('終點時間nf= ');	% 要求鍵入終點時間
n = min(nsx,nsy):nf;		% 生成總時間數組
% 求無限序列yi和直接序列yd
yi=(r(1)*p(1).^(n-nsy)+r(2)*p(2).^(n-nsy)).*stepseq(nsy,n(1),nf);
yd =k(1)*impseq(nsy,n(1),nf)+k(2)*impseq(-1-nsy,n(1),nf);
y=yi+yd;										 	% 合成輸出
xe = zeros(1,length(n));							% 初始化,將x延拓為xe,
xe(find((n>=nsx)&(n<=nfx)==1))=x;					% 在對應的n處把xe賦值x
subplot(2,1,1),stem(n,xe,'.'),line([min(n(1),0),nf],[0,0]) % 繪圖
subplot(2,1,2),stem(n,y,'.'),line([min(n(1),0),nf],[0,0])
subplot(2,1,1),stem(n,xe,'.'),
line([min(n(1),0),nf],[0,0]) % 繪圖
subplot(2,1,2),stem(n,y,'.'),
line([min(n(1),0),nf],[0,0])