• 大小: 1KB
    文件类型: .rar
    金币: 1
    下载: 0 次
    发布日期: 2021-10-11
  • 语言: 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


评论

共有 条评论

相关资源