f1.m

利用PSO算法优化求解PID参数,已发表文章,效果很好

%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;