%《MATLAB及在電子信息課程中的應用》第八章例8.14程序q814
% 四階連續和離散系統開閉環bode圖
% 電子工業出版社出版  陳懷琛 吳大正 高西全合著 2001年10月初版，2003年7月第二版
%
clear,
Ts=0.1
s=zpk(-6,[0,-1,-10,-10],200)	    % 生成四階連續系統s
sd=c2d(s,Ts)						% 變換為四階離散系統sd
sb=feedback(s,1)					% 把連續系統閉合，生成閉環系統sb
sbd=feedback(sd,1)				    % 把離散系統閉合，生成閉環系統sbd
figure(1),bode(s,'--',sb,'.-')		% 在圖1中同時繪出連續系統的開閉環頻率特性
set(gcf,'color','w')                % 設置圖形背景色為白色
figure(2),bode(sd,'--',sbd,'.-')	% 在圖2中同時繪出離散系統的開閉環頻率特性
set(gcf,'color','w')                % 設置圖形背景色為白色
damp(sb)							% 閉環sb的根的固有頻率和阻尼系數
damp(sbd)							% 判斷閉環sbd的根的固有頻率和阻尼系數
[Gm,Pm,wcg,wcp]=margin(s)		    % 判斷閉環sb的穩定裕度
[Gmd,Pmd,wcgd,wcpd]=margin(sd)	    % 判斷閉環sb的穩定裕度