%《MATLAB在電子信息課程中的應用》第八章例8.15程序q815
% nyquist圖判穩
% 電子工業出版社出版  陳懷琛 吳大正 高西全合著 2001年10月

clear,
s1=zpk([ ],[-6,-1,1.2],50 );		% 生成連續系統s1
figure(1)
subplot(2,2,1),nyquist(s1)	,grid			% 畫nyquist圖
sb1=feedback(s1,1)
subplot(2,2,2),impulse(s1),grid		% 畫開環脈沖響應圖	
subplot(2,2,3),margin(s1) ,grid			% 畫nichols圖
subplot(2,2,4),impulse(sb1),grid		% 畫閉環脈沖響應圖
s2=zpk([-.5],[-6,-1,1.2],50 );			% 生成加微分的連續系統s2
figure(2)
subplot(2,2,1),nyquist(s2),grid			% 畫nyquist圖
sb2=feedback(s2,1)
subplot(2,2,2),impulse(s2),grid		% 畫開環脈沖響應圖	
subplot(2,2,3),margin(s2), grid		% 畫nichols圖
subplot(2,2,4),impulse(sb2),grid		% 畫閉環脈沖響應圖