PQ分解法計(jì)算大電網(wǎng)潮流程序

1、PQ分解法計(jì)算大電網(wǎng) 潮流程序 -CAL-FENGHAI-(2020YEAR-YICAI) JINGBIAN function PQ %用卩(1分解法計(jì)算大電網(wǎng)潮流 % %bus數(shù)組1.節(jié)點(diǎn)編號(hào)2.節(jié)點(diǎn)電壓3.節(jié)點(diǎn)電壓角度4.注入有功5.注入無(wú)功6.節(jié)點(diǎn)類(lèi)型(1PQ 2PV 3 平衡) %line數(shù)組1.始端節(jié)點(diǎn)編號(hào)2.末端節(jié)點(diǎn)編號(hào)3.電阻4電抗5電導(dǎo)G 6電納B 7.變比 %打開(kāi)數(shù)據(jù)文件 clear clc bus=load(,!); line=load(); linenum(:/l/2)=line(:,1,2); (nb/v=size(bus); (nl/=size(line); % no

2、denum=(l:nb) bus(: J); %帶入子函數(shù)數(shù)據(jù)處理 busJinePCbnPVSWodenum =changel_busline( busjine );%對(duì)盯點(diǎn)晅新編匕 Y = admittance(busjine/l);%生成節(jié)點(diǎn)導(dǎo)納矩陣 Yl= admittancefbusjine );%生成化簡(jiǎn)條件3的矩陣B1 Y2=admittance(bus/line/3 );%生成化簡(jiǎn)條件 3 的矩陣 B2 %臨時(shí)添加的測(cè)試數(shù)據(jù) % nPQ二4; nPV=0;nSW=l;nb=5; % Y= +5i+5i+-5+15i % +5i 10+30i 0 + % +5i 10+30i +

3、 0 % + 0 + 0 % -5+15i+ 0 0 % %Yl=+5i+5i + -5+15i % +5i 10+30i 0 + % +5i 10+30i + 0 % + 0 + 0 % -5+15i+ 0 0 % % bus=l 101 %21 0 1 %31 01 %41 0 1 %5 0 0 0 3; % % line=5 20 0 %23 10 -30 0 0 0 %34 0 0 0 %41 0 0 0 %12 5 0 0 0 %13 5 0 0 0 %15 5 15 0 0 0； % bus_PV0=bus(nPQ+l):end/2),;%*ones(l/nPV+nSW); bus

4、_U=ones(l/nPQ) bus_PVO;%電JE幅值 bus_e=zeros(nb,l); %電壓角度 delta_P=zeros(nPQ+nPV/l); delta_Q=zeros(nPCbl)； % delta_e=zerosnb-l/l); %delta_U=zeros( nPQl); c=O；KP=l;KQ=l；%KP KQ用來(lái)判斷有功、無(wú)功是否收斂 G=real(Y);B=imag(Y);B10=imag(Yl);B20=imag(Y2);% BO 是進(jìn)化簡(jiǎn)三后的譏導(dǎo)納矩陣虛部 %形成解耦潮流的系數(shù)矩陣B1和B2 Bl=B10(l:nb-l,l:nb-l); B2=B2O(l:

5、nPCbl：nPQ); while c delta_e=-(UP#Bl)error_P; bus_e=bus_e+delta_e;O); c=c+l; KQ 二 1; else KP=O; if KCT二0 else break end end for ii=l:nPQ delta_Q(ii)=bus(ii/5); for jj=l:nb delta_Q(ii)=delta_Q(ii)-bus_U(ii)*bus_U(jj)#(G(ii/jj)*sin(bus_e(ii)-bus_e(jj)-B(ii/jj)#cos(bus_e(ii)- bus_e(jj); end end UQ=diag(

6、bus_U(l:( nbnPV nSW); erro r_Q 二 U Qd e 11 a_Q; if max(abs(error_Q) delta_U=-B2error_Q; bus_U=bus_U+delta_U;zeros( nPV+nSW)，l); c=c+l; KP=1; else KQ=O; if KP JO else break end end end %至此得到收斂的節(jié)點(diǎn)電壓值 % % %對(duì)計(jì)算結(jié)果進(jìn)行數(shù)據(jù)處理 %將節(jié)點(diǎn)結(jié)果用原節(jié)點(diǎn)編號(hào)農(nóng)示 bus_Ue=zeros(nb/3)； bus_Ue(:/l/2,3)=nodenum(:/2) bus_U bus_e/pi*180; f

7、or ii=l:nb for jj=ii+l:nb if bus_Ue(ii/l)bus_Ue(jj/l) t=bus_Ue(ii,:); bus_Ue(ii,:)二 bus_UG(jj,:); bus_Ue(jj，:)=t; end end end %r_U是收斂的電壓衣達(dá)成復(fù)數(shù)的形式 r_U=zeros(nbzl); for k=l:nb r_U(k)=bus_U(k)#(cos(bus_e(k)+li#sin(bus_e(k); end %計(jì)算平衡節(jié)點(diǎn)功率 sw_s=o； SW_S 二 SW_S+r_U(nb 廣 conj(Y( nb,:)廣 c on j(r_U); %計(jì)算各支路功率S

8、ij line_S=zeros( nb,nb);line_SO 二 zeros( nb，nb); for ii=l:nb for jj=l:nb line_S(ii/jj)=r_U(ii)*(conj(r_U(ii)*conj(Yii/ii)+(conj(r_U(ii)-conj(r_U(jj)*conj(Y(ii/jj); end end % %把線路結(jié)果還原成原節(jié)點(diǎn)編號(hào)對(duì)應(yīng)的結(jié)果 for ii=l:nb for jj=l:nb line_S0(nodenum(ii/2)/nodenum(jj/2)=line_S(ii/jj); end end line_P=real(line_SO);li

9、ne_Q=imag(line_SO); %計(jì)算各支路損耗 delta_S=zeros(nhl); for k=l:nl a=linenum(k,l);b 二 line nu m(k,2); delta_S(k)=line_SO(a/b)+line_S(b,a); end %計(jì)算網(wǎng)絡(luò)總損耗 SO=sum(delta_S); % %將計(jì)算結(jié)果輸入指定文件 fid=fopen(C:UserslrDesktopmatlab 練習(xí)訓(xùn)練題大電網(wǎng)潮流計(jì)算twt1)； fprintf(fid； ij點(diǎn)號(hào)t節(jié)點(diǎn)電Jk幅值t節(jié)點(diǎn)電壓角度n; for k=l:nb fprintf(fid/!%dt%ft%fn,/k

10、,bus_Ue(kzl)/bus_Ue(k/2); end fprintf(fid|支路首端t支路末端t支路仃功t支路無(wú)功t支路損耗W); for k=l:nl fprintf(fid/%dtt%dtt%At%ft%fn,Jinenum(k/l)/linenum(k/2)/line_P(linenum(k/l) ,linenum(k/2)/line_Q (linGnum(kJ)inenum(k,2),delta_S(k); end fprintf(fid；平衡節(jié)點(diǎn)功率=%fn：SW_S); fprintf(fid；網(wǎng)絡(luò)總損耗=%fn : SO); fclose(fid); end functi

11、on busJinePCbnPVSWodenum =changel_busline( busjine ) %此函數(shù)用來(lái)對(duì)原始輸入節(jié)點(diǎn)、線路數(shù)據(jù)進(jìn)行重新編號(hào) % %bus數(shù)組1.節(jié)點(diǎn)編號(hào)2.節(jié)點(diǎn)電壓3.節(jié)點(diǎn)電壓角度4.注入有功5.注入無(wú)功6.節(jié)點(diǎn)類(lèi)型(1PQ2PV3 平衡) %line數(shù)組1 始端節(jié)點(diǎn)編號(hào)2末端節(jié)點(diǎn)編號(hào)3電阻4電抗5電導(dǎo)G 6電納B 7變比 nb/sizbus); nl/=size(line); % nodenum=(l:nb) bus nPQ=O; nPV=0;nSW=0;%PQ=;PQ=zeros(nb/6);PZ=zeros(nb/6);SW=zeros(nb/6);%PQ

12、 PV 平衡節(jié)丄 的個(gè)數(shù) for k=l:nb switch bus(k,6) case 1 nPQ=nPQ+l; PQ(nPQ:)=bus(匕:); case 2 nPV=nPV+l; PV(nPV/:)=bus(k/:); case 3 nSW=nSW+l; SW(nSW,:)=bus(k/:); otherwise dispC節(jié)點(diǎn)數(shù)據(jù)類(lèi)型出錯(cuò)！ 7； end end %生成重新編號(hào)后的節(jié)點(diǎn)數(shù)據(jù)矩陣 bus=PQ;PV;SW; nodenum=(l:nb)bus(:zl);%第列為新的節(jié)點(diǎn)編號(hào)，第二列為對(duì)應(yīng)的舊節(jié)點(diǎn)編號(hào) bus(:,l)=(l:n; %至此實(shí)現(xiàn)廣節(jié)點(diǎn)數(shù)據(jù)的重新編號(hào) % %對(duì)

13、線路數(shù)據(jù)重新編號(hào) % nodenum=(l:nb)1 bus(:zl);%第列為新的節(jié)點(diǎn)編號(hào)，第二列為對(duì)應(yīng)的I口節(jié)點(diǎn)編號(hào) for ii=l:nl rl/v=find(nodenum(:/2)=line(ii/l); line(ii/l)=nodenum(rl/l); r2/v=find(nodenum(:/2)=line(ii/2); Iine(ii/2)=nodenum(r2/l); end end function Y = admittance(busjin已c ) %此函數(shù)用來(lái)形成節(jié)點(diǎn)導(dǎo)納矩陣 % %bus數(shù)組1 節(jié)點(diǎn)編號(hào)2節(jié)點(diǎn)電壓3節(jié)點(diǎn)電壓角度4注入有功5注入無(wú)功6節(jié)點(diǎn)類(lèi)型(1PQ 2

14、PV 3 平衡) %line數(shù)組1.始端節(jié)點(diǎn)編號(hào)2.末端節(jié)點(diǎn)編號(hào)3.電阻4電抗5電導(dǎo)6電納B/2 7.變比 %c是用來(lái)控制形成節(jié)點(diǎn)導(dǎo)納矩陣的方式的，c=l形成般的節(jié)點(diǎn)導(dǎo)納矩陣，可以用來(lái)確定B2, c=2,形成化簡(jiǎn)條件3的節(jié)點(diǎn)導(dǎo)納矩陣，確定B1 nbz=size(bus); (nl/AZ=size(line); Y=zeros (n b,nb); zt=zeros( nl,l);yt二zeros( nl，l);ym 二 ZGros(nl,l);l=Zros(nlJ);J=zeros( nl);K二zeros(nl,l); switch c case 1 for k=l:nl zt(k)=line

15、(kz3)+i 產(chǎn) line(匕 4); yt(k)=l/zt(k); ym(k)=line(k/5)+li*line(k/6); l(k)=line(k/l);J(k)=line(kz2);K(k)=line(k/7); end case 2 for k=l:nl yt(k)=l/(line(k/3)+li*line(k/4); ym=zeros(nbl); l(k)二 line(kJ);J(k)=line(k，2); if line(k,7)性 1 K(k)=l; else K(k)=line(k/7); end end case 3 for k=l:nl yt(k)=l/(li*line(k/4); ym(k)=line(kz5)+li*line(k/6); l(k)=line(k/l);J(k)=line(kz2);K(k)=line(k/7); end %針對(duì)不同線路設(shè)置節(jié)點(diǎn)導(dǎo)納值 for k=l:nl if (K(k)=l) Y(J(k),J(k)=Y(J(k),J(k)+yt(k)+ym(k); Y(l(k),J(k)=Y(l(k),J(k)-yt(k); Y(J(k),l(k)=Y(l(k),J(k); end if