【文章內(nèi)容簡介】 — ； ? )0(3? =— ； ? )0(4? =— ； ? )0(5? =— ； 所以 )1(2? =— ； )1(3? =— ； )1(4? =— ； )1(5? =— ； 內(nèi)蒙古科技大學(xué)電力系統(tǒng)穩(wěn)態(tài)分析課程設(shè)計說明 書 15 iQ? 按下式計算各節(jié)點(diǎn)無功功率不平衡量 ??0iQ? = iQ ? ????njj iU10 ??0jU ? ? ? ?? ?11 c o ss in ijijijij BG ?? ? ??02Q? =— [— （ — 0） — （ — 0） ]— [+]— [（ +）— （ + ） ] — [ （ —+0,091810） — （ 0,036952+） ]— [（ — +） — （ +） ]= 相似地可得 ??03Q? =— ； ??04Q? =— ； ??05Q? =— ； iU 由方程式 — 139。39。)( ?B ( UQ? )= ? U 可得 ??02U? =； ??03U? =； ??04U? =； ??05U? =； 所以 ??12U =； ??13U =； ??14U =； ??15U =； 求得各節(jié)點(diǎn)電壓的新值后，就開始第二次迭代。每次迭代所得適于表 1~表 3。由表 1 可見經(jīng) 6 次迭代就可滿足 ε ≤ 510? 的要求。 表 1 迭代過程中各節(jié)點(diǎn)功率的不平衡量 K )(2kP? ??kQ2? )(3kP? ??kQ3? )(4kP? ??kQ4? )(5kP? ??kQ5? 0 1 2 內(nèi)蒙古科技大學(xué)電力系統(tǒng)穩(wěn)態(tài)分析課程設(shè)計說明 書 16 3 4 0,000024 5 6 表 2 迭代過程中各節(jié)點(diǎn)電壓的修正量 k ?)(2k? ??kU2? ? )(3k? ??kU3? ? )(4k? ??kU4? ? )(5k? ??kU5? 0 1 2 3 4 5 6 表 3 迭代過程中各節(jié)點(diǎn)電壓 k )(2k? ??kU2 )(3k? ??kU3 )(4k? ??kU4 )(5k? ??kU5 0 1 0,108571 2 0,089348 0,103974 3 4 5 6 0,083906 2U = ?? ； 3U = ?? 4U = ?? ； 5U = ?? 內(nèi)蒙古科技大學(xué)電力系統(tǒng)穩(wěn)態(tài)分析課程設(shè)計說明 書 17 ~1S 和線路功率 ~ijS 迭代收斂后，就可以計算平衡節(jié)點(diǎn)功率和線路功率。結(jié)果如下 ~1S =+ 各線路功率如表 4 表 4 各線路功率 ~ijS j i 1 2 3 4 5 1 + + 2 + + + 3 4 + + 5 計算完畢。 第三章 程序設(shè)計 設(shè)計流程 內(nèi)蒙古科技大學(xué)電力系統(tǒng)穩(wěn)態(tài)分析課程設(shè)計說明 書 18 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 disp(39。平衡節(jié)點(diǎn)為： 39。)。 1 Y=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。 內(nèi)蒙古科技大學(xué)電力系統(tǒng)穩(wěn)態(tài)分析課程設(shè)計說明 書 19 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:5 dt(m)=0。 %電壓虛部初值 u(m)=。 %電壓實(shí)部初值 end for m=2:5 p(m)=input(39。請輸入 PQ 節(jié)點(diǎn)的有功功率： p=39。)。 q(m)=input(39。請輸入 PQ 節(jié)點(diǎn)的無功功率： q=39。)。 end k=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:5 pt(n)=u(m)*u(n)*(G(m,n)*cos(dt(m)dt(n))+B(m,n)*sin(dt(m)dt(n)))。 % ? ????njj iU10 ??0jU? ? ? ?? ?00 s inc o s ijijijij BG ?? ? end disp(39。各節(jié)點(diǎn)有功率不平衡量為： 39。) dp(m)=p(m)sum(pt) % ??0iP? = ip ? ????njj iU10 ??0jU ? ? ? ?? ?00 s inc o s ijijijij BG ?? ? end for m=2:5 u1=diag(u)。 %對角矩陣 內(nèi)蒙古科技大學(xué)電力系統(tǒng)穩(wěn)態(tài)分析課程設(shè)計說明 書 20 u2=inv(u1(2:5,2:5))。 %取逆 u3=flipud(rot90(u(2:5)))。 %逆時針 90176。上下對稱 u4=u3(1:51,:)。 dp1=flipud(rot90(dp(2:5)))。 %— 139。)( ?B ( UP? )=U ??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) end for m=2:5 for 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 % ? ????njj iU10 ??0jU ? ? ? ?? ?11 c o ss in ijijijij BG ?? ? disp(39。各節(jié)點(diǎn)無功率不平衡量為： 39。) qq(m)=q(m)sum(qt) % ??0iQ? = iQ ? ????njj iU10 ??0jU ? ? ? ?? ?11 c o ss in ijijijij BG ?? ? end for m=2:5 qq1=flipud(rot90(qq(2:5)))。 qq2=qq1(1:51,:)。 ut=b2*qq2/u4。 disp( 39。各節(jié)點(diǎn)電壓的修正量： 39。) ut disp( 39。各節(jié)點(diǎn)電壓的大小為： 39。) u(m)=u(m)+ut(m1) 內(nèi)蒙古科技大學(xué)電力系統(tǒng)穩(wěn)態(tài)分析課程設(shè)計說明 書 21 end for n=1:5 U(n)=u(n)*(cos(dt(n))+j*sin(dt(n)))。 %電壓 end for n=1:5 I(n)=Y(1,n)*U(n)。 end disp(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))。 %功率 end end disp(39。各節(jié)點(diǎn)功率 sij 為： 39。) S wucha=max(abs(ut)) k=k+1 %迭代 +1 end 程序運(yùn)行結(jié)果 內(nèi)蒙古科技大學(xué)電力系統(tǒng)穩(wěn)態(tài)分析課程設(shè)計說明 書 22 節(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= 請輸入 PQ 節(jié)點(diǎn)的無功功率： q= 各節(jié)點(diǎn)有功率不平衡量為： dp = 0 各節(jié)點(diǎn)有功率不平衡量為： dp = 內(nèi)蒙古科技大學(xué)電力系統(tǒng)穩(wěn)態(tài)分析課程設(shè)計說明 書 23 0 各節(jié)點(diǎn)有功率不平衡量為： dp = 0 各節(jié)點(diǎn)有功率不平衡量為： dp = 0 dtt = 各節(jié)點(diǎn)電壓的相位角 dt 為： dt = 0 0 0 0 dtt = 各節(jié)點(diǎn)電壓的相位角 dt 為： dt = 0 0 0 dtt = 各節(jié)點(diǎn)電壓的相位角 dt 為： dt = 0 0 dtt = 各節(jié)點(diǎn)電壓的相位角 dt 為： dt = 0 各節(jié)點(diǎn)無功率不平衡量為： 內(nèi)蒙古科技大學(xué)電力系統(tǒng)穩(wěn)態(tài)分析課程設(shè)計說明 書 24 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 = 0 0 0 0 0 0 0 0 0