⑶ 第三类平衡节点。潮流计算时一般只设一个平衡节点。等值负荷功率、是给定的，节点电压的大小和相位也是给定的。担负调整系统频率任务的发电厂母线往往被选作为平衡节点。 二.用牛顿法计算潮流时，有以下的步骤： ⑴输入线路，电气元件参数，形成节点导纳矩阵。 ⑵给这各节点电压初始值。
⑶将以上电压初始值代入式（4—38a）~式（4—38c）或式（4—45c）、（4—45a），求出修正方程式中的不平衡量。
⑷将各节点电压的初值代入式（4-41a）、式（4-41b）或式（4-49a）~式（4-49d）,求修正方程式的系数矩阵——雅克比矩阵的各个元素
?0??0??0??0??0?。 Hij、Nij、Jij、L?ij0?以及Rij、Sij⑸解修正方程式，求各节点电压的变化量，即修正量。 ⑹计算各节点电压的新值，即修正后值
ei?1??ei?0???ei?0?;fi?1??fi?0???fi?0?
⑺运用各节点电压的新值自第三步开始进入下一次迭代。 ⑻计算平衡节点功率和线路功率。 其中，平衡节点功率为
线路功率为
Sij?UiIij?Ui[Uiyi0?(Ui?Uj)]?Pij?jQijSji?UjIji?Uj[Ujyj0?(Uj?Ui)?Pji?jQji从而，线路上损耗的功率为
~???????~???????
?Sij?Sij?Sji??Pij?j?Qij
三.MATLAB程序画牛顿拉夫逊法的流程图
~~~问题分析
一.节点设置及分类
根据系统图可知此系统为两端供电网路，将母线1，2设为节点1，6，将变电所1、2、3、4的高低压侧分别设为节点2，3，4，5，7，8，9，10。并且，将节点1设为平衡节点，将节点6设为PV节点，其余节点设为PQ节点。 二.参数求取
设定基准值，，所以根据题目原始资料，计算发电厂、变压器及线路的参数。
（1）运用下列公式计算变压器参数（折算到低压测）：
计算变压器分接头变比（折算到低压测）：
,
（2）计算线路参数
单位阻抗值为，长度为，则线路总阻抗为：
单位电纳，总导纳为：
（3）变电所负荷分别为：
变电所1 =60+j37.18 变电所2 =40+j24.79 变电所3 =40+j24.79 变电所4 =50+j30.98
表1：发电厂、变压器及线路的参数
首端号 1 1 1 6 7 2 4 7 9 末端号 2 4 7 9 9 3 5 8 10 阻抗有名值 8.5+j20.1 4.25+j10.05 7.65+j18.09 8.5+j20.1 6.8+j16.08 0.05455+j1.65024 0.003735+j0.10083 0.05455+j1.65024 0.003735+j0.10083 电纳有名值 j0.000556 j0.000278 J0.000556 J0.000448 0 0 0 0
三.等值电路图
四.MATLAB潮流计算 1. B1、B2矩阵的形成：
B1矩阵：1、支路首端号；2、末端号；3、支路阻抗4、线路对地电纳 (或变压器导纳)；5、支路的变比；6、支路首端处于K侧为1，1侧为0；7、线路变压器标识（01）。
根据所求参数，以及B1矩阵的含义，列写B1矩阵如下:
B1=[1 2 8.5+20.1i 5.56e-4i 1 0 0;
1 4 4.25+10.05i 2.78e-4i 1 0 0; 1 7 7.65+18.09i 5.004e-4i 1 0 0; 6 9 8.5+20.1i 5.56e-4i 1 0 0;
7 9 6.8+16.08i 4.448e-4i 1 0 0; 2 3 5.455e-2+1.65024i 0 22038.5 1 1; 4 5 3.735e-3+0.10083i 0 22011 1 1; 7 8 5.455e-2+1.65024i 0 22038.5 1 1; 9 10 3.735e-3+0.10083i 0 22011 1 1;] B2矩阵：1、该节点发电机功率；2、该节点负荷功率； 3、PQ节点电压初始值，或平衡节点及PV节点电压的给定值。4、节点所接无功补偿并联电容（感）的电纳。5、节点分类标号：1为平衡节点（应为1号节点）；2为PQ节点；3为PV节点。
由题目只两个电厂共同供给功率，结合发电机效率，功率损耗等条件选择合理的分配方式，即：发电厂二供给功率为60MW，其余由发电厂一供给。
B2=[0 0 231 0 1;
0 0 220 0 2; 0 60+37.18i 35 0 2; 0 0 220 0 2; 0 40+24.79i 10 0 2; 60 0 231 0 3 0 0 220 0 2; 0 40+24.79i 35 0 2; 0 0 220 0 2; 0 50+30.98i 0 0 2;] 2.编写程序并运行 (1) 编写程序见附录1