因为没完整的程序，我没看清楚，你的程序种a(k)的意义，我猜想你说的支路类型，可能是k与节点N的关系类型，一般是支路是n指向K，a(k)是1，否则是-1，这么一来，(a(k)^2=1。
后面支路潮流计算中，需要用到支路的方向。但导纳矩阵，不需要计算y(k)的方向。
关于潮流计算的基础理论，你看看：
1,王锡凡，现代电力系统分析，科学出版社。
2.张伯明，高等电力网络分析，清华大学出版社。
3.于尔铿，能量管理系统，科学出版社。
内容全部来源于网络收集，如有侵权，请联系网站删除：QQ：24596024