opt_fun.m

用MATLAB的脚本语言M语言编写的

function y=opt_fun(x,G,key,nn,nd,dc)
ff0=1e10; alph=[1,x(nn+1:nn+nd)]; 
beta=x(1:nn+1); beta(end)=alph(end)*dc; 
g=tf(beta,alph); 
if key==1, 
   tau=x(end); if tau<=0, tau=eps; end
   [nP,dP]=pade(tau,3); gP=tf(nP,dP);
else, gP=1; end
G_err=G-g*gP; 
G_err.num{1}=[0,G_err.num{1}(1:end-1)];
[y,ierr]=geth2(G_err);
if ierr==1, y=10*ff0; else, ff0=y; end