function DC = PSO(V,I)
...
DC=0;D=zeros(1,3);P=zeros(1,3);
...
**********************************
iter = 0 ;
while ( iter < bird_step )
iter = iter + 1;
for i = 1:3,
P(i)=V*I;
D(i)=current_dutycycl(i);
DC=D(i);
pause(0.1)
end
current_power=P;
for i = 1 : 3
if current_power(i) > local_best_power(i)
local_best_power(i) = current_power(i) ;
local_best_dutycycl(:,i) = current_dutycycl(:,i) ;
end
end
[current_global_best_power,g] = max(local_best_power);
if current_global_best_power > global_best_power
global_best_power = current_global_best_power;
for i=1:3
globl_best_dutycycl(:,i) = local_best_dutycycl(:,g);
end
end
velocity = w *velocity + c1*(R1.*(local_best_dutycycl-current_dutycycl)) + c2*(R2.*(globl_best_dutycycl-current_dutycycl));
current_dutycycl = current_dutycycl + velocity ;
end
[Peak_power,I] = max(current_power)
current_dutycycl(:,I)
DC=current_dutycycl(:,I);
end
