%《MATLAB在電子信息課程中的應(yīng)用》第八章例8.14程序q814
% 四階連續(xù)和離散系統(tǒng)開(kāi)閉環(huán)bode圖
% 電子工業(yè)出版社出版  陳懷琛 吳大正 高西全合著 2001年10月

clear,
Ts=0.1
s=zpk(-6,[0,-1,-10,-10],200)	% 生成四階連續(xù)系統(tǒng)s
sd=c2d(s,Ts)						% 變換為四階離散系統(tǒng)sd
sb=feedback(s,1)					% 把連續(xù)系統(tǒng)閉合，生成閉環(huán)系統(tǒng)sb
sbd=feedback(sd,1)				% 把離散系統(tǒng)閉合，生成閉環(huán)系統(tǒng)sbd
figure(1),bode(s,'--',sb,'.-')		% 在圖1中同時(shí)繪出連續(xù)系統(tǒng)的開(kāi)閉環(huán)頻率特性
figure(2),bode(sd,'--',sbd,'.-')	% 在圖2中同時(shí)繪出離散系統(tǒng)的開(kāi)閉環(huán)頻率特性
damp(sb)							% 閉環(huán)sb的根的固有頻率和阻尼系數(shù)
damp(sbd)							% 判斷閉環(huán)sbd的根的固有頻率和阻尼系數(shù)
[Gm,Pm,wcg,wcp]=margin(s)		% 判斷閉環(huán)sb的穩(wěn)定裕度
[Gmd,Pmd,wcgd,wcpd]=margin(sd)	 % 判斷閉環(huán)sb的穩(wěn)定裕度