Shri, I think you are right. There may be some problems with case3375wp.
sprank(J) = 6355 < 6357, the size of J. So the Jacobian matrix is surely singular. As far as I know, Matpower doesn’t apply a topological check before PF calculation. So if some parts of the system are actually isolated but not assigned the bus type “4” (none) to, Jacobian matrix will have all-zero rows and columns, then singular. For case3375wp, bus 591, 3574 (inner index) are actually isolated. You could find their external index; assigned them bus type 4; then run the case again. I think Dr. Zimmerman may correct this mistake in the next version. Shiyang Li -----邮件原件----- 发件人: [email protected] [mailto:[email protected]] 代表 Shrirang Abhyankar 发送时间: 2012年5月16日 23:57 收件人: MATPOWER discussion forum 主题: case3375wp There seems to be something amiss with the case3375wp data. Matlab complains of singular Jacobian matrix. mpopt= mpoption; mpoption(31) = 2; runpf('case3375wp',mpopt,'pfresults'); (printpf is commented out) it max P & Q mismatch (p.u.) ---- --------------------------- 0 2.442e+01Warning: Matrix is singular to working precision. > In newtonpf at 125 In runpf at 224 1 NaNWarning: Matrix is singular to working precision. > In newtonpf at 125 In runpf at 224 2 NaNWarning: Matrix is singular to working precision. > In newtonpf at 125 In runpf at 224 3 NaNWarning: Matrix is singular to working precision. > In newtonpf at 125 In runpf at 224 4 NaNWarning: Matrix is singular to working precision. > In newtonpf at 125 In runpf at 224 5 NaNWarning: Matrix is singular to working precision. > In newtonpf at 125 In runpf at 224 6 NaNWarning: Matrix is singular to working precision. > In newtonpf at 125 In runpf at 224 7 NaNWarning: Matrix is singular to working precision. > In newtonpf at 125 In runpf at 224 8 NaNWarning: Matrix is singular to working precision. > In newtonpf at 125 In runpf at 224 9 NaNWarning: Matrix is singular to working precision. > In newtonpf at 125 In runpf at 224 10 NaN Newton's method power flow did not converge in 10 iterations. Thanks, Shri
