关于电力系统潮流计算牛顿拉夫逊算法的MATLAB程序中 形成节点导纳矩阵时有个语句理解不了,Ybus(n,n) = Ybus(n,n)+y(k)/(a(k)^2) + Bc(k);一句中y(k)是线路导纳a(k)为线路类型 他们的商又是什么呢 源程

Question

关于电力系统潮流计算牛顿拉夫逊算法的MATLAB程序中 形成节点导纳矩阵时有个语句理解不了,
Ybus(n,n) = Ybus(n,n)+y(k)/(a(k)^2) + Bc(k);一句中y(k)是线路导纳a(k)为线路类型 他们的商又是什么呢 源程序相关内容 % This program obtains th Bus Admittance Matrix for power flow solution j=sqrt(-1); i = sqrt(-1); nl = linedata(:,1); nr = linedata(:,2); R = linedata(:,3); X = linedata(:,4); Bc = j*linedata(:,5); a = linedata(:,6); nbr=length(linedata(:,1)); nbus = max(max(nl),max(nr)); Z = R

Answer 1

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