含风光发电的电力系统概率潮流计算，考虑负荷波动，风力和光伏出力不确定性，一个概率潮流计算的程序，用于考虑分布式电源、发电机和负荷随机波动的情况

含风光发电的电力系统概率潮流计算，考虑负荷波动，风力和光伏出力不确定性，算法方面：基于蒙特卡洛法和半不变量法（gram-charlier和corn-fisher级数）。
这是一个概率潮流计算的程序，用于考虑分布式电源、发电机和负荷随机波动的情况。它应用在电力系统领域，用于分析电力系统中节点电压和支路功率的概率分布情况。

使用了半不变量法来计算电力系统的概率潮流。它首先读取输入数据，包括节点和线路的信息，然后形成交流系统节点导纳矩阵。接下来进行潮流计算，得到节点电压、功率损耗等信息。然后，根据发电机和负荷的随机参数输入，计算发电机和负荷的八阶半不变量。接着，根据光伏随机特性建模，计算光伏的八阶半不变量。最后，计算电压幅值及相角的八阶半不变量、支路潮流的灵敏度矩阵以及支路有功和无功功率的八阶半不变量。
通过概率分析来评估电力系统中节点电压和支路功率的概率分布情况。输入是电力系统的节点和线路信息，输出是节点电压和支路功率的概率分布情况。代码中使用了半不变量法来计算概率分布，这种方法可以更准确地估计电力系统的概率分布情况。

代码用于分析电力系统中发电机和负荷的随机波动对潮流的影响。代码的输入是电力系统的基础参数，包括节点、线路、发电机和负荷的信息。代码的输出是电力系统中各节点的电压幅值和相角，以及支路的有功和无功功率。

YID:84150675311825751

主要工作包括以下几个步骤：

1. 读取输入数据：从文件中读取电力系统的基础参数，包括节点、线路、发电机和负荷的信息。

2. 形成节点导纳矩阵：根据输入的节点和线路信息，计算得到电力系统的节点导纳矩阵。

3. 进行潮流计算：使用牛顿-拉夫逊法进行潮流计算，得到电力系统中各节点的电压幅值和相角。

4. 计算发电机和负荷的随机参数：根据输入的发电机和负荷信息，计算得到发电机和负荷的八阶半不变量。

5. 计算电压和功率的八阶半不变量：根据电压幅值和相角的八阶半不变量，计算得到随机变量电压幅值和功率的八阶半不变量。

6. 计算电压和功率的概率分布：根据八阶半不变量的数据，使用Cornish-Fisher级数展开法，计算得到电压和功率的概率分布。

7. 绘制概率分布曲线：将计算得到的概率分布数据绘制成曲线，用于分析电压和功率的概率分布情况。

代码中使用的算法主要是牛顿-拉夫逊法进行潮流计算，以及Cornish-Fisher级数展开法计算概率分布。牛顿-拉夫逊法是一种迭代算法，用于求解非线性方程组，适用于潮流计算。Cornish-Fisher级数展开法是一种近似计算概率分布的方法，通过展开级数来逼近实际的概率分布。

这段代码的目的是分析电力系统中发电机和负荷的随机波