不对称短路电流计算机算法实现.docVIP

电力系统分析实验三 不对称短路电流计算机算法实现 电力系统的简化等值电路如图，以此为模型应用计算机算法计算计算当节点3发生单相短路接地、两相短路的瞬时，1）节点1和2的电压；（2）线路1-2、1-3和2-3的电流； 计算的源程序如下： clear ZZ1(1,2)=j*0.1; ZZ1(1,3)=j*0.1; ZZ1(2,3)=j*0.1;%节点m,n之间的正序阻抗（mn） ZZ2(1,2)=j*0.1; ZZ2(1,3)=j*0.1; ZZ2(2,3)=j*0.1;%节点m,n之间的负序阻抗（mn） ZZ0(1,2)=j*0.2; ZZ0(1,3)=j*0.2; ZZ0(2,3)=j*0.2;%节点m,n之间的零序阻抗（mn） Y1=[-j*26.6266 j*10 j*10 j*10 -j*33.2933 j*10 j*10 j*10 -j*19.96];%输入正序网络节点导纳矩阵 Y2=[-j*26.6266 j*10 j*10 j*10 -j*33.2933 j*10 j*10 j*10 -j*19.96];%输入负序网络节点导纳矩阵 Y0=[-j*30 j*5 j*5 j*5 -j*50 j*5 j*5 j*5 -j*10];%输入零序网络节点导纳矩阵 YY1=[-j*39.96 j*10 j*10 j*20 0 j*10 -j*59.96 j*10 0 j*40 j*10 j*10 -j*19.96 0 0 j*20 0 0 -j*30 0 0 j*40 0 0 -j*60]; YY2=YY1;%输入包括发电机机端电压节点的正，负序网络节点导纳矩阵 N1=3;%输入网络的节点数 N2=5;%输入包括所有发电机节点的网络的节点数 k=3;%输入短路点的节点号 fault=1;%输入短路类型f(3)=3;f(1)=1;f(2)=2;f(1,1)=4 %第一部分：计算所有节点的a,b,c三相电压 for p=1:N1 if p==k I(p)=1; else I(p)=0; end end Z1(:,k)=Y1\I;Zk1=Z1(:,k);%正序网络中节点m的自阻抗和互阻抗 Z2(:,k)=Y2\I;Zk2=Z2(:,k);%负序网络中节点m的自阻抗和互阻抗 Z0(:,k)=Y0\I;Zk0=Z0(:,k);%零序网络中节点m的自阻抗和互阻抗 if fault==1%根据故障类型选择不同的计算公式 Ik1=1/(Z1(k,k)+Z2(k,k)+Z0(k,k)); Ik2=Ik1;Ik0=Ik1; else if fault==2 Ik1=1/(Z1(k,k)+Z2(k,k)); Ik2=-Ik1;Ik0=0; else if fault==3 Ik1=1/Z1(k,k);Ik2=0;Ik0=0; else if fault==4 Ik1=1/(Z1(k,k)+Z2(k,k)*Z0(k,k)/(Z2(k,k)+Z0(k,k))); Ik2=-Ik1*Z0(k,k)/(Z2(k,k)+Z0(k,k)); Ik0=-Ik1*Z2(k,k)/(Z2(k,k)+Z0(k,k)); end end end end Ik1 %计算短路节点的正序电流 for p=1:N1 if p==k I1(p)=-Ik1; I2(p)=-Ik2; I0(p)=-Ik0; else I1(p)=0; I2(p)=0; I0(p)=0; end end uu1(:,k)=Y1\I1.; uu2(:,k)=Y2\I2.; uu0(:,k)=Y0\I0.; for p=1:N1