Aha. Yes, that is (or was) a bug and you found the solution. Sorry, I should have remembered this issue and guessed the cause immediately.
I assume you are using MATPOWER 4.0. This bug was fixed in 4.1. Just download the latest version. -- Ray Zimmerman Senior Research Associate 419A Warren Hall, Cornell University, Ithaca, NY 14853 phone: (607) 255-9645 On Feb 16, 2012, at 12:21 PM, Santiago Chamba wrote: > Prof. Zimmerman > > Whit respect H, I think that the error is in get_cost_params.m, specifically > in line 61: > > cp.H = cp.H(idx,:); > > > I think that the cp.H = cp.H(idx,idx), what do you think about this?? This > change generates some errors??? > > Cordially, > > Santiago > De: [email protected] > [mailto:[email protected]] En nombre de Santiago > Chamba > Enviado el: jueves, 16 de febrero de 2012 14:02 > Para: 'MATPOWER discussion forum' > Asunto: RE: Modified toogle- reserves > > I want to coordinate a decentralized DC-OPF with reserves for interconnected > power system. The problem: > > <image001.jpg> > > I add the new variables Tet, what is equal to Va, then add cost alfa for Tet. > > The error is respect to the dimmension of H. For example, for a 4 buses > system and 3 generators, the new reserves variables are 3 and the new > variables for penalizing the angles are 4. In the first iteration for > calculate de cost, H is 7x3, because added the 4 new variables for penalizing > angles. > > I am going to search the error inside, cp = get_cost_params(om, name); > because cp.H=7x3. > > Best Wishes > > Santiago Chamba > > De: [email protected] > [mailto:[email protected]] En nombre de Ray Zimmerman > Enviado el: jueves, 16 de febrero de 2012 13:02 > Para: MATPOWER discussion forum > Asunto: Re: Modified toogle- reserves > > Type ... > > dbstop if error > > ... then re-run your code. When you get the error, check the dimensions of H, > w and Cw to find out where the mismatch is. There is no way for me to suggest > the cause of the error without more details. > > What changes did you make to toggle_reserves? What inputs are you using? ... > etc. > > -- > Ray Zimmerman > Senior Research Associate > 419A Warren Hall, Cornell University, Ithaca, NY 14853 > phone: (607) 255-9645 > > > > > > On Feb 16, 2012, at 8:45 AM, Santiago Chamba wrote: > > > Dear Professor Zimmerman, > > I want to add new variable and linear constraints to the OPF with reserves. > For this cause, I modified toogle_ reserves , but I have problems when the > program calculated to the new cost (compute_cost.m): > > Inner matrix dimensions must agree. > > Error in ==> opf_model.compute_cost at 105 > f = full((w' * H * w) / 2 + Cw' * w); > > The matrix H is define with others dimensions and this case no add the new > cost. > > My new variables are the nodals angles, that are penalizing when not meet > inside specific range > > Please, you explain me the error. > > Best Wishes > > Santiago Chamba > >
