潮流计算报告

一、系统结构图：

二、网络参数：

四、matlab程序：

clear;

Un=input('请输入Un:'); %输入所需的额定电压PQ=[

%节点电压有功无功

Un 0 0

Un 4 2

Un 6 3.2

Un 3 1.44

Un 4 3.2

Un 2 1.1

];

FT=[

%首端末端

4 3

3 2

6 5

5 2

2 1

];

RX=[

% R X

4 8

3 6

4 4

1 2

2 4

];

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

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

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

maxd=1

k=1

while maxd>0.0001

k=k+1;

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

losss(i,1)=RX(i,1)*x; %计算线路有功损耗 /MW

losss(i,2)=RX(i,2)*x; %计算线路无功损耗/MW

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; %计算沿支路的无功/MW 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

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

fullloss(1,1)=0;%计算线路总损耗

fullloss(1,2)=0;

finalPQ=max(PQ1);

for i=1:NB

fullloss(1,1)=fullloss(1,1)+losss(i,1);

fullloss(1,2)=fullloss(1,2)+losss(i,2);

end

disp('辐射网迭代次数：')

k

disp('辐射网系统电压差精度：')

maxd

disp('辐射网系统末端节点有功和无功：')

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

disp('辐射网系统总功率损耗：')

fullloss %线路总损耗MVA

disp('辐射网系统各支路功率损耗：')

losss %各支路损耗MVA

disp('辐射网系统各节点电压幅值：')

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

disp('辐射网系统各节点电压相角：')

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

end

clc

disp('辐射网迭代次数：')

k

disp('辐射网系统电压差精度：')

maxd

disp('辐射网系统末端节点有功和无功/MVA：')

FinPQ=finalPQ(1,1)+finalPQ(1,2)*j %潮流分布即支路首端潮流MVA

disp('辐射网系统总功率损耗/MVA：')

Fulloss=fullloss (1,1)+fullloss(1,2)*j %线路总损耗MVA

disp('辐射网系统各支路功率损耗/MVA：')

for(a=1:5)

LOSS=losss (a,1)+losss(a,2)*j %各支路损耗MVA

end

disp('辐射网系统各节点电压幅值/KV：')

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

disp('辐射网系统各节点电压相角：')

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

n=5; %input('节点数');

nl=6; %input('支路数');

isb=1; %input('平衡母线节点号');

pr=0.000001; %input('误差精度：pr=');

B1=[1,2,13.6+125.5i,0.00006785i,1,0;

1,3,8.321+130.5i,0.00005224i,1,0;

3,5,10.2+128.8i,0.00007499i,1,0;

2,3,8.5+105.4i,0.00002836i,1,0;

1,4,7.579+129.6i,0.00005145i,1,0;

4,5,13.84+125.31i,0.0000278i,1,0]; %input('由支路参数形成的矩阵');

B2=[-FinPQ,0,Un,0,0,1;

100,0,Un,Un,0,3;

0,15+9.4i,Un,0,0,2;

0,27+6i,Un,0,0,2;

0,35.5+25.5i,Un,0,0,2]; %input('各节点参数形成的矩阵');

Y=zeros(n);e=zeros(1,n);f=zeros(1,n);V=zeros(1,n);sida=zeros(1,n);S1= zeros(nl);%对各矩阵置零

%－－－－－－－修改部分－－－－－－－－－－－－

ym=1;

SB=100;UB=Un; %定义视在功率和电压基值

if ym~=0 %若不是标幺值

YB=SB./UB./UB; %定义导纳标幺值

BB1=B1;

BB2=B2;

for i=1:nl

B1(i,3)=B1(i,3)*YB; %切换为阻抗标幺值

B1(i,4)=B1(i,4)./YB; %切换为导纳标幺值

end

disp('支路矩阵B1=');

sparseB1=sparse(B1);

disp(sparseB1) %输出标幺值稀疏矩阵B1

disp('-----------------------------------------------------');

for i=1:n

B2(i,1)=B2(i,1)./SB; %切换为视在功率标幺值

B2(i,2)=B2(i,2)./SB; %切换为视在功率标幺值

B2(i,3)=B2(i,3)./UB; %切换为电压标幺值

B2(i,4)=B2(i,4)./UB; %切换为电压标幺值

B2(i,5)=B2(i,5)./SB; %切换为视在功率标幺值

end