【正文】 ⑵對(duì)第i列的元素進(jìn)行消去運(yùn)算：（p=k+1，k+2,…n）然后按照步驟⑴、 ⑵對(duì)其進(jìn)行矩陣的消去運(yùn)算，直到對(duì)所有的列進(jìn)行消去運(yùn)算為止。 直接法常用于解線性方程組，隨著計(jì)算技術(shù)的發(fā)展，特別是應(yīng)用電力網(wǎng)絡(luò)的線性方程組的系數(shù)矩陣的稀疏性和對(duì)稱性處理技巧后，直接法解線性方程組更顯得比迭代法優(yōu)越。一般的發(fā)電機(jī)模型為PQ、PV和平衡節(jié)點(diǎn)。 與牛頓拉夫遜相比，PQ分解法的具有如下的優(yōu)點(diǎn)： ①系數(shù)矩陣的階數(shù)很大的降低，降低了儲(chǔ)存容量的要求，加快了計(jì)算的速度； ②系數(shù)矩陣是保持不變的，而牛頓拉夫遜法的系數(shù)矩陣在每一次迭代的過(guò)程中都會(huì)發(fā)生改變，從而大大的加快了計(jì)算速度，使潮流計(jì)算更加的快速； ③分解法中的系數(shù)矩陣是對(duì)稱的，在求逆等運(yùn)算量和所需的儲(chǔ)存量都大為的減少。②將設(shè)定的初值代入到式（23）、（24）中，求取修正方程式中的功率不平衡量和電壓不平衡量；③利用電壓的初值代入式（25）中，求取雅可比矩陣中的各個(gè)元素；④解修正方程式，求得電壓的變化量 ，計(jì)算出新的各節(jié)點(diǎn)電壓值；⑤比較求取的不平衡量是否滿足題目所要求的精度，不滿足則從第二部開(kāi)始進(jìn)行新的一次迭代，滿足則計(jì)算平衡節(jié)點(diǎn)的功率和線路功率。前推回代法的計(jì)算前提條件是已知配電網(wǎng)的開(kāi)始端電壓和線路的末端負(fù)荷其計(jì)算單位是以饋線為基本。⑵熟練的掌握了有關(guān)PQ算法的原理，熟悉了該算法的計(jì)算步驟，了解了其在運(yùn)算過(guò)程中的優(yōu)缺點(diǎn)。 文獻(xiàn)【13】主要介紹了針對(duì)配網(wǎng)的潮流計(jì)算結(jié)果的不收斂提出了母線法和支路法來(lái)解決該問(wèn)題，并取得了較好的結(jié)果。文獻(xiàn)【6】介紹了因子表的形成的原理以及在解方程時(shí)的應(yīng)用過(guò)程。為了能夠解決實(shí)際的潮流計(jì)算問(wèn)題，Dariush提出了一種初等的配電系統(tǒng)的三相潮流計(jì)算的實(shí)時(shí)解析性此方法對(duì)弱環(huán)網(wǎng)潮流計(jì)算進(jìn)行了擴(kuò)展，加強(qiáng)了對(duì)分布式電源的控制。 本次設(shè)計(jì)主要用的潮流算法是快速分解法，即PQ算法。⑷可靠性很高?，F(xiàn)如今，分布式電源已經(jīng)從原來(lái)的柴油發(fā)電機(jī)組的緊急備用電源或者是燃煤自備的小熱電廠發(fā)展成為了燃料電池、微型的燃?xì)廨啓C(jī)或者是包括太陽(yáng)能、風(fēng)力、生物質(zhì)能在內(nèi)的可再生能源發(fā)電系統(tǒng)?？偟膩?lái)說(shuō)BX是能夠改進(jìn)PQ算法的，但是還存在不足，需要對(duì)這個(gè)算法進(jìn)行改善，一定會(huì)改善PQ算法的，使該算法能夠完全的適用于所有的配電網(wǎng)。分布式電源在并入電網(wǎng)以后，會(huì)對(duì)配電網(wǎng)的電壓、網(wǎng)絡(luò)損耗等產(chǎn)生一定的影響，所以對(duì)分布式電源并網(wǎng)后的潮流計(jì)算研究顯得尤為重要。然而，每個(gè)國(guó)家的政策不同，對(duì)其的理解也就不同，所以至今為止也沒(méi)有一個(gè)統(tǒng)一的定義。⑵有利于環(huán)境。⑸由于分布式電源的規(guī)模較小，因此在投資方面也遠(yuǎn)遠(yuǎn)小于傳統(tǒng)大電廠的投資費(fèi)用，風(fēng)險(xiǎn)也就小很多了，而且它的建設(shè)周期也很短相對(duì)于傳統(tǒng)的大電廠而言，因此能夠比較快的解決某些地區(qū)的電力短缺的問(wèn)題，給其供電。本次設(shè)計(jì)的任務(wù)就是對(duì)其計(jì)算結(jié)果不收斂的問(wèn)題進(jìn)行分析，并加以改善。文獻(xiàn)【2】主要是對(duì)PQ分解法在配電網(wǎng)潮流計(jì)算中結(jié)果不收斂的的影響因素進(jìn)行了介紹，并對(duì)不收斂的問(wèn)題提出了解決方案。 文獻(xiàn)【9】主要介紹了分布式電源模型的建立，其節(jié)點(diǎn)類型有PV、PQ、PI、PQ(V)并對(duì)這幾種節(jié)點(diǎn)類型的算法進(jìn)行了詳細(xì)的介紹。這就是所謂的分解算法對(duì)比值的敏感性【2】。⑸最后總結(jié)本學(xué)期的畢業(yè)設(shè)計(jì)并撰寫畢業(yè)設(shè)計(jì)論文。 (21) (22) 以上就是前推回代法的全部過(guò)程，可以根據(jù)公式（21）、（22）來(lái)進(jìn)行潮流計(jì)算、并且通過(guò)循環(huán)此過(guò)程最終達(dá)到所需的精度，求得線路的參數(shù)。PQ分解法的基礎(chǔ)依然是牛頓拉夫遜法，與牛頓拉夫遜法相比，其最大的優(yōu)點(diǎn)是將純數(shù)學(xué)的牛頓拉夫遜法與電力系統(tǒng)網(wǎng)絡(luò)的具體特點(diǎn)相結(jié)合起來(lái)，對(duì)牛頓拉夫遜法進(jìn)行了了合理的簡(jiǎn)化和改進(jìn)，大大的提高了計(jì)算速度并且節(jié)省了內(nèi)存。所以，總的來(lái)收PQ分解法的計(jì)算速度要比牛頓拉夫遜法快得多，所以在潮流計(jì)算中只要能夠成功的解決它的收斂性問(wèn)題，那么利用該方法來(lái)進(jìn)行潮流計(jì)算將會(huì)使潮流計(jì)算更加的快速。在考慮到分布式電源的運(yùn)行在額定情況下的時(shí)候可以認(rèn)為其為PQ節(jié)點(diǎn)【4】。因子表適用于在求解方程組的時(shí)候，當(dāng)每次的運(yùn)算過(guò)程僅改變常數(shù)項(xiàng)矩陣，但其系數(shù)矩陣不改變，適用因子表法運(yùn)算的速度則會(huì)比較快。當(dāng)和一個(gè)算法相關(guān)的方程組需要多次迭代，而它的系數(shù)矩陣又不發(fā)生改變的時(shí)候，就利用因子表來(lái)求解該方程組，也能夠改善迭代所導(dǎo)致的結(jié)果不收斂的問(wèn)題【5】。則此公式就可以簡(jiǎn)化成的形式，其中皆為矩陣。 在向量范數(shù)中，描述了向量與實(shí)數(shù)之間的關(guān)系。下來(lái)的具體問(wèn)題就是如何解決PQ算法對(duì)比值的敏感性的問(wèn)題。其中分布式電源加在節(jié)點(diǎn)2上面。因此可以得出在形成系數(shù)矩陣的時(shí)候的時(shí)候，忽略支路的電阻對(duì)于PQ算法有較好的改善。表413 十四節(jié)點(diǎn)網(wǎng)絡(luò)的運(yùn)行結(jié)果算法類型是否收斂迭代次數(shù)時(shí)間(s)一般PQ法否含BX法否BX法+二范數(shù)否 從表413可以看出，對(duì)于十四節(jié)點(diǎn)的配電網(wǎng)絡(luò)，在加入了BX法和二范數(shù)的情況下結(jié)果均不收斂，說(shuō)明改進(jìn)之后的PQ算法不適合14節(jié)點(diǎn)的配電網(wǎng)絡(luò)。本次設(shè)計(jì)用了BX法對(duì)PQ算法進(jìn)行了改進(jìn)，又利用二范數(shù)來(lái)處理數(shù)據(jù)收斂條件的，加快了計(jì)算速度，減少了迭代次數(shù)。,39。%nl為支路數(shù)[isb]=xlsread(39。B4:B439。,39。%終止節(jié)點(diǎn)編號(hào)(tbus)x(:,3)=xlsread(39。G3:G100039。39。)。39。)。,39。%PV節(jié)點(diǎn)電壓給定值y(:,7)=xlsread(39。,39。%對(duì)地阻抗值%X=[z(:,1),z(:,2)]。ZB=UB^2/SB。f=zeros(1,n)。else p=B1(i,2)。%對(duì)變壓器支路求自導(dǎo)納 Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2。B=imag(YI)。%設(shè)V(i)為PV節(jié)點(diǎn)電壓的給定值幅值endfor i=1:nif B2(i,5)~=3 %當(dāng)節(jié)點(diǎn)為PQ節(jié)點(diǎn)時(shí)V(i)=sqrt(e(i)^2+f(i)^2)。q=0。else k=i+1。K=1。%進(jìn)行下一次迭代 for i=1:n if i~=isb %非平衡節(jié)點(diǎn)行 C(i)=0。 %本次迭代有功不滿足的節(jié)點(diǎn)個(gè)數(shù) end%=============== 用因子表求解修正方程的修正量角度sida =============== Np(K)=ICT2。 DP(I)=DP(I)B(I,L)*DP(L)。 end DQ1(i)=Q(i)V(i)*C(i)。 if ICT3~=0 L=0。 %前代消去 end end end %=================== 對(duì)PQ節(jié)點(diǎn)電壓模進(jìn)行修正 ======================= L=0。t2=clock。)。=====================================