直角坐標(biāo)系牛頓拉夫遜算法--畢業(yè)設(shè)計(jì)53頁

1、洛陽理工學(xué)院畢業(yè)設(shè)計(jì)（論文）基于直角坐標(biāo)系的牛頓拉夫法潮流計(jì)算摘 要潮流計(jì)算是電力網(wǎng)絡(luò)設(shè)計(jì)及運(yùn)行中最基本的計(jì)算，對電力網(wǎng)絡(luò)的各種設(shè)計(jì)方案及各種運(yùn)行方式進(jìn)行潮流計(jì)算，可以得到各種電網(wǎng)各節(jié)點(diǎn)的電壓，并求得網(wǎng)絡(luò)的潮流及網(wǎng)絡(luò)中各元件的電力損耗，進(jìn)而求得電能損耗。因此，網(wǎng)絡(luò)通過潮流計(jì)算可以分析該網(wǎng)絡(luò)的電壓水平高低，功率分布與電力損耗的合理性及經(jīng)濟(jì)性等，從而對該網(wǎng)絡(luò)的設(shè)計(jì)及運(yùn)行做出評價(jià)。本文主要章節(jié)內(nèi)容包括節(jié)點(diǎn)導(dǎo)納矩陣，電力系統(tǒng)潮流計(jì)算數(shù)學(xué)模型，直角坐標(biāo)的牛頓拉夫遜法，以及程序設(shè)計(jì)和算例結(jié)果分析。牛頓拉夫遜法潮流計(jì)算主要采用C語言編寫潮流計(jì)算程序，首先編寫原始數(shù)據(jù)的輸入程序，將原始數(shù)據(jù)形成一個(gè)目標(biāo)文件，打

2、開文件取出數(shù)據(jù)，根據(jù)電力網(wǎng)絡(luò)連接圖列寫節(jié)點(diǎn)電壓方程，編寫計(jì)算節(jié)點(diǎn)導(dǎo)納矩陣各個(gè)元素的程序，計(jì)算各個(gè)節(jié)點(diǎn)的不平衡量，進(jìn)行校驗(yàn)，如不收斂則編寫計(jì)算雅克比矩陣元素的程序，列寫修正方程，利用高斯消去法，對增廣矩陣進(jìn)行消去，從而解出修正量的值。經(jīng)過多次迭代，如果經(jīng)校驗(yàn)節(jié)點(diǎn)電壓的不平衡量已收斂，則結(jié)束迭代，轉(zhuǎn)入編寫計(jì)算平衡節(jié)點(diǎn)的功率和線路的潮流分布程序。最后，編寫潮流計(jì)算的輸出程序，輸出節(jié)點(diǎn)電壓的近似值。關(guān)鍵詞：潮流計(jì)算，牛頓拉夫遜法，C語言，節(jié)點(diǎn)導(dǎo)納矩陣，雅克比矩陣 Newton Laphson flow calculation based on Cartesian coordinate systemAB

3、STRACTFlow calculation is the most basic calculation in the design of the electricity network and the various ways of the operation. Flow calculation is used to calculation of all kinds of designs and the various ways of operation. So we can acquire a variety of nodal voltage and achieve the power l

4、oss in the network and the trend and various components in the network. Then can seek power loss. Therefore, the network can be calculated through the analysis of the level of voltage, power distribution and power losses and the reasonableness, the economy in the network.Node admittance matrix，Digit

5、al Analog Function Table ，Newton Raphson of Rectangular coordinates , formula design and results analysis are the main chapter of this article Newton- Raphson flow calculation method is main program through C-language. First of all, write input the program of the raw data and make raw data form a ta

6、rget file. Open the file and retrieve data, according to the electricity network plan to write out he raw data to form a target file, according to the power network, plan to write out nodal voltage equation, computing various elements of nodal Admittance matrix, calculated the imbalance of each node

7、, to check, if the imbalance isnt convergence, then calculate various elements of Jacques matrix and write that equation, use Gaussian elimination method to eliminate broaden matrix, answer the equation. After several rounds iteration, if the imbalance of nodal voltage is convergence, then end the i

8、teration, turn into calculate the power of balance nodes and the distribution of power in various branches. Finally, write the output program, output approximation of nodal voltage. KEY WORDS: Flow calculation；Newton - Laphson method；C language；Node admittance matrix；Jacques matrix4目錄前言1第一章 電力系統(tǒng)潮流計(jì)算

9、的基本知識21.1潮流計(jì)算的定義及應(yīng)用21.2 潮流計(jì)算方法的發(fā)展與前景3第二章 電力網(wǎng)絡(luò)的數(shù)學(xué)模型52.1節(jié)點(diǎn)導(dǎo)納矩陣的形成及修改52.1.1節(jié)點(diǎn)導(dǎo)納矩陣的形成52.1.2節(jié)點(diǎn)導(dǎo)納矩陣的修改72.2節(jié)點(diǎn)導(dǎo)納矩陣元素的物理意義8第三章 牛頓拉夫遜法潮流分布計(jì)算113.1牛頓拉夫遜法簡介113.1.1牛頓拉夫遜法簡介113.1.2牛頓拉夫遜法的幾何意義153.2牛頓拉夫遜法計(jì)算潮流分布16第四章 C程序設(shè)計(jì)的基本知識214.1 C語言的特點(diǎn)214.2 C語言的發(fā)展過程22第五章 計(jì)算實(shí)例235.1算例235.2程序計(jì)算結(jié)果245.3平衡節(jié)點(diǎn)功率和線路功率265.4數(shù)據(jù)分析26結(jié)論28謝 辭29參

10、考文獻(xiàn)30附錄31外文資料翻譯45前言如今，電能已成為我國國民經(jīng)濟(jì)緊密聯(lián)系的，不能儲(chǔ)存的能量，它的生產(chǎn)、輸送、分配和使用同時(shí)完成，暫態(tài)過程非常迅速，以電磁波的形式傳播。電力系統(tǒng)對電能生產(chǎn)、傳輸和管理實(shí)現(xiàn)自動(dòng)控制、自動(dòng)調(diào)度和自動(dòng)化管理，主要目標(biāo)是保證供電的電能質(zhì)量、頻率、電壓，保證系統(tǒng)運(yùn)行的安全可靠，提高經(jīng)濟(jì)效益和管理效能。電力系統(tǒng)潮流計(jì)算是研究電力系統(tǒng)穩(wěn)態(tài)運(yùn)行情況的一種計(jì)算，是電力系統(tǒng)安全運(yùn)行的基礎(chǔ)。它是根據(jù)給定的運(yùn)行條件和網(wǎng)絡(luò)結(jié)構(gòu)確定整個(gè)系統(tǒng)的運(yùn)行狀態(tài)，如各母線上的電壓、各原件中流過的功率、系統(tǒng)的功率損耗等等。電力系統(tǒng)潮流計(jì)算是計(jì)算系統(tǒng)動(dòng)態(tài)穩(wěn)定和靜態(tài)穩(wěn)定的基礎(chǔ)。在電力系統(tǒng)規(guī)劃設(shè)計(jì)和現(xiàn)有電力系

11、統(tǒng)運(yùn)行方式的研究中，都需要運(yùn)用電力系統(tǒng)潮流計(jì)算來定量的比較供電方案或運(yùn)行方式的穩(wěn)定性、可靠性和經(jīng)濟(jì)性。第一章 電力系統(tǒng)潮流計(jì)算的基本知識電力系統(tǒng)潮流計(jì)算是研究電力系統(tǒng)穩(wěn)態(tài)運(yùn)行情況的一種計(jì)算，它根據(jù)給定的運(yùn)行條件及系統(tǒng)接線情況確定整個(gè)電力系統(tǒng)各部分的運(yùn)行狀態(tài)：各母線的電壓，各元件中流過的功率，系統(tǒng)的功率損耗等等。在電力系統(tǒng)規(guī)劃設(shè)計(jì)和現(xiàn)有電力系統(tǒng)運(yùn)行方式的研究中，都需要利用潮流計(jì)算來定量地分析比較供電方案或運(yùn)行方式的合理性.可靠性和經(jīng)濟(jì)性。此外，電力系統(tǒng)潮流計(jì)算也是計(jì)算系統(tǒng)動(dòng)態(tài)穩(wěn)定和靜態(tài)穩(wěn)定的基礎(chǔ)。電力系統(tǒng)潮流計(jì)算也分為離線計(jì)算和在線計(jì)算兩種，前者主要用于系統(tǒng)規(guī)劃設(shè)計(jì)和安排系統(tǒng)的運(yùn)行方式，后者則用

12、于正在運(yùn)行系統(tǒng)的經(jīng)常監(jiān)視及實(shí)時(shí)控制【2】。1.1潮流計(jì)算的定義及應(yīng)用潮流計(jì)算是研究電力系統(tǒng)穩(wěn)態(tài)運(yùn)行情況的一種基本電氣計(jì)算，常規(guī)潮流計(jì)算的任務(wù)是根據(jù)給定的運(yùn)行條件和網(wǎng)路結(jié)構(gòu)確定整個(gè)系統(tǒng)的運(yùn)行狀態(tài)，如各母線上的電壓（幅值及相角）、網(wǎng)絡(luò)中的功率分布以及功率損耗等。潮流計(jì)算的結(jié)果是電力系統(tǒng)穩(wěn)定計(jì)算和故障分析的基礎(chǔ)，具體表現(xiàn)在以下方面： (1)在電網(wǎng)規(guī)劃階段,通過潮流計(jì)算,合理規(guī)劃電源容量及接入點(diǎn),合理規(guī)劃網(wǎng)架,選擇無功補(bǔ)償方案,滿足規(guī)劃水平的大、小方式下潮流交換控制、調(diào)峰、調(diào)相、調(diào)壓的要求。 (2)在編制年運(yùn)行方式時(shí),在預(yù)計(jì)負(fù)荷增長及新設(shè)備投運(yùn)基礎(chǔ)上,選擇典型方式進(jìn)行潮流計(jì)算,發(fā)現(xiàn)電網(wǎng)中薄弱環(huán)節(jié),供調(diào)

13、度員日常調(diào)度控制參考,并對規(guī)劃、基建部門提出改進(jìn)網(wǎng)架結(jié)構(gòu),加快基建進(jìn)度的建議。(3)正常檢修及特殊運(yùn)行方式下的潮流計(jì)算,用于日運(yùn)行方式的編制,指導(dǎo)發(fā)電廠開機(jī)方式,有功、無功調(diào)整方案及負(fù)荷調(diào)整方案,滿足線路、變壓器熱穩(wěn)定要求及電壓質(zhì)量要求。 (4)預(yù)想事故、設(shè)備退出運(yùn)行對靜態(tài)安全的影響分析及作出預(yù)想的運(yùn)行方式調(diào)整方案。 總結(jié)為在電力系統(tǒng)運(yùn)行方式和規(guī)劃方案的研究中，都需要進(jìn)行潮流計(jì)算以比較運(yùn)行方式或規(guī)劃供電方案的可行性、可靠性和經(jīng)濟(jì)性。同時(shí)，為了實(shí)時(shí)監(jiān)控電力系統(tǒng)的運(yùn)行狀態(tài)，也需要進(jìn)行大量而快速的潮流計(jì)算。因此，潮流計(jì)算是電力系統(tǒng)中應(yīng)用最廣泛、最基本和最重要的一種電氣運(yùn)算。在系統(tǒng)規(guī)劃設(shè)計(jì)和安排系統(tǒng)的

14、運(yùn)行方式時(shí)，采用離線潮流計(jì)算；在電力系統(tǒng)運(yùn)行狀態(tài)的實(shí)時(shí)監(jiān)控中，則采用在線潮流計(jì)算【3】。1.2 潮流計(jì)算方法的發(fā)展與前景電力系統(tǒng)潮流計(jì)算屬于穩(wěn)態(tài)分析范疇，不涉及系統(tǒng)元件的動(dòng)態(tài)特性和過渡過程。因此其數(shù)學(xué)模型不包含微分方程，是一組高階非線性方程。非線性代數(shù)方程組的解法離不開迭代，因此，潮流計(jì)算方法首先要求它是能可靠的收斂，并給出正確答案。隨著電力系統(tǒng)規(guī)模的不斷擴(kuò)大，潮流問題的方程式階數(shù)越來越高，目前已達(dá)到幾千階甚至上萬階，對這樣規(guī)模的方程式并不是采用任何數(shù)學(xué)方法都能保證給出正確答案的。這種情況促使電力系統(tǒng)的研究人員不斷尋求新的更可靠的計(jì)算方法。在用數(shù)字計(jì)算機(jī)求解電力系統(tǒng)潮流問題的開始階段，人們普遍

15、采用以節(jié)點(diǎn)導(dǎo)納矩陣為基礎(chǔ)的高斯-賽德爾迭代法（以下簡稱導(dǎo)納法）。這個(gè)方法的原理比較簡單，要求的數(shù)字計(jì)算機(jī)的內(nèi)存量也比較小，適應(yīng)當(dāng)時(shí)的電子數(shù)字計(jì)算機(jī)制作水平和電力系統(tǒng)理論水平，于是電力系統(tǒng)計(jì)算人員轉(zhuǎn)向以阻抗矩陣為主的逐次代入法（以下簡稱阻抗法）。 利用電子計(jì)算機(jī)進(jìn)行潮流計(jì)算從20世紀(jì)50年代中期就已經(jīng)開始。此后，潮流計(jì)算曾采用了各種不同的方法，這些方法的發(fā)展主要是圍繞著對潮流計(jì)算的一些基本要求進(jìn)行的。對潮流計(jì)算的要求可以歸納為下面幾點(diǎn)： （1）算法的可靠性或收斂性 （2）計(jì)算速度和內(nèi)存占用量 （3）計(jì)算的方便性和靈活性 20世紀(jì)60年代初，數(shù)字計(jì)算機(jī)已經(jīng)發(fā)展到第二代，計(jì)算機(jī)的內(nèi)存和計(jì)算速度發(fā)生了

16、很大的飛躍，從而為阻抗法的采用創(chuàng)造了條件。阻抗矩陣是滿矩陣，阻抗法要求計(jì)算機(jī)儲(chǔ)存表征系統(tǒng)接線和參數(shù)的阻抗矩陣。這就 需要較大的內(nèi)存量。阻抗法改善了電力系統(tǒng)潮流計(jì)算問題的收斂性，解決了導(dǎo)納法無法解決的一些系統(tǒng)的潮流計(jì)算，但是，阻抗法的主要缺點(diǎn)就是占用計(jì)算機(jī)的內(nèi)存很大，每迭代的計(jì)算量很大。當(dāng)系統(tǒng)不斷擴(kuò)大時(shí)，這些缺點(diǎn)就更加突出。近20多年來，潮流算法的研究仍然非?；钴S，但是大多數(shù)研究都是圍繞改進(jìn)牛頓法和P-Q分解法進(jìn)行的。此外，隨著人工智能理論的發(fā)展，遺傳算法、人工神經(jīng)網(wǎng)絡(luò)、模糊算法也逐漸被引入潮流計(jì)算。但是，到目前為止這些新的模型和算法還不能取代牛頓法和P-Q分解法的地位。由于電力系統(tǒng)規(guī)模的不斷擴(kuò)

17、大，對計(jì)算速度的要求不斷提高，計(jì)算機(jī)的并行計(jì)算技術(shù)也將在潮流計(jì)算中得到廣泛的應(yīng)用，成為重要的研究領(lǐng)域。通過幾十年的發(fā)展，潮流算法日趨成熟。近幾年，對潮流算法的研究仍然是如何改善傳統(tǒng)的潮流算法，即高斯-塞德爾法、牛頓法和快速解耦法。牛頓法，由于其在求解非線性潮流方程時(shí)采用的是逐次線性化的方法，為了進(jìn)一步提高算法的收斂性和計(jì)算速度，人們考慮采用將泰勒級數(shù)的高階項(xiàng)或非線性項(xiàng)也考慮進(jìn)來，于是產(chǎn)生了二階潮流算法。后來又提出了根據(jù)直角坐標(biāo)形式的潮流方程是一個(gè)二次代數(shù)方程的特點(diǎn)，提出了采用直角坐標(biāo)的保留非線性快速潮流算法【6】。第二章 電力網(wǎng)絡(luò)的數(shù)學(xué)模型 線性網(wǎng)絡(luò)的常用解法有節(jié)點(diǎn)電壓法和回路法，前者須列寫節(jié)

18、點(diǎn)電流平衡方程，后者則須列寫回路方程。本章重點(diǎn)介紹節(jié)點(diǎn)方程，以及節(jié)點(diǎn)導(dǎo)納矩陣【1】。2.1節(jié)點(diǎn)導(dǎo)納矩陣的形成及修改2.1.1節(jié)點(diǎn)導(dǎo)納矩陣的形成在圖2-1（a）的簡單電力系統(tǒng)中，若略去變壓器的勵(lì)磁功率和線路電容，負(fù)荷用阻抗表示，便可以得到一個(gè)有5個(gè)節(jié)點(diǎn)（包括零電位點(diǎn)）和7條支路的等值網(wǎng)絡(luò)，如圖2-1（b）所示。將接于節(jié)點(diǎn)1和4的電勢源和阻抗的串聯(lián)組合變換成等值的電流源和導(dǎo)納的并聯(lián)組合，變得到圖（c）的等值網(wǎng)絡(luò)，其中和分別稱為節(jié)點(diǎn)1和4的注入電流源。 圖2-1 電力系統(tǒng)及其網(wǎng)絡(luò)以零電位點(diǎn)作為計(jì)算節(jié)點(diǎn)電壓的參考點(diǎn)，根據(jù)基爾霍夫定律，可以寫出4個(gè)獨(dú)立節(jié)點(diǎn)的電流平衡方程如下 （2-1）上述方程組經(jīng)過整理

19、可以寫成 （2-2）式中，；。一般的，對于有個(gè)獨(dú)立節(jié)點(diǎn)的網(wǎng)絡(luò)，可以列寫個(gè)節(jié)點(diǎn)方程 （2-3）也可以用矩陣寫成 （2-4）或縮寫為 （2-5）矩陣稱為節(jié)點(diǎn)導(dǎo)納矩陣。它的對角線元素稱為節(jié)點(diǎn)的自導(dǎo)納，其值等于接于節(jié)點(diǎn)的所有支路導(dǎo)納之和。非對角線元素稱為節(jié)點(diǎn)、 間的互導(dǎo)納，它等于直接接于節(jié)點(diǎn)、間的支路導(dǎo)納的負(fù)值。若節(jié)點(diǎn)、間不存在直接支路，則有。由此可知節(jié)點(diǎn)導(dǎo)納矩陣是一個(gè)稀疏的對稱矩陣。2.1.2節(jié)點(diǎn)導(dǎo)納矩陣的修改在電力系統(tǒng)中，接線方式或運(yùn)行狀態(tài)等均會(huì)發(fā)生變化，從而使網(wǎng)絡(luò)接線改變。比如一臺變壓器支路的投入或切除，均會(huì)使與之相連的節(jié)點(diǎn)的自導(dǎo)納或互導(dǎo)納發(fā)生變化，而網(wǎng)絡(luò)中其它部分結(jié)構(gòu)并沒有改變，因此不必從新形

20、成節(jié)點(diǎn)導(dǎo)納矩陣，而只需對原有的矩陣作必要的修改就可以了。現(xiàn)在幾種典型的接線變化說明具體的修改方法。 圖2-2 電力接線的改變（a）增加支路和節(jié)點(diǎn)；（b）增加支路；（c）切除支路；（d）改變支路參數(shù)；（e）改變變壓器變比 （1）從原有網(wǎng)絡(luò)的節(jié)點(diǎn)引出一條導(dǎo)納為的支路，為新增加的節(jié)點(diǎn)，如圖2-2(a)所示。由于新增加了一個(gè)節(jié)點(diǎn)，所以節(jié)點(diǎn)導(dǎo)納矩陣增加一階，矩陣作如下修改： 1）原有節(jié)點(diǎn)的自導(dǎo)納的增量； 2）新增節(jié)點(diǎn)的自導(dǎo)納； 3）新增的非對角元素；其它新增的非對角元均為零。（2）在原有網(wǎng)絡(luò)的節(jié)點(diǎn)與j之間增加一條導(dǎo)納為的支路，如圖2-2(b)所示。則與、有關(guān)的元素應(yīng)作如下修改： 1）節(jié)點(diǎn)、的自導(dǎo)納增量；

21、 2）節(jié)點(diǎn)、的互導(dǎo)納增量。 （3）在網(wǎng)絡(luò)的原有節(jié)點(diǎn)、之間切除一條導(dǎo)納為的支路，如圖2-2（c）所示，其相當(dāng)在、之間增加一條導(dǎo)納為的支路，因此與、有關(guān)的元素應(yīng)作以下修改：1）節(jié)點(diǎn)、的自導(dǎo)納增量；2）節(jié)點(diǎn)、j之間的互導(dǎo)納增量；（4）原有網(wǎng)絡(luò)節(jié)點(diǎn)、之間的導(dǎo)納由變成，相當(dāng)于在節(jié)點(diǎn)、之間切除一條導(dǎo)納為的支路，再增加一條導(dǎo)納為的支路，如圖2-2（d）所示。則與、有關(guān)的元素應(yīng)作如下修改： 1）節(jié)點(diǎn)、的自導(dǎo)納增量； 2）節(jié)點(diǎn)、的互導(dǎo)納增量。（5）原有網(wǎng)絡(luò)節(jié)點(diǎn)、之間變壓器的變比由變?yōu)?，即相?dāng)于切除一臺變比為的變壓器，再投入一臺變比為的變壓器，如圖2-2（e）變壓器型等值電路，圖中為與變壓器原邊基準(zhǔn)電壓對應(yīng)的變壓

22、器導(dǎo)納標(biāo)幺值，則與、有關(guān)的元素應(yīng)作如下修改：1）節(jié)點(diǎn)的自導(dǎo)納增量；節(jié)點(diǎn)的自導(dǎo)納增量；2）節(jié)點(diǎn)與之間的互導(dǎo)納增量。2.2節(jié)點(diǎn)導(dǎo)納矩陣元素的物理意義節(jié)點(diǎn)導(dǎo)納矩陣的元素已在上一節(jié)作了說明，現(xiàn)在進(jìn)一步討論這些元素的物理意義。如果令 代入2-3的各式，可得 或 （2-6）當(dāng)時(shí)，公式2-6說明，當(dāng)網(wǎng)絡(luò)中除節(jié)點(diǎn)以外所有節(jié)點(diǎn)都接地時(shí)，從節(jié)點(diǎn)注入網(wǎng)絡(luò)的電流同施加于節(jié)點(diǎn)的電壓之比，即等于節(jié)點(diǎn)的自導(dǎo)納。換句話說，自導(dǎo)納是節(jié)點(diǎn)以外的所有節(jié)點(diǎn)都接地時(shí)節(jié)點(diǎn)對地的總導(dǎo)納。顯然，應(yīng)等于與節(jié)點(diǎn)相接的各支路導(dǎo)納之和，即 （2-7）式中 ，為節(jié)點(diǎn)與零電位節(jié)點(diǎn)之間的支路導(dǎo)納；為節(jié)點(diǎn)與節(jié)點(diǎn)之間的支路導(dǎo)納。當(dāng)時(shí)，公式2-6說明，當(dāng)網(wǎng)絡(luò)中

23、除節(jié)點(diǎn)以外所有節(jié)點(diǎn)都接地時(shí)，從節(jié)點(diǎn)注入網(wǎng)絡(luò)的電流同施加于節(jié)點(diǎn)的電壓之比，即等于節(jié)點(diǎn)、的互導(dǎo)納。在這種情況下，節(jié)點(diǎn)的電流實(shí)際上是自網(wǎng)絡(luò)流出并進(jìn)入地中的電流，所以應(yīng)等于與節(jié)點(diǎn)、之間的支路導(dǎo)納的負(fù)值，即 （2-8）不難理解。若節(jié)點(diǎn)和沒有支路直接相聯(lián)時(shí)，便有。在圖2-2 所示的網(wǎng)絡(luò)中，單獨(dú)在節(jié)點(diǎn)2接上電源，而將其余節(jié)點(diǎn)都接地。圖2-3 自導(dǎo)納和互導(dǎo)納的確定根據(jù)上述節(jié)點(diǎn)自導(dǎo)納和互導(dǎo)納的定義，可得因，故。從圖中也可以清楚地看到，節(jié)點(diǎn)4、5和6同節(jié)點(diǎn)2都沒有直接的支路關(guān)系。導(dǎo)納矩陣元素的其它元素也可以用類似方法確定。節(jié)點(diǎn)導(dǎo)納矩陣的主要特點(diǎn)是：（1）節(jié)點(diǎn)導(dǎo)納矩陣的元素很容易根據(jù)網(wǎng)絡(luò)接線圖和支路參數(shù)直觀的求得，

24、形成節(jié)點(diǎn)導(dǎo)納矩陣的程序比較簡單。（2）節(jié)導(dǎo)納矩陣是稀疏矩陣。它的對角線元素一般不為零，但在非對角線元素中則存在不少零元素。在電力系統(tǒng)的接線圖中，一般每個(gè)節(jié)點(diǎn)同平均不超過個(gè)其它節(jié)點(diǎn)有直接的支路聯(lián)接，因此在導(dǎo)納矩陣的非對角線元素中每行平均僅有個(gè)非零元素，其余的元素都為零。如果在程序設(shè)計(jì)中設(shè)法排除零元素的貯存和運(yùn)算，就可以大大地節(jié)省貯存單元和提高計(jì)算速度。10洛陽理工學(xué)院畢業(yè)設(shè)計(jì)（論文）第三章 牛頓拉夫遜法潮流分布計(jì)算 牛頓拉夫遜法(簡稱牛頓法)在數(shù)學(xué)上是求解非線性代數(shù)方程式的有效方法，其要點(diǎn)是把非線性方程式的求解過程變成反復(fù)地對相應(yīng)的線性方程式進(jìn)行求解的過程，即通常所稱的逐次線性化過程【1】。3.

25、1牛頓拉夫遜法簡介3.1.1牛頓拉夫遜法簡介牛頓拉夫遜法（NewtonRaphson法）是求解非線性方程代數(shù)方程組的有效迭代計(jì)算方法。在牛頓拉夫遜法的每一次迭代過程中，對非線性方程通過線性化處理逐步近似。下面以單變量加以說明。設(shè)有單變量非線性方程 (3-1)求解此方程時(shí)。先給出解的近似值它與真解的誤差為，則將滿足方程，即 (3-2)將(3-8)式左邊的函數(shù)在附近展成泰勒級數(shù)，于是便得 （3-3）式中,分別為函數(shù)在處的一階導(dǎo)數(shù)，.，n階導(dǎo)數(shù)。如果差值很小，3-9式右端的二次及以上階次的各項(xiàng)均可略去。于是，3-9便簡化為 0 (3-4)這是對于變量的修正量的現(xiàn)行方程式，亦稱修正方程式。解此方程可得

26、修正量 (3-5)用所求的去修正近似解，變得 (3-6)由于3-10是略去高次項(xiàng)的簡化式，因此所解出的修正量也只是近似值。修正后的近似解同真解仍然有誤差。但是，這樣的迭代計(jì)算可以反復(fù)進(jìn)行下去，迭代計(jì)算的通式是 (3-7)迭代過程的收斂判據(jù)為 (3-8)或 (3-9)式中，為預(yù)先給定的小正數(shù)。 這種解法的幾何意義可以從圖31得到說明。函數(shù)yf(x)為圖中的曲線。f(x)0的解相當(dāng)于曲線與x軸的交點(diǎn)。如果第k次迭代中得到，則過點(diǎn)作一切線，此切線同x軸的交點(diǎn)便確定了下一個(gè)近似值。由此可見，牛頓拉夫遜法實(shí)質(zhì)上就是切線法，是一種逐步線性化的方法。 應(yīng)用牛頓法求解多變量非線性方程組3-1時(shí)，假定已給出各變

27、量的初值，. ，令，. 分別為各變量的修正量，使其滿足方程3-2即 (3-10)將上式中的n個(gè)多元函數(shù)在初始值附近分別展成泰勒級數(shù),并略去含有，二次及以上階次的各項(xiàng)，便得. (3-11)方程式3-17也可以寫成矩陣形式 (3-12)方程式3-18是對于修正量， 的線性方程組,稱為牛頓法的修正方程式.利用高斯消去法或三角分解法可以解出修正量，。然后對初始近似值進(jìn)行修正 (i=1,2,.,n) (3-13)如此反復(fù)迭代，在進(jìn)行k1次迭代時(shí)，從求解修正方程式 (3-14)得到修正量，并對各變量進(jìn)行修正 (i=1,2,n) (3-15)式3-20和3-21也可以縮寫為 (3-16)和 (3-17) 式

28、中的X和分別是由n個(gè)變量和修正量組成的n維列向量；F(X)是由n個(gè)多元函數(shù)組成的n維列項(xiàng)量；J是n階方陣，稱為雅可比矩陣，它的第i、j個(gè)元素是第n個(gè)函數(shù)對第j個(gè)變量的偏導(dǎo)數(shù)；上角標(biāo)(k)表示陣的每一個(gè)元素都在點(diǎn)處取值。迭代過程一直到滿足收斂判據(jù) (3-18)或 (3-19)為止。和為預(yù)先給定的小正數(shù)。 將牛頓拉夫遜法用于潮流計(jì)算，要求將潮流方程寫成形如方程式3-1的形式。由于節(jié)點(diǎn)電壓可以采用不同的坐標(biāo)系表示，牛頓拉夫遜法潮流計(jì)算也將相應(yīng)的采用不同的計(jì)算公式。3.1.2牛頓拉夫遜法的幾何意義從幾何意義上，牛頓拉夫遜法實(shí)質(zhì)上就是切線法，是一種逐步線性化的方法。圖3-1 牛頓拉夫遜方法的幾何意義27

29、3.2牛頓拉夫遜法計(jì)算潮流分布以下討論的是用直角坐標(biāo)形式的牛頓拉夫遜法潮流的求解過程。當(dāng)采用直角坐標(biāo)時(shí)，潮流問題的待求量為各節(jié)點(diǎn)電壓的實(shí)部和虛部兩個(gè)分量，由于平衡節(jié)點(diǎn)的電壓向量是給定的，因此待求量共2(n-1)需要2(n-1)個(gè)方程式。事實(shí)上，除了平衡節(jié)點(diǎn)的功率方程式在迭代過程中沒有約束作用以外，其余每個(gè)節(jié)點(diǎn)都可以列出兩個(gè)方程式。功率方程可以寫為 （3-20）也可以寫成 （3-21）根據(jù)節(jié)點(diǎn)電壓的兩種不同表示方法，可以得到兩種不同的牛頓拉夫遜法潮流計(jì)算方法。鑒于本次的設(shè)計(jì)題目，現(xiàn)只對節(jié)點(diǎn)電壓以直角坐標(biāo)表示時(shí)的牛頓拉夫遜法潮流計(jì)算做以下介紹。節(jié)點(diǎn)電壓以直角坐標(biāo)表示，即，為節(jié)點(diǎn)電壓實(shí)部，為節(jié)點(diǎn)電壓

30、虛部，功率方程可寫為 （3-22） （3-23）對于節(jié)點(diǎn)，電壓有效值為設(shè)定值，而實(shí)部和虛部的比例是可變的，它們之間的關(guān)系為 （3-24）等式425、式427構(gòu)成關(guān)于節(jié)點(diǎn)的兩個(gè)約束條件。對于第次迭代，可寫出 （3-25） （3-26） （3-27）可列出修正方程。對于除平衡節(jié)點(diǎn)以外的所有節(jié)點(diǎn) （3-28）對于節(jié)點(diǎn) （3-29）對于節(jié)點(diǎn) （3-30）如系統(tǒng)具有個(gè)節(jié)點(diǎn)，其中個(gè)為節(jié)點(diǎn)，則式4-284-29式可寫成如下矩陣形式 （3-31）式中、和為向量，為向量，為向量。、為矩陣，、為矩陣，、為矩陣。其雅可比矩陣元素為 （3-32）上式中和分別是節(jié)點(diǎn)注入電流的實(shí)部和虛部，即 （3-33） （3-34）由

31、式434可求得第次迭代的修正量和，從而可得到新的解 （3-35）這樣反復(fù)計(jì)算，直到收斂至要求的精度。收斂指標(biāo)一般取所有節(jié)點(diǎn)的和或。3、牛頓拉夫遜法潮流計(jì)算的步驟求取雅可比矩陣是牛頓拉夫遜法的一項(xiàng)重要的工作。電力系統(tǒng)潮流計(jì)算的雅可比矩陣具有以下特點(diǎn)：（1）雅可比矩陣為一非奇異方陣。當(dāng)節(jié)點(diǎn)電壓以極坐標(biāo)表式時(shí)，該矩陣為階方陣；節(jié)點(diǎn)電壓以直角坐標(biāo)表式時(shí)，為階。（2）矩陣元素與節(jié)點(diǎn)電壓有關(guān)，故每次迭代時(shí)都要重新計(jì)算。（3）與導(dǎo)納矩陣具有相似的結(jié)構(gòu)，當(dāng)時(shí)，、均為0，因此，也是高度稀疏的矩陣。針對這一特點(diǎn)利用求解稀疏矩陣技巧，對減少計(jì)算所需的內(nèi)存和時(shí)間是很有好處的。（4）具有結(jié)構(gòu)對稱性，但數(shù)值不對稱。例如，

32、由于各節(jié)點(diǎn)電壓不同，因而。需要指出的是，當(dāng)計(jì)算過程中發(fā)生節(jié)點(diǎn)的無功功率越限時(shí)，節(jié)點(diǎn)要轉(zhuǎn)化成節(jié)點(diǎn)。此時(shí)，對節(jié)點(diǎn)電壓以極坐標(biāo)形式表示的修正方程，需增加一對應(yīng)于該節(jié)點(diǎn)無功功率不平衡量的關(guān)系式，因而式333中的誤差向量、電壓向量和雅可比矩陣都有相應(yīng)變動(dòng)。當(dāng)采用直角坐標(biāo)表式時(shí)，要增加一個(gè)對應(yīng)于該節(jié)點(diǎn)的無功功率不平衡量的關(guān)系式，同時(shí)要減少一個(gè)對應(yīng)于設(shè)定節(jié)點(diǎn)電壓約束條件的形式。用牛頓拉夫遜法計(jì)算電力系統(tǒng)潮流的基本步驟：求解過程大致可以分為以下步驟：(1)形成節(jié)點(diǎn)導(dǎo)納矩陣；(2)將各節(jié)點(diǎn)電壓設(shè)初值 ；(3)將節(jié)點(diǎn)初值代入相關(guān)求式，求出修正方程式的常數(shù)項(xiàng)向量；(4)將節(jié)點(diǎn)電壓初值代入求式，求出雅可比矩陣元素；(

33、5)求解修正方程，求修正向量；(6)求取節(jié)點(diǎn)電壓的新值；(7)檢查是否收斂，如不收斂，則以各節(jié)點(diǎn)電壓的新值作為初值自第3步重新開始進(jìn)行狹義次迭代，否則轉(zhuǎn)入下一步；(8)計(jì)算支路功率分布，PV節(jié)點(diǎn)無功功率和平衡節(jié)點(diǎn)功率。圖中32中示出牛頓拉夫遜法計(jì)算潮流的程序框圖。圖中為實(shí)現(xiàn)給定的最大迭代次數(shù)，當(dāng)實(shí)際迭代次數(shù)時(shí)，即認(rèn)為計(jì)算不收斂。圖32 牛頓拉夫遜法計(jì)算潮流程序框圖 第四章 C程序設(shè)計(jì)的基本知識C語言是一種計(jì)算機(jī)程序設(shè)計(jì)語言。它既有高級語言的特點(diǎn)，又具有匯編語言的特點(diǎn)。它可以作為系統(tǒng)設(shè)計(jì)語言，編寫工作系統(tǒng)應(yīng)用程序，也可以作為應(yīng)用程序設(shè)計(jì)語言，編寫不依賴計(jì)算機(jī)硬件的應(yīng)用程序。因此，它的應(yīng)用范圍廣泛

34、。本設(shè)計(jì)就是就是依據(jù)C程序進(jìn)行編程運(yùn)算【9】。4.1 C語言的特點(diǎn)1. C是中級語言它把高級語言的基本結(jié)構(gòu)和語句與低級語言的實(shí)用性結(jié)合起來。C語言可以象匯編語言一樣對位、字節(jié)和地址進(jìn)行操作, 而這三者是計(jì)算機(jī)最基本的工作單元。 2. C是結(jié)構(gòu)式語言結(jié)構(gòu)式語言的顯著特點(diǎn)是代碼及數(shù)據(jù)的分隔化,即程序的各個(gè)部分除了必要的信息交流外彼此獨(dú)立。這種結(jié)構(gòu)化方式可使程序?qū)哟吻逦? 便于使用、維護(hù)以及調(diào)試。C 語言是以函數(shù)形式提供給用戶的,這些函數(shù)可方便的調(diào)用, 并具有多種循環(huán)、條件語句控制程序流向, 從而使程序完全結(jié)構(gòu)化。 3. C語言功能齊全C 語言具有各種各樣的數(shù)據(jù)類型, 并引入了指針概念,可使程序效率

35、更高。另外C 語言也具有強(qiáng)大的圖形功能, 支持多種顯示器和驅(qū)動(dòng)器。而且計(jì)算功能、邏輯判斷功能也比較強(qiáng)大,可以實(shí)現(xiàn)決策目的編游戲，編3D游戲，做數(shù)據(jù)庫，做聯(lián)眾世界，做聊天室，做PHOTOSHOP做FLASH，做3DMAX。 4. C語言適用范圍大C 語言還有一個(gè)突出的優(yōu)點(diǎn)就是適合于多種操作系統(tǒng), 如DOS、UNIX,也適用于多種機(jī)型。 C語言對操作系統(tǒng)和系統(tǒng)使用程序以及需要對硬件進(jìn)行操作的場合，用C語言明顯優(yōu)于其它解釋型高級語言，有一些大型應(yīng)用軟件也是用C語言編寫的。 C語言具有繪圖能力強(qiáng)，可移植性，并具備很強(qiáng)的數(shù)據(jù)處理能力，因此適于編寫系統(tǒng)軟件，三維，二維圖形和動(dòng)畫。它是數(shù)值計(jì)算的高級語言。4

36、.2 C語言的發(fā)展過程 C語言的原型ALGOL 60語言（也稱為A語言） 。 1963年，劍橋大學(xué)將ALGOL 60語言發(fā)展成為CPL(Combined Programming Language)語言。 1967年，劍橋大學(xué)的Matin Richards 對CPL語言進(jìn)行了簡化，于是產(chǎn)生了BCPL語言。 1970年，美國貝爾實(shí)驗(yàn)室的Ken Thompson將BCPL進(jìn)行了修改，并為它起了一個(gè)有趣的名字“B語言”。意思是將CPL語言煮干，提煉出它的精華。并且他用B語言寫了第一個(gè)UNIX操作系統(tǒng)。 而在1973年，B語言也給人“煮”了一下，美國貝爾實(shí)驗(yàn)室的D.M.RITCHIE在B語言的基礎(chǔ)上最終

37、設(shè)計(jì)出了一種新的語言，他取了BCPL的第二個(gè)字母作為這種語言的名字，這就是C語言。 為了使UNIX操作系統(tǒng)推廣，1977年Dennis M. Ritchie 發(fā)表了不依賴于具體機(jī)器系統(tǒng)的C語言編譯文本可移植的C語言編譯程序。1978年Brian W. Kernighian和Dennis M. Ritchie出版了名著The C Programming Language，從而使C語言成為目前世界上流行最廣泛的高級程序設(shè)計(jì)語言。 1987年，隨著微型計(jì)算機(jī)的日益普及,出現(xiàn)了許多C語言版本。由于沒有統(tǒng)一的標(biāo)準(zhǔn),使得這些C語言之間出現(xiàn)了一些不一致的地方。為了改變這種情況,美國國家標(biāo)準(zhǔn)研究所(ANSI)

38、為C語言制定了一套ANSI標(biāo)準(zhǔn), 成為現(xiàn)行的C語言標(biāo)準(zhǔn) 3.C語言的主要特點(diǎn) 。C語言發(fā)展迅速, 而且成為最受歡迎的語言之一,主要因?yàn)樗哂袕?qiáng)大的功能。許多著名的系統(tǒng)軟件, 如DBASE PLUS、DBASE 都是由C語言編寫的。用C語言加上一些匯編語言子程序, 就更能顯示C語言的優(yōu)勢了,象PC- DOS 、WORDSTAR等就是用這種方法編寫的。 1990年，國際化標(biāo)準(zhǔn)組織ISO接受了87 ANSI C為ISO C的標(biāo)準(zhǔn)（ISO9899-1990）。1994年，ISO修訂了C語言的標(biāo)準(zhǔn)。目前流行的C語言編譯系統(tǒng)大多是以ANSI C為基礎(chǔ)進(jìn)行開發(fā)的，但不同版本的C編譯系統(tǒng)所實(shí)現(xiàn)的語言功能和語法

39、規(guī)則略有差別。第五章 計(jì)算實(shí)例5.1算例圖5-1為一五結(jié)點(diǎn)系統(tǒng)，各支路參數(shù)均為標(biāo)幺值。假定結(jié)點(diǎn)1、2、3為PQ節(jié)點(diǎn)，結(jié)點(diǎn)4為PV節(jié)點(diǎn)、結(jié)點(diǎn)5為平衡結(jié)點(diǎn)，試分別用直角坐標(biāo)和極坐標(biāo)牛頓拉夫遜法計(jì)算其潮流。取收斂判據(jù)為|DPi|£10-5和|DQi(DVi2)|£10-5。給定：S1=-1.6-j0.8 S2=-2.0-j1.0 S3=-3.7-j1.3 P4=5.0|V1(0)|=|V2(0)|=|V3(0)|=1.00 |V4(0)|=|V5(0)|=1.05 圖5-1 節(jié)點(diǎn)系統(tǒng)5.2程序計(jì)算結(jié)果節(jié)點(diǎn)原始數(shù)據(jù)：scanf the 1th numbers:2 -4 0 0.01

40、5 1.05scanf the 2th numbers:4 3 0.08 0.30 0.25scanf the 3th numbers:4 5 0.04 0.25 0.25scanf the 4th numbers:5 3 0.1 0.35 0scanf the 5th numbers:-3 1 0 0.03 1.05導(dǎo)納矩陣：Y11=0.000000-j33.33333Y12=0.000000+j0.000000Y13=0.000000+j1.796945 Y14=0.000000+j0.000000Y15=0.000000+j0.000000Y21=0.000000+j0.000000Y2

41、2=0.000000-j66.666672Y23=0.000000+j0.000000Y24=0.000000+j63.492069Y25=0.000000+j0.000000Y31=0.000000+j1.796945Y32=0.000000+j0.000000Y33=1.584592-j7.214919Y34=-0.829875+j3.112033Y35=-0.754717+j2.641510Y41=0.000000+j0.000000Y42=0.000000+j63.492069Y43=-0.829875+j3.112033Y44=1.453900-j66.980827Y45=-0.62

42、4025+j3.900156Y51=0.000000+j0.000000Y52=0.000000+j0.000000Y53=-0.754717+j2.641510Y54=-0.624025+j3.900156Y55=1.378742-j6.291666 5.3平衡節(jié)點(diǎn)功率和線路功率5.4數(shù)據(jù)分析有題意可得：收斂判據(jù)取|DPi|£10-5和|DQi(DVi2)|£10-5，程序計(jì)算結(jié)果滿足要求。該數(shù)據(jù)是采用的軟件Microsoft visual C+ 6.0編譯得到的結(jié)果，與同組采用軟件Matlab6.5編譯得到的結(jié)果相同。牛頓拉夫遜法潮流計(jì)算的核心問題是修正方程式的建立和求

43、解。運(yùn)用Matlab6.5進(jìn)行編譯時(shí)，每次迭代都需要先解修正方程式，然后用解得的各節(jié)點(diǎn)電壓變量（修正量）求各節(jié)點(diǎn)電壓的新值（修正后值），在得到迭代收斂后才能進(jìn)行計(jì)算平衡節(jié)點(diǎn)功率和線路功率。結(jié)論本設(shè)計(jì)是在全面理解用牛頓法計(jì)算潮流原理的基礎(chǔ)上，利用程序設(shè)計(jì)語言C語言編制用計(jì)算機(jī)計(jì)算電力系統(tǒng)潮流的程序設(shè)計(jì)方案，此設(shè)計(jì)方案的任務(wù)是完成牛頓迭代法進(jìn)行潮流計(jì)算。本設(shè)計(jì)采用直角坐標(biāo)形式的牛頓拉夫遜法作常規(guī)潮流計(jì)算。P-Q分解法利用了電力系統(tǒng)的一些特有的運(yùn)行特點(diǎn)，對牛頓拉夫遜法作了簡化，可提高計(jì)算速度，但較難理解，牛頓拉夫遜法的雅克比矩陣在每次迭代過程中都有變化，需要重新形成和求解，這占用了計(jì)算的大部分時(shí)間，

44、成為牛頓拉夫遜計(jì)算速度不能提高的主要因素，但收斂性好，物理概念也較為清晰。應(yīng)用C語言進(jìn)行編程層次清晰，便于按模塊化方式組織，易于調(diào)試和維護(hù)，C語言的表現(xiàn)能力和處理能力極強(qiáng)，極大的提高了計(jì)算效率，同時(shí)也培養(yǎng)了我們快速的掌握應(yīng)用程序的能力和解決實(shí)際問題的能力。潮流計(jì)算本身，實(shí)際上就是求解一組非線性方程的代數(shù)方程，鑒于非線性系統(tǒng)所固有的特性，完全有可能出現(xiàn)重解，而其機(jī)理是亟待深入探討的課題之一；潮流計(jì)算的電壓不穩(wěn)定性，而無論怎樣調(diào)整系統(tǒng)中的無功電源都不能促使遠(yuǎn)遠(yuǎn)低于要求值的受端電壓回復(fù)正常，或者缺乏抗干擾的阻尼能力而導(dǎo)致系統(tǒng)運(yùn)行狀態(tài)的持續(xù)振蕩。通過本次設(shè)計(jì)，我們更好的理解和掌握電力系統(tǒng)潮流計(jì)算的原理

45、和方法，也更加系統(tǒng)的理解C程序設(shè)計(jì)語言的編程技巧。今后的研究方向仍圍繞著如何改善傳統(tǒng)潮流計(jì)算方法以及提高計(jì)算機(jī)潮流計(jì)算的可靠性、收斂性以及速度和靈活性。謝 辭本次畢業(yè)設(shè)計(jì)是在我的老師吳茜瓊的親切關(guān)懷和悉心指導(dǎo)下完成的。她嚴(yán)肅的科學(xué)態(tài)度，嚴(yán)謹(jǐn)?shù)闹螌W(xué)精神，精益求精的工作作風(fēng)，深深地感染和激勵(lì)著我。從課題的選擇到項(xiàng)目的最終完成，吳老師都始終給予我細(xì)心的指導(dǎo)和不懈的支持。兩年多來，吳老師不僅在學(xué)業(yè)上給我以精心指導(dǎo)，同時(shí)還在思想、生活上給我以無微不至的關(guān)懷，在此謹(jǐn)向吳老師致以誠摯的謝意和崇高的敬意。    在此，我還要感謝我們畢業(yè)設(shè)計(jì)組的其他三位同學(xué)，正是由于你們的幫助和支持，我才能

46、克服一個(gè)一個(gè)的困難和疑惑，直至本文的順利完成。特別感謝我的同學(xué)張童年，他對本課題做了不少工作，給予我不少的幫助。    在論文即將完成之際，我的心情無法平靜，從開始進(jìn)入課題到論文的順利完成，有多少可敬的師長、同學(xué)、朋友給了我無言的幫助，在這里請接受我誠摯的謝意!最后我還要感謝培養(yǎng)我長大含辛茹苦的父母，謝謝你們!29參考文獻(xiàn)1 何仰贊等.電力系統(tǒng)分析上冊M武漢：華中理工大學(xué)出版社.2 何仰贊等.電力系統(tǒng)分析下冊M武漢：華中理工大學(xué)出版社.3 諸俊偉等.電力系統(tǒng)分析M.北京：中國電力出版社，1995.4 周全仁等.電網(wǎng)計(jì)算與程序設(shè)計(jì)M.長沙：湖南科學(xué)技術(shù)出版社,1983.5丁化

47、成.單片機(jī)應(yīng)用技術(shù)A.北京:北京航空航天大學(xué)出版社,2000.6張伯明,陳壽孫.高等電力網(wǎng)絡(luò)分析M.北京：清華大學(xué)出版社,1996.7陳珩.電力系統(tǒng)穩(wěn)態(tài)分析. 水利電力出版社,19948程衛(wèi)國,馮峰,王雪梅,劉藝品. MATLAB5.3 精要、編程及高級應(yīng)用. 機(jī)械工業(yè)出版社,20009Steven Holzner(美)著. 詳實(shí)翻譯組譯. Visual Basic 6 技術(shù)內(nèi)幕. 機(jī)械工業(yè)出版社,200010楊少兵,駱平. 電力系統(tǒng)分析的教學(xué)軟件開發(fā)，電力系統(tǒng)潮流分析. 華北電力技術(shù),2000(10)11劉支斂等. 應(yīng)用程序接口用戶指南. 科學(xué)出版社,200012張晗,潘正運(yùn),金曉燕.智能“

48、旅游電子超市”系統(tǒng)的研究與設(shè)計(jì)13微計(jì)算機(jī)信息,2005,12- 3:13- 15.14趙茂泰.智能儀器原理及應(yīng)用. 電子工業(yè)出版社，2002.15 熊靜琪. 計(jì)算機(jī)控制技術(shù). 電子工業(yè)出版社，ISBN 7- 5053-8706- 5.16王幸之.單片機(jī)應(yīng)用系統(tǒng)抗干擾技術(shù)A.北京:北京航空航天大學(xué)出版社,2002.附錄 計(jì)算機(jī)程序#include <stdio.h>#include <math.h>#include<stdlib.h> #include<conio.h> #define NULL 0#define N 30#define LEN

49、sizeof(struct ZL)struct ZL int I; int J; float R; float X; float K; float C; numN=0;struct Nodeint type; float P; float Q_V; nodeN=0,balance=0;struct Node_Cfloat e; float f;CSN=0;float b2*N=0;float YJ2*N2*N=0;int m;int L;int MARK=1;float x2*N=0;float max=0;struct YYfloat G; float B; YNN=0;FILE *fp;c

50、har filename20;ZDJ(struct ZL A) int i; i=A.J; Yi-1i-1.G+=A.R/(A.K*A.K*(A.R*A.R+A.X*A.X); Yi-1i-1.B+=-A.X/(A.K*A.K*(A.R*A.R+A.X*A.X);ZDI(struct ZL A) int i; i=A.I; Yi-1i-1.G+=A.R/(A.R*A.R+A.X*A.X); Yi-1i-1.B+=-A.X/(A.R*A.R+A.X*A.X);YL(struct ZL A)int i,j; j=A.J;i=A.I; Yj-1j-1.B+=A.C; Yi-1i-1.B+=A.C;H

51、D(struct ZL A)int i,j; j=A.J;i=A.I; Yi-1j-1.B=Yj-1i-1.B=A.X/(A.K*(A.R*A.R+A.X*A.X); Yi-1j-1.G=Yj-1i-1.G=-A.R/(A.K*(A.R*A.R+A.X*A.X);XG()int i,j,r=1,x,y;char chs;struct ZL p; for(;r;) system("cls"); printf("1.修改數(shù)據(jù)n2.添加數(shù)據(jù)n3.刪除數(shù)據(jù)n4.退出n請選擇:"); scanf("%d",&j); getchar(); switch(j) case 1:printf("請輸入要修改的首末結(jié)點(diǎn)標(biāo)號I,J:"); scanf("%d,%d",&x,&y); getchar(); for(i=0;i<L;i+) if(x=numi.I&&y=numi.J)|(x=numi.J&&y=numi.I) printf("這條支路參數(shù):I:%d,J:%d,R:%f,X:%f,K:%f,C:%fn",numi.I,numi.J,numi.R,numi.X,nu