Dear Mr. Abhyankar .G., I'm not sure about it and I'm not getting correct way how to extract the elements from solution vector you are talking about. But what I've done is- 1) I didn't put tap ratio and shift values in branch matrix. I have created another matrix data for transformer like mpc.tapdata in matpower case file (case9.m) and assigned values with lower and upper bounds for both ratio and shift. Also created idx_tapdata where I have defined constants for named column. 2) Added these extra variables and their constraints in opf_setup by add _vars and add_constraints. (Total variables were 24 and now they are 42 in Jacobian matrix verified by spy function.) 3) Created new function where I get ratio and shift values.
Problem what I'm facing is how to extract values of ratio and shift from solver to update the Y bus in opf_hessfcn and opf_consfcn ? But facing the problem in opf_args. Thanks . Mirish On Fri, Jan 9, 2015 at 3:05 PM, Abhyankar, Shrirang G. <[email protected]> wrote: > Are you setting the solution vector (x) size and the indexing > information for extracting the elements from it correctly? > > On Jan 8, 2015, at 4:50 AM, "Mirish Thakur" <[email protected]> > wrote: > > Dear Prof. R. Zimmerman, > > I'm master student working on a project dealing with power system > optimization with additional variables 'tap ratio' and 'shift' on 9 bus > system by *mipsopf_solver*. I have made necessary changes in respective > function files. Also I concern with only current flow constraints in > *opf_hessfcn*. For that I used command like *>>* *mpopt = > mpoption('OPF_FLOW_LIM', 2);* and* >> runopf(case9,mpopt);*. But I'm > getting errors - > > Did not converge in 17 iterations. > Error using opf_args (line 315) > opf_args.m: Incorrect input parameter order, number or type > > Error in opf (line 199) > [mpc, mpopt] = opf_args(varargin{:}); > > Error in opf (line 262) > opf(varargin) > > Error in runopf (line 96) > [r, success] = opf(casedata, mpopt); > > As I'm not much faliliar with MATPOWER so I could not understand what > exactly I need to rectify in *opf_args.m* file. I need your help .Thanks > in advance. > > > Sincerely > Mirish > >
