%《MATLAB在電子信息課程中的應用》第六章例6.3程序q603
% 求任意線性系統的沖擊函數
% 電子工業出版社出版  陳懷琛 吳大正 高西全合著  2001年10月

a=input('多項式分母系數向量a=(書上取 poly([0,-1+2j,-1-2j,-2,-5]) ）  ');
b=input('多項式分子系數向量b=(書上取 [8,3,1] ）  ');
[r,p] = residue(b, a), 	 % 求留數
disp('解析式h(t)=Σr(i)*exp(p(i)*t)')
disp('給出時間數組t=[0:dt:tf]')
dt=input('dt= ');					% 給定時間數組
tf=input('tf= ');	
t=0:dt:tf;h=zeros(1,length(t)); 	%  h的初始化
for i=1:length(a)-1 				% 根數為a的長度減1
   h = h+ r(i)*exp(p(i)*t); 	% 疊加各根分量
end
plot(t,h),grid