Hello!
Dear Creators and Users of Matpower
so apologize ! cause my English language is weak!
so !, this is very good that ,I see you are .I am very glad to meet you!
I am student ! and put in the begin of way ! you are professional !
let me ! that give help me?
pleas write comment for any line of your program .
I am so sorry ! that I can not help to you ! but ! if learn as you ! i am
very glad help to you !
Ok ! so ! I beginning
nt = 6;%I don't know why?
mpc = loadcase('MMcase118v5');% load case to mpc I know it
Pg = mpc.gen(:,PG);% put generator from mpc.gen on the Pg
for i = 1:4 %make loop for I think 4 iteration is not it ?
result = MM_118bus_UC(nt,Xwind,Xload,cw,cl,Pg); % I think put many
generator on the network that has 118 bus but I don't know how ?!
if result.results.f>0 % if "f" results grate than zero
res(j) = result; % results put in the J (th) res
Pg = result.results.Pc([1:19],nt);% I dont know result.Pc
cw = cw+nt; % ok !!cw and cl defined ! but I know nt
cl = cl+nt; % cl defined !
j= j+1; % any time doing Loop .so one add to the "j"
else
break
end
end
is possible for you that help to me ?!
thank s alot
my subject study is :
Distributed generation planning .and use TLBO , PSO algorithm for answers
Optimization .
ok !
so tow variables 1- where must be generation on the buses( place of
generator)
2- how many product energy (capacity of generator)
I have tow scope 1- the network parameters be better than last.
2-The who's have power plane ( DG ,PV, Wind, etc ) have Max benefit
so I do this
I study TLBO and PSO .then give ones exercise of them.
In the MATPOWER ues
The create M-file flow the below
clc;
clear all;
%% identify system
mpc = loadcase('case9');% load case9 in the mpc
sizbus = size(mpc.bus);
numbus =sizbus(1,1); %get number fo buses
[idx] = (1:9)';% make index for branch
%% get information of system
% for i=1:1:numbus
%
% end
results = runuopf(mpc); %optimization power flow
[Ploss] = real(get_losses(results));%sperate losses active
indPloss = [idx Ploss];% the matrix that it has index colume and active
losses
[maxploss,idx1] = max(Ploss);% maximon loss active power with index
Qloss = imag(get_losses(results));%sperate losses diactive
indQloss = [idx Qloss];%the matrix that it has index colume and deactive
losses
[maxqloss,idx2] = max(Qloss);%maximon loss deactive power with index
%% chenge system
mpc.gen(4,:) = [4,5,0,5,-5,1,7,1,7,0.5,0,0,0,0,0,0,0,0,0,0,0]% add generator
mpc.gencost(4,:) = [2,300,0,3,0.122500000000000,1,335]% add generator cost
mpc.bus(4,2) = 2; % change bus type
%% get results twovice
results2 = runuopf(mpc);
On Mon, Sep 10, 2018 at 5:55 PM Michal Polecki <[email protected]>
wrote:
> Dear Creators and Users of Matpower,
>
> I have a problem with Unit Commiment problem. I've studied all examples
> attached to Matpower and also I've created some more cases before going to
> my large problem. It seemed to give me correct results. I'm using Gurobi.
> Now I'm trying to run 118 bus case with 13 additional wind farms. What I
> have observed is, that MATPOWER is sticking to the first plan of work (no
> startups or shutdowns though cost of it is 0). What is more is some cases
> it is impossible for MATPOWER to calculate long period e.i nt = 24. I have
> written simple script to cut this long period into 4 6-hour periods. Then
> Matpower was able to calculate each of them. Plot_UC of each 6-hour period
> is attached.
>
> Do you have an idea what could be problem?
>
>
> nt = 6;
> mpc = loadcase('MMcase118v5');
> Pg = mpc.gen(:,PG);
>
> for i = 1:4
> result = MM_118bus_UC(nt,Xwind,Xload,cw,cl,Pg);
> if result.results.f>0
> res(j) = result;
> Pg = result.results.Pc([1:19],nt);
> cw = cw+nt;
> cl = cl+nt;
> j= j+1;
> else
> break
> end
> end
>
> Michał Połecki
>
> --
> Michał Połecki
clc;
clear all;
%% identify system
mpc = loadcase('case9');% load case9 in the mpc
sizbus = size(mpc.bus);
numbus =sizbus(1,1); %get number fo buses
[idx] = (1:9)';% make index for branch
%% get information of system
% for i=1:1:numbus
%
% end
results = runuopf(mpc); %optimization power flow
[Ploss] = real(get_losses(results));%sperate losses active
indPloss = [idx Ploss];% the matrix that it has index colume and active losses
[maxploss,idx1] = max(Ploss);% maximon loss active power with index
Qloss = imag(get_losses(results));%sperate losses diactive
indQloss = [idx Qloss];%the matrix that it has index colume and deactive losses
[maxqloss,idx2] = max(Qloss);%maximon loss deactive power with index
%% chenge system
mpc.gen(4,:) = [4,5,0,5,-5,1,7,1,7,0.5,0,0,0,0,0,0,0,0,0,0,0]% add generator
mpc.gencost(4,:) = [2,300,0,3,0.122500000000000,1,335]% add generator cost
mpc.bus(4,2) = 2; % change bus type
%% get results twovice
results2 = runuopf(mpc);
