节点电压用直角坐标表示时的牛顿拉夫逊法潮流计算（电力系统分析华中科技大学出版社）...

题目描述：

在图所示的电力系统中，网络各元件参数的标幺值如下.......试用牛顿法计算潮流分布。 

 

程序代码

 

clc %清空command windows
clear all %清空workspace
N=4;
NPQ=2;
NPV=1;
K=1;
Kmax=1000 ;
s=10^(-5);
G=[1.042093 -0.588235 0 -0.453858;
-0.588235 1.069005 0 -0.480769;
0 0 0 0 ;
-0.453858 -0.480769 0 0.934627 ];
B=[ -8.242876 2.352941 3.666667 1.891074 ;
2.352941 -4.727377 0 2.403846 ;
3.666667 0 -3.3333333 0 ;
1.891074 2.403846 0 -4.261590 ];
Y=complex(G,B);
Pnode=[ -0.3 -0.55 0.5 ];
Qnode=[ -0.18 -0.13 0 ];
Vnode=[ 0 0 1.10 ];
e=[ 1.0 1.0 1.1 1.05 ];
f=[ 0 0 0 0 ];
%利用for循环来实现多次迭代。
for K=1:Kmax,
%计算△W。
for i=1:2, %PQ