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;