资源描述：
《潮流程序的开发》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、课题组前期培训阶段二—潮流程序的开发1、数据结构的设计1.1输入数据输入的数据为线路，双绕组变压器，三绕组变压器的铭牌值及母线的参数。具体如下：线路参数存储到line矩阵中，格式如下：首端节点序号末端节点序号单位长度电阻单位长度电抗单位长度电纳线路长度双绕组变压器参数存储到trans矩阵中，格式如下：高压侧节点序号低压侧节点序号容量高压侧额定电压低压侧额定电压短路损耗短路电压百分数空载损耗空载电流百分数三绕组变压器参数存储到tw_trans矩阵中，格式如下：高压侧节点序号中压侧节点序号低压侧节点序号高压侧容量中压侧容量低压侧容量高压侧额定电压中压侧额定电压低压侧额定电压短路损耗P12短路损耗
2、P23短路损耗P31短路电压百分数V12短路电压百分数V23短路电压百分数V31空载电流百分数I0空载损耗P0母线参数存储到bus矩阵中，格式如下：节点序号基准电压节点类型注入有功注入无功电压幅值电压相角节点类型中：1代表PQ节点2代表PV节点3代表平衡节点1.2输出数据输出数据存储到poweFlow矩阵中格式如下节点序号注入有功注入无功电压幅值电压相角2、程序流程框图启动生成节点导纳矩阵Y设置节点初始电压ei(0),fi(0)或Ui(0),δi(0)置迭代次数k=1输入原始数据bus,gen，transform,Line计算各个节点的不平衡量:∆Pi(k),∆Qi(k),∆Ui2(k)计算
3、雅可比矩阵J(k)计算平衡节点功率、支路功率和支路损耗计算ei(k+1),fi(k+1)或Ui(k+1),δi(k+1)max{
4、∆ei(k),∆fi(k)
5、}<ε或max{
6、∆Ui(k),∆δi(k)
7、}<εk=k+1解修正方程求取∆ei(k),∆fi(k)或∆Ui(k),∆δi(k)是否停止潮流流程图1、程序功能模块介绍3.1节点导纳矩阵的形成函数名称：daonaFormfunctiondaona=daonaForm(Sb,line,bus,trans,tw_trans)输入参数：Sb(基准功率),line(线路参数),bus(母线参数),trans(双绕组变压器参数)，tw_trans
8、(三绕组变压器参数)输出参数:daona(节点导纳矩阵)程序步骤：1、根据线路参数通过公式计算出线路的的阻抗（取倒数用导纳表示），对地导纳。把上述数据按如下格式存储到line_admi矩阵中。如果为对地支路，首末端节点序号取同一个。支路首端节点序号支路末端节点序号支路导纳数值2、根据双绕组铭牌参数计算双绕组变压器的阻抗支路和励磁支路，均用导纳表示。把以上数据按如下格式存储到trans_admi矩阵中。如果为对地支路，首末端节点序号取相同值。支路首端节点序号支路末端节点序号支路导纳数值3、根据公式计算三绕组变压器的三侧的阻抗及励磁支路参数，注意三绕组变压器会增加一个节点。把以上数据按如下格式存
9、储到twtrans_admi矩阵中。如果为对地支路，首末端节点序号取同一个。支路首端节点序号支路末端节点序号支路导纳数值4、根据line_admi、trans_admi、twtrans_admi合成节点导纳矩阵。自导纳矩阵按如下方法形成：如Y11，查询上述矩阵的节点序号把含节点1的所有导纳加起来即可。互导纳按如下方法形成：支路的首末端序号不同的即为互导纳，首末端序号就是其在节点导纳矩阵的行、列坐标，注意取负号。3.2重新对母线矩阵进行排列处理函数名：busRankfunction[busRanking,busGain,numPQ,numPV]=busRank(bus,tw_trans)该函数
10、的功能是对bus矩阵进行处理。包括一下两方面：1、三绕组变压器的模型会增加PQ节点,这些节点要在不平衡功率方程及雅克比矩阵和求解修正量中反应出来。2、对节点矩阵进行重新排列，按照PQ节点、PV节点、平衡节点的顺序来排序，以方便后面形成雅克比矩阵。输入参数：bus原始的母线参数，tw_trans三绕组变压器参数矩阵输出参数：busRanking经过按PQ节点、PV节点、平衡节点的顺序排列过的母线参数；busGain仍然按照节点顺序排列的（但考虑三绕组增加节点后）母线参数；numPQ、numPV分别为PQ、PV节点的个数程序的步骤1、busRanking的形成Bus矩阵的第三列放置的是其节点类型
11、，通过查询第三列即可确定节点类型。PQ节点从第一行开始放，PV节点从倒数第二行开始放，平衡节点放最后一行。存放顺序依次为：节点序号，基准电压，节点类型，有功，无功，电压，相角。2、busGain的形成在原来bus的基础上增加N个节点。N为三绕组变压器的个数。增加的节点均为PQ节点，且P、Q都为0。3.3不平衡功率的形成函数名：unbalancedPower该函数的功能是形成不平衡功率。输入参数：daona,节
