《電力系統(tǒng)分析》(附錄Matlab程序)公開課_第1頁
《電力系統(tǒng)分析》(附錄Matlab程序)公開課_第2頁
《電力系統(tǒng)分析》(附錄Matlab程序)公開課_第3頁
《電力系統(tǒng)分析》(附錄Matlab程序)公開課_第4頁
《電力系統(tǒng)分析》(附錄Matlab程序)公開課_第5頁
全文預(yù)覽已結(jié)束

下載本文檔

版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請進行舉報或認領(lǐng)

文檔簡介

附錄Matlab程序

第四章牛頓?拉夫遜法潮流計算:

%牛頓-拉夫遜法潮流計算

yO=O.25j;%接地導納都為此值

z24=O.OI5j;

zl2=0.04+0.25j;

z23=0.()8+0.3j;

zl3=O.l+O.35j;

z35=0.03j;

k=1.05;

%求導納矩陣

Yll=yO+l/zl3+l/zl2;

Y12=-l/zl2;

Y21=Y12;

Y13=-l/zl3;

Y31=Y13;

Y14=0;

Y41=0;

Y15=0;

Y51=0;

Y22=y()+I/ZI2+1/z23+y()+1/z24/kA2;

Y23=-l/z23;

Y32=Y23;

Y24=-l/z24/k;

Y42=Y24;

Y25=0;

Y52=0;

Y33=yO+I/ZI3+1/z23+1Zz35/kA2;

Y35=-l/z35/k;

Y53=Y35;

Y34=0;

Y43=0;

Y44=1/Z24;

Y45=0;

Y54=0;

Y55=l/z35;

YY=[Y!1,Y12,Y13,YI4,Y15;Y21,Y22,Y23,Y24,Y25;Y31,Y32,Y33,Y34,Y35;Y41,Y42,Y43,Y4

4.Y45;Y51,Y52,Y53,Y54,Y55]

%形成Bl矩陣(導納矩隴去掉平衡節(jié)點形成)

%B1=[YI1,Y12,YI3,Y14;Y21,Y22,Y23,Y24;Y31,Y32,Y33,Y34;Y41,Y42,Y43,Y44];

%B1l=imag(Bl)

%形成B2矩陣(導納矩除去掉平衡節(jié)點和PV節(jié)點形成)

%B2=[Y1LY12,Y13;Y21,Y22,Y23;Y31,Y32,Y33];

%B22=imag(B2)

G=real(YY);

B=imag(YY);

%賦電壓初值

U=[1.0,1.0,1.Oj.05,1.05]';

Dlt=[0,0,0,0,0],;

P=[-1.6,-2,-375]';

Q=[-0.8,-L-1.3]*;

Ep=le-4;%迭代精度為0.0001

fork=1:10

forii=l:4

aa=0;

forjj=l:5

aa=aa+U(ii)*U(jj)*(G(ii,jj)*cos(Dlt(ii)-Dlt(jj))+B(ii,jj)*sin(Dlt(ii)-Dlt(jj)));

end

dP(ii)=P(ii)-aa;

end

forii=l:3

bb=0;

forjj=l:5

bb=bb+U(ii)*U(jj)*(G(ii,jj)*sin(Dlt(ii)-Dlt(jj))-B(ii,jj)*cos(Dlt(ii)-Dlt(jj)));

end

dQ(ii)=Q(ii)-bb;

end

%求雅可比矩陣

fori=1:4

forj=1:4

if(j==i)

bb=0;

forjj=l:5

ifg)

bb=bb;

else

bb=bb+U(i)*UGj)*(G(ijj)*sin(Dlt(i)-Dlt(jj))-B(i,jj)*cos(Dlt(i)-Dlt(jj)));

end

end

H(iJ)=bb;

N(ij)=-U(i)*U(j)*G(ij)-P(i);

else

H(ij)=-U(i)*U(j)*(G(ij)*sin(Dlt(i)-Dlt(j))-B(ij)*cos(Dlt(i)-Dlt(j)));

N(i,j)=-U(i)s)Ud)*(G(i,j)*cos(Dlt(i)-Dk(j))+B(ij)*sin(Dlt(i)-Dlt(j)));

end

J(i,j)=H(i,j);

ifj<=3

J(i,j+4尸N(i,j);

end

end

end

fori=l:3

forj=1:4

ifj==i

K(i,j)=U(i)*U(i)*G(ij)日(i);

L(iJ)=U(i)*U(i)*B(ij)-Q(i);

else

K(iJ)=-N(i,j);

L(ij)=H(i.j);

end

J(i+4j)=K(ij);

ifj<=3

J(i+4j+4)=L(i,j);

end

end

end

%解修正方程

XX=inv(J)*[dRdQ],;

fori=l:4

dDlt(i)=XX(i);

ifi<=3

dU(i)=XX(i+4)*U(i);

U(i)=U(i)-dU(i);

end

end

dP

dQ

Dl(

U

C=max(abs(dP));

CC=max(abs(dQ));

if(C<Ep&&CC<Ep)

break

end

end

%求平衡節(jié)點的輸入功率

Pn=0;

forjj=l:5

YY=IY11,Y12,Y13,Y14,Y15;Y21,Y22,Y23,Y24,Y25;Y31,Y32,Y33,Y34,Y35;Y41,Y42,Y43,Y4

4.Y45;Y5LY52,Y53,Y54,Y55]

%形成Bl矩陣(導納矩隴去掉平衡節(jié)點形成)

B1=[Y11,Y12,Y13,Y14;Y21,Y22,Y23,Y24;Y31,Y32,Y33,Y34;Y41,Y42,Y43,Y44];

Bll=imag(Bl)

%形成B2矩陣(導納矩I笆去抻平衡節(jié)點和PV節(jié)點形成)

B2=[Y11,Yl2,Y13;Y21,Y22,Y23;Y31,Y32,Y33];

B22=imag(B2)

G=reaI(YY);

B=imag(YY);

%賦電壓初值

U=[1.0,l.oJ.oJ.05J.05]';

Dk=[0,0,0,0,0]';

P=[-1.6,-2,-3.7,5T;

,

Q=[-0.8rl,-1.3J;

Ep=le-4%迭代精度為0.0001

fork=1:10()

forii=l:4

aa=0;

forjj=l:5

aa=aa+U(ii)*U(jj)*(G(ii,jj)*cos(Dlt(ii)-Dlt(jj))+B(ii,jj)*sin(Dlt(ii)-Dlt(jj)));

end

dP(ii)=P(ii)-aa;

end

forii=l:4

UP(ii)=dP(ii)/U(ii);

end

dDlt=(inv(BU))*UP';

forii=l:4

dDlt(ii)=dDlt(ii)/U(ii);%求出相位角修正值

Dlt(ii)=Dlt(ii)-dDlt(ii);

end

forii=l:3

bb=0;

forjj=l:5

bb=bb+U(ii)*U(jj)*(G(iijj)*sin(Dlt(ii)-Dlt(jj))-B(ii,jj)*cos(Dlt(ii)-Dlt(jj)));

end

dQ(ii)=Q(ii)-bb;

溫馨提示

  • 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
  • 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會有圖紙預(yù)覽,若沒有圖紙預(yù)覽就沒有圖紙。
  • 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
  • 5. 人人文庫網(wǎng)僅提供信息存儲空間,僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護處理,對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對任何下載內(nèi)容負責。
  • 6. 下載文件中如有侵權(quán)或不適當內(nèi)容,請與我們聯(lián)系,我們立即糾正。
  • 7. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

評論

0/150

提交評論