配电网辐射网前推回代潮流计算程序matlab

配电网辐射网前推回代潮流计算程序（matlab）配电网辐射网前推回代潮流计算程序（matlab）

clear
n=40;%n=input（‘请输入节点数‘）;
b=39;%b=input（‘请输入支路数‘）;
%disp（‘请输入支路阻抗和节点功率矩阵，第一列存支路号，第二列存首节点号，第三列存尾节点号，第四列存支路自阻抗第五列存尾节点给定功率‘）;
Sb=15;
Ub=23;
Zb=Ub^2/Sb;
Z=[1010.169/Zb+i*0.655/Zb0
2120.864/Zb+i*0.751/Zb0
3130.196/Zb+i*0.655/Zb0
4241.374/Zb+i*0.774/Zb0.522/Sb+i*0.174/Sb
5250.864/Zb+i*0.751/Zb0
6360.444/Zb+i*0.439/Zb0
7370.196/Zb+i*0.655/Zb0
8580.864/Zb+i*0.751/Zb0
9590.864/Zb+i*0.751/Zb0
106101.374/Zb+i*0.774/Zb1.08/Sb+i*0.36/Sb
116110.864/Zb+i*0.751/Zb0
127120.279/Zb+i*0.051/Zb0
137130.279/Zb+i*0.051/Zb0
148141.374/Zb+i*0.774/Zb0.36/Sb+i*0.12/Sb
158151.374/Zb+i*0.774/Zb0.63/Sb+i*0.21/Sb
169161.374/Zb+i*0.774/Zb0.45/Sb+i*0.15/Sb
179170.864/Zb+i*0.751/Zb0
1811181.374/Zb+i*0.774/Zb0.54/Sb+i*0.18/Sb
1911191.374/Zb+i*0.774/Zb0.675/Sb+i*0.225/Sb
2012201.374/Zb+i*0.774/Zb0.45/Sb+i*0.15/Sb
2112210.444/Zb+i*0.439/Zb0
2213220.444/Zb+i*0.439/Zb1.8/Sb+i*0.6/Sb
2313230.444/Zb+i*0.439/Zb0
2417241.374/Zb+i*0.774/Zb0.675/Sb+i*0.225/Sb
2517251.374/Zb+i*0.774/Zb0.63/Sb+i*0.21/Sb
2621261.374/Zb+i*0.774/Zb0.54/Sb+i*0.18/Sb
2721270.444/Zb+i*0.439/Zb0
2823281.374/Zb+i*0.774/Zb0.855/Sb+i*0.285/Sb
2923290.864/Zb+i*0.751/Zb0
3027301.374/Zb+i*0.774/Zb0.63/Sb+i*0.21/Sb
3127310.864/Zb+i*0.751/Zb0
3229321.374/Zb+i*0.774/Zb0.54/Sb+i*0.18/Sb
3327331.374/Zb+i*0.774/Zb0.72/Sb+i*0.24/Sb
3431341.374/Zb+i*0.774/Zb0.765/Sb+i*0.225/Sb
3531350.864/Zb+i*0.751/Zb0
3635361.374/Zb+i*0.774/Zb0.54/Sb+i*0.18/Sb
3735370.864/Zb+i*0.751/Zb0
3837381.374/Zb+i*0.774/Zb0.675/Sb+i*0.225/Sb
3937391.374/Zb+i*0.774/Zb0.522/Sb+i*0.174/Sb];
k=0;
V=ones（n1）;
t=0;
%迭代开始处
while t%算节点注入电流
 x1=Z（b3）;x=x1-n;
   for l=1:b
       j=Z（l3）;
       ua=V（j+11）;
     I（j1）=conj（Z（j5）/ua）;
   end
%回推算支路电流
J=zeros（b1）;
l=b;
J（l）=J（l）+I（l）;
for jj=1:b-1
    l=l-1;
 for m=l+1:b
     if Z（m2）==Z（l3） 
     J（l）=J（l）+J（m）;
    end 
end
J（l）=J（l）+I（l）;
end
%前推算节点电压
   for l=1:b
    j=Z（l3）+1;
    i=Z（l2）+1;
    V（j1）=V（i1）-Z（l4）*J（l1）;
  end
%收敛判定
t=0;
for j=2:n
 SS=V（j1）*conj（I（j-11））;
dp=real（SS-Z（j-15））;
dq=imag（SS-Z（j-15））;
S（j-11）=SS;
ddp=abs（dp）;
ddq=abs（dq）;
L1=（ddp<0.000001）&（ddq<0.000001）;
F（j-11）=L1;
if L1==1
 t=t+1;
end
end
 k=k+1;
end

%输出结果与否

 disp（‘输出直角坐标各节点电压‘）;
 disp（V）;
 disp（‘显示迭代次数‘）;
 disp（k）;
disp（‘显示收敛节点情况“1“表示收敛，“0“表示不收敛‘）;
 disp（F）;
 for j=1:b
     if F（j1）==0
         disp（‘显示不收敛节点号、计算功率‘）;
disp（j）;disp（S（j1））;
     end
 end
for j=1:n
Vm（j1）=abs（V（j1））;Va（j1）=angle（V（j1））;
end
disp（‘输出各节点电压幅值‘）;
disp（Vm）;
disp（‘输出各节点电压相角‘）;
disp（Va）
 %返回开始处

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       2924  2008-06-19 09:15  cthd39-40.m

----------- ---------  ---------- -----  ----

                 2924                    1