% Matlab code of the PSO Algorithm applied to RIP
% Corresponding Contact: Nguyen Duc Hoang (ndhoang@hcmut.edu.vn),
% (C) 2008 Ho Chi Minh City of University Technology
function RIP_Fuzzy_PSO
clc; close all; clear all;
max_iteration = 50;
n = 50; % bird population size
rang=[ 0 0 0 0 0.00 0.50 0.1 0.01 0.1 0.01 1.0;... % range of the variables
1 1 1 1 0.50 1.00 2.0 1.00 2.0 1.00 100.0];
index = 1:11;
C1=2.5;
C2=1.8;
if((C1+C2)<4)
C1 = C1*4.1/(C1+C2);
C2 = C2*4.1/(C1+C2);
w = 2/(abs(2-(C1+C2)-sqrt((C1+C2)*(C1+C2)-4*(C1+C2))));
else
w = 2/(abs(2-(C1+C2)-sqrt((C1+C2)*(C1+C2)-4*(C1+C2))));
end
U=X_random(n,rang);
for i=1:n
Var = U(i,:);
evaluate;
J(i) = JJ;
end
Par_N=sortrows([U,J'],12)
Par_O=Par_N;
Pbest=Par_N;
Gbest=Pbest(1,:);
V_N=zeros(n,11);
V_O=zeros(n,11);
paras;
clear U; l=0;
%Iterations of the algorithm
for i=1:max_iteration % iterations
pause(1)
disp(sprintf('Interation Number: %02.0f',i));
[fit,ind]=min(Pbest(:,12));
Gbest=Pbest(ind,:);
Best(i,:)=Gbest;
save BestVal Best;
for j = 1 : n,
% Update bird
V_N(j,:) = w*V_O(j,:)+C1*rand*(Pbest(j,index)-Par_O(j,index))+C2*rand*(Gbest(index)-Par_O(j,index));
Par_N(j,index) = Par_O(j,index)+V_N(j,:);
Var = Par_N(j,index);
evaluate;
if JJ < Pbest(j,12),
Pbest(j,:)=[Par_N(j,index),JJ];
end % first if
end % for j=1:n
Par_O = Par_N;
V_O=V_N;
end % for i=1:iter
plot(Best(:,12));
Gbest
end% End NDH 2008
评论0