电力系统不对称故障计算的Matlab算法程序文档格式.docx
- 文档编号:16699285
- 上传时间:2022-11-25
- 格式:DOCX
- 页数:10
- 大小:15.57KB
电力系统不对称故障计算的Matlab算法程序文档格式.docx
《电力系统不对称故障计算的Matlab算法程序文档格式.docx》由会员分享,可在线阅读,更多相关《电力系统不对称故障计算的Matlab算法程序文档格式.docx(10页珍藏版)》请在冰豆网上搜索。
05010;
06010;
2700;
08010;
3900;
45011;
46011;
57011;
69011;
78011;
89011];
B2=[01010;
B0=[01010;
45*3011;
46*3011;
57*3011;
69*3011;
78*3011;
89*3011];
A1=[111;
1-1./2-i*sqrt(3)/2-1./2+i*sqrt(3)/2;
1-1./2+i*sqrt(3)/2-1./2-i*sqrt(3)/2];
Vcs=[;
;
];
Z1=zeros(n1);
Z2=zeros(n2);
Z0=zeros(n0);
Y1=zeros(n1);
Y2=zeros(n2);
Y0=zeros(n0);
V1=zeros(n1,1);
V2=zeros(n2,1);
V0=zeros(n0,1);
I1=zeros(nl1,1);
I2=zeros(nl2,1);
I0=zeros(nl0,1);
form1=1:
3
m=0;
ifm1==1
nl=nl1;
B=B1;
Z=Z1;
n=n1;
elseifm1==2
nl=nl2;
B=B2;
Z=Z2;
n=n2;
elseifm1==3
nl=nl0;
B=B0;
Z=Z0;
n=n0;
end
Z=zeros(n);
V=zeros(n);
I=zeros(nl);
fork1=1:
nl
p=B(k1,1);
q=B(k1,2);
ifB(k1,6)==0
k=1./B(k1,5);
elsek=B(k1,5);
ifp==0
ifq>
m
Z(q,q)=B(k1,3);
m=m+1;
else
fori=1:
Z(i,m+1)=-Z(i,q);
Z(m+1,i)=-Z(q,i);
Z(m+1,m+1)=Z(q,q)+B(k1,3);
forj=1:
Z(i,j)=Z(i,j)-Z(i,m+1)*Z(m+1,j)./Z(m+1,m+1);
Z(i,m+1)=0;
m+1
Z(m+1,i)=0;
end
elseifq>
Z(i,q)=Z(i,p)*k;
Z(q,i)=Z(p,i)*k;
Z(q,q)=k^2*Z(p,p)+k^2*B(k1,3);
m=m+1;
else
Z(i,m+1)=k*Z(i,p)-Z(i,q);
Z(m+1,i)=k*Z(p,i)-Z(q,i);
Z(m+1,m+1)=k^2*Z(p,p)+Z(q,q)-2*k*Z(p,q)+k^2*B(k1,3);
Z(i,j)=Z(i,j)-Z(i,m+1)*Z(m+1,j)./Z(m+1,m+1);
ifm1==1
Z1=Z;
elseifm1==2
Z2=Z;
Z0=Z;
ifIf==0
Z1(f,f)=Z1(f,f);
Z2(f,f)=Z2(f,f);
Z0(f,f)=Z0(f,f);
Z1(f,f)=Z1(f,f)+Z1(1,1)-2*Z1(f,1);
Z2(f,f)=Z2(f,f)+Z2(1,1)-2*Z2(f,1);
Z0(f,f)=Z0(f,f)+Z0(1,1)-2*Z0(f,1);
ifLf==1
Z4=Z2(f,f)+Z0(f,f)+3*zf;
k2=1;
k0=1;
I1(f)=Vcs(f)./(Z1(f,f)+Z4);
I2(f)=k2*I1(f);
I0(f)=k0*I1(f);
a1=I1(f);
a2=I2(f);
a0=I0(f);
elseifLf==2
Z4=Z2(f,f)+2*zf;
k2=-1;
k0=0;
elseifLf==3
Z4=zf+(Z2(f,f)+zf)*(Z0(f,f)+zf+3*zg)./(Z2(f,f)+Z0(f,f)+2*zf+3*zg);
k2=-(Z0(f,f)+zf+3*zg)./(Z2(f,f)+Z0(f,f)+2*zf+3*zg);
k0=-(Z2(f,f)+zf)./(Z2(f,f)+Z0(f,f)+2*zf+3*zg);
a1=I1(f);
elseifLf==4
Z4=zf+(Z2(f,f)+zf)*(Z0(f,f)+zf)./(Z2(f,f)+Z0(f,f)+2*zf);
k2=-(Z0(f,f)+zf)./(Z2(f,f)+Z0(f,f)+2*zf);
k0=-(Z2(f,f)+zf)./(Z2(f,f)+Z0(f,f)+2*zf);
elseifLf==5
k0=1;
disp('
故障点处的正序电流标幺值为:
'
disp(a1);
故障点处的负序电流标幺值为:
disp(a2);
故障点处的零序电流标幺值为:
disp(a0);
n1
V1(i)=Vcs(i)-I1(f)*Z1(i,f);
n2
V2(i)=-I2(f)*Z2(i,f);
n0
V0(i)=-I0(f)*Z0(i,f);
end
各节点的正序电压标幺值为(节点号从小到大排):
disp(V1);
各节点的负序电压标幺值为(节点号从小到大排):
disp(V2);
各节点的零序电压标幺值为(节点号从小到大排):
disp(V0);
ifi==1
B=B1;
V=V1;
nl=nl1;
elseifi==2
B=B2;
V=V2;
nl=nl2;
elseifi==3
B=B0;
V=V0;
nl=nl0;
ifB(j,6)==0
k=B(j,5);
elsek=1./B(j,5);
p=B(j,1);
q=B(j,2);
e=0;
b=B(j,3);
I(j)=(e-V(q)./k)./b;
I(j)=(V(p)-V(q)./k)./B(j,3);
I1=I;
I2=I;
I0=I;
各支路的正序电流标幺值为(顺序同您输入B时一样);
disp(I1(i));
各支路的负序电流标幺值为(顺序同您输入B时一样);
disp(I2(i));
各支路的零序电流标幺值为(顺序同您输入B时一样);
disp(I0(i));
s0=a0;
s1=a1;
s2=a2;
s0=V0;
s1=V1;
s2=V2;
else
s0=I0;
s1=I1;
s2=I2;
CA=s0*A1(j,1)+s1*A1(j,2)+s2*A1(j,3);
ifj==1
A=CA;
elseifj==2
B=CA;
elseifj==3
C=CA;
Iadl=A;
Ibdl=B;
Icdl=C;
Vajy=A;
Vbjy=B;
Vcjy=C;
Iazl=A;
Ibzl=B;
Iczl=C;
故障点处A相电流Iadl为:
disp(Iadl);
故障点处B相电流Ibdl为:
disp(Ibdl);
故障点处C相电流Icdl为:
disp(Icdl);
各节点处A相电压Vajy分别为(节点号从小到大排):
disp(Vajy);
各节点处B相电压Vbjy分别为(节点号从小到大排):
disp(Vbjy);
各节点处C相电压Vcjy分别为(节点号从小到大排):
disp(Vcjy);
各支路A相电流Iazl分别为(顺序同您输入B时一样'
disp(Iazl(i));
各支路B相电流Ibzl分别为(顺序同您输入B时一样'
disp(Ibzl(i));
各支路C相电流Iczl分别为(顺序同您输入B时一样'
disp(Iczl(i));
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 电力系统 不对称 故障 计算 Matlab 算法 程序