3.6.2支路15断开
此时支路的数目为9，相应地， B1矩阵变为：
?1?1??4??2B1??5??3?4??5?6?34566789108.5+20.1i0.000556i13.6+32.16i0.0002224i10.2+24.12i0.0001668i8.5+20.1i0.0004448i6.8+16.08i0.000556i1.494+40.333i01.781+53.885i01.494+40.333i01.781+53.885i01111111110?0??0??0?0? ?1?1??1?1??将程序中的B1矩阵替换后，运行并进行调节，此系统经过一次调节后达到
要求，对应的调节方法和结果见下表：
调节方法：
表22 调节电压方法（支路15断开）
调节方法 电厂1电压 电厂2电压 分接头1 分接头2 分接头3 分接头4 未调整 第一次 1.05 1.05 1.05 1.05 1 1 1 0.95 1 1 1 0.95 调节结果： 表23 调节结果（支路15断开）
节点电压 未调整 第一次 1 2 3 4 5 6 7 8 9 10 1.05 1.05 1.03 1.00 1.00 1.01 1.01 0.93 0.98 0.96 1.05 1.05 1.03 1.00 1.00 1.01 1.01 0.98 0.98 1.02 由电压的标幺值可以观察到，变比调到极限后仍然不符合题目要求，但我们可以调节无功功率来达到题目要求。
此时各节点的功率S为(节点号从小到大排列)：
1.0e+002 *
0.3440 + 0.8892i 2.0000 + 0.0026i -0.0000 - 0.0000i -0.0000 - 0.0000i -0.0000 - 0.0000i -0.0000 - 0.0000i -0.5000 - 0.3099i -0.8000 - 0.4958i -0.3000 - 0.1859i -0.6000 - 0.3718i 得到的图像如下：
1.1520电压标幺值1.11.0510.95电压角度05节点号10100-1005节点号1031节点注入有功支路首端无功12345678910210-10.50-0.5123456789图8 支路15断开时MATLAB图
3.6.2支路45断开
此时支路的数目为9，相应地， B1矩阵变为：
?1?1??1??2B1??5??3?4??5?6?34566789108.5+20.1i0.000556i13.6+32.16i0.0002224i13.6+32.16i0.0002224i8.5+20.1i0.0004448i6.8+16.08i0.000556i1.494+40.333i01.781+53.885i01.494+40.333i01.781+53.885i01111111110?0??0??0?0? ?1?1??1?1??将程序中的B1矩阵替换后，运行并进行调节，此系统经过一次调节后达到
要求，对应的调节方法和结果见下表：
调节方法：
表24 调节电压方法（支路45断开）
调节方法 电厂1电压 电厂2电压 分接头1 分接头2 分接头3 分接头4 未调整 第一次 1.05 1.05 1.05 1.05 1 1 1 0.95 1 1 1 0.975 调节结果： 表25 调节结果（支路45断开）
节点电压 未调整 第一次 1 2 3 4 5 6 7 8 9 10 1.05 1.05 1.03 1.00 1.03 1.03 1.01 0.92 1.02 0.98 1.05 1.05 1.03 1.00 1.03 1.03 1.01 0.98 1.02 1.01 由电压的标幺值可以观察到，变比调到极限后仍然不符合题目要求，但我们可以调节无功功率来达到题目要求。
此时各节点的功率S为(节点号从小到大排列)：
1.0e+002 *
0.3691 + 1.2863i 2.0000 - 0.3831i -0.0000 - 0.0000i -0.0000 - 0.0000i -0.0000 + 0.0000i -0.0000 - 0.0000i -0.5000 - 0.3099i -0.8000 - 0.4958i -0.3000 - 0.1859i -0.6000 - 0.3718i 得到的图像如下：
1.1520电压标幺值1.11.0510.95电压角度05节点号10100-1005节点号1031节点注入有功支路首端无功12345678910210-10.50-0.5123456789
图9 支路45断开时MATLAB图
4.7DDRTS仿真
利用DDRTS软件，进行绘制系统图进行各种情况的潮流的分析。
4.7.1变电所负荷不变时的求解
当4个变电所为正常负荷情况时 ，进行潮流计算，得到潮流图如下：
图10 负荷正常时的仿真图
将DDRTS仿真得到的潮流图中个点电压的标幺值与MATLAB计算得到的各节点电压进行比较：
表26：MATLAB与DDRTS结果比较（正常负荷） 比较 1 2 3 4 5 6 7 8 9 10 MAT1.05 1.05 1.03 1.01 1.03 1.03 1.01 1.00 1.01 1.01 LAB DDR1.05 1.05 1.036 1.018 1.032 1.035 1.109 1.05 1.014 1.08 ST 误差 0 0 0.006 0.008 0.002 0.005 0.101 0.007 0.004 0.08 通过上表我们将MATLAB所得结果与DDRTS潮流模拟结果进行比较，发现两组计算结果大部分相差为1%左右，可知此时计算的计算结果较准确。
由上述比较可观察到: 应用DDRTS仿真和MATLAB计算的到的结果大致相同，可以认为此潮流计算是正确的。
4.7.2 4个变电所负荷同时以2%的比例增大
当4个变电所的负荷同时以2%的比例上升时 ，进行潮流计算。 得到潮流图如下：
图11 负荷增大2%时的仿真图
比较结果如下表：
表27：MATLAB与DDRTS结果比较（负荷增大2%）
比较 1 MAT1.05 LAB DDR1.05 TS 误差 0 2 1.05 3 1.03 4 1.01 5 1.02 6 1.03 7 1.00 1.10 0.10 8 1.00 1.05 0.05 9 1.01 1.017 0.007 10 1.01 1.09 0.09 1.05 1.036 1.018 1.03 1.034 0 0.006 0.008 0.01 0.004 通过上表我们将MATLAB所得结果与DDRTS潮流模拟结果进行比较，发现两组计算结果大部分相差为1%左右，可知此时计算的计算结果较准确。
4.7.3 4个变电所负荷同时以2%的比例下降
当4个变电所的负荷同时以2%的比例下降时 ，进行潮流计算。 得到潮流图如下：