optimal power flow by using pso

hi this is my program . but i want to change my program for power flow .i want to compare ieee 30 bus result ,from newton raphson method

please help me ......

i am waiting your result....

clc;
clear all;
close all;

ipf=fopen('opres75n.in','r');
opf=fopen('opres75.out','w+');

%reading data 
no_part=fscanf(ipf,'%f',1);
no_var=fscanf(ipf,'%f',1);
no_units=fscanf(ipf,'%f',1);
itermax=fscanf(ipf,'%f',1);
temp=fscanf(ipf,'%f',[5,no_units]);
temp=temp';
c=temp(:,1);
b=temp(:,2);
a=temp(:,3);
pmin=temp(:,5);
pmax=temp(:,4);

bng=fscanf(ipf,'%f',no_var);
pg=fscanf(ipf,'%f',no_units);
pbest=zeros(no_part,1);
gbest=0;
part_v=zeros(no_part,no_var);
wmin=0.4; wmax=0.9;
c1=2; c2=2;
vmax=(pmax-pmin)/50;

% Random Generation of particles
% for i=1:no_part
%     for j=1:no_var
%         ng=bng(j);
%         part(i,j)=pmin(ng)+(pmax(ng)-pmin(ng))*rand-pg(ng);
%     end
% end

i=1;
while i~=no_part+1
    for j=1:no_var
        ng=bng(j);
        part(i,j)=pmin(ng)+(pmax(ng)-pmin(ng))*rand-pg(ng);
    end
    sum_part=sum(part(i,:));
    i=i+1;
    if abs(sum_part)>=0.1
        i=i-1;
    end
    if i==0
        i=1;
    end
end

for iter=1:itermax
    i=1;
    tic;
    while i~=no_part+1
        sum_part=sum(part(i,:));
        if abs(sum_part)>=0.1
            for j=1:no_var
                ng=bng(j);
                part(i,j)=pmin(ng)+(pmax(ng)-pmin(ng))*rand-pg(ng);
            end
            i=i-1;
        else
            i=i+1;
        end
        if i==0
            i=1;
        end
    end

      % Binding Fitness calculation
      bf=zeros(no_part,no_var);
      dpmin=pg-pmin;
      dpmax=pmax-pg;
      i=1;
      while i~=no_part+1
          for j=1:no_var
              ng=bng(j);
              if part(i,j)<dpmin(ng)
                  bf(i,j)=dpmin(ng)-part(i,j);
              end
              if part(i,j)>dpmax(ng)
                  bf(i,j)=part(i,j)-dpmax(ng);
              end
          end
          k=0;
          for j=1:no_var
              if bf(i,j)==0
                  part(i,j)=pmin(ng)+(pmax(ng)-pmin(ng))*rand-pg(ng);
                  k=1;
              end
          end
          if k==1
              i=i-1;
          end
          i=i+1;
      end

      % Optimal Fitness Calculation
      cost=zeros(no_part,1);
      for i=1:no_part
          for j=1:no_var
              ng=bng(j);
              cost(i)=cost(i)+(a(ng)+b(ng)*part(i,j)+c(ng)*part(i,j)*part(i,j))*part(i,j);
          end
      end

      % PBEST
      if iter==1
          pbest=cost;
          part_pb=part;
      else
          for i=1:no_part
              if pbest(i)>cost(i)
                  pbest(i)=cost(i);
                  part_pb(i,:)=part(i,:);
              end
          end
      end

      %GBEST
      gbest=min(pbest);
      for i=1:no_part
          if gbest==cost(i)
              part_gb=part(i,:);
          end
      end

      % Velocity
      w=wmax-((wmax-wmin)/itermax)*iter;
      for i=1:no_part
          for j=1:no_var
              ng=bng(j);
              part_v(i,j)=w*part_v(i,j)+c1*rand*(part_pb(i,j)-part(i,j))+c2*rand*(part_gb(j)-part(i,j));    
              if part_v(i,j)>vmax(ng)
                  part_v(i,j)=vmax(ng);
              end
              if part_v(i,j)<-vmax(ng)
                  part_v(i,j)=-vmax(ng);
              end
          end
      end

      % Position Update
      for i=1:no_part
          for j=1:no_var
              ng=bng(j);
              part(i,j)=part(i,j)+part_v(i,j);
              if part(i,j)>dpmax(ng)
                  part(i,j)=dpmax(ng);
              end
              if part(i,j)<dpmin(ng)
                  part(i,j)=dpmin(ng);
              end
          end
      end
      time(iter)=toc;
      fitmin(iter)=min(cost);
      costa(iter)=mean(cost);
      iter
  %     max(max(bf))
  %     min(min(bf))
  end
  subplot(2,1,1);
  plot(1:iter,costa);
  xlabel('Iterations');
  ylabel('Cost');
  subplot(2,1,2);
  plot(1:iter,time);
  xlabel('Iterations');
  ylabel('Time');

for j=1:no_var
    for k=1:no_units
        ng=bng(j);
        if ng==k
            pg(ng)=pg(ng)+part_gb(j);
            if pg(ng)<pmin(ng)
                pg(ng)=pmin(ng);
            end
            if pg(ng)>pmax(ng)
                pg(ng)=pmax(ng);
            end
        end
    end
end
fc=0;
for k=1:no_units
    fc=fc+a(k)+b(k)*pg(k)+c(k)*pg(k)*pg(k);
end
s=0;
for i=1:6
    s=s+pg(i);
end

[Walter Roberson]

You did not ask a question. You presented code (that appears to implement Particle Swarm) and you indicated a general desire ("I want to change my program for power flow"), but you did not ask a question. Please be clearer as to what you are asking.

[Walter Roberson]

duplicate is at http://www.mathworks.com/matlabcentral/answers/38336-power-flow

[Chirag Tanti]

Dear sir, I want to the PSO code for optimal location of STATCOM in IEEE30 bus system

[Stephen Ndubuisi]

Pls, I need a matlab code I can use to simulate the transmission loss equation using PSO. How can I apply the constraints to the PSO programming Syntex

[Akhil Madarapu]

Hi Manoj sir, You had given a code for optimal power flow by using pso & I tried a lot to execute it in Matlab but it always shows an error... could u help me please

[Muhammad Rizaldy]

Where is your data input sir?

[Mehmet Cumali Köylü]

Dear sir, I want to the PSO code for optimal location of STATCOM in IEEE30 bus system

Answer 1

[MANOJ MAHAJAN]

sir, i apply this program for economic load dispatch ,it is run but the problem is that i want result for power flow on bus and real power loss for different ieee bus

[Walter Roberson]

This is a statement:

"I found instructions for trimming the edge of my garden with electric clippers.I want to cut my lawn with an gasoline lawnmower."

And *this* is a question:

"How far out should I pull the choke when I am starting a Banana Junior 9000 gasoline lawnmower?"

Answer 2

[MANOJ MAHAJAN]

sir ,there is any tool box to solve the problem of optimal power flow by using pso like a economic load dispatch

[Walter Roberson]

Not that I know of.

http://www.mathworks.com/matlabcentral/fileexchange/28291-particle-swarm-optimization-research-toolbox

Answer 3

[kishan bhayani]

have you get your answer of economic load dispatch of ieee 30 bus system using pso by using above program????

Answer 4

[Stephen Ndubuisi]

Sir pls I need a MATLAB code that will help me minimise transmission loss using PSO.

[Walter Roberson]

Look in the File Exchange to see what is available. If you do not see what you need there, then you will have to write the code yourself (or hire someone to write it.)

[Akhil Madarapu]

Sir please could you give me the code for " transmission loss minimization using PSO" It's my Btech final year project & I have to submit the code and execute it within 2 days. Help me please........

[mohammed hamouda]

please sir

if you finished your project can you send me the code in my email

[e-mail]

[TEBBAKH NOUR]

hello please if you have got this code send it to me

[e-mail]

Answer 5

[Nandhini Sankar]

sir please i need a MATLAB code that will help me Rician noise removal method using PSO

Answer 6

[mohammed hamouda]

pleas sir

i nead code for optimal planning to penetrate of renewable resource in electrical network (sizing & type) dealing with reliability?

[TEBBAKH NOUR]

Hello sir please if you have got this code send it to me

[e-mail]ù