電力系統(tǒng)穩(wěn)態(tài)分析基于matlab的潮流計算快速分解法-資料下載頁
2025-01-13 14:14本頁面
【正文】 er=x(:,2)。 EndPoint=x(:,3)。 EndNumber=x(:,4)。 R=x(:,5)。 X=x(:,6)。 B=x(:,7)。 k=x(:,8)。 L=length(HeadNumber)。 Y=zeros(L,L)。for i=1:L。 m=HeadNumber(i)。n=EndNumber(i)。 if n~=0。 Y(m,m)=Y(m,m)+1j*B(i)+1/(R(i)+1j*X(i))。 Y(n,n)=Y(n,n)+1j*B(i)+1/(R(i)+1j*X(i))。 Y(m,n)=Y(m,n)1/(R(i)+1j*X(i))。 Y(n,m)=Y(n,m)1/(R(i)+1j*X(i))。 else Y(m,m)=Y(m,m)+R(i)+1j*X(i)。 end。end。maxm=HeadNumber(1)。for i=1:L。 if maxm=HeadNumber(i)。 maxm=HeadNumber(i)。 end。end。maxn=EndNumber(1)。for i=1:L。 if maxn=EndNumber(i)。 maxn=EndNumber(i)。 end。end。Y=Y(1:max(maxm,maxn),1:max(maxm,maxn))。形成矩陣function Ba=B1(x,y)x=[1,1,2,2,0。 1,1,3,3,0。 2,2,3,3,0。 2,2,4,4,0,。 3,3,5,5,0,]。y=[ 0,5,0,0。 1,4,0。 2,3,1,0。 2,2,1,0。 2,1,1,0]。X=x(:,6)。for i=1:4Ba(i)=1/X(i)。i=i+1。if i==4 Ba(i)=1/X(i+1)。end。end。Bc=zeros(4,4)。for m=1:4。 for n=1:4。 if m==n。 Bc(m,m)=Ba(m)。 else Bc(m,n)=0。 end。 end。end。Ba=Bc。形成矩陣function Bb=B11(x,y)x=[1,1,2,2,0。 1,1,3,3,0。 2,2,3,3,0。 2,2,4,4,0,。 3,3,5,5,0,]。y=[ 0,5,0,0。 1,4,0。 2,3,1,0。 2,2,1,0。 2,1,1,0]。Y=Yn(x,y)。B=imag(Y)。Pointstyle=y(:,1)。i=1。j=1。while Pointstyle(i)=1。 i=i+1。 j=j+1。end。Bb=B(j:end,j:end)。程序2運行結(jié)果：運行結(jié)果分析：通過另外一種算法，我們同樣得到一樣的運行結(jié)果，可見快速分解法的潮流計算還是較為可靠的。5．總結(jié)快速分解潮流算法，實際上是“定雅可比”潮流計算的另一種簡化迭代格式。本文介紹了如何從“定雅可比”潮流迭代公式推導(dǎo)快速分解潮流算法的迭代公式。實際上，有功功率和無功功率的耦合作用已隱含在快速分解法的特殊迭代過程中，因此它具有良好的收斂性。需要注意的是，快速分解法不是PQ分解法，它是分別求兩個低階的修正方程，系數(shù)矩陣、是常數(shù)矩陣，在迭代過程中保持不變，所以這種算法速度快，在電網(wǎng)實時計算中得到廣泛應(yīng)用。[1] 張伯明,陳壽孫,:清華大學(xué)出版社,2000 [2] Stott B,Alsac decoupled load ,(3):859869[3] Monticelli A,et decoupled load flow:hypothesis,derivations and power system,(4):14251431[4] Van Amerongen R A generalpurpose version of the fast decoupled load power system,(2):760770[5] (I)——(自然科學(xué)報),1986,26(4):7685[6] (II)——在線潮流計算中的應(yīng)用,1986,26(5):2231[7] 柳春光,李光華,2010
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1