电力系统的潮流计算

电力系统的潮流计算,以下程序参考文献电力系统毕业设计中国水利电力出版社该文献用极坐标下的牛顿拉夫逊方法实现，在此为了与课本一致做了修改为了计算方便将原来的下标做以下修改S2S3S4S5U2U3U4U5改为S1S2S3S4U1U2U3U4,即原题的平衡点1就变为现在的平衡点51形成节点导纳矩阵,YB55625018750JYB51500015000JYB5212503750JYB5300000000JYB5400000000JYYB111083432500JYB1216675000JYB1316675000JYB1425007500JYB2512503750JYB2116675000JYB221291738750JYB231000030000JYB2400000000JYB3500000000JYB3116675000JYB321000030000JYB331291738750JYB3412503750JYB4500000000JYB4125007500JYB4200000000JYB4312503750JYB44375011250JYBYB11YB12YB13YB14YB15YB21YB22YB23YB24YB25YB31YB32YB33YB34YB35YB41YB42YB43YB44YB45YB51YB52YB53YB54YB55计算各节点功率的不平衡量设UEJF；YGBJ；E1100E2100E3100E4100F10F20F30F40GREALYBBIMAGYB设SPBJ；S1020020IS2045015IS3040005IS4060010IPREALSQIMAGSK0PRECISION1N14WHILEPRECISION000000001E5106F50FORM1N1FORN1N11计算PI,QI,设PIPTQIQTPTNEMGM,NENBM,NFNFMGM,NFNBM,NENQTNFMGM,NENBM,NFNEMGM,NFNBM,NENEND设P,Q的改变量为DP,DQDPMPMSUMPTDQMQMSUMQTENDFORM1N1FORN1N11计算HIJNIJJIJLIJHM,NBM,NEMGM,NFMNM,NGM,NEMBM,NFMJM,NBM,NFMGM,NEMLM,NGM,NFMBM,NEMENDENDFORM1N1FORN1N11BINGM,NFNBM,NENAINGM,NENBM,NFNEND计算HII,NII,JII,LII,由公式444B左侧公式实现,SUMAI，SUMBI用于实现公式中的SIGERMA从J到N的求和HM,MSUMBIBM,MEMGM,MFM2GM,MFMNM,MSUMAIGM,MEMBM,MFM2GM,MEMJM,M2BM,MFMSUMAIGM,MEMBM,MFMLM,M2BM,MEMSUMBIBM,MEMGM,MFMEND设雅可比矩阵为JJ，以下语句用来实现雅可比矩阵中对角线上元素HNJL的排列FORM1N1JJ2M1,2M1HM,MJJ2M1,2MNM,MJJ2M,2M1JM,MJJ2M,2MLM,MEND以下语句用于实现雅可比矩阵非对角线上元素的排列FORM1N1FORN1N1IFMNELSEHM,NBM,NEMGM,NFMNM,NGM,NEMBM,NFMJM,NBM,NFMGM,NEMLM,NGM,NFMBM,NEMJJ2M1,2N1HM,NJJ2M1,2NNM,NJJ2M,2N1JM,NJJ2M,2NLM,NENDENDEND设由P,Q的改变量组成的81矩阵为PQ，由E,F的改变量组成的81矩阵为DUFORM1N1PQ2M1DPMPQ2MDQMENDDUINVJJPQPRECISIONMAXABSDUFORN1N1FNFNDU2N1ENENDU2NENDFORN1N11UNENFNJENDKK1K1,DUDU,PQ,UEND计算S5,也就是题目中的S1，即平衡节点功率FORM1N11IMYB5,MUMENDS5U5SUMCONJI设网络总损耗为SS,计算输电效率EFFICIENCYFORM1N11S0MSMPMREALSMENDSSSUMS0EFFICIENCYA