【正文】 。②基于MATLAB的編程。ji)稱互導(dǎo)納，由此可得互導(dǎo)納數(shù)值上就等于在節(jié)點(diǎn)i施加單位電壓，其他節(jié)點(diǎn)全部接地時(shí)，經(jīng)節(jié)點(diǎn)j注入網(wǎng)絡(luò)的電流，因此可定義為： (24)節(jié)點(diǎn)j，i之間的互導(dǎo)納數(shù)值上就等于連接節(jié)點(diǎn)j，i支路到導(dǎo)納的負(fù)值。輸電線路的這些參數(shù)通?？梢哉J(rèn)為是沿全長均勻分布的，每單位長度的參數(shù)為電阻r0。對(duì)現(xiàn)有電力系統(tǒng)的運(yùn)行和擴(kuò)建，對(duì)新的電力系統(tǒng)進(jìn)行規(guī)劃設(shè)計(jì)以及對(duì)電力系統(tǒng)進(jìn)行靜態(tài)和暫態(tài)穩(wěn)定分析都是以潮流計(jì)算為基礎(chǔ)。可以選擇出線數(shù)多或者靠近電網(wǎng)中心的發(fā)電廠母線作平衡節(jié)點(diǎn)。其要點(diǎn)是把非線性方程式的求解過程變成反復(fù)地對(duì)相應(yīng)的線性方程式進(jìn)行求解的過程。解決這個(gè)問題的辦法可以用高斯法迭代1~2次，以此迭代結(jié)果作為牛頓拉夫遜法的初值。但是， H、L 在迭代過程中仍然不斷變化，而且又都是不對(duì)稱矩陣，因此，牛頓法的第二個(gè)簡(jiǎn)化，就是把式(224)中的系統(tǒng)矩陣簡(jiǎn)化為在迭代過程中不變的對(duì)稱矩陣。⑤折算到哪一側(cè)的標(biāo)志一欄，表示如果系統(tǒng)簡(jiǎn)圖的支路起始節(jié)點(diǎn)處于高壓側(cè)則請(qǐng)輸入“1”，否則輸入“0”。)。I=zeros(y,1)。P=point(:,4)。 Y(J(m),I(m))=1/Zij(m)。BB=B。 enddeltaPi(m)=P(notph(m))sum1。 pqnum=pqnum(1)。 程序運(yùn)行結(jié)果輸出程序段fid=fopen(39。總結(jié)整個(gè)程序設(shè)計(jì)過程，作者得出幾下結(jié)論：面對(duì)當(dāng)今復(fù)雜的電力網(wǎng)絡(luò)，運(yùn)用計(jì)算機(jī)完成潮流計(jì)算是可行的，并且具有速度快，誤差小等優(yōu)點(diǎn)。與同組同學(xué)在課題工作中共同奮斗的時(shí)光給作者留下了難以磨滅的深刻記憶，遇到困難時(shí)的相互幫助，成功時(shí)的喜悅都值得作者終生回味。)。I=zeros(y,1)。P=point(:,4)。 Y(J(m),I(m))=1/(KT(m)*ZT(m))。 %%求B39。 deltaPi=zeros(x1,1)。 end else Uq=U。 end Uq=U。 Sji(m)=(U(J(m))*cos(a(J(m)))+i*U(J(m))*sin(a(J(m))))*((U(J(m))*cos(a(J(m)))i*U(J(m))*sin(a(J(m))))*((1KT(m))/(KT(m))^2*conj(1/ZT(m)))+(U(J(m))*cos(a(J(m)))i*U(J(m))*sin(a(J(m)))U(I(m))*cos(a(I(m)))+i*U(I(m))*sin(a(I(m))))*conj(Y(I(m),J(m))))。)。fprintf(fid,39。 ================================================\n39。,m,P(m))。for m=1:y fprintf(fid,39。 線路功率損耗deltaSij為: \n39。fprintf(fid,39。)。 節(jié)點(diǎn)有功功率P為: \n39。for m=1:x fprintf(fid,39。,ph,real(Sph),imag(Sph))。wt39。 Sji(m)=(U(J(m))*cos(a(J(m)))+i*U(J(m))*sin(a(J(m))))*((U(J(m))*cos(a(J(m)))i*U(J(m))*sin(a(J(m))))*(i*B0(m)/2)+(U(J(m))*cos(a(J(m)))i*U(J(m))*sin(a(J(m)))U(I(m))*cos(a(I(m)))+i*U(I(m))*sin(a(I(m))))*conj(Y(I(m),J(m))))。 end kq=1。 end deltaQi(m)=Q(pq(m))sum2。amp。 BB(:,ph)=[]。W=zhilu(:,9)。KT=zeros(y,1)。U=zeros(x,1)。e=xlsread(39。衷心感謝付老師無微不至的教導(dǎo)和關(guān)心。圖42 算例原始數(shù)據(jù)Excel輸入界面 算例運(yùn)行結(jié)果在程序運(yùn)行后新生成的TXT文檔中輸出如下結(jié)果：*************電力系統(tǒng)06白云霄 0603010213**************************潮流上機(jī)輸出結(jié)果*************迭代次數(shù)k為: 12 ================================================平衡節(jié)點(diǎn)1的復(fù)功率Sph為: +j*() ================================================節(jié)點(diǎn)電壓U為: 第1個(gè)節(jié)點(diǎn)電壓: 第2個(gè)節(jié)點(diǎn)電壓: 第3個(gè)節(jié)點(diǎn)電壓: 第4個(gè)節(jié)點(diǎn)電壓: 第5個(gè)節(jié)點(diǎn)電壓: ================================================節(jié)點(diǎn)相角a為: 第1個(gè)節(jié)點(diǎn)相角: 第2個(gè)節(jié)點(diǎn)相角: 第3個(gè)節(jié)點(diǎn)相角: 第4個(gè)節(jié)點(diǎn)相角: 第5個(gè)節(jié)點(diǎn)相角: ================================================節(jié)點(diǎn)復(fù)功率S為: 第1個(gè)節(jié)點(diǎn)復(fù)功率: +i*()第2個(gè)節(jié)點(diǎn)復(fù)功率: +i*()第3個(gè)節(jié)點(diǎn)復(fù)功率: +i*()第4個(gè)節(jié)點(diǎn)復(fù)功率: +i*()第5個(gè)節(jié)點(diǎn)復(fù)功率: +i*() ================================================節(jié)點(diǎn)有功功率P為: 第1個(gè)節(jié)點(diǎn)有功功率P: 第2個(gè)節(jié)點(diǎn)有功功率P: 第3個(gè)節(jié)點(diǎn)有功功率P: 第4個(gè)節(jié)點(diǎn)有功功率P: 第5個(gè)節(jié)點(diǎn)有功功率P: ================================================節(jié)點(diǎn)無功功率Q為: 第1個(gè)節(jié)點(diǎn)無功功率Q: 第2個(gè)節(jié)點(diǎn)無功功率Q: 第3個(gè)節(jié)點(diǎn)無功功率Q: 第4個(gè)節(jié)點(diǎn)無功功率Q: 第5個(gè)節(jié)點(diǎn)無功功率Q: ================================================線路功率Sij和Sji為: 節(jié)點(diǎn)1到節(jié)點(diǎn)2的功率為: +i*() 節(jié)點(diǎn)2到節(jié)點(diǎn)1的功率為: +i*()節(jié)點(diǎn)2到節(jié)點(diǎn)3的功率為: +i*() 節(jié)點(diǎn)3到節(jié)點(diǎn)2的功率為: +i*() 節(jié)點(diǎn)2到節(jié)點(diǎn)4的功率為: +i*() 節(jié)點(diǎn)4到節(jié)點(diǎn)2的功率為: +i*() 節(jié)點(diǎn)3到節(jié)點(diǎn)4的功率為: +i*() 節(jié)點(diǎn)4到節(jié)點(diǎn)3的功率為: +i*() 節(jié)點(diǎn)3到節(jié)點(diǎn)5的功率為: +i*() 節(jié)點(diǎn)5到節(jié)點(diǎn)3的功率為: +i*() ================================================網(wǎng)絡(luò)總損耗sumdeltaS為: +i*() ================================================線路功率損耗deltaSij為: 12線路的功率損耗為: +i*() 23線路的功率損耗為: +i*() 24線路的功率損耗為: +i*() 34線路的功率損耗為: +i*() 35線路的功率損耗為: +i*()千萬不要?jiǎng)h除行尾的分節(jié)符，此行不會(huì)被打印。Upq=Uq。②設(shè)置pqnum1階空行矩陣deltaQi。 計(jì)算不平衡功率ΔPi程序段notph=find(TYPE(:,1)3)。 ④在B矩陣中找出平衡節(jié)點(diǎn)和PV節(jié)點(diǎn)，并將其所在的行和列置零。W=zhilu(:,9)。 KT=zeros(y,1)。 U=zeros(x,1)。e=xlsread(39。圖32 Excel輸入界面在此，對(duì)數(shù)據(jù)的輸入有以下幾點(diǎn)說明：①在節(jié)點(diǎn)信息里，節(jié)點(diǎn)電壓及相角為迭代計(jì)算時(shí)所設(shè)的初值。 PQ分解法潮流計(jì)算PQ分解法是從簡(jiǎn)化牛頓法極坐標(biāo)形式計(jì)算潮流程序的基礎(chǔ)上提出來的。牛頓拉夫遜法的可靠收斂取決于有一個(gè)良好的啟動(dòng)初值。 牛頓拉夫遜法求解潮流計(jì)算 牛頓拉夫遜法原理潮流計(jì)算的結(jié)果可用如電力系統(tǒng)穩(wěn)態(tài)研究，安全估計(jì)或最優(yōu)潮流等對(duì)潮流計(jì)算的模型和方法有直接影響的分析。③平衡節(jié)點(diǎn)在潮流計(jì)算中，這類節(jié)點(diǎn)一般只設(shè)一個(gè)。圖210 阻抗歸算到高壓側(cè)變壓器等值電路 電力系統(tǒng)潮流計(jì)算 電力系統(tǒng)潮流計(jì)算數(shù)學(xué)模型電力系統(tǒng)潮流計(jì)算是電力系統(tǒng)分析中的一種最基本的計(jì)算，是對(duì)復(fù)雜電力系統(tǒng)正常和故障條件下穩(wěn)態(tài)運(yùn)行狀態(tài)的計(jì)算。因此，在一般情況下，節(jié)點(diǎn)導(dǎo)納矩陣的對(duì)角元素往往大于非對(duì)角元素的負(fù)值。 導(dǎo)納矩陣的形成節(jié)點(diǎn)導(dǎo)納矩陣的對(duì)角元素Yii (i=1，2，n)成為自導(dǎo)納。①探討學(xué)習(xí)潮流計(jì)算的基本原理。為了克服阻抗法在內(nèi)存和速度方面的缺點(diǎn)，60年代中期發(fā)展了以阻抗矩陣為基礎(chǔ)的分塊阻抗法。潮流計(jì)算的目標(biāo)是求取電力系統(tǒng)在給定運(yùn)行狀態(tài)的計(jì)算。 減少內(nèi)存占用量，用于牛頓潮流計(jì)算中，不僅減少了內(nèi)存占用量，還大大提高了潮流計(jì)算的速度。為實(shí)現(xiàn)這一功能，電力系統(tǒng)在各個(gè)環(huán)節(jié)和不同層次還具有相應(yīng)的信息與控制系統(tǒng)，對(duì)電能的生產(chǎn)過程進(jìn)行測(cè)量、調(diào)節(jié)、控制、保護(hù)、通信和調(diào)度，以保證用戶獲得安全、經(jīng)濟(jì)、優(yōu)質(zhì)的電能如圖11所示。在電力系統(tǒng)規(guī)劃的設(shè)計(jì)和現(xiàn)有電力系統(tǒng)運(yùn)行方式的研究中，都需要利用潮流計(jì)算來定量地分析比較供電方案或運(yùn)行方式的合理性，可靠性和經(jīng)濟(jì)性。 計(jì)算的方便性和靈活性完善的潮流程序應(yīng)該可以重復(fù)應(yīng)用，可移植性高。實(shí)際電力系統(tǒng)的潮流技術(shù)主要采用牛頓拉夫遜法及其衍生的PQ分解法。自從60年代中期，在牛頓法中利用了最佳順序消去法以后，牛頓法在收斂性。本文對(duì)MATLAB軟件編程中所用的M語言有一定的介紹，說明了常用函數(shù)的作用，并對(duì)M語言中所涉及到的邏輯關(guān)系符，運(yùn)算符，矩陣的正確輸入以及簡(jiǎn)單的人機(jī)對(duì)話功能有了初步的解釋。顯然，恒等于。電感L0，電導(dǎo)g0及電容C0，其一項(xiàng)等值電路如圖23所示。采用導(dǎo)納矩陣時(shí)，節(jié)點(diǎn)注入電流和節(jié)點(diǎn)電壓構(gòu)成如式(21)所示線性方程組可展開如下形式： (210)由于實(shí)際電網(wǎng)中測(cè)量的節(jié)點(diǎn)注入量一般不是電流而是功率，因此，必須將式中的注入電流用節(jié)點(diǎn)注入功率來表示。以上三類節(jié)點(diǎn)4個(gè)運(yùn)行參數(shù)P，Q，V，中，已知量都是兩個(gè)，待求量也是兩個(gè)，只是類型不同而已。即通常所稱的逐次線性化過程。也可以先用直流法潮流求解一次以求得一個(gè)較好的角度初值，然后轉(zhuǎn)入牛頓拉夫遜法迭代。一般情況下，線路兩端電壓相角差是不大的（通常不超過10176。⑥本條程序默認(rèn)節(jié)點(diǎn)數(shù)，支路數(shù)均在在100以下，可以解決絕大部分電力系統(tǒng)的潮流問題，若遇到超大系統(tǒng)，可對(duì)程序做稍加調(diào)整，仍然使用。[point]=xlsread(39。 J=zeros(y,1)。Q=point(:,5)。 else Y(I(m),J(m))=1/(KT(m)*ZT(m))。BB(:,ph)=[]。end 計(jì)算修正相角Δθi①將之前求得的求逆，形成新的矩陣。deltaQi=zeros(pqnum,1)。39。本設(shè)計(jì)采用的PQ分解法利用了電力系統(tǒng)的一些運(yùn)行特性，對(duì)牛頓拉夫遜法做了簡(jiǎn)化，收斂性較好，可提高計(jì)算速度，物理概念較為清晰。在此對(duì)他們表示衷心的感謝!感謝父母對(duì)我的教誨與鼓勵(lì)，沒有他們的支持和犧牲，論文的完成是難以想象的。[point]=xlsread(39。J=zeros(y,1)。Q=point(:,5)。 endendfor m=1:x for n=1:y if KT(n)==0 if(I(n)==m|J(n)==m) Y(m,m)=Y(m,m)Y(I(n),J(n))+j*B0(n)/2。39。 pq=find(TYPE(:,1)==1)。 Uq(phpv)=[]。 Uq(phpv)=[]。 end end enddeltaSij=Sij+Sji。fprintf(fid,39。 節(jié)點(diǎn)電壓U為: \n39。)。endfprintf(fid,39。 節(jié)點(diǎn)%d到節(jié)點(diǎn)%d的功率為: %f+i*(%f) \n39。)。,real(sumdeltaS),imag(sumdeltaS))。 =======