電力系統(tǒng)穩(wěn)態(tài)分析基于matlab的潮流計算快速分解法-文庫吧資料
2025-01-19 14:14本頁面
【正文】 (p))+j*U(p)*sin(O(p)))*sum3。 end end sum3=0+j*0。 %求dQ/U絕對值的最大值 if max2=e %如果最大值滿足要求,則eq置為1,表示收斂 eq=1。%求得包括PQ節(jié)點電壓的電壓列向量 dU=B2*(dQi./Upq)。%Uq為中間變量 Uq(npq)=[]。 end dQi(m)=Q(pq(m))sum2。 end for m=1:pqn %求dQi sum2=0。%求dP/U絕對值的最大值 for m=1:(x2) if abs(dPi(m)/U(np(m)))abs(dPi(m+1)/U(np(m+1))) max1=abs(dPi(m+1)/U(np(m+1)))。%求相角O的不平衡量 for m=1:(x1) %求相角O的新迭代值矩陣 O(np(m))=O(np(m))+dO(m)。%將平衡節(jié)點所在行置空 Unp=Up。 end Up=U。 for n=1:x sum1=sum1+U(np(m))*U(n)*(G(np(m),n)*cos(O(np(m))O(n))+B(np(m),n)*sin(O(np(m))O(n)))。(k100)) k=k+1。%求PQ節(jié)點的個數(shù)(因B1矩陣的維數(shù)等于PQ節(jié)點數(shù)) dQi=zeros(pqn,1)。%找出PQ節(jié)點編號 pqn=size(B2)。%找出非平衡節(jié)點編號 dPi=zeros(x1,1)。 %計算P不平衡量dPi的收斂標志(0表示不收斂,1表示收斂) eq=0。%求得B2矩陣 %計算各節(jié)點有功功率不平衡量dPi k=0。%非PQ節(jié)點編號對應行置空 B22 B2=B22。 %BB矩陣為中間變量 B22(:,npq)=[]。及其逆矩陣B2 npq=find(TYPE(:,1)1)。%B1求逆后得B1矩陣 %%求B39。%平衡節(jié)點編號對應列置空 B11 B1=B11。%找出平衡節(jié)點編號 B11(:,p)=[]。 end。 end。 end。 for m=1:y B11(I(m),J(m))=1/(Xij(m))。矩陣及其逆矩陣B1B=imag(Y)。YG=real(Y)。 end。for m=1:x %求Y中的Yii for n=1:y if(I(n)==m|J(n)==m)Y(m,m)=Y(m,m)Y(I(n),J(n))+j*B0(n)/2。 Y(J(m),I(m))=1/Zij(m)。KT=branch(:,6)。Zij=Rij+j*Xij。Rij=branch(:,3)。I=branch(:,1)。P=point(:,4)。%將point矩陣的第一列賦給TYPE,以下類似U=point(:,2)。%B0為n*1階線路對地電納值K=zeros(y,1)。%G為n階節(jié)點電導方陣B=zeros(x)。%Yij為線路阻抗Y=zeros(x)。%R為線路電阻Xij=zeros(y,1)。%I為起始節(jié)點編號矩陣J=zeros(y,1)。%P為節(jié)點有功功率Q=zeros(x,1)。%U為節(jié)點電壓矩陣a=zeros(x,1)。%支路矩陣TYPE=zeros(x,1)。 3 5 0 。 2 3 0。%節(jié)點矩陣branch=[1 2 0。 3 0 。 1 1 0 。%支路數(shù)ye=。function [Y,U,P,Q,dSij, O,S,Sij,Sji,sumdS]=kuaisufenjie1776_1()%輸出：U——節(jié)點電壓，P節(jié)點有功，Q節(jié)點無功，dSij支路功率損耗，O為相角；Sij從節(jié)點i流向節(jié)點j的功率,S節(jié)點復功率,sumdS網(wǎng)絡總損耗；point為節(jié)點矩陣，branch為支路矩陣；[x]=5。利用快速分解法，當?shù)β势顦藴蕿?，迭代結束。 3 快速分解法在電力系統(tǒng)安全分析中的應用某地區(qū)電力供應系統(tǒng)由2個發(fā)電廠，4個變電站組成，系統(tǒng)模型見圖2，其中節(jié)點3為PQ節(jié)點，4為PV節(jié)點，5為平衡節(jié)點。該例中，和分別為 可見，比更接近于，而用代替即得到XB型快速分解法。這可以用下面的例子來說明。對于環(huán)形網(wǎng)絡，如果電網(wǎng)是均一網(wǎng)，即對任一條支路有，則得并有所以，故有 （18）如果電網(wǎng)不是均一網(wǎng)，上述
環(huán)評公示相關推薦
鄂ICP備17016276號-1