12节点matlab潮流计算例题,电力系统分析大作业-matlab三机九节点潮流计算报告.docx...

电力系统分析大作业-matlab三机九节点潮流计算报告.docx

电力系统分析大作业一、设计题目本次设计题目选自课本第五章例58,美国西部联合电网WSCC系统的简化三机九节点系统，例题中已经给出了潮流结果，计算结果可以与之对照。取Ε000001。二、计算步骤第一步，为了方便编程，修改节点的序号，将平衡节点放在最后。如下图92132745683第二步，这样得出的系统参数如下表所示第三步，形成节点导纳矩阵。第四步，设定初值；01605403201UU，。87Q87第五步，计算失配功率0，125，09，0，1，0，163，01P0203P0405P0607P085；808614，02590，00420，06275，01710，01Q0203Q0405Q07101。6显然，。5108640|},MAX{|IIP第六步，形成雅克比矩阵(阶数为1414)第七步，解修正方程，得到00371，00668，01020300628，00732，00191，00422，01726，009456070808；00334，00084，00223，00372，002601U0203U0405U6，00400。从而00371，00668，00628，00732，00191，11213141500422，01726，00908；167810334，10084，10223，10372，10266，1041U11U11U100。然后转入下一次迭代。经三次迭代后。55101840|},MAX{|IIQP迭代过程中失配功率的变化情况如下表K0123ΔP10001060000100000000845ΔP2125003790000500000000896ΔP309004390000500000000778ΔP40004210001200000003421ΔP5100610000900000001845ΔP60002690000700000001631ΔP7163005790000400000000278ΔP8085003360000200000000103ΔQ108614005010000400000000561ΔQ20259007140001200000002774ΔQ30042004240000600000001236ΔQ406275018750002100000003279ΔQ50171002410000400000000805ΔQ607101008280000700000000799MAX16300610000900000001845迭代过程中节点电压变化情况如下表KU1U2U3U4U5U60111111110334100841022310372102661040021025909958101281025910160103243102580995610127102581015910324迭代收敛后各节点的电压和功率KUΘPQ110258221680000000000209956398881250005000310127368740900003000410258371970000000000510159072751000003500610324196670000000000710250928001630000665810250466480850001086910400000000716402705同课本上给出的潮流相比较，结果完全一致，证明计算过程与程序编写正确。最后得出迭代收敛后各支路的功率和功率损耗三、源程序及注释由于计算流程比较简单，所以编写程序过程中没有采用模块化的形式，直接按顺序一步步进行。DISP 【节点数】 N1XLSREAD XLS , A3A3 节点数DISP 【支路数】 NXLSREAD XLS , B3B3 支路数DISP 【精度】 ACCURACYXLSREAD XLS , B4B4 精度BRANCHXLSREAD XLS , E4K12 NODEXLSREAD XLS , M4S12 DATA_B1BRANCH支路参数DATA_B2NODE节点参数T1ZEROSN,2T2ZEROSN1,3ISQRT1ATSHORTFORJ1NT1J,1DATA_B1J,3DATA_B1J,41IT1J,2DATA_B1J,51IIJPIJQIJIIJPJIQJIIJIPLQL120409402289045720406803869056390002601579130307000103029950305401654034300001701551240843201131085450866200838084840023001969360594601346060200608201807061460013503153450763800080074470759001070075460004801150560241002430033680241800312023620000902118910716402705073630716402392073630000000312741630000665159161630000918159160000001583860850001086083600850001496083600000000410ENDFORJ1N1T2J,1DATA_B2J,1DATA_B2J,21IT2J,2DATA_B2J,3DATA_B2J,41IENDB1ZEROSN,6B2ZEROSN1,5FORJ1NB1J,1DATA_B1J,1B1J,2DATA_B1J,2B1J,3T1J,1B1J,4T1J,2B1J,5DATA_B1J,6B1J,6DATA_B1J,7ENDFORJ1N1B2J,1T2J,1B2J,2T2J,2B2J,3DATA_B2J,5B2J,4DATA_B2J,6B2J,5DATA_B2J,7ENDDISP 【支路参数矩阵】 B1显示支路参数矩阵DISP 【节点参数矩阵】 B2显示节点参数矩阵以上为从EXCEL中导入初值的程序YZEROSN1FORI1NIFB1I,60不含变压器的支路PB1I,1QB1I,2YP,QYP,Q1/B1I,3YQ,PYP,QYP,PYP,P1/B1I,305B1I,4YQ,QYQ,Q1/B1I,305B1I,4ELSE含有变压器的支路PB1I,1QB1I,2YP,QYP,Q1/B1I,3B1I,5YQ,PYP,QYP,PYP,P1/B1I,3YQ,QYQ,Q1/B1I,52B1I,3ENDENDDISP 【导纳矩阵】 Y显示导纳矩阵M0FORI1N1IFB2I,52MM1ENDENDMPQ节点个数L0FORI1N1IFB2I,51LL1ENDENDLPV节点个数MISMATCH_POWERZEROSLM2,1FORI1N11PJ0FORJ1N1PJPJB2I,3B2J,3REALYI,JCOSB2I,4B2J,4IMAGYI,JSINB2I,4B2J,4ENDMISMATCH_POWERI,1REALB2I,1REALB2I,2PJENDFORKN1LM2QJ0FORJ1