MATLAB牛拉法计算潮流,matlab潮流计算

《matlab潮流计算》由会员分享，可在线阅读，更多相关《matlab潮流计算(14页珍藏版)》请在装配图网上搜索。

1、附录1使用牛顿拉夫逊法进行潮流计算的Matlab程序代码 牛拉法计算潮流程序- B1矩阵1、支路首端号；2、末端号；3、支路阻抗；4、支路对地电纳 5、支路的变比；6、支路首端处于K侧为1，1侧为0 B2矩阵1、该节点发电机功率；2、该节点负荷功率；3、节点电压初始值 4、PV节点电压V的给定值；5、节点所接的无功补偿设备的容量 6、节点分类标号1为平衡节点(应为1号节点)；2为PQ节点；3为PV节点；clear all;at long;n请输入节点数nodes;nl请输入支路数lines;isb请输入平衡母线节点号balance;pr请输入误差精度precision;B1请输入由各支路参数形。

2、成的矩阵B1;B2请输入各节点参数形成的矩阵B2;Yzerosn;ezeros1,n;fzeros1,n;Vzeros1,n;sidazeros1,n;S1zerosnl;for i1nl支路数 if B1i,60左节点处于1侧 pB1i,1;qB1i,2; else 左节点处于K侧 pB1i,2;qB1i,1;end Yp,qYp,q-1./B1i,3*B1i,5;非对角元 Yq,pYp,q; 非对角元 Yq,qYq,q1./B1i,3*B1i,52B1i,4;对角元K侧 Yp,pYp,p1./B1i,3B1i,4;对角元1侧end求导纳矩阵disp导纳矩阵 Y;dispY-GrealY;B。

3、imagY;分解出导纳阵的实部和虚部for i1n给定各节点初始电压的实部和虚部eirealB2i,3;fiimagB2i,3; ViB2i,4;PV节点电压给定模值endfor i1n给定各节点注入功率 SiB2i,1-B2i,2; i节点注入功率SG-SL Bi,iBi,iB2i,5;i节点无功补偿量end-PrealS;QimagS; 分解出各节点注入的有功和无功功率ICT10;IT21;N02*n;NN01;a0; 迭代次数ICT1、a；不满足收敛要求的节点数IT2while IT20 N02*n 雅可比矩阵的阶数；NN01扩展列 IT20;aa1;for i1n if iisb非平。

4、衡节点Ci0;Di0;for j11n CiCiGi,j1*ej1-Bi,j1*fj1;Gij*ej-Bij*fj DiDiGi,j1*fj1Bi,j1*ej1;Gij*fjBij*ejend P1Ci*eifi*Di;节点功率P计算eiGij*ej-Bij*fjfiGij*fjBij*ej Q1Ci*fi-ei*Di;节点功率Q计算fiGij*ej-Bij*fj-eiGij*fjBij*ej求i节点有功和无功功率P,Q的计算值 V2ei2fi2;电压模平方以下针对非PV节点来求取功率差及Jacobi矩阵元素- if B2i,63非PV节点 DPPi-P1;节点有功功率差 DQQi-Q1; 节。

5、点无功功率差以上为除平衡节点外其它节点的功率计算求取Jacobi矩阵for j11n if j1isb dP/de-dQ/df X2Bi,j1*ei-Gi,j1*fi; dP/dfdQ/de X3X2; X2dp/df X3dQ/de X4-X1; X1dP/de X4dQ/df p2*i-1;q2*j1-1; Jp,qX3;Jp,NDQ;mp1; X3dQ/de Jp,NDQ节点无功功率差 Jm,qX1;Jm,NDP;qq1; X1dP/de Jm,NDP节点