基于IEEE33节点的配电网重构，采用最优流法开展了配电网重构工作，得到重构方案

基于IEEE33节点的配电网重构，采用最优流法开展了配电网重构工作，得到重构方案，应打开的开关数等，同时对比了重构前后的网损和电压结果
–以下内容来源于第三方解读，仅供参考
这段代码是一个用于电力系统潮流计算的程序。它使用了牛顿-拉夫逊法（Newton-Raphson method）来迭代计算节点电压和功率的平衡。下面我会逐步解释代码的功能和算法。

首先，代码定义了一些变量和数据，包括节点个数、支路参数、节点参数等。其中，节点类型分为平衡节点、PQ节点、PV节点和PI节点。支路参数包括支路电阻、电抗、变压器变比和电纳。节点参数包括注入有功、注入无功、电压幅值和电压相位。

代码的主要部分是一个while循环，用于进行潮流计算的迭代。循环中的代码包括以下几个步骤：

计算节点导纳矩阵Y：根据支路参数计算节点之间的导纳矩阵Y。导纳矩阵是一个复数矩阵，表示节点之间的导纳关系。

初始化节点功率参数OrgS：根据节点类型和参数，计算初始的有功功率和无功功率。

创建DetaS：根据节点类型和参数，计算初始的有功功率、无功功率和电压幅值的不平衡量。

创建Jacbi矩阵：根据节点参数和导纳矩阵，计算雅可比矩阵。雅可比矩阵是一个复数矩阵，表示节点电压和功率之间的关系。

求解修正方程：使用雅可比矩阵和DetaS，通过求解线性方程组，得到节点电压的修正量DetaU。

修正节点电压：根据修正量DetaU，更新节点电压。

判断是否达到收敛条件：通过比较修正量DetaU的绝对值与误差精度pr，判断是否达到收敛条件。如果未达到，则返回第3步继续迭代。

计算系统网损：根据节点电压和导纳矩阵，