【文章內(nèi)容簡介】 內(nèi)蒙古科技大學(xué)電力系統(tǒng)穩(wěn)態(tài)分析課程設(shè)計說明書14 0 ] 由此可見，網(wǎng)絡(luò)的節(jié)點(diǎn)導(dǎo)納矩陣雖為奇異矩陣，但它的虛數(shù)部分的子陣 Bˊ或 B〞則是非奇異矩陣，可以求逆，其逆陣為= =139。)(?B139。[ 0。057143 ] iP? 取 =， =0； = = = =； = = = =0，按1U1?)0(2U)(3)0(4)(5)0(2?)(3)0(4)(5?下式計算各節(jié)點(diǎn)有功功率不平衡量 = ??0iP?ip????njji10j????00sincosjjijijBG??=—（ +）—??02（—）—（+）—1.（+）—(+)=相似地可得 =； =； =??03P???04P???05P? （弧度）i? 由方程式 — ( )=U139。)?B??0可得=—； =—； =—； =—?)0(2??)0(3??)0(4??)0(5?內(nèi)蒙古科技大學(xué)電力系統(tǒng)穩(wěn)態(tài)分析課程設(shè)計說明書15；所以=—； =—； =—； =—)1(2?)1(3?)1(4?)1(5?； iQ?按下式計算各節(jié)點(diǎn)無功功率不平衡量 = ??0ii????njjiU10j????11cossinijijjjBG??=—[—（—0）—（—0）]??02Q?—[+]—[（+）—（+）] —[（—+0,091810）—（0,036952+）]—[（—+）—（+）]=相似地可得 =—； =—； =—；??03Q???04Q???05Q? iU 由方程式 — ( )= U139。)?B?可得=； =； =； =；??02U???03U??04??05U?所以=； =； =； =；??12??13??14U??15求得各節(jié)點(diǎn)電壓的新值后，就開始第二次迭代。每次迭代所得適于表 1~表3。由表 1 可見經(jīng) 6 次迭代就可滿足 ε≤ 的要求。50?內(nèi)蒙古科技大學(xué)電力系統(tǒng)穩(wěn)態(tài)分析課程設(shè)計說明書16表 1 迭代過程中各節(jié)點(diǎn)功率的不平衡量K )(2kP???kQ2)(3kP???kQ3)(4kP???kQ4)(5kP???kQ50 1 2 3 4 0,000024 5 6 表 2 迭代過程中各節(jié)點(diǎn)電壓的修正量k ?)(k???kU2?)(3k???kU3?)(4k???kU4?)(5k???kU50 1 2 3 4 5 6表 3 迭代過程中各節(jié)點(diǎn)電壓k )(2k???kU2)(3k???kU3)(4k???kU4)(5k???kU50 1 0,108571 2 0,089348 0,103974 3 4 內(nèi)蒙古科技大學(xué)電力系統(tǒng)穩(wěn)態(tài)分析課程設(shè)計說明書175 6 0,083906 = ； =U??U?? = ； = 455 9265 和線路功率~1S~ijS 迭代收斂后，就可以計算平衡節(jié)點(diǎn)功率和線路功率。結(jié)果如下 =+~1各線路功率如表 4表 4 各線路功率~ijSj i1 2 3 4 51 ++2 +++3 4 ++5 計算完畢。內(nèi)蒙古科技大學(xué)電力系統(tǒng)穩(wěn)態(tài)分析課程設(shè)計說明書18第三章 程序設(shè)計 設(shè)計流程3．2 程序設(shè)計針對兩機(jī)五節(jié)點(diǎn)網(wǎng)絡(luò)潮流計算基于 MATLAB 的程序如下：Clc %清除當(dāng)前屏幕上顯示的內(nèi)容，但并不清楚工作空間的數(shù)據(jù)Clear %清除工作空間內(nèi)的所有變量disp(39。節(jié)點(diǎn)總數(shù)為：39。)。 %顯示內(nèi)容N=5 內(nèi)蒙古科技大學(xué)電力系統(tǒng)穩(wěn)態(tài)分析課程設(shè)計說明書19disp(39。平衡節(jié)點(diǎn)為：39。)。1Y=input(39。請輸入導(dǎo)納矩陣：Y=39。)。pr=input(39。請輸入誤差精度：pr=39。)。G=real(Y)。B=imag(Y)。 %real 實(shí)部，image 虛部B1=B(2:5, 2:5)。 %去掉平衡節(jié)點(diǎn)B2=B1。 b1=inv(B1)。 %取逆b2=inv(B2)。ure=input(39。請輸入平衡節(jié)點(diǎn) 1 電壓初值的實(shí)部，Re=39。)。dim=input(39。請輸入平衡節(jié)點(diǎn) 1 電壓初值的虛部，Im=39。)。for m=2:5dt(m)=0。 %電壓虛部初值u(m)=。 %電壓實(shí)部初值endfor m=2:5p(m)=input(39。請輸入 PQ 節(jié)點(diǎn)的有功功率： p=39。)。q(m)=input(39。請輸入 PQ 節(jié)點(diǎn)的無功功率： q=39。)。endk=0。wucha=1。 %迭代次數(shù)，while wucha(1)pr。u(1)=ure。dt(1)=dim。 %平衡節(jié)點(diǎn)電壓實(shí)部for m=2:5 %平衡節(jié)點(diǎn)電壓虛部for n=1:5pt(n)=u(m)*u(n)*(G(m,n)*cos(dt(m)dt(n))+B(m,n)*sin(dt(m)dt(n)))。 % ????njjiU10j????00sincosjjijijBG??內(nèi)蒙古科技大學(xué)電力系統(tǒng)穩(wěn)態(tài)分析課程設(shè)計說明書20enddisp(39。各節(jié)點(diǎn)有功率不平衡量為：39。) dp(m)=p(m)sum(pt) % = ??0iP?ip????njjiU10j????00sincosjjijijBG??endfor m=2:5u1=diag(u)。 %對角矩陣u2=inv(u1(2:5,2:5))。 %取逆u3=flipud(rot90(u(2:5)))。 %逆時針 90176。上下對稱u4=u3(1:51,:)。 dp1=flipud(rot90(dp(2:5)))。 %— ( )=U139。)?BP???0?dp2=dp1(1:51,:)。dtt=(b1*dp2/u4)*u2。dtt=dtt(:,1:1)。dtt=flipud(rot90(dtt))disp(39。各節(jié)點(diǎn)電壓的相位角 dt 為：39。) dt(m)=dt(m)+dtt(m1)endfor m=2:5for n=1:5 qt(n)=u(m)*u(n)*(G(m,n)*sin(dt(m)dt(n))B(m,n)*cos(dt(m)dt(n)))。end % ????njjiU10j????11cossinijijjjBG??disp(39。各節(jié)點(diǎn)無功率不平衡量為：39。)qq(m)=q(m)sum(qt) % = ??0iQ?i???njji10j????11cossinijijjj ??endfor m=2:5內(nèi)蒙古科技大學(xué)電力系統(tǒng)穩(wěn)態(tài)分析課程設(shè)計說明書21qq1=flipud(rot90(qq(2:5)))。 qq2=qq1(1:51,:)。ut=b2*qq2/u4。disp( 39。各節(jié)點(diǎn)電壓的修正量：39。) utdisp( 39。各節(jié)點(diǎn)電壓的大小為：39。) u(m)=u(m)+ut(m1) endfor n=1:5 U(n)=u(n)*(cos(dt(n))+j*sin(dt(n)))。 %電壓endfor n=1:5 I(n)=Y(1,n)*U(n)。enddisp(39。各節(jié)點(diǎn)功率 s 為：39。)S1=U(1)*sum(conj(I)) %電流for m=1:5 for n=1:5 S(m,n)=U(m)*(conj(U(m))conj(U(n)))*conj(Y(m,n))。 %功率endenddisp(39。各節(jié)點(diǎn)功率 sij 為：39。)Swucha=max(abs(ut))k=k+1 %迭代+1end內(nèi)蒙古科技大學(xué)電力系統(tǒng)穩(wěn)態(tài)分析課程設(shè)計說明書22 程序運(yùn)行結(jié)果節(jié)點(diǎn)總數(shù)為：N = 5平衡節(jié)點(diǎn)為：ans = 1請輸入導(dǎo)納矩陣：Y=[ 5+15i + 0 0 5+15i +5i +5i + + +5i 10+30i 0 0 +5i 10+30i +0 + 0 + ]請輸入誤差精度：pr=請輸入平衡節(jié)點(diǎn) 1 電壓初值的實(shí)部，Re=請輸入平衡節(jié)點(diǎn) 1 電壓初值的虛部，Im=0請輸入 PQ 節(jié)點(diǎn)的有功功率：p=請輸入 PQ 節(jié)點(diǎn)的無功功率：q=請輸入 PQ 節(jié)點(diǎn)的有功功率：p=請輸入 PQ 節(jié)點(diǎn)的無功功率：q=請輸入 PQ 節(jié)點(diǎn)的有功功率：p=請輸入 PQ 節(jié)點(diǎn)的無功功率：q=請輸入 PQ 節(jié)點(diǎn)的有功功率：p=內(nèi)蒙古科技大學(xué)電力系統(tǒng)穩(wěn)態(tài)分析課程設(shè)計說明書23請輸入 PQ 節(jié)點(diǎn)的無功功率：q=各節(jié)點(diǎn)有功率不平衡量為：dp = 0 各節(jié)點(diǎn)有功率不平衡量為：dp = 0 各節(jié)點(diǎn)有功率不平衡量為：dp = 0 各節(jié)點(diǎn)有功率不平衡量為：dp = 0 dtt = 各節(jié)點(diǎn)電壓的相位角 dt 為：dt = 0 0 0 0dtt = 各節(jié)點(diǎn)電壓的相位角 dt 為：dt = 0 0 0dtt = 各節(jié)點(diǎn)電壓的相位角 dt 為：dt = 0 0內(nèi)蒙古科技大學(xué)電力系統(tǒng)穩(wěn)態(tài)分析課程設(shè)計說明書24dtt = 各節(jié)點(diǎn)電壓的相位角 dt 為：dt = 0 各節(jié)點(diǎn)無功率不平衡量為：qq = 0 各節(jié)點(diǎn)無功率不平衡量為：qq = 0 各節(jié)點(diǎn)無功率不平衡量為：qq = 0 各節(jié)點(diǎn)無功率不平衡量為：qq = 0 各節(jié)點(diǎn)電壓的修正量：ut = 0 0 0 0 0 0 0 0 0 0 0 0各節(jié)點(diǎn)電壓的大小為：u = 各節(jié)點(diǎn)電壓的修正量：ut =內(nèi)蒙古科技大學(xué)電力系統(tǒng)穩(wěn)態(tài)分析課程設(shè)計說明