matlab潮流计算仿真.doc

附录AMATLAB程序本程序的功能是用牛顿拉夫逊法进行潮流计算B1矩阵1、支路首端号；2、末端号；3、支路阻抗；4、支路对地电纳5、支路的变比；6、支路首端处于K侧为1，1侧为0B2矩阵1、该节点发电机功率；2、该节点负荷功率；3、节点电压初始值4、PV节点电压V的给定值；5、节点所接的无功补偿设备的容量6、节点分类标号CLEARN10 请输入节点数N NL11 请输入支路数NL ISB1 请输入平衡母线节点号ISB PR000001 请输入误差精度PR B1121755E24155E2I026I10143159E27479E2I01215I10163159E27479E2I01215I1023368E3011135I00909145368E3011135I009091462808E26648E2I0108I106730865E300833I009091683159E27479E2I01215I106102457E25817E2I00945I108930865E300833I0090918102808E26648E2I0108I10 请输入由支路参数形成的矩阵B1 B2001051050100100200603718I10020010020040247I1002001002003502169I100200100200503099I100208010510503 请输入各节点参数形成的矩阵B2 YZEROSNEZEROS1,NFZEROS1,NVZEROS1,NSIDAZEROS1,NS1ZEROSNL－－－－－－－修改部分－－－－－－－－－－－－YM0SB100UB220YM 您输入的参数是标么值(若不是则输入一个不为零的数值) IFYM0SB 请输入功率基准值SB UB 请输入电压基准值UB YBSB/UB/UBBB1B1BB2B2FORI1NLB1I,3B1I,3YBB1I,4B1I,4/YBENDDISP B1矩阵B1 DISPB1FORI1NB2I,1B2I,1/SBB2I,2B2I,2/SBB2I,3B2I,3/UBB2I,4B2I,4/UBB2I,5B2I,5/SBENDDISP B2矩阵B2 DISPB2ENDFORI1NL支路数IFB1I,60左节点处于低压侧PB1I,1QB1I,2ELSEPB1I,2QB1I,1ENDYP,QYP,Q1/B1I,3B1I,5非对角元YQ,PYP,QYQ,QYQ,Q1/B1I,3B1I,52B1I,4/2对角元K侧YP,PYP,P1/B1I,3B1I,4/2对角元1侧END求导纳矩阵DISP 导纳矩阵Y DISPYGREALYBIMAGY分解出导纳阵的实部和虚部FORI1N给定各节点初始电压的实部和虚部EIREALB2I,3FIIMAGB2I,3VIB2I,4PV节点电压给定模值ENDFORI1N给定各节点注入功率SIB2I,1B2I,2I节点注入功率SGSLBI,IBI,IB2I,5I节点无功补偿量ENDPREALSQIMAGSICT10IT21N02NNN01A0WHILEIT20IT20AA1FORI1NIFIISB非平衡节点CI0DI0FORJ11NCICIGI,J1EJ1BI,J1FJ1ΣGIJEJBIJFJDIDIGI,J1FJ1BI,J1EJ1ΣGIJFJBIJEJENDP1CIEIFIDI节点功率P计算EIΣGIJEJBIJFJFIΣGIJFJBIJEJQ1CIFIEIDI节点功率Q计算FIΣGIJEJBIJFJEIΣGIJFJBIJEJ求P ,Q V2EI2FI2电压模平方以下针对非PV节点来求取功率差及JACOBI矩阵元素IFB2I,63非PV节点DPPIP1节点有功功率差DQQIQ1节点无功功率差以上为除平衡节点外其它节点的功率计算求取JACOBI矩阵FORJ11NIFJ1ISBDP/DEDQ/DFX2BI,J1EIGI,J1FIDP/DFDQ/DEX3X2X2DP/DFX3DQ/DEX4X1X1DP/DEX4DQ/DFP2I1Q2J11JP,QX3JP,NDQMP1JM,QX1JM,NDPQQ1JP,QX4JM,QX2ELSEIFJ1IDP/DEX2DIBI,IEIGI,IFIDP/DFX3DIBI,IEIGI,IFIDQ/DEX4CIGI,IEIBI,IFIDQ/DFP2I1Q2J11JP,QX3JP,NDQ扩展列△QMP1JM,QX1QQ1JP,QX4JM,NDP扩展列△PJM,QX2ENDENDELSE下面是针对PV节点来求取JACOBI矩阵的元素DPPIP1PV节点有功误差DVVI2V2PV节点电压误差FORJ11NIFJ1ISBDP/DEX2BI,J1EIGI,J1FIDP/DFX50X60P2I1Q2J11JP,QX5JP,NDVMP1JM,QX1JM,NDPQQ1JP,QX6JM,QX2ELSEIFJ1IDP/DEX2DIBI,IEIGI,IFIDP/DFX52EIX62FIP2I1Q2J11JP,QX5JP,NDVMP1JM,QX1JM,NDPQQ1JP,QX6JM,QX2ENDENDENDENDEND以上为求雅可比矩阵的各个元素FORK3N0N02N(从第三行开始，第一、二行是平衡节点)K1K1N1NNN01即N2N1扩展列△P、△QFORK2K1N1扩展列△P、△QJK,K2JK,K2/JK,K非对角元规格化ENDJK,K1对角元规格化IFK3不是第三行K4K1FORK33K4用K3行从第三行开始到当前行前的K4行消去FORK2K1N1K3行后各行下三角元素JK3,K2JK3,K2JK3,KJK,K2消去运算ENDJK3,K0ENDIFKN0BREAKENDFORK3K1N0FORK2K1N1JK3,K2JK3,K2JK3,KJK,K2消去运算ENDJK3,K0ENDELSEFORK3K1N0FORK2K1N1JK3,K2JK3,K2JK3,KJK,K2消去运算ENDJK3,K0ENDENDEND上面是用线性变换方式将JACOBI矩阵化成单位矩阵FORK32N01LK1/2ELELJK,N修改节点电压实部K1K1FLFLJK1,N修改节点电压虚部END修改节点电压FORK3N0DETABSJK,NIFDETPR电压偏差量是否满足要求IT2IT21不满足要求的节点数加1ENDENDICT2AIT2ICT1ICT11END用高斯消去法解“WJV“DISP 迭代次数 DISPICT1DISP 没有达到精度要求的个数 DISPICT2FORK1NVKSQRTEK2FK2SIDAKATANFK/EK180/PIEKEKFKJEND计算各输出量DISP 各节点的实际电压标幺值E为节点号从小到大排列 DISPEEEEUBDISPEEDISP DISP 各节点的电压大小V为节点号从小到大排列 DISPVVVVUBDISPVVDISP DISP 各节点的电压相角SIDA为节点号从小到大排列 DISPSIDAFORP1NCP0FORQ1NCPCPCONJYP,QCONJEQENDSPEPCPENDDISP 各节点的功率S为