电力系统随机潮流概率潮流计算MATLAB程序包含蒙特卡洛模拟法、半不变量法＋级数展开（Gram-Charlie，Cornish-Fisher）；
考虑光伏不确定性（Beta分布），以IEEE34节点为例，计算节点电压、支路潮流概率密度、累计概率并绘制曲线。
有注释，附带参考文献，不代做。
缺点是该节点系统不是很具有代表性，因为仅存在一个发电机节点，附带该节点系统的拓扑数据供参考。
这个程序主要是进行概率潮流计算，并考虑了分布式电源、发电机和负荷的随机波动。它应用在电力系统领域，用于解决概率潮流计算问题。程序的主要思路是先进行基础潮流计算，形成雅克比矩阵，然后计算输入的半不变量，包括发电机和负荷的八阶半不变量。接下来，根据光伏随机特性建模，计算太阳能光伏随机输出的八阶半不变量。然后，计算电压幅值及相角的八阶半不变量，并计算随机变量电压幅值的八阶半不变量。接着，计算支路潮流的灵敏度矩阵和八阶半不变量。最后，使用Monte Carlo模拟方法和Cornish-Fisher级数展开方法进行概率分析，并计算各节点电压越限概率。
程序中使用了一些函数和数据文件，例如dataIn函数用于读取数据文件，formACY函数用于形成交流系统节点导纳矩阵，NR_main函数用于进行潮流计算，NcalGCum函数用于计算发电机的八阶半不变量，NcalPLCum函数用于计算负荷的八阶半不变量，ProbCMGC函数用于计算电压越限概率（基于CMGC方法），ProbMC函数用于计算电压越限概率（基于Monte Carlo模拟方法），ProbCMCF函数用于计算电压越限概率（基于Cornish-Fisher级数展开方法）。
程序中的代码逻辑结构比较清晰，分为多个部分，每个部分都有相应的注释说明。程序运行过程中，首先读取输入数据，并进行基础潮流计算，形成雅克比矩阵。然后，根据输入的半不变量和随机特性数据，计算发电机和负荷的八阶半不变量，以及光伏随机输出的八阶半不变量。接着，计算电压幅值及相角的八阶半不变量，以及随机变量电压幅值的八阶半不变量。然后，计算支路潮流的灵敏度矩阵和八阶半不变量。最后，使用Monte Carlo模拟方法和Cornish-Fisher级数展开方法进行概率分析，并计算各节点电压越限概率。
整个程序的结构比较清晰，代码注释详细，适合零基础的程序员阅读和理解。程序主要涉及到电力系统概率分析、潮流计算、随机变量建模等知识点。
YID:2749653475056796
