文档介绍：5电力系统分析
集中了以上各种类型元件的简单网络如图
节点注入电流和节点电压构成以下线性方程组
其中
可展开为如下形式：
（）
若
可展开如下形式：
（）
式中n为网络节点数
中系数矩阵（雅可比矩阵）的各元素
解修正方程式（），求出修正量
修正各节点电压
将 再代入（），（）式
求
校验是否收敛，即
如果收敛，迭代到此结束，进一步计算各线路潮流和平衡节点功率，并打印输出结果。如果不收敛， 转回第（2）步进行下一次迭代计算，直到收敛为止。
牛顿法的框图及求解过程
试用牛顿－。
网络图
解：
⑴ n－节点数、nl－支路数、isb－平衡母线节点号（固定为1）、pr－误差精度、如果要输入 则输入eps即可。
⑵ 请输入由支路参数形成的矩阵B1
矩阵B1的每行是由下列参数构成的：
某支路的首端号P；
某支路末端号Q，且P<Q；
支路的阻抗（R+jX）；
支路的对地容抗；
支路的变比K；
折算到哪一侧的标志（如果支路的首端P处于高压侧则请输入“1”，否则请输入“0”）。
⑶ 请输入各节点参数形成的矩阵B2
矩阵B2的每行是由下列参数构成的：
节点所接发电机的功率SG；
节点负荷的功率SL；
节点电压的初始值。
PU节点电压U的给定值。
节点所接的无功补偿设备的容量。
节点分类标号igl。
1 －－平衡节点
igl= 2－－PQ节点
3－－PU节点
（4）请输入由节点号及其容抗形成的矩阵X
。
（）和（）求出修正方程式的常数项向量
（）和（）求出雅可比矩阵各元素值，即可得到第一次迭代时的修正方程式。
（），求
，即得出第一次迭代后各节点的电压值。
，当收敛精度取
（即PR=）时，需要进行四次迭代。求出了各节点电压后，即可求各支路的潮流分布。
：
请输入节点数:n=5
请输入支路数:nl=5
请输入平衡母线节点号:isb=1
请输入误差精度:pr=
请输入由支路参数形成的矩阵:B1=[1 2 0 0;
2 3 + 1 0;2 4 + 0 1 0;
3 4 + 1 0;3 5 0 1]
请输入各节点参数形成的矩阵:B2=[0 0 0 1;
0 + 1 0 0 2;0 2+1i 1 0 0 2;0 + 1 0 0 2;
5 0 0 3]
请输入由节点号及其对地阻抗形成的矩阵:X=[1 0;2 0;3 0;
4 0;5 0]
结果如下：
迭代次数 5
没有达到精度要求的个数 7 8 8 6 0
各节点的实际电压标么值E为(节点号从小到大排列):
- +
- +
各节点的电压大小U为(节点号从小到大排列):
各节点的电压相角为θ(节点号从小到大排列):
0 - -
各节点的功率S为(节点号从小到大排列):
+ - - - -
- - +
各条支路的首端功率Si为(顺序同您输入B1时一样):
+ - +
+ +
+