std_tf.m

反馈控制系统设计与分析（matlab程序）

function G=stdtf(id,wn,n)
TypeISys=[1,1,0,0,0,0,0; 
          1,1.41,1,0,0,0,0; 
          1,1.75,2.15,1,0,0,0; 
          1,2.1,3.4,2.7,1,0,0; 
          1,2.8,5.0,5.5,3.4,1,0; 
          1,3.25,6.6,8.6,7.45,3.95,1];
TypeIISys=[1,3.2,1,0,0,0,0;
           1,1.75,3.25,1,0,0,0; 
           1,2.41,4.93,5.14,1,0,0;
           1,2.19,6.5,6.3,5.24,1,0; 
           1,6.12,13.42,17.16,14.14,6.76,1];
TypeBSys=[1,1,0,0,0,0,0; 
          1,1.4,1,0,0,0,0;
          1,2,2,1,0,0,0;
          1,2.6,3.4,2.6,1,0,0; 
          1,3.24,5.24,5.24,3.24,1,0; 
          1,3.86,7.46,9.13,7.46,3.86,1];
if n>6, 
   error('the order is too high'); return; 
end
if id==1, mats=TypeISys; 
elseif id==2, mats=TypeIISys; n=n-1;
else, mats=TypeBSys; end
dd=mats(n,1:n+1).*(wn*ones(1,n+1)).^[0:n]; 
nn=dd(n+1); if id==2, nn=dd(n:n+1); end
G=tf(nn,dd);