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