Or ... 3. You can call the makePTDF function which will compute the shift factors for you directly (see section 4.4 in the manual). It essentially does what you suggest in number 1, but using makeBdc (which does not ignore tap) rather than makeB.
-- Ray Zimmerman Senior Research Associate 419A Warren Hall, Cornell University, Ithaca, NY 14853 phone: (607) 255-9645 On Sep 25, 2011, at 9:22 PM, Cuong P Nguyen wrote: > Ray, > > I understand there are two ways of calculating generation shift factor, > > 1. From makeB function, one could calculate the B prime matrix( Bp), then one > can remove the row/column associated with the reference bus. THE TRANSFORMER > TAP IS IGNORED while Bp is formed. Therefore, we can compute the flow on the > lines by injecting 1 MW at a bus i by Bp\dP where dP(i)=1 only at bus i, > zeros everywhere else, this flow is called generation shift factor. > > 2. We repeat the procedure above by setting 1MW injection at bus i, the > injection at all other buses set to 0. Then call rundcpf, we could end up > with the line flow which we is essentially generation shift factor. But the > RUNDCPF DOES CONSIDER TAP. > > Which one is correct in calculating generation shift factor? > > Thanks > Cuong
