Respuesta autom�tica de Yahoo!
Por este medio quiero informarte que estoy de vacaciones familiares, dejame tu mensaje que pronto te devolveré. Gracias
Re: I want to add a 10MW generator to bus 7
Enviado desde mi dispositivo BlackBerry® de Claro Dominicana -Original Message- From: iman Sender: bounce-38276187-9651...@list.cornell.edu Date: Wed, 16 Nov 2011 19:16:48 To: Reply-To: "MATPOWER discussion forum" Subject: I want to add a 10MW generator to bus 7 Dear All, I have a quick question about MATPOWER. From MATPOWER documentation we know if we want to load the 30-bus system data from case30.m,increase its real power demand at bus 2 to 30 MW, then run an AC optimal power flow with default options, this could be accomplished as follows: >> define_constants; >> mpc = loadcase('case30'); >> mpc.bus(2, PD) = 30; >> runopf(mpc); Now I want to add a 10MW generator to bus 7 and run pf (not opf) for 14 bus system so I write: >> define_constants; >> mpc = loadcase('case14'); >> mpc.gen (7, PG) = 10; //add 10MW to bus 7 >>mpc.bus (7, BUS_TYPE) = 2;// changes the type of bus 7 to PV >> runpf(mpc); but I recieve this error: ??? Subscript indices must either be real positive integers or logicals. Error in ==> ext2int at 230 gs = ( mpc.gen(:, GEN_STATUS) > 0 & ... %% gen status Error in ==> runpf at 131 mpc = ext2int(mpc); Do u know what is wrong with my code? Thanks --part2265-boundary-861242848-228417209 Content-Transfer-Encoding: base64 Content-Type: text/html; charset="utf-8" Enviado desde mi dispositivo BlackBerry® de Claro DominicanaFrom: imanSender: bounce-38276187-9651...@list.cornell.edu Date: Wed, 16 Nov 2011 19:16:48 +To: ReplyTo: "MATPOWER discussion forum" Subject: I want to add a 10MW generator to bus 7Dear All,I have a quick question about MATPOWER.From MATPOWER documentation we know if we want to load the 30-bus system data from case30.m,increase its real power demand at bus 2 to 30 MW, then run an AC optimal power flow with default options, this could be accomplished as follows: >> define_constants;>> mpc = loadcase('case30');>> mpc.bus(2, PD) = 30;>> runopf(mpc); Now I want to add a 10MW generator to bus 7 and run pf (not opf) for 14 bus system so I write: >> define_constants; >> mpc = loadcase('case14'); >> mpc.gen (7, PG) = 10; //add 10MW to bus 7>>mpc.bus (7, BUS_TYPE) = 2; // changes the type of bus 7 to PV >> runpf(mpc);but I recieve this error:??? Subscript indices must either be real positive integers or logicals.Error in ==> ext2int at 230 gs = ( mpc.gen(:, GEN_STATUS) > 0 & ... %% gen status Error in ==> runpf at 131mpc = ext2int(mpc);Do u know what is wrong with my code?Thanks --part
Re: Power flow results question
Bar Enviado desde mi dispositivo BlackBerry® de Claro Dominicana -Original Message- From: Sam Hazim Sender: bounce-49184047-9651...@list.cornell.edu Date: Fri, 20 Apr 2012 13:07:04 To: MATPOWER discussion forum Reply-To: "MATPOWER discussion forum" Subject: Re: Power flow results question Thank you so much for your reply, If i put them under BS column. What about the cap banks cost function, i can't put them under the gencost,, where can i put cap cost function? Thank you, From: Ray Zimmerman To: MATPOWER Discussion List Sent: Saturday, 21 April 2012 6:26 AM Subject: Re: Power flow results question Capacitors are constant impedance elements that should be added to the BS column of the bus matrix, not to the gen matrix. Ray On Apr 20, 2012, at 1:33 AM, Sam Hazim wrote: > > > > > > > > >Hi everyone, >I am new to Matpower, i have some few questions about power flow . > >I have been given in my assignment some possible locations for capacitor banks > to be placed on the system. > >1: I want to make sure, to include any cap bank on any bus, is it correct be >placed within generator data matrix and set Pg,Qg,Qmin =0, vg=1, Qgmax= 300 >(value given in my assignment). >2: from my knowledge, we use cap banks order to minimize losses or support >voltage.i tried to connect them to the high voltage bus or to the main >transform, or other locations and then runpf. the answer doesn't change at all >( Why is that). i have got min voltage at one bus is about 0.831 and it should >be between a recommend range 0.95 to 1.05. > >I have included the case file for two modified cases.and runpf , the results >are the same for both of them. > >Please, any suggestions will be very helpful for me to tackle this problem. > >Raad > >functionmpc = case18%%- Power Flow Data -%% >%% system MVA basempc.baseMVA = 500;%% bus data >% bus_i type Pd Qd Gs Bs area Vm Va baseKV zone Vmax Vminmpc.bus = [ >1 3 0 0 0 0 1 1 0 220 1 1.05 0.95; >2 1 0 0 0 0 1 1 0 220 1 1.05 0.95; >3 2 0 0 0 0 1 1 0 220 1 1.05 0.95; >4 1 0 0 0 0 1 1 0 220 1 1.05 0.95; >5 2 0 0 0 0 1 1 0 220 1 1.05 0.95; >6 1 200 50 0 0 1 1 0 220 1 1.05 0.95; >7 1 450 115 0 0 1 1 0 220 1 1.05 0.95; >8 2 0 0 0 0 1 1 0 220 1 1.05 0.95; >9 1 0 0 0 0 1 1 0 220 1 1.05 0.95; >10 1 0 0 0 0 1 1 0 220 1 1.05 0.95; >11 1 0 0 0 0 1 1 0 220 1 1.05 0.95; >12 1 200 50 0 0 1 1 0 110 1 1.05 0.95; >13 1 150 40 0 0 2 1 0 110 1 1.05 0.95; >14 1 250 65 0 0 2 1 0 110 1 1.05 0.95; >15 1 750 190 0 0 2 1 0 110 1 1.05 0.95; >16 2 0 0 0 0 2 1 0 33 1 1.05 0.95; >17 1 0 0 0 0 2 1 0 220 1 1.05 0.95; >18 2 0 0 0 0 2 1 0 110 1 1.05 0.95; >];%% generator data >% bus Pg Qg Qmax Qmin Vg mBase status Pmax Pmin Pc1 Pc2 Qc1min Qc1max Qc2min >Qc2max ramp_agc ramp_10 ramp_30 ramp_q apfmpc.gen = [ >1 500 0 750 -200 1.05 500 1 1250 0 0 0 0 0 0 0 0 0 0 0 0; >3 450 0 50 -150 1.05 500 1 500 0 0 0 0 0 0 0 0 0 0 0 0; >5 550 0 200 -150 1.028 500 1 600 0 0 0 0 0 0 0 0 0 0 0 0; >8 50 0 65 -200 0.95 500 1 120 0 0 0 0 0 0 0 0 0 0 0 0; >16 400 0 300 -40 1.022 500 1 500 0 0 0 0 0 0 0 0 0 0 0 0; >18 250 0 220 -100 0.99 500 1 350 0 0 0 0 0 0 0 0 0 0 0 0; >9 0 0 300 0 1 500 1 0 0 0 0 0 0 0 0 0 0 0 0 0; >12 0 0 300 0 1 500 1 0 0 0 0 0 0 0 0 0 0 0 0 0; >13 0 0 300 0 1 500 1 0 0 0 0 0 0 0 0 0 0 0 0 0; >14 0 0 300 0 1 500 1 0 0 0 0 0 0 0 0 0 0 0 0 0; >15 0 0 300 0 1 500 1 0 0 0 0 0 0 0 0 0 0 0 0 0; >];%% branch data >% fbus tbus r x b rateA rateB rateC ratio angle status mpc.branch = [ >1 3 0.0712 0.3595 0.0982 524 0 0 0 0 1; >4 5 0.0978 0.5855 0.0773 455 0 0 0 0 1; >4 6 0.023 0.1615 0.01 671 0 0 0 0 1; >4 7 0.0201 0.1548 0.0395 766 0 0 0 0 1; >4 2 0.0112 0.1136 0.0071 766 0 0 0 0 1; >7 2 0.0112 0.1136 0.0071 766 0 0 0 0 1; >7 9 0.01 0.0954 0.0058 915 0 0 0 0 1; >1 6 0.042 0.2715 0.017 493 0 0 0 0 1; >1 7 0.088 0.5675 0.0356 323 0 0 0 0 1; >1 17 0.0444 0.2104 0.0247 323 0 0 0 0 1; >7 17 0.0444 0.2104 0.0247 323 0 0 0 0 1; >7 18 0. 0.08 0.00 2000 0 0 1 0 1; >7 8 0.002 0.0175 0.0011 915 0 0 0 0 1; >7 15 0. 0.05 0.00 2000 0 0 1 0 1; >8 9 0.009 0.0775 0.0047 915 0 0 0 0 1; >9 10 0.0095 0.078 0.005 827 0 0 0 0 1; >9 11 0.079 0.4815 0.0191 408 0 0 0 0 1; >9 14 0.00 0.15 0.00 2000 0 0 1 0 1; >10 11 0.044 0.3585 0.0233 456 0 0 0 0 1; >10 12 0. 0.15 0.00 2000 0 0 1 0 1; >11 13 0.0 0.30 0.00 2000 0 0 1 0 1; >4 16 0.00 0.131 0.00 2000 0 0 1 0 1; >3 5 0.036 0.2396 0.0667 870 0 0 0 0 1; >];%%- OPF Data -%% >%% generator cost data >% 1 startup shutdown n x1 y1 ... xn yn >% 2 startup shutdown n c(n-1) ... c0mpc.gencost = [ >2 0 0 3 0.039 20 0; >2 0 0 3 0.029 20 0; >2 0 0 3 0.031 20 0; >2 0 0 3 0.12 20 0; >2 0 0 3 0.25 20 0; >2 0 0 3 0.1 20 0; >2 0.01 0 3 0.01 40 0; >2 0.01 0 3 0.01 40 0; >2 0.01 0 3 0.01 40 0; >2 0.01 0 3 0.01 40 0; >2 0.01 0 3 0.01 40 0; >]; > >This is the modifed case > >functionmpc = case18%%- Power Flow Data -%% >%% system MVA basempc.baseMVA = 500;%% bus data >% bus_i type Pd Qd Gs Bs area Vm Va baseKV zone Vmax Vm
Re: Power flow results question
F Enviado desde mi dispositivo BlackBerry® de Claro Dominicana -Original Message- From: Sam Hazim Sender: bounce-49184047-9651...@list.cornell.edu Date: Fri, 20 Apr 2012 13:07:04 To: MATPOWER discussion forum Reply-To: "MATPOWER discussion forum" Subject: Re: Power flow results question Thank you so much for your reply, If i put them under BS column. What about the cap banks cost function, i can't put them under the gencost,, where can i put cap cost function? Thank you, From: Ray Zimmerman To: MATPOWER Discussion List Sent: Saturday, 21 April 2012 6:26 AM Subject: Re: Power flow results question Capacitors are constant impedance elements that should be added to the BS column of the bus matrix, not to the gen matrix. Ray On Apr 20, 2012, at 1:33 AM, Sam Hazim wrote: > > > > > > > > >Hi everyone, >I am new to Matpower, i have some few questions about power flow . > >I have been given in my assignment some possible locations for capacitor banks > to be placed on the system. > >1: I want to make sure, to include any cap bank on any bus, is it correct be >placed within generator data matrix and set Pg,Qg,Qmin =0, vg=1, Qgmax= 300 >(value given in my assignment). >2: from my knowledge, we use cap banks order to minimize losses or support >voltage.i tried to connect them to the high voltage bus or to the main >transform, or other locations and then runpf. the answer doesn't change at all >( Why is that). i have got min voltage at one bus is about 0.831 and it should >be between a recommend range 0.95 to 1.05. > >I have included the case file for two modified cases.and runpf , the results >are the same for both of them. > >Please, any suggestions will be very helpful for me to tackle this problem. > >Raad > >functionmpc = case18%%- Power Flow Data -%% >%% system MVA basempc.baseMVA = 500;%% bus data >% bus_i type Pd Qd Gs Bs area Vm Va baseKV zone Vmax Vminmpc.bus = [ >1 3 0 0 0 0 1 1 0 220 1 1.05 0.95; >2 1 0 0 0 0 1 1 0 220 1 1.05 0.95; >3 2 0 0 0 0 1 1 0 220 1 1.05 0.95; >4 1 0 0 0 0 1 1 0 220 1 1.05 0.95; >5 2 0 0 0 0 1 1 0 220 1 1.05 0.95; >6 1 200 50 0 0 1 1 0 220 1 1.05 0.95; >7 1 450 115 0 0 1 1 0 220 1 1.05 0.95; >8 2 0 0 0 0 1 1 0 220 1 1.05 0.95; >9 1 0 0 0 0 1 1 0 220 1 1.05 0.95; >10 1 0 0 0 0 1 1 0 220 1 1.05 0.95; >11 1 0 0 0 0 1 1 0 220 1 1.05 0.95; >12 1 200 50 0 0 1 1 0 110 1 1.05 0.95; >13 1 150 40 0 0 2 1 0 110 1 1.05 0.95; >14 1 250 65 0 0 2 1 0 110 1 1.05 0.95; >15 1 750 190 0 0 2 1 0 110 1 1.05 0.95; >16 2 0 0 0 0 2 1 0 33 1 1.05 0.95; >17 1 0 0 0 0 2 1 0 220 1 1.05 0.95; >18 2 0 0 0 0 2 1 0 110 1 1.05 0.95; >];%% generator data >% bus Pg Qg Qmax Qmin Vg mBase status Pmax Pmin Pc1 Pc2 Qc1min Qc1max Qc2min >Qc2max ramp_agc ramp_10 ramp_30 ramp_q apfmpc.gen = [ >1 500 0 750 -200 1.05 500 1 1250 0 0 0 0 0 0 0 0 0 0 0 0; >3 450 0 50 -150 1.05 500 1 500 0 0 0 0 0 0 0 0 0 0 0 0; >5 550 0 200 -150 1.028 500 1 600 0 0 0 0 0 0 0 0 0 0 0 0; >8 50 0 65 -200 0.95 500 1 120 0 0 0 0 0 0 0 0 0 0 0 0; >16 400 0 300 -40 1.022 500 1 500 0 0 0 0 0 0 0 0 0 0 0 0; >18 250 0 220 -100 0.99 500 1 350 0 0 0 0 0 0 0 0 0 0 0 0; >9 0 0 300 0 1 500 1 0 0 0 0 0 0 0 0 0 0 0 0 0; >12 0 0 300 0 1 500 1 0 0 0 0 0 0 0 0 0 0 0 0 0; >13 0 0 300 0 1 500 1 0 0 0 0 0 0 0 0 0 0 0 0 0; >14 0 0 300 0 1 500 1 0 0 0 0 0 0 0 0 0 0 0 0 0; >15 0 0 300 0 1 500 1 0 0 0 0 0 0 0 0 0 0 0 0 0; >];%% branch data >% fbus tbus r x b rateA rateB rateC ratio angle status mpc.branch = [ >1 3 0.0712 0.3595 0.0982 524 0 0 0 0 1; >4 5 0.0978 0.5855 0.0773 455 0 0 0 0 1; >4 6 0.023 0.1615 0.01 671 0 0 0 0 1; >4 7 0.0201 0.1548 0.0395 766 0 0 0 0 1; >4 2 0.0112 0.1136 0.0071 766 0 0 0 0 1; >7 2 0.0112 0.1136 0.0071 766 0 0 0 0 1; >7 9 0.01 0.0954 0.0058 915 0 0 0 0 1; >1 6 0.042 0.2715 0.017 493 0 0 0 0 1; >1 7 0.088 0.5675 0.0356 323 0 0 0 0 1; >1 17 0.0444 0.2104 0.0247 323 0 0 0 0 1; >7 17 0.0444 0.2104 0.0247 323 0 0 0 0 1; >7 18 0. 0.08 0.00 2000 0 0 1 0 1; >7 8 0.002 0.0175 0.0011 915 0 0 0 0 1; >7 15 0. 0.05 0.00 2000 0 0 1 0 1; >8 9 0.009 0.0775 0.0047 915 0 0 0 0 1; >9 10 0.0095 0.078 0.005 827 0 0 0 0 1; >9 11 0.079 0.4815 0.0191 408 0 0 0 0 1; >9 14 0.00 0.15 0.00 2000 0 0 1 0 1; >10 11 0.044 0.3585 0.0233 456 0 0 0 0 1; >10 12 0. 0.15 0.00 2000 0 0 1 0 1; >11 13 0.0 0.30 0.00 2000 0 0 1 0 1; >4 16 0.00 0.131 0.00 2000 0 0 1 0 1; >3 5 0.036 0.2396 0.0667 870 0 0 0 0 1; >];%%- OPF Data -%% >%% generator cost data >% 1 startup shutdown n x1 y1 ... xn yn >% 2 startup shutdown n c(n-1) ... c0mpc.gencost = [ >2 0 0 3 0.039 20 0; >2 0 0 3 0.029 20 0; >2 0 0 3 0.031 20 0; >2 0 0 3 0.12 20 0; >2 0 0 3 0.25 20 0; >2 0 0 3 0.1 20 0; >2 0.01 0 3 0.01 40 0; >2 0.01 0 3 0.01 40 0; >2 0.01 0 3 0.01 40 0; >2 0.01 0 3 0.01 40 0; >2 0.01 0 3 0.01 40 0; >]; > >This is the modifed case > >functionmpc = case18%%- Power Flow Data -%% >%% system MVA basempc.baseMVA = 500;%% bus data >% bus_i type Pd Qd Gs Bs area Vm Va baseKV zone Vmax Vmi