基于前推回代潮流计算(MATLAB).docx

本文主要运用前推回代潮流法计算菲宝线配电网络的线损，该网络图如下所示：

图一菲宝线配网图及原始节点编号
所给的物理参数如下：

在本文中，主要的解答过程如下：
一、节点和支路编号
为了编程的方便，在此源点编号为1，其他节点随机编号，而支路的编号与该支路终点的编号取值一样。菲宝线配网的编号图如下所示：

图二节点及支路编号
二、给节点分层
在给节点分层之前，举例说明节点分层的方法。
1
3
2
4
5
6
7
8
9
1000
图三

在图三中，先找出网络中的一颗树干（树干是随机的），同时去掉目标树干上的末梢节点。先将树干上的节点分层（树干上有n个节点就分成年n层），再找树干上各节点对应的子节点，将这些子节点分在同一层，依此类推，最后将末梢点分在同一层。在图一中，可将该网络分成六层，分层情况如下：
第一层：1
第二层：2
第三层：3
第四层：4
第五层：6，8
第六层：5，7，9，10
在菲宝线潮流计算中，先建立原始数据矩阵Raw_data,第1列、第二列分别为支路始节点和末节点；第3-5列为该支路对应的线路参数，分别对应线路的长度、线路的单位电阻和单位阻抗（由于线路小于100km，线路的导纳约等于零）；第6列表示支路的末节点是否装有变压器，若无置为0，若有且配网抄表终端安装在变压器的低压侧置为1，若配网抄表终端安装在变压器的高压侧置为2；第7-12列为对应变压器参数，分别为该变压器的额定容量、额定电压、短路损耗、短路电压百分比、空载损耗和空载电流百分比；第13-14列分别为该变压器所接负额的有功功率和无功功率。
原始数据矩阵Raw_data，如图四所示

图四原始数据矩阵Raw_data
由原始数据矩阵Raw_data生成弧结构矩阵C和末梢节点矩阵O。按照上述的节点分层方法将该网络节点分成17层并存入节点分层矩阵layer中，如图五所示：

图五节点分层矩阵layer
三、计算过程
第一步：给节点电压赋初值，计算出变压器的总注入功率。
第二步：计算末梢点的注入电流，由外层向根节点推进，计算各节点的注入电流Ii（该电流就是支路电流），也就是回推求解各支路电流。
第三步：前推求解各节点电压Ui。
Ui(c+1)=Uj（c+1）-Iij(c)Zij
第四步：求出各节点电压的修正量△Ui,同时计算出max(△Ui).。
△Ui(c+1)=Ui（c+1）-Ui(c)其中c为循环次数
第五步：判定收敛性
max(△Ui)<
若满足，则计算终止。否则继续计算。
四、计算结果
1.变压器总损耗：
T_total=
6.7097+4.4818i（KVA）
2.支路总损耗：
L_total=
0.2461+0.0407i(KVA)
3.线路总损耗：
S_total=
6.9559+4.5224i(KVA)