Have you looked at the examples in Appendix A in the MATPOWER User's Manual<http://www.pserc.cornell.edu/matpower/manual.pdf>? You can also find a few more examples in t/t_mips.m. It sounds like you have already implemented the objective function and constraint evaluation function for you problem, which includes the gradients. So, the only thing missing is a function for efficiently evaluating the Hessian of the Lagrangian function for your problem. The opf_hessian function computes this for MATPOWER's OPF and the derivation of the formulas used can be found in the MATPOWER Technical Note 2, AC Power Flows, Generalized OPF Costs and their Derivatives using Complex Matrix Notation<http://www.pserc.cornell.edu/matpower/TN2-OPF-Derivatives.pdf>. You may find it helpful in the derivation of the Hessian for your problem.
Feel free to ask if you have specific implementation questions. -- Ray Zimmerman Senior Research Associate 211 Warren Hall, Cornell University, Ithaca, NY 14853 phone: (607) 255-9645 On May 31, 2010, at 5:45 PM, b v wrote: Hello all, I have a project to do OPF for distribution systems. The objective function is to minimizing the Power loss. I have done it already by MATLAB using fmincon function. However, for using mips solver in Matpower and taking the hessian matrix from the opf_hessfcn function I would like to do it in Matpower too. But I do not have any idea how to do it and if it is possible at all. If it works then we can use the matpower for doing OPF in distribution systems as well. Can anyone please explain me simply how I can use Matpower to optimize my own objective function respect to my own constraint. I should add that my data file is exactly same as casedata format in Matpower and just the obj. function and constraints are different. Many Thanks in advance. Best Regards
