%《MATLAB及在電子信息課程中的應用》第八章例8.15程序q815
% 乃奎斯特(nyquist)圖判穩
% 電子工業出版社出版  陳懷琛 吳大正 高西全合著 2001年10月初版，2003年7月第二版
%
clear,
s1=zpk([ ],[-6,-1,1.2],50 );		    % 生成連續系統s1
figure(1)
subplot(2,2,1),nyquist(s1)	,grid		% 畫nyquist圖
set(gcf,'color','w')                    % 設置圖形背景色為白色
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圖
set(gcf,'color','w')                    % 設置圖形背景色為白色
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		% 畫閉環脈沖響應圖