Info
이 질문은 마감되었습니다. 편집하거나 답변을 올리려면 질문을 다시 여십시오.
optimal power flow by using pso
조회 수: 5 (최근 30일)
이전 댓글 표시
이 질문에 2명의 기여자가 플래그를 지정함
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
댓글 수: 7
Mehmet Cumali Köylü
2021년 1월 4일
Dear sir, I want to the PSO code for optimal location of STATCOM in IEEE30 bus system
답변 (6개)
MANOJ MAHAJAN
2012년 6월 3일
댓글 수: 1
Walter Roberson
2012년 6월 3일
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?"
kishan bhayani
2013년 3월 4일
have you get your answer of economic load dispatch of ieee 30 bus system using pso by using above program????
댓글 수: 0
Stephen Ndubuisi
2017년 10월 6일
Sir pls I need a MATLAB code that will help me minimise transmission loss using PSO.
댓글 수: 4
mohammed hamouda
2019년 6월 26일
please sir
if you finished your project can you send me the code in my email
eng_mohammedhamouda@azhar.edu.eg
TEBBAKH NOUR
2021년 9월 26일
hello please if you have got this code send it to me
nouropport2010@yahoo.com
Nandhini Sankar
2018년 1월 14일
sir please i need a MATLAB code that will help me Rician noise removal method using PSO
댓글 수: 0
mohammed hamouda
2019년 6월 26일
pleas sir
i nead code for optimal planning to penetrate of renewable resource in electrical network (sizing & type) dealing with reliability?
댓글 수: 1
TEBBAKH NOUR
2021년 9월 26일
Hello sir please if you have got this code send it to me
nouropport2010@yahoo.comù
이 질문은 마감되었습니다.
참고 항목
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!