You'll have to write the code yourself, but a description is found in section 8.1.2 of my dissertation, where you can find several references as well.
-- Ray Zimmerman Senior Research Associate 419A Warren Hall, Cornell University, Ithaca, NY 14853 phone: (607) 255-9645 On Oct 13, 2012, at 3:13 AM, Zeinab Ghofrani <[email protected]> wrote: > Dear Dr Zimmerman, > > Thanks for your reply. I would appreciate it if you let me know how to find > the code for running Implicit Zbus Gauss method. > > > > Regards > Zeinab Ghofrani > > From: Ray Zimmerman <[email protected]> > To: MATPOWER discussion forum <[email protected]> > Sent: Friday, October 12, 2012 7:28 AM > Subject: Re: convergence of 32-bus distribution network PF > > That is not the Implicit Zbus Gauss method described in my dissertation. That > is the standard Gauss-Seidel method that is already implemented in MATPOWER. > And by the way, the easy way to run that one is … > > opt = mpoption('PF_ALG', 4); > runpf(case_name, opt); > > -- > Ray Zimmerman > Senior Research Associate > 419A Warren Hall, Cornell University, Ithaca, NY 14853 > phone: (607) 255-9645 > > > > > On Oct 10, 2012, at 3:26 AM, Zeinab Ghofrani <[email protected]> wrote: > >> Dear Dr Zimmerman >> >> I implemented the Gauss method using the following codes, but still there is >> convergence problem: >> >> >> case_name='case_dist_32bus'; >> mpc = loadcase(case_name); >> >> ref=[33]; pv=[]; >> pq=1:32; pq=pq'; >> V0=ones(size(mpc.bus,1),1); >> Sbus=makeSbus(mpc.baseMVA, mpc.bus, mpc.gen); >> Ybus=makeYbus(mpc); >> results=gausspf(Ybus, Sbus, V0, ref, pv, pq); >> >> The system matpower data is attached to this email. >> >> Thanks >> Zeinab Ghofrani >> >> From: Ray Zimmerman <[email protected]> >> To: MATPOWER discussion forum <[email protected]> >> Sent: Monday, October 8, 2012 7:49 AM >> Subject: Re: convergence of 32-bus distribution network PF >> >> I wish I could. Unfortunately, the code I used for my thesis was tightly >> tied to some proprietary data formats and it also does not run on current >> versions of Matlab. I you already have the data in MATPOWER format, it would >> probably be rather straightforward to implement the Implicit Zbus Gauss >> method, since MATPOWER can easily build the Ybus matrix for you. That's >> probably the approach I would take if I were you. >> >> -- >> Ray Zimmerman >> Senior Research Associate >> 419A Warren Hall, Cornell University, Ithaca, NY 14853 >> phone: (607) 255-9645 >> >> >> >> >> On Oct 6, 2012, at 1:27 AM, Zeinab Ghofrani <[email protected]> wrote: >> >>> Dear Dr Zimmerman; >>> >>> Thanks for sending me your thesis, which was very useful. Could you please >>> provide me with a reliable MATLAB code for running backward-forward sweep >>> load flow on a radial distribution matpower case? >>> >>> >>> >>> >>> Zeinab Ghofrani >>> >>> From: Ray Zimmerman <[email protected]> >>> To: MATPOWER discussion forum <[email protected]> >>> Sent: Friday, October 5, 2012 6:21 AM >>> Subject: Re: convergence of 32-bus distribution network PF >>> >>> My guess is that it may require a different power flow formulation, namely >>> one that is specifically designed for radial distribution systems. Feel >>> free to look at my dissertation and the references in it for some ideas. >>> >>> -- >>> Ray Zimmerman >>> Senior Research Associate >>> 419A Warren Hall, Cornell University, Ithaca, NY 14853 >>> phone: (607) 255-9645 >>> >>> >>> >>> On Oct 5, 2012, at 2:37 AM, Zeinab Ghofrani <[email protected]> wrote: >>> >>>> Dear all; >>>> >>>> I'm tryning to run power flow on the 32 bus distribution system described >>>> in the following paper: >>>> M. E. Baran and F. F. Wu, “Network Reconfiguration in Distribution Systems >>>> for Loss Reduction and Load Balancing”, IEEE Trans on Power Delivery, >>>> 1989, Vol 4(2), pp >>>> 1401-1407. >>>> I have chnaged the place of the tie switches so that the modified system >>>> has a long latteral with high voltage drop. The power flow doesn't >>>> converge, unless we set Vg of mpc.gen matrix to 1.3. How can I have a >>>> converged load flow when Vg is set to 1? >>>> The data I use is attached to this email. Waiting for your reply. >>>> >>>> Zeinab Ghofrani >>>> <case_dist_32bus.m> >>> >>> >>> >> >> >> >> <case_dist_32bus.m> > > >
