前推回代法计算流程

前推回代法计算流程

要看懂前推回代法计算程序，报告叙述计算原理及计算流程。绘制计算流程框图。确定前推回代支路次序（广度优先，或深度优先），编写前推回代计算输入文件。进行潮流计算。

下列为节点配电网结构图及系统支路参数和系统负荷参数表。

图1-2 节点配电网结构图

表1 系统支路参数

表2 系统负荷参数

主程序清单：[PQ,FT,RX]=case114(); %调用数据文件

NN=size(PQ,1); %节点数

NB=size(FT,1); %支路数数

V=PQ(:,1); %V初始电压相量

maxd=1

k=1

while maxd>0.0001

PQ2=PQ; %每一次迭代各节点的注入有功和无功相同PL=0.0;

for i=1:NB

kf=FT(i,1); %前推始节点号

kt=FT(i,2); %前推终节点号

x=(PQ2(kf,2)^2+PQ2(kf,3)^2)/V(kf)/V(kf); %计算沿线电流平方A

PQ1(i,1)=PQ2(kf,2)+RX(i,1)*x; %计算支路首端有功/MW RX(i,1)~R

PQ1(i,2)=PQ2(kf,3)+RX(i,2)*x; %计算沿支路的无功损耗/Mvar RX(i,2)~X

PQ2(kt,2)= PQ2(kt,2)+PQ1(i,1); %用PQ1去修正支路末端节点的有功P 单位MW PQ2(kt,3)= PQ2(kt,3)+PQ1(i,2); %用PQ1去修正支路末端节点的有功Q 单位Mvar PL=PL+RX(i,1)*x;

end

angle(1)=0.0;

for i=NB:-1:1

kf=FT(i,2); %回代始节点号

kt=FT(i,1); %回代终节点号

dv1=(PQ1(i,1)*RX(i,1)+PQ1(i,2)*RX(i,2))/V(kf); %计算支路电压损耗的纵分量dv1

dv2=(PQ1(i,1)*RX(i,2)-PQ1(i,2)*RX(i,1))/V(kf); %计算支路电压损耗的横分量dv2

V2(kt)=sqrt((V(kf)-dv1)^2+dv2^2); %计算支路末端电压/kV

angle(kt)=angle(kf)+atand(dv2/(V(kf)-dv1)); %计算支路

end

maxd=abs(V2(2)-V(2));

V2(1)=V(1);

for i=3:1:NN

if abs(V2(i)-V(i))>maxd;

maxd=abs(V2(i)-V(i));

end

end

maxd

k=k+1

PQ1 %潮流分布即支路首端潮流MVA

V=V2 %节点电压模计算结果kV

angle %节点电压角度计算结果单位度

PL %网损单位MW

end

clear

输入文件清单：function [PQ,FT,RX]=case114()

PQ=[

%节点电压有功无功

10.4 0 0

10.0 0.0342 0.0301

10.0 0.0693 0.0642

10.0 0.0845 0.0763

10.0 0.0295 0.0261

10.0 0.0474 0.0409

10.0 0.1176 0.0957

10.0 0.0946 0.0857

10.0 0.0916 0.0859

10.0 0.0271 0.0229

10.0 0.0696 0.0643

10.0 0.0676 0.0579

10.0 0.0298 0.0242

];

FT=[

%首端末端5 4

13 4

4 3

10 3

12 11 11 3

7 6

6 2

9 8

8 2

3 2

2 1 ];

RX=[

% R X 4.524 5.04 3.521 3.966 1.145 1.28 4.14 4.696 2.436 2.866

1.328 1.763

2.745 2.965 0.856 1.14

2.237 2.756

3.743

4.251

2.356 2.541

3.367 3.685 ];

计算过程

maxd =

1

k =

1

maxd =

0.1780

k =

2

PQ1 =

0.0296 0.0262

0.0299 0.0243

0.1443 0.1272

0.0272 0.0230

0.0678 0.0581

0.1378 0.1230

0.1182 0.0964

0.1660 0.1378

0.0920 0.0863

0.1890 0.1748

0.3847 0.3439

0.8099 0.7260

V =

Columns 1 through 8

10.4000 9.8807 9.8220 9.9672 9.9734 9.9701 9.9390 9.8550 Columns 9 through 13

9.9556 9.9780 9.9600 9.9668 9.9799

angle =

Columns 1 through 8

0 0.3011 0.3986 0.4211 0.4387 0.3421 0.3916 0.3878 Columns 9 through 13

0.4225 0.4173 0.4444 0.4747 0.4400

PL =

0.0471

maxd =

0.1787

k =

3