Dear Dr. Ray, Good morning. I found that Y matrix from mathpower and it manual: Y = Cf'*Yf + Ct'*Yt+Ysh values are not same! They are same when their is no no shunt admittance, i.e, Y = Cf'*Yf + Ct'*Yt. Why this happens?
As a result, the active/reactive power injection: conj(diag((Y) * v*v')) is different from Cf'*Pf+Ct'*Pt+Ysh (active power for example, Pf and Pt are power flow from and to). When their is no shunt admittance both are equal but the problem comes when there is shunt admittance. Please help: how can I get theoretical active power injection real(conj(diag((Y) * v*v'))) is equal to calculated Cf'*Pf+Ct'*Pt+Ysh values when there is shunt admittance? Thanks. Regards, Rana ________________________________ From: Md Masud Rana Sent: Tuesday, 12 November 2019 4:06 PM To: [email protected] <[email protected]>; Ray Zimmerman <[email protected]> Subject: About Y-bus values from Matpower and theoretical values. Thanks. Dear Dr. Ray, Good morning. I found that Y matrix from mathpower and it manual: Y = Cf'*Yf + Ct'*Yt+Ysh values are not same! They are same when their is no no shunt admittance, i.e, Y = Cf'*Yf + Ct'*Yt. Why this happens? As a result, the active/reactive power injection: conj(diag((Y) * v*v')) is different from Cf'*Pf+Ct'*Pt+Ysh (active power for example, Pf and Pt are power flow from and to). When their is no shunt admittance both are equal but the problem comes when there is shunt admittance. Please help: how can I get theoretical active power injection real(conj(diag((Y) * v*v'))) is equal to calculated Cf'*Pf+Ct'*Pt+Ysh values when there is shunt admittance? Thanks. Regards, Rana
