sensfun.m

文件是编辑好的m文件

function [S,T,senfun,tfun]=sensfun(G1,H1,kc,ti);
% Sensfun calculate sensitivity of system with PI controller.
% G1 The Process transfer function,H1 The feedback transfer function.
% kc The gain of controller, ti The integral time of controller.
%kc=5;ti=5/6;G=tf([2 1],[1 3 2 0]);H=tf([1],[.01 1]);

% Copyright 2001-2002 ECUST. 
%   $Revision: 1.2 $  $Date: 2002/02/13 18:48:23 $

syms dc nc nG dG nH dH num den kcs tis dd senfun tfun sn sd tn td s ;
nc=[kcs*tis kcs];dc=[tis 0];G=tf(G1);H=tf(H1);
dc=poly2sym(dc,'s');nc=poly2sym(nc,'s');
nG=poly2sym(G.num{1},'s');dG=poly2sym(G.den{1},'s');
nH=poly2sym(H.num{1},'s');dH=poly2sym(H.den{1},'s');
num=nc*nG*nH;den=dc*dG*dH+nc*nG*nH;dd=dc*dG*dH/den;
senfun=subs(dd,{kcs, tis},{kc,ti});tfun=subs(num/den,{kcs, tis},{kc,ti});
[sn sd]=numden(senfun);[tn td]=numden(tfun);
snp=sym2poly(sn);sdp=sym2poly(sd);tnp=sym2poly(tn);tdp=sym2poly(td);
disp('灵敏度传递函数：');S=tf(snp,sdp)
disp('补灵敏度传递函数：');T=tf(tnp,tdp)