www.gusucode.com > matlab编程PSO算法优化求解PID参数源码程序 > gusucode/code/f1.m

    %global y t
function value=f1(x);
s2=tf(1,[1 6 5 0]);
s1=tf([x(1) x(2) x(3)],[1 0]);
sope=s1*s2;
sys=feedback(sope,1);
%step(sys);
[y,t]=step(sys);
%nuber1=length(y)
%nuber2=length(t)
[errs,sigma,tp,ts]=ste(y,t);
value=0.9*ts+0.1*errs;