%《MATLAB在電子信息課程中的應用》第八章例8.10程序q810
% 帶時延環節系統的分析
% 電子工業出版社出版  陳懷琛 吳大正 高西全合著  2001年10月

clear, 
K=0.1;Tm=5;Td=1;
s =  tf(K*[Tm,1],[1,0,0],'Td',Td)		% 建立系統在開環狀態下的模型
spd=pade(s,4)						% 以四階多項式來近似代替延時環節的模型
sbpd=feedback(spd,1)				% 求閉環模型
damp(sbpd),pause					% 求閉環模型的阻尼系數
subplot(2,2,1),pzmap(s),pause			% 畫出原始開環模型s的零極點
subplot(2,2,2),pzmap(spd),pause		% 畫出近似開環模型的零極點
subplot(2,2,3),pzmap(sbpd),pause		% 畫出閉環模型的零極點
subplot(2,2,4),impulse(sbpd) 			% 畫出閉環環系統的脈沖響應
pause,subplot(1,1,1)