程序说明书(4)

z1=inv(y1);

%*******************************************************************计算短路点电流 Z1=diag(z1);%正（负）序网节点入端阻抗 Z0=diag(z0);%零序网节点入端阻抗

%********************* %给If1，If2，If0赋值 %*********************

if malfunction_type==3%发生三相短路时短路点的正序短路电流 If1_three=U./Z1; If2_three=zeros(n,1); If0_three=zeros(n,1);

If1=If1_three; If2=zeros(n,1); If0=zeros(n,1); end

if malfunction_type==1%发生单相短路接地时短路点的正负零序短路电流 If1_one0=U./(Z1+Z1+Z0); If2_one0=If1_one0; If0_one0=If1_one0;

If1=If1_one0; If2=If2_one0; If0=If0_one0; end

if malfunction_type==2%发生两相短路时短路点的正负序短路电流 If1_two=U./(Z1+Z1); If2_two=-If1_two; If0_two=zeros(n,1);

If1=If1_two; If2=If2_two; If0=zeros(n,1);

end

if malfunction_type==11%发生两相短路接地时短路点的正负零序短路电流 If1_two0=U./(Z1+(Z1.*Z0./(Z1+Z0))); If2_two0=-If1_two0.*(Z0./(Z1+Z0)); If0_two0=-If1_two0.*(Z1./(Z1+Z0));

If1=If1_two0; If2=If2_two0; If0=If0_two0; end

If1_malfunction=If1(malfunction_nodenumber) If2_malfunction=If2(malfunction_nodenumber) If0_malfunction=If0(malfunction_nodenumber)

If1_abs=abs(If1);

If1_angle=angle(If1)*180/3.1416; If2_abs=abs(If2);

If2_angle=angle(If2)*180/3.1416; If0_abs=abs(If0);

If0_angle=angle(If0)*180/3.1416;

If1_abs_malfunction=If1_abs(malfunction_nodenumber) If1_angle_malfunction=If1_angle(malfunction_nodenumber) If2_abs_malfunction=If2_abs(malfunction_nodenumber) If2_angle_malfunction=If2_angle(malfunction_nodenumber) If0_abs_malfunction=If0_abs(malfunction_nodenumber) If0_angle_malfunction=If0_angle(malfunction_nodenumber) %a相

Ifa=If1+If2+If0; Ifa_abs=abs(Ifa)

Ifa_angle=angle(Ifa)*180/3.1416 %b相

Ifb=a*a*If1+a*If2+If0;

Ifb_abs=abs(Ifb)

Ifb_angle=angle(Ifb)*180/3.1416 %c相

Ifc=a*If1+a*a*If2+If0; Ifc_abs=abs(Ifc)

Ifc_angle=angle(Ifc)*180/3.1416

%*******************************************************************计算发生短路时各节点正负零序电压

I=zeros(n,1);

I(malfunction_nodenumber)=-If1(malfunction_nodenumber); U1=z1*I; U1=U1+U

I=zeros(n,1);

I(malfunction_nodenumber)=-If2(malfunction_nodenumber); U2=z1*I

I=zeros(n,1);

I(malfunction_nodenumber)=-If0(malfunction_nodenumber); U0=z0*I

U1_abs=abs(U1)

U1_angle=angle(U1)*180/3.1416 U2_abs=abs(U2)

U2_angle=angle(U2)*180/3.1416 U0_abs=abs(U0)

U0_angle=angle(U0)*180/3.1416 %a相 Ua=U1+U2+U0; Ua_abs=abs(Ua)

Ua_angle=angle(Ua)*180/3.1416

%b相

Ub=a*a*U1+a*U2+U0; Ub_abs=abs(Ub)

Ub_angle=angle(Ub)*180/3.1416 %c相

Uc=a*U1+a*a*U2+U0; Uc_abs=abs(Uc)

Uc_angle=angle(Uc)*180/3.1416

%*******************************************************************计算支路电流 %求各序支路电流

%******************************************************************* branch=size(x1);

I_branch1=zeros(branch,1); I_branch2=zeros(branch,1); I_branch0=zeros(branch,1);

for k=1:branch

if Fr_i(k)~=0&To_j(k)~=0

I_branch1(k)=(U1(Fr_i(k))-U1(To_j(k)))*b1(k); I_branch2(k)=(U2(Fr_i(k))-U2(To_j(k)))*b1(k); if Fr_i(k)<=3&To_j(k)<=3

I_branch0(k)=(U0(Fr_i(k))-U0(To_j(k)))*b0(k); end end end

I_branchmatrix1=[Fr_i,To_j,I_branch1] I_branchmatrix2=[Fr_i,To_j,I_branch2] I_branchmatrix0=[Fr_i,To_j,I_branch0]

I_branch_1_abs=abs(I_branch1)

I_branch_1_angle=angle(I_branch1)*180/3.1416 I_branch_2_abs=abs(I_branch2)

I_branch_2_angle=angle(I_branch2)*180/3.1416 I_branch_0_abs=abs(I_branch0)

I_branch_0_angle=angle(I_branch0)*180/3.1416

%a相

I_branch_a=[Fr_i,To_j,I_branch1+I_branch2+I_branch0] I_branch_a_abs=abs(I_branch1+I_branch2+I_branch0)

I_branch_a_angle=angle(I_branch1+I_branch2+I_branch0)*180/3.1416 %b相

I_branch_b=[Fr_i,To_j,a*a*I_branch1+a*I_branch2+I_branch0] I_branch_b_abs=abs(a*a*I_branch1+a*I_branch2+I_branch0)

I_branch_b_angle=angle(a*a*I_branch1+a*I_branch2+I_branch0)*180/3.1416 %c相

I_branch_c=[Fr_i,To_j,a*I_branch1+a*a*I_branch2+I_branch0] I_branch_c_abs=abs(a*I_branch1+a*a*I_branch2+I_branch0)

I_branch_c_angle=angle(a*I_branch1+a*a*I_branch2+I_branch0)*180/3.1416

百度搜索“77cn”或“免费范文网”即可找到本站免费阅读全部范文。收藏本站方便下次阅读，免费范文网，提供经典小说教育文库程序说明书(4)在线全文阅读。