q802.m

<MATLAB及在电子信息工程中的应用>原程序,

%《MATLAB在电子信息课程中的应用》第八章例8.2程序q802
% 含并联、串联和反馈环节的系统传递函数
% 电子工业出版社出版  陈怀琛 吴大正 高西全合著 2001年10月

format compact
% 方法一  用旧的控制系统工具箱命令
K2=0.08
for K1=[250,1000]
   [fb1,gb1]=feedback(0.1*K1,[0.01,1,1],4*K2,[0.2,1]);
   [fb,gb]=cloop(fb1,conv(gb1,[1,0]));
   K1,K2,printsys(fb,gb)
   p=roots(gb)'
end
% 方法二  用新的控制系统工具箱命令
K1=250
   s1=tf(0.1*K1,[0.01,1,1]);
   s2=tf(4*K2,[0.2,1]);
   s3=zpk([ ], 0, 1);
   sb1= feedback (s1,s2);
   s = feedback(series(sb1,s3),1)