%《MATLAB及其在理工課程中的應用指南》第九章例9－3－1程序
% 傳遞函數和零極點的關系
% 西安電子科技大學出版社出版  陳懷琛編著 1999年10月初版，2004年11月第二版
%
f=[8,2,1],g=[1.0000, 3.3000, 9.9000, 2.7000],
z=roots(f), p=roots(g),k=f(1)/g(1)		                % 由傳遞函數求零極增益
[r,p,h] = residue(f,g);						            % 由零極增益求極點留數
t=0:0.1:5;
x = r(1)*exp(p(1)*t) + r(2)*exp(p(2)*t) + r(3)* exp(p(2)*t);	% 求過渡函數
plot(t,x),grid
set(gcf,'color','w')                                    % 置圖形背景色為白色