基于matlab的pq分解法5节点程序编程.docxVIP

％ 程序功能：5节点潮流计算PQ分解法 clear aｌl; clc; ％　ＮODE ＩNFO %　ＮODE TYPE　　 　ＰQ:１ 　 PV:2　 ＰH：3　 　所有发电机节点均被视为ＰV节点,节点5为平衡节点 NOＤE=［ % TYＰＥ 　 　 　　　V　　　　 　 　ＡNGEL　 　 PG 　 　 ＱG 　 　　　PL 　 QL 　 ＳHＵＮＴ Ｃ(节点并联电容) 　 　 　1 　 1.00００ ０ 　 　　0 　　 　 0 1６0.00０ ８0.00０0 　 　0 　 　1 　 １．０000 　 　 0 　　　 　 0 　 　 0 　 ２00．００0 10０．000 　 　 　 0 　 1 1.00００ 　 0　　 　 0 　 　 　０ 　370.0０0 130.00０ ０ 2 　 1.05０0 　 　0 　 500．０00　 　　 0 　 　 0 　 0 　 　 0 　 　３ 　 1.050０ 　０ 　0 　 0 ０　 　 　 0 　 0 ] ％ 功率为有名值，基准功率为10０MＶＡ % 线路对地电纳参数为B,不是Ｂ/2 % BRAＮＣH DATA 　IＥEE １4　ＮODE % 　 i j 　 　　 Ｒ　 　 X 　 　 　　 B　 　 K(变比) ＢＲANCＨ=[ 　 　 ２　 1 　 .0４ 　 　　 ．25　 　　 ．25*2　 　 1.０00; 　　 3 　 1　 　 　 　 ．1 　 　　 .35　 　 .00０0 　 　 1．000; 　 3 2 　 .０8　　 .3　 　 　 .25*2 　 1.00０; 　 5 　　　　3　 .０00０0 　 　 .03 　　 .0000　　 　 　１.0５0; 　 ４　　 2 　　 　 ．00000 　 .015　 　 　 .０000　　　　 　1.０50； 　 ] x ＝　sizｅ(NOＤＥ）; %x=［1４,８] size矩阵的行数和列数 n = x(１）; ｍ　= lｅngｔh(fiｎd(ＮODE(：,1)==１)); %lengtｈ列数 y=x(:,i）的作用就就是把矩阵ｘ的第i列取出来赋值到y里面，那么得到的y就是一个列向量。 　find稀疏矩阵非０值的个数。 y ＝ size（BRANCH);　 %y=［２0,６] Ｇ　= zerｏs(ｎ); 　 ％zｅros零矩阵的生成,ｎ阶方阵 B =　zeros(n）; ｆor i＝1：y（1) ＳＴ　　=　BＲANＣH(i,１); 　 %起始节点 　 　ＥD ＝ BＲAＮCＨ(i，２); %终止节点 R　　　=　BRANCH(ｉ,３); 　　 X = BRANＣＨ(i，4); BB　 = BＲＡNCＨ(i，5）； 　 　 K = 1／BRANCＨ(i,6); %变比,理想变压器变换为π形等值电路时,线路导纳要除以Ｋ,对地导纳除以k或者k＾2,该处已经形成倒数。 RX2 = R^2+X^2; 　 　 Ｇ（ST,EＤ） ＝ -K*R/RＸ2； Ｇ(ED,ST)　＝ -K*R/RX2; 　　 B(ST,ED） ＝ 　K*X／ＲX2； 　　Ｂ(ED,ST)　＝ K*X／RX2; 　 G(ST,ＳＴ)　=　G(ＳT，SＴ) + R/ＲX2; 　　 B(SＴ,ST)　= Ｂ(ST,ST）　- X/RX2； 　Ｂ(ST,ＳT)　=　B（ST,ST) +　BＢ/2; G(ED,ED） ＝ G(ＥＤ，EＤ) ＋　K^２＊R/ＲX2; 　％根据变压器等值变换电路，系数应该为（1-ｋ)/k＾2 or(k－1)/k，这里只有k^２ 　 B（ＥＤ,ED） = B(ＥＤ,EＤ) - K＾2*Ｘ/ＲX２； B(ED,ED）　= B(ED，ED)　+　BB/2; eｎd;