matlab牛頓拉夫遜法算潮流分析

《matlab牛頓拉夫遜法算潮流分析》由會員上傳分享，免費在線閱讀，更多相關(guān)內(nèi)容在行業(yè)資料-天天文庫。

1、MATLAB牛頓拉夫遜法算潮流分析%本程序的功能是用牛頓——拉夫遜法進行潮流計算%B1矩陣：1、支路首端號；2、末端號；3、支路阻抗；4、線路對地電納(或變壓器導(dǎo)納)；%5、支路的變比；6、支路首端處于K側(cè)為1，1側(cè)為0；%7、線路/變壓器標識（0/1）變壓器參數(shù)當(dāng)支路首端處于K側(cè)標識為1時歸算至末端側(cè)，0歸算至首端側(cè)%B2矩陣：1、該節(jié)點發(fā)電機功率；2、該節(jié)點負荷功率；%3、PQ節(jié)點電壓初始值，或平衡節(jié)點及PV節(jié)點電壓的給定值%4、節(jié)點所接無功補償并聯(lián)電容（感）的電納%5、節(jié)點分類標號：1為平衡節(jié)點（應(yīng)為1號節(jié)點）；2為PQ節(jié)點；3為PV節(jié)點；clear;isb=1;%input('請輸入

2、平衡母線節(jié)點號：isb=');pr=1e-5;%input('請輸入誤差精度：pr=');%---------------------------------------------------n=10;%input('請輸入節(jié)點數(shù)：n=');nl=10;%input('請輸入支路數(shù)：nl=');B1=[123.4+12.8i1.4e-4i100;145.1+19.2i2.1e-4i100;234.25+16i1.75e-4i100;454.25+16i7e-4i100;135.1+19.2i2.1e-4i100;645.95+22.4i9.8e-4i100;271.78+53.89i038

3、.5/23101;381.49+48.02i011/23101;491.49+48.02i011/23101;5102.46+70.17i038.5/23101]%input('請輸入由支路參數(shù)形成的矩陣：B1=');B2=[00225.501;0022002;0022002;0022002;0022002;120023103;061.11+37.87i3502;047.53+29.46i1002;054.32+33.66i1002;040.74+25.25i3502]%input('請輸入各節(jié)點參數(shù)形成的矩陣：B2=');%---------------------------------

4、----------------------------%n=4;%input('請輸入節(jié)點數(shù)：n=');%nl=4;%input('請輸入支路數(shù)：nl=');%B1=[124+16i0100;%134+16i0100;%232+8i0100;%241.49+48.02i011/11001]%input('請輸入由支路參數(shù)形成的矩陣：B1=');%B2=[0011501;%0011002;%020+4i11002;%010+6i1002]%input('請輸入各節(jié)點參數(shù)形成的矩陣：B2=');%-------------------------------------------------

5、------------Y=zeros(n);e=zeros(1,n);f=zeros(1,n);V=zeros(1,n);sida=zeros(1,n);S1=zeros(nl);%%%-----------求導(dǎo)納矩陣------------------------fori=1:nl%從1到n1（總支路數(shù)）ifB1(i,7)==1%-----------如果是變壓器支路--------ifB1(i,6)==0%左節(jié)點(首端)處于1側(cè)p=B1(i,1);q=B1(i,2);else%左節(jié)點(首端)處于K側(cè)p=B1(i,2);q=B1(i,1);endY(p,q)=Y(p,q)-1./(B1(

6、i,3)*B1(i,5));%非對角元Y(q,p)=Y(p,q);%非對角元Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)^2);%對角元K側(cè)Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4);%對角元1側(cè)+勵磁導(dǎo)納else%------------否則為線路支路--------------------p=B1(i,1);q=B1(i,2);Y(p,q)=Y(p,q)-1./B1(i,3);%非對角元Y(q,p)=Y(p,q);%非對角元Y(q,q)=Y(q,q)+1./B1(i,3)+B1(i,4)./2.0;%對角元j側(cè)+線路電納的一半Y(p,p)=Y(

7、p,p)+1./B1(i,3)+B1(i,4)./2.0;%對角元i側(cè)+線路電納的一半endenddisp('導(dǎo)納矩陣Y=');disp(Y);%-----------給定各節(jié)點初始電壓及給定各節(jié)點注入功率--------------------------G=real(Y);B=imag(Y);%分解出導(dǎo)納陣的實部和虛部fori=1:n%給定各節(jié)點初始電壓的實部和虛部e(i)=real(B2(i,3));f