【正文】 else p=B1(s,1)。for s=1:m if B1(s,6)==0 %不含變壓器的支路 p=B1(s,1)。 Jacobi(h,r+k)=G(i,j)*f(i)+B(i,j)*e(i)。 Jacobi(h,r+k)=G(i,j)*f(i)+B(i,j)*e(i)D(h,1)。 end end for s=1:n if s~=isbamp。 for j=1:n C(s,1)=C(s,1)+G(s,j)*e(j)B(s,j)*f(j)。while abs(max(DetaU))pr disp(39。 end end end end end k=0。B2(s,3)==3 h=h+1。 OrgS(h,1)=B2(s,4)。 h=0。 else %含有變壓器的支路 p=B1(s,1)。誤差精度:pr=39。 再者我要感謝我的家人多年來對我的關(guān)懷和支持。稀疏存儲(chǔ)技術(shù)源于數(shù)據(jù)結(jié)構(gòu)中的數(shù)據(jù)存儲(chǔ)方式。 二層鏈表 二層鏈表的改進(jìn)用于牛頓拉夫遜潮流計(jì)算的二層鏈表存儲(chǔ)結(jié)構(gòu)中，節(jié)點(diǎn)導(dǎo)納矩陣與節(jié)點(diǎn)分塊雅可比矩陣中的每一個(gè)節(jié)點(diǎn)單元均對應(yīng)一個(gè)二層鏈表單元。(a)所示，單鏈表結(jié)構(gòu)中每個(gè)節(jié)點(diǎn)包含非零元素的元素值、編號(hào)等信息以及與該節(jié)點(diǎn)相連的下一個(gè)非零元素地址(next)，另外，還需要有一個(gè)指針 (head) 指向第一個(gè)非零元素存儲(chǔ)單元。其存儲(chǔ)格式是：U ——存儲(chǔ)矩陣A的上三角部分的非零元素的值，按行依次存儲(chǔ)JU——存儲(chǔ)矩陣A的上三角部分的非零元素的列號(hào)IU——存儲(chǔ)矩陣A中上三角部分每行第一個(gè)非零元素在U中的位置(首地址)L ——存儲(chǔ)矩陣A中下三角部分的非零元素的值，按列依次存儲(chǔ)IL——存儲(chǔ)矩陣A的下三角部分的非零元素的行號(hào)JL——存儲(chǔ)矩陣A中下三角部分每列第一個(gè)非零元素在L中的位置(首地址)D ——存儲(chǔ)矩陣A的對角元素的值，其檢索下標(biāo)不需要存儲(chǔ)三角檢索存儲(chǔ)格式示例：則矩陣A的存儲(chǔ)可表示為：IU（3）為4，表明矩陣A上三角部分第3行的第1個(gè)非零元如果有的話應(yīng)在U的第4個(gè)位置，而U表中第4個(gè)位置沒有非零元素，為了檢索方便，IU(3)仍應(yīng)賦值4。所以極坐標(biāo)形式的牛頓拉夫遜潮流計(jì)算的未知變量少了個(gè)，方程數(shù)也相應(yīng)少了個(gè)。除了一個(gè)平衡節(jié)點(diǎn)的電壓已知外，共有個(gè)未知數(shù)待求。上述情況表明無功功率已達(dá)到極限值，不能再用調(diào)節(jié)無功功率來保持PV節(jié)點(diǎn)的電壓值?？焖俳怦罘ㄊ怯膳ｎD拉夫遜法極坐標(biāo)形式簡化而來的，它以節(jié)點(diǎn)有功功率偏差作為電壓相角的修正量，以節(jié)點(diǎn)無功功率偏差作為電壓幅值的修正量，有功功率和無功功率迭代是分開進(jìn)行的。 功率方程電力網(wǎng)絡(luò)的節(jié)點(diǎn)電壓方程就是潮流計(jì)算時(shí)的數(shù)學(xué)模型。作者在文章指出，采用VEPD法存儲(chǔ)矩陣所占用的內(nèi)存空間較小，且元素的查找、修改、刪除等操作方便快捷。關(guān)鍵詞：潮流計(jì)算，稀疏技術(shù)，二層鏈表ABSTRCTPower flow calculation is one of most basic eletrical operation applied in the electrical system analysis,its efficiency has been a matter of concern to ,the sparse technique has been widely used in the areas of engineering calculations,scientific ,to improve the efficiency of the sparse technique by bing the power flow calculation features could enhance the efficiency of the electrical proposed a special kind of twolayer linked list to make the bus admittance matrix and the bus block jacobian matrix stored,and the function of location searching between the two matrix implemented under the same the efficiency of the foemation and modification of jacobian matrix has benn largely research then used the symmetry of the symmetric matrix during the factorization,halved the amount of the layer two linked list storage cell,and formed the lower triangular and twolayer linked the storage space could be better the mean time,through exacting the elements from the being revised jacobian matrix and the unbalanced power to reduce the operation so as to enhance the efficiency.Keywords:Power Flow Calculation,Sparse Technique,Linked ListII1 緒論 課題研究目的和意義潮流計(jì)算是電力系統(tǒng)分析中最根本的電氣運(yùn)算。一般而言，數(shù)據(jù)的有效運(yùn)算是非零元素之間的運(yùn)算，但稀疏數(shù)據(jù)中大量的零元素將會(huì)給計(jì)算帶來大量的冗余度，降低計(jì)算機(jī)算法的速度和效率。文獻(xiàn)[8]提出了一種在輻射網(wǎng)絡(luò)的導(dǎo)納矩陣消去分解過程中不出現(xiàn)注入元的節(jié)點(diǎn)優(yōu)化編號(hào)方法，即逆流編號(hào)法。功率方程隨著節(jié)點(diǎn)電壓向量表示形式不同有兩種不同的形式。但是目前為止，無論新的或是改進(jìn)的算法，仍不能取代牛頓拉夫遜法和快速解耦法的地位。求出節(jié)點(diǎn)電壓后，即可計(jì)算各條線路的潮流、平衡節(jié)點(diǎn)功率以及各支路的功率損耗。因而其收斂判據(jù)為其中為預(yù)先給定的小數(shù)。具體的操作方式有很多，如果按照存儲(chǔ)原理可以分為兩大類：一是靜態(tài)數(shù)組結(jié)構(gòu)存儲(chǔ)方法，二是動(dòng)態(tài)鏈?zhǔn)浇Y(jié)構(gòu)存儲(chǔ)方法。但是如果在計(jì)算過程中，矩陣A的稀疏結(jié)構(gòu)發(fā)生了改變，即其中的非零元素的分布位置發(fā)生變化，相應(yīng)的檢索信息也要隨之變化，這時(shí)用三角檢索存儲(chǔ)格式很不方便。對于NN 的二維稀疏數(shù)組，如果將每一行看作是1個(gè)一維稀疏數(shù)組，則該NN的二維數(shù)組可以由N 個(gè)一維稀疏數(shù)組組成。二層鏈表下層為二叉單元鏈表層，只存儲(chǔ)節(jié)點(diǎn)導(dǎo)納矩陣下三角非零元素。用于牛頓拉夫遜潮流計(jì)算的二層鏈表中，節(jié)點(diǎn)導(dǎo)納矩陣與節(jié)點(diǎn)分塊雅可比矩陣中的每一個(gè)節(jié)點(diǎn)單元均對應(yīng)了一個(gè)鏈表單元。clc。B2=input(39。 Y(p,p)=Y(p,p)+B1(s,5)^2/B1(s,3)。 end OrgS(h,1)=C(s,1)*e(s)+f(s)*D(s,1)。B2(s,3)==2 h=h+1。 %對P的處理 for i=1:n if B2(i,3)==2 h=h+1。 Jacobi(r+h,k)=2*e(h)。 end %創(chuàng)建S，用于存儲(chǔ)參數(shù) OrgS=zeros(2*n2,1)。 OrgS(h,1)=f(s)*C(s,1)D(s,1)*e(s)。B2(s,3)==2 h=h+1。 Jacobi(r+h,k)=G(i,j)*f(i)+B(i,j)*e(i)。 %電壓偏移量 U=zeros(1,n)。 q=B1(s,2)。各條支路的功率為：39。endS=zeros(n)。 end end %節(jié)點(diǎn)電壓 for j=1:n U(1,j)=B2(j,2)。 if i==j Jacobi(h,k)=G(i,j)*e(i)+B(i,j)*f(i)C(h,1)。 h=0。 r=0。 for j=1:n C(s,1)=C(s,1)+G(s,j)*e(j)B(s,j)*f(j)。 end if k==(n1) k=0。 else Jacobi(h,k)=G(i,j)*e(i)B(i,j)*f(i)。 DetaS(h,1)=real(B2(s,1))OrgS(h,1)。 end end %Q的處理 for s=1:n if s~=isbamp。 for s=1:n e(s)=real(B2(s,2))。 q=B1(s,2)。)。是仰老師的耐心指導(dǎo)使我可以克服論文寫作過程中的困難和難題。*i 1 2 0。如將傳統(tǒng)的雅可比矩陣分塊，將每一個(gè)節(jié)點(diǎn)或節(jié)點(diǎn)與節(jié)點(diǎn)之間的鏈接看成一個(gè)整體單元，即可形成以階子陣、構(gòu)成的分塊雅可比矩陣。十字鏈表實(shí)質(zhì)上是動(dòng)態(tài)鏈表，主要用于存儲(chǔ)稀疏矩陣。 按行（列）存儲(chǔ)格式按行（列）存儲(chǔ)格式的原理是按行(列)順序依次存儲(chǔ)階矩陣A中的非零元，同一行(列)元素依次排在一起。因此，牛頓拉夫遜潮流計(jì)算時(shí)的修正方程還有另一種表達(dá)形式。?直角坐標(biāo)形式的牛頓拉夫遜潮流計(jì)算、無功功率與電壓幅值方程。對于PV節(jié)點(diǎn)，因其無功功率是未知量，只能在迭代開始時(shí)給定初始值，此后的迭代值必須在逐次的迭代過程中計(jì)算得出。高斯賽德爾迭代法是數(shù)學(xué)上求解線性或非線性方程組的一種常用的迭代方法。本文的章節(jié)具體安排如下：第一章 緒論第二章 電力系統(tǒng)潮流計(jì)算第三章 稀疏技術(shù)第四章 用于牛頓拉夫遜潮流計(jì)算的稀疏技術(shù)第五章 結(jié)論注：1. 第一章寫出本文的研究目的和意義，并對目前稀疏技術(shù)在電力系統(tǒng)潮流計(jì)算中的研究現(xiàn)狀加以介紹；2. 第二章寫出了電力系統(tǒng)潮流計(jì)算的原理，介紹了兩種潮流計(jì)算的方法，其中作者詳細(xì)的介紹了牛頓拉夫遜潮流計(jì)算法；3. 第三章作者介紹了稀疏技術(shù)內(nèi)容，著重對各類稀疏存儲(chǔ)技術(shù)進(jìn)行了詳細(xì)的講解與說明；4. 第四章作者將稀疏技術(shù)運(yùn)用到電力系統(tǒng)的潮流計(jì)算中，旨在提高計(jì)算的速度和效率，并且給出Matlab編寫的算法實(shí)例；5. 第五章總結(jié)。文獻(xiàn)[2]通過研究數(shù)組存儲(chǔ)方式和鏈表存儲(chǔ)方式在牛頓拉夫遜潮流計(jì)算中的效率問題、以及內(nèi)存占用情況，證實(shí)鏈表存儲(chǔ)方式在內(nèi)存空間占用、計(jì)算效率上有明顯的優(yōu)勢。 畢 業(yè) 設(shè) 計(jì)（論 文）設(shè)計(jì)(論文)題目：基于潮流計(jì)算的稀疏技術(shù)研究 目 錄摘要 IABSTRCT II1 緒論 1 課題研究目的和意義 1 稀疏技術(shù)在電力系統(tǒng)潮流計(jì)算中的研究 1 本文完成的主要工作和章節(jié)安排 22 電力系統(tǒng)潮流計(jì)算 4 電力系統(tǒng)潮流計(jì)算的數(shù)學(xué)模型 4 節(jié)點(diǎn)電壓方程 4 功率方程 5 電力系統(tǒng)潮流計(jì)算的方法 5 高斯賽德爾潮流計(jì)算 6 牛頓拉夫遜潮流計(jì)算 83 稀疏存儲(chǔ)技術(shù) 15 靜態(tài)數(shù)組結(jié)構(gòu)存儲(chǔ)方法 15 散居格式 15 按行（列）存儲(chǔ)格式 15 三角檢索存儲(chǔ)格式 16 17 帶行指針數(shù)組的單鏈表表示法 17 十字鏈表表示法 18 二叉單元鏈表數(shù)組 194 用于牛頓拉夫遜潮流計(jì)算的稀疏技術(shù) 21 節(jié)點(diǎn)分塊雅可比矩陣 21 二層鏈