loadflowcalculations.m

配电网潮流计算 配电网潮流计算 配电网潮流计算 配电网潮流计算

clear
format long;
n=10;b=9;
Z=[1,0,1,0.0233+i*0.078,0.0184+i*0.00460
   2,1,2,0.0026+i*0.1144,0.0098+i*0.0034
   3,2,3,0.1411+i*0.2278,0.0179+i*0.00446
   4,3,4,0.132+i*0.115,0.01598+i*0.0184
   5,4,5,0.3749+i*0.3266,0.0161+i*0.006
   6,5,6,0.1711+i*0.1491,0.0078+i*0.0011
   7,6,7,0.3885+i*0.22,0.0115+i*0.0006
   8,7,8,0.9065+i*0.5134,0.0098+i*0.0013
   9,8,9,1.0101+i*0.5721,0.0164+i*0.002];
%支路损耗设为零，求解各节点（1-9）功率
%gonglv=zeros(b,1);
for t=1:b
    sum=0;
    for x=t:b
        sum=sum+Z(x,5);
    end
    gonglv(t,1)=sum;
end
%取出电阻和节点功率的实部虚部
%dianya=zeros(b+1,1);
%dianzuShiBu=zeros(b,1);dianzuXuBu=zeros(b,1);
%gonglvShiBu=zeros(b,1);gonglvXuBu=zeros(b,1);
for t=1:b
    dianzuShiBu(t,1)=real(Z(t,4));dianzuXuBu(t,1)=imag(Z(t,4));
    gonglvShiBu(t,1)=real(gonglv(t,1));gonglvXuBu(t,1)=imag(gonglv(t,1));
end
%初始化节点电压值
dianya1(1,1)=1;
e=0.0001;k=1;%k2=1;
while (k>e)
for t=1:b
    dianya1(t+1,1)=sqrt(1/4*(dianya1(t,1))^2-gonglvShiBu(t,1)*dianzuShiBu(t,1)-gonglvXuBu(t,1)*dianzuXuBu(t,1))+1/2*dianya1(t,1);
end
%计算第i条线段的线路损耗
LP=zeros(b,1);
LQ=zeros(b,1);
for t=1:b
    LP(t,1)=dianzuShiBu(t,1)*(gonglvShiBu(t,1)^2+gonglvXuBu(t,1)^2)/dianya1(t+1,1);
    LQ(t,1)=dianzuXuBu(t,1)*(gonglvShiBu(t,1)^2+gonglvXuBu(t,1)^2)/dianya1(t+1,1);
end
%计算各节点功率
for t=1:b
    gonglvShiBu(t,1)=0;
    gonglvXuBu(t,1)=0;
    for x=t:b
        gonglvShiBu(t,1)=gonglvShiBu(t,1)+real(Z(x,5))+LP(x,1);
        gonglvXuBu(t,1)=gonglvXuBu(t,1)+imag(Z(x,5))+LQ(x,1);
    end
end
%第二次计算的节点电压值
dianya2(1,1)=1;
for t=1:b
    dianya2(t+1,1)=sqrt(1/4*(dianya2(t,1))^2-gonglvShiBu(t,1)*dianzuShiBu(t,1)-gonglvXuBu(t,1)*dianzuXuBu(t,1))+1/2*dianya2(t,1);
end
k=max(abs(dianya1-dianya2));
%k2=min(dianya1-dianya2);
dianya11=dianya1;  %dianya11记录dianya1的值
dianya1=dianya2;
end
gonglv=gonglvShiBu+i*gonglvXuBu;
for t=1:b
    dianya22(t,1)=dianya2(t+1,1);
end
%F1='各节点(1-9)注入功率';disp(F1);disp(gonglv);
F2='各节点(1-9)基波电压';disp(F2);disp(dianya22);