电力系统潮流及短路电流计算程序.doc

word班级：：学号：一、 作业要求编写程序计算图1所示算例系统的潮流与三相短路电流潮流计算：方法不限，计算系统的节点电压和相角短路电流：4号母线发生金属性三相短路时〔zf=0〕，分别按照准确算法和近似算法计算短路电流、系统中各节点电压以与网络中各支路的电流分布，并对两种情况下的计算结果进展比拟二、 电路图与参数图1 3机9节点系统表1 9节点系统支路参数支路R〔p.u.〕X〔p.u.〕B/2(TK) 〔p.u.〕1~402~703~904~54~65~76~97~88~9表2 9节点系统发电机参数发电机编号节点类型PG〔p.u.〕VG〔p.u.〕〔p.u.〕〔p.u.〕1平衡2PV3PV表3 9节点系统负荷参数节点编号节点类型Pi〔p.u.〕Qi〔p.u.〕4PQ005PQ6PQ7PQ008PQ19PQ00三、 计算步骤(1) 进展系统正常运行状态的潮流计算，求得(2) 形成不含发电机和负荷的节点导纳矩阵YN；(3) 将发电机表示为电流源()和导纳()的并联组合；节点负荷用恒阻抗的接地支路表示；形成包括所有发电机支路和负荷支路的节点导纳矩阵Y，即在YN中的发电机节点和负荷节点的自导纳上分别增加发电机导纳和负荷导纳〔〕；(4) 利用，计算节点阻抗矩阵，从而得到阻抗矩阵中的第f列；(5) 利用公式(6-7)或(6-10)计算短路电流；(6) 利用公式(6-8)或(6-11)计算系统中各节点电压；(7) 利用公式(6-9)计算变压器支路的电流；对输电线路利用П型等值电路计算支路电流。

四、 计算结果节点导纳矩阵Yn： Columns 1 through 5 0 -17.3611i 0 0 0 +17.3611i 0 0 0 -16.0000i 0 0 0 0 0 0 -17.0648i 0 0 0 0 0 -1.9422 +10.5107i 0 0 0 0 0 0 0 0 0 +17.0648i 0 0 Columns 6 through 9 0 0 0 0 0 0 +16.0000i 0 0 -1.9422 +10.5107i 0 0 0 0 -1.1876 + 5.9751i 0 0 0 2.8047 -35.4456i -1.6171 +13.6980i 0 电压幅值：电压相角：节点有功：节点无功：修正后的节点导纳矩阵Y： Columns 1 through 5 0 -20.6944i 0 0 0 +17.3611i 0 0 0 -19.3333i 0 0 0 0 0 0 -20.3982i 0 0 0 0 0 -1.9422 +10.5107i 0 0 0 0 0 0 0 0 0 +17.0648i 0 0 Columns 6 through 9 0 0 0 0 0 0 +16.0000i 0 0 -1.9422 +10.5107i 0 0 0 0 -1.1876 + 5.9751i 0 0 0 2.8047 -35.4456i -1.6171 +13.6980i 0 节点阻抗矩阵Z的第4列：准确计算结果：短路电流：模值：相角：节点电压模值：支路电流： i j Iij近似计算结果：短路电流：模值：相角：节点电压模值：五、程序流程图YNNY结 束输出结果计算潮流分布输出潮流不收敛信息k←k+1k>Kmax?Max|fi〔x[k]〕|<ɛ?求解修正方程获得∆x[k]x[k+1] ←x[k]+∆x[k]计算节点导纳矩阵参数设置节点电压初值x[0]设置k←0与最大迭代次数Kmax计算f〔x[k]〕与雅克比矩阵开 始输入电网与节点注入数据六、程序与输入文件 文件：Sheet19节点系统支路参数首节点编号i尾节点编号jR〔p.u.〕X〔p.u.〕B/2〔p.u.〕140027003900454657697889Sheet29节点系统发电机与负荷参数节点编号节点类型P〔p.u.〕Q〔p.u.〕U〔p.u.〕Xd〔p.u.〕E〔p.u.〕1平衡2PV3PV4PQ005PQ6PQ7PQ008PQ-19PQ00powerflow_cal.m 文件：l=9;%支路数n=9;%节点数m=6;%PQ节点数Yn=zeros(n);%初始化节点导纳矩阵YDATA1=xlsread('input_data.xls',1);%计算节点导纳矩阵Yfor k=1:l i=DATA1(k,1); j=DATA1(k,2); R=DATA1(k,3); X=DATA1(k,4); B2=DATA1(k,5); Yn(i,i)=Yn(i,i)+1i*B2+1/(R+1i*X); Yn(j,j)=Yn(j,j)+1i*B2+1/(R+1i*X); Yn(i,j)=Yn(i,j)-1/(R+1i*X); Yn(j,i)=Yn(j,i)-1/(R+1i*X);enddisp('节点导纳矩阵Yn：');disp(Yn);G=real(Yn);B=imag(Yn);DATA2=xlsread('input_data.xls',2);P=zeros(1,n);Q=zeros(1,n);U=ones(1,n);P(2:n)=DATA2(2:n,3);Q(4:n)=DATA2(4:n,4);U(1:3)=DATA2(1:3,5);%设置节点电压初值e(1)=DATA2(1,5);e(2:n)=1.0;f(1:n)=0.0;%设置迭代次数t=0;tmax=10;while t<=tmax%计算f(x) a(1:n)=0.0; c(1:n)=0.0;for i=2:nfor j=1:n a(i)=a(i)+G(i,j)*e(j)-B(i,j)*f(j); c(i)=c(i)+G(i,j)*f(j)+B(i,j)*e(j);endendfor i=2:n deltaP(i)=P(i)-e(i)*a(i)-f(i)*c(i);endfor j=4:n deltaQ(j)=Q(j)-f(j)*a(j)+e(j)*c(j);endfor k=2:3 deltaU2(k)=U(k)*U(k)-e(k)*e(k)-f(k)*f(k);end fx=[deltaP(2:n) delt。