matlab连续潮流程序绘制PV曲线
静态电压稳定
该程序为连续潮流IEEE14节点和33节点的程序
运行出来有分岔点和鼻点
可移植性强,注释详细
这段程序主要是用来计算电力系统中的潮流分布,并绘制PV曲线。下面我将对程序进行详细的分析。
首先,程序开始时使用`clc`、`clear`和`close all`清除命令窗口、清除工作区变量和关闭所有图形窗口。
接下来,程序定义了一些基准值,包括电压基准值`Vbase`、功率基准值`Sbase`和阻抗基准值`Zbase`。
然后,程序通过`xlsread`函数从Excel文件中读取节点数据和支路数据,并将其存储在`BusData`和`BranchData`中。
接下来,程序对读取的数据进行标幺化处理,将功率和阻抗转换为标幺值。
然后,程序调用`Calculate_Ybus`函数计算节点导纳矩阵`Ybus`。
接着,程序记录下Slack、PV和PQ节点的索引,以及节点数量`N`。
然后,程序设置初始值,包括电压和相角。电压初始值为1,相角初始值为0。
接下来,程序计算节点的注入功率,即发电减去负荷。
然后,程序定义了一些参数,包括矫正环节的收敛精度和最大迭代次数。
接下来,程序对PV曲线的计算进行了分步处理。
第一步是以lambda为连续参数进行计算。程序使用一个while循环,不断进行预测和校正环节,直到潮流计算收敛或达到最大迭代次数。在预测环节,程序根据当前的相角、电压和lambda值计算雅可比矩阵,并根据雅可比矩阵计算修正量。然后,程序根据修正量更新相角、电压和lambda值。在校正环节,程序调用`NRLF`函数进行潮流计算,并判断是否收敛。如果潮流计算收敛,则将当前节点的电压值存储在`Y_Vph1`中,将lambda值存储在`X_LamdaPh1`中,并记录上一次的相角、电压和lambda值。如果潮流计算不收敛,则跳出第一步,进行第二步。
第二步是以V为连续参数进行计算。程序使用一个while循环,不断进行预测和校正环节,直到潮流计算收敛或达到最大迭代次数。在预测环节,程序根据当前的相角、电压和lambda值计算雅可比矩阵,并根据雅可比矩阵计算修正量。然后,程序根据修正量更新相角、电压和lambda值。在校正环节,程序计算修正值,并根据修正值更新相角、电压和lambda值。如果lambda值大于上一次的lambda值,则更新上一次的lambda值和节点的电压值。如果lambda值小于等于上一次lambda值的一定比例(由`changefactor`定义),则跳出循环,避免在拐点附近不收敛。
第三步是重新将连续参数变为lambda。程序使用一个while循环,不断进行预测和校正环节,直到潮流计算不收敛或lambda减小到0。在预测环节,程序根据当前的相角、电压和lambda值计算雅可比矩阵,并根据雅可比矩阵计算修正量。然后,程序根据修正量更新相角、电压和lambda值。在校正环节,程序调用`NRLF`函数进行潮流计算,并判断是否收敛。如果潮流计算收敛,则将当前节点的电压值存储在`Y_Vph3`中,将lambda值存储在`X_LamdaPh3`中。
最后,程序绘制PV曲线,将`X_LamdaPh1`、`X_LamdaPh2`和`X_LamdaPh3`以及对应的`Y_Vph1`、`Y_Vph2`和`Y_Vph3`作为横坐标和纵坐标进行绘制。同时,程序设置横坐标和纵坐标的标签和标题。
总结来说,这段程序主要是用来计算电力系统中的潮流分布,并绘制PV曲线。它涉及到了潮流计算、雅可比矩阵的计算、修正量的计算和更新、收敛判断等知识点。通过对节点数据和支路数据的读取和处理,以及对节点导纳矩阵的计算,可以得到电力系统中各节点的电压和相角。然后,通过预测和校正的过程,可以得到PV曲线。这个程序可以应用在电力系统的潮流计算和分析中,帮助工程师了解电力系统的潮流分布情况,优化电力系统的运行和规划。
YID:53120
详询客服 微信shujuqudong1 或shujuqudong6 或 qq68823886 或 27699885
图文详情请查看: http://matup.cn/676973227785.html