電力系統(tǒng)分析課程設(shè)計(jì)報(bào)告-4

PAGE26-電力系統(tǒng)分析課程設(shè)計(jì)報(bào)告題目：電力系統(tǒng)三相對(duì)稱短路計(jì)算專業(yè)：電氣工程及其自動(dòng)化班級(jí)：姓名：學(xué)號(hào)：指導(dǎo)教師：平時(shí)表現(xiàn)(10%)驗(yàn)收答辯(50%)設(shè)計(jì)報(bào)告(40%)總評(píng)成績(jī)目錄電力系統(tǒng)分析 -0-第一章設(shè)計(jì)目的與任務(wù) -2-1.1設(shè)計(jì)目的 -2-1.2設(shè)計(jì)任務(wù) -2-第二章基礎(chǔ)理論與原理 -2-2.1對(duì)稱短路計(jì)算的基本方法 -2-2.2用節(jié)點(diǎn)阻抗矩陣的計(jì)算方法 -4-2.3用節(jié)點(diǎn)導(dǎo)納矩陣的計(jì)算方法 -6-2.4用三角分解法求解節(jié)點(diǎn)阻抗矩陣 -7-2.5短路發(fā)生在線路上任意處的計(jì)算方法 -8-第三章程序設(shè)計(jì) -9-3.1變量說(shuō)明 -9-3.2程序流程圖 -10-3.2.1主程序流程圖 -11-3.2.2導(dǎo)納矩陣流程圖 -12-3.2.3三角分解法流程圖 -13-3.3程序源代碼見(jiàn)附錄1 -14-第四章結(jié)果分析 -14-第五章收獲與建議 -15-參考文獻(xiàn) -17-附錄 -18-附錄1：程序源代碼 -18-附錄2：測(cè)試系統(tǒng)數(shù)據(jù)與系統(tǒng)圖 -23-附錄3：測(cè)試系統(tǒng)的運(yùn)行結(jié)果 -25-第一章設(shè)計(jì)目的與任務(wù)1.1設(shè)計(jì)目的加深理解并鞏固電力系統(tǒng)發(fā)生短路的基本知識(shí)。掌握電力輸電、配電系統(tǒng)發(fā)生三相對(duì)稱短路計(jì)算的計(jì)算方法和計(jì)算機(jī)計(jì)算的方法。熟悉C語(yǔ)言或Matlab計(jì)算程序的基本編寫(xiě)過(guò)程。為今后從事與電力系統(tǒng)相關(guān)的設(shè)計(jì)、運(yùn)行、管理及科學(xué)研究工作打下良好的基礎(chǔ)。1.2設(shè)計(jì)任務(wù)在兩周的時(shí)間內(nèi)，利用C語(yǔ)言或Matlab編寫(xiě)電力系統(tǒng)三相對(duì)稱短路計(jì)算的程序，并調(diào)試通過(guò)。利用教材《電力系統(tǒng)暫態(tài)分析》例題3-4測(cè)試程序是否正確，并與數(shù)個(gè)標(biāo)注測(cè)試系統(tǒng)（IEEE節(jié)點(diǎn)系統(tǒng)）解出算例結(jié)果并分析。第二章基礎(chǔ)理論與原理2.1對(duì)稱短路計(jì)算的基本方法在正常運(yùn)行情況下，三相電力系統(tǒng)可以認(rèn)為是三相對(duì)稱的，只要知道某一相的值就可以推算出其他兩相的值，所以可以用單相電路來(lái)求解三相電力系統(tǒng)，即將所有Δ形接法的三相電力系統(tǒng)元件變換成等值的Υ形接法后，求其一臂的單相值。對(duì)稱短路的一種計(jì)算方法就是在電力系統(tǒng)三相的同一點(diǎn)D上直接短接或經(jīng)一很小的阻抗（如電弧電阻）三相短接（圖3-1（a）、（b））。在單相電路圖上就是在短路點(diǎn)D直接接地或經(jīng)一電阻R接地（圖3-1(c)、（d））。但這種短路故障并不破壞整個(gè)系統(tǒng)的對(duì)稱性，所以對(duì)稱短路的計(jì)算與一般的對(duì)稱三相電路的計(jì)算沒(méi)有什么本質(zhì)差別。但是，在計(jì)算短路電流時(shí)，應(yīng)計(jì)及由于短路而發(fā)生系統(tǒng)參數(shù)突變時(shí)，因磁鏈?zhǔn)睾愕脑?，發(fā)電機(jī)的次暫態(tài)電勢(shì)（或在無(wú)阻尼繞組發(fā)電機(jī)中的暫態(tài)電勢(shì)'）保持不變。這個(gè)在突變前后保持不變的虛擬電勢(shì)，相當(dāng)于在各發(fā)電機(jī)次暫態(tài)電抗Xd"（或暫態(tài)電抗Xd'）后的電勢(shì)，可以根據(jù)電力系統(tǒng)正常運(yùn)行方式的計(jì)算（潮流計(jì)算）結(jié)果，按各發(fā)電機(jī)的負(fù)荷電流和端電壓算出。然后，利用這電勢(shì)進(jìn)行t=0秒時(shí)的短路電流計(jì)算。在對(duì)稱短路計(jì)算時(shí)，上述次暫態(tài)電勢(shì)經(jīng)過(guò)發(fā)電機(jī)阻抗加于電網(wǎng)中響應(yīng)的發(fā)電機(jī)母線節(jié)點(diǎn)i。但是，用這種方法將使一臺(tái)發(fā)電機(jī)增加一個(gè)節(jié)點(diǎn)。由于一般用節(jié)點(diǎn)注入電流作為已知條件，所以可將電壓源等值電路用電流源等值電路來(lái)代替。這樣發(fā)電機(jī)電勢(shì)的作用就相當(dāng)于在其母線的節(jié)點(diǎn)上注入一個(gè)數(shù)值等于的電流，同時(shí)在該節(jié)點(diǎn)上對(duì)地并聯(lián)一個(gè)阻抗，其數(shù)值等于發(fā)電機(jī)的阻抗這樣就可使網(wǎng)絡(luò)節(jié)點(diǎn)數(shù)不改變。在進(jìn)行短路電流計(jì)算時(shí)，一般用恒定阻抗表示各節(jié)點(diǎn)的負(fù)荷，負(fù)荷的阻抗值ZH可用下式求得(2-1)式中為故障前負(fù)荷節(jié)點(diǎn)的電壓；為故障前負(fù)荷電流；PH、QH為負(fù)荷的有功功率和無(wú)功功率。因?yàn)橐话阖?fù)荷電流較短路電流要小得多，所以在簡(jiǎn)化的短路電流計(jì)算中可以忽略負(fù)荷電流的影響，這時(shí)就可以不計(jì)算節(jié)點(diǎn)的負(fù)荷阻抗。另一種計(jì)算方法是應(yīng)用重迭原理，即將故障后的網(wǎng)絡(luò)狀態(tài)看成是兩種情況的迭加，一種是故障前的狀態(tài)，見(jiàn)圖3-2（a），即正常運(yùn)行方式的計(jì)算結(jié)果；還有一種情況是各發(fā)電機(jī)電勢(shì)均等于零，而僅在故障點(diǎn)D加一電勢(shì)，該電勢(shì)值剛好與第一種情況下故障點(diǎn)D的電壓值Ud相等，但極性相反，如圖3-2（b）所示。這樣，當(dāng)兩種情況迭加時(shí)，剛好使故障點(diǎn)D的電壓在單相圖中等于零，即相當(dāng)于三相短路。按第二種情況計(jì)算所得的電流即短路電流中的故障分量，與第一種情況的計(jì)算結(jié)果迭加起來(lái)，就可得到總的短路電流。2.2用節(jié)點(diǎn)阻抗矩陣的計(jì)算方法任一網(wǎng)絡(luò)用節(jié)點(diǎn)阻抗矩陣表示的節(jié)點(diǎn)電壓方程為（2-2）當(dāng)某一節(jié)點(diǎn)發(fā)生三相對(duì)稱故障時(shí)，相當(dāng)于在此接點(diǎn)接上故障阻抗Zf，流過(guò)Zf的故障電流If，其參考方向（流向）地，其他節(jié)點(diǎn)沒(méi)有外接電路，所以其注入電流為零。設(shè)在節(jié)點(diǎn)f發(fā)生三相對(duì)稱故障，故障阻抗為Zf時(shí)，其邊界條件為（2-3）可求得故障點(diǎn)電流為（2-4）各點(diǎn)電壓計(jì)算公式為（2-5）各支路電流計(jì)算公式（2-6）Zij為連接節(jié)點(diǎn)i與j的支路阻抗。在略去輸電線電容的條件下，支路電流也就是輸電線電流。2.3用節(jié)點(diǎn)導(dǎo)納矩陣的計(jì)算方法任一網(wǎng)絡(luò)用節(jié)點(diǎn)導(dǎo)納矩陣表示的網(wǎng)絡(luò)節(jié)點(diǎn)方程為（2-7）節(jié)點(diǎn)導(dǎo)納矩陣ＹB的特點(diǎn)是易于形成，當(dāng)網(wǎng)絡(luò)結(jié)構(gòu)變化時(shí)也容易修改，而且矩陣本身是很稀疏的，但是應(yīng)用它計(jì)算短路電流不如用節(jié)點(diǎn)阻抗矩陣ＺB那樣直接。一般是利用節(jié)點(diǎn)導(dǎo)納矩陣ＹB求逆的方法求得節(jié)點(diǎn)阻抗矩陣ＺB，或者利用三角分解法求得節(jié)點(diǎn)阻抗矩陣ＺB，再根據(jù)用節(jié)點(diǎn)阻抗矩陣的計(jì)算方法求得短路點(diǎn)電流，以及各節(jié)點(diǎn)電壓和線路電流。應(yīng)用節(jié)點(diǎn)導(dǎo)納矩陣計(jì)算短路電流，實(shí)質(zhì)上是先用它計(jì)算與短路點(diǎn)f相關(guān)的節(jié)點(diǎn)阻抗矩陣的第f列元素：Z1fZffZnf，然后根據(jù)式（2-4）~式（2-6）進(jìn)行短路電流的相關(guān)計(jì)算。Z1f~Znf是在f點(diǎn)通以單位電流（其他節(jié)點(diǎn)電流均為零）是1-n點(diǎn)的電壓，故可用式（2-7）求解下列方程(2-8)求得的~即為Z1f~Znf。2.4用三角分解法求解節(jié)點(diǎn)阻抗矩陣節(jié)點(diǎn)導(dǎo)納矩陣的因子表可以采用三角分解法來(lái)形成。將方程（2-8）簡(jiǎn)寫(xiě)為：YU=I(2-9)將Y按照矩陣的三角分解法，把Y分解為：Y=LDR(2-10)式中，D為對(duì)角陣，L為下三角陣，R為上三角陣，L和R互為轉(zhuǎn)置陣。節(jié)點(diǎn)方程變?yōu)椋篖DRU=I(2-11)這個(gè)方程又可以分為三個(gè)方程：(2-12)依次求解這三個(gè)方程就能得到U1~Uf，即一列阻抗元素值。2.5短路發(fā)生在線路上任意處的計(jì)算方法若短路不是發(fā)生在網(wǎng)絡(luò)原有節(jié)點(diǎn)上，而是如圖2-3所示，發(fā)生在線路的任意點(diǎn)上，則網(wǎng)絡(luò)增加了一個(gè)節(jié)點(diǎn)，其阻抗矩陣（和導(dǎo)納矩陣）增加了一階，即與f點(diǎn)有關(guān)的一列和一行元素。顯然，采取重新形成網(wǎng)絡(luò)矩陣的方法是不可取的，以下將介紹利用原網(wǎng)絡(luò)阻抗矩陣中j和k兩列元素直接計(jì)算與f點(diǎn)有關(guān)的一列阻抗元素（Z1f,，Zif,，Zff）的方法。（1）互阻抗Zfi（=Zif）根據(jù)節(jié)點(diǎn)阻抗矩陣元素的物理意義，當(dāng)網(wǎng)絡(luò)中任意節(jié)點(diǎn)i注入單位電流，而其余節(jié)點(diǎn)注入電流均為零時(shí)，f點(diǎn)的對(duì)地電壓即為Zfi，故：(2-13)式中，Zji和Zki為已知的原網(wǎng)絡(luò)的j、k對(duì)i的互阻抗元素。（2）自阻抗Zff當(dāng)f點(diǎn)注入單位電流時(shí)，f點(diǎn)的對(duì)地電壓即為Zff。則有：(2-14)將電壓用相應(yīng)的阻抗元素表示，則得：(2-15)化簡(jiǎn)后得：(2-16)式中Zjf和Zkf用式（4-13）代入，則：(2-17)式中Zjj、Zkk、Zjk和zjk均為已知。由上述算式即可求得f列的阻抗元素,從而可用式（2-4）~式（2-6）作短路電流的相關(guān)計(jì)算。第三章程序設(shè)計(jì)3.1變量說(shuō)明structdata_type{intnode1;//首節(jié)點(diǎn)號(hào)intnode2;//末節(jié)點(diǎn)號(hào)floatx;//電抗值}data[SIZE];//結(jié)構(gòu)體變量，用于導(dǎo)入數(shù)據(jù) inti,j,k;//循環(huán)變量 intchoose;//選擇短路方式變量intf,p,q;//短路位置變量 floatl;//短路點(diǎn)在線路上的位置變量floatIf;//短路電流 floatx[n][n]={0};//原始電抗矩陣 floaty[n][n]={0};//網(wǎng)絡(luò)節(jié)點(diǎn)導(dǎo)納矩陣 floatR[n][n]={0};//R矩陣 floatD[n][n]={0};//D矩陣 floatL[n][n]={0};//L矩陣 floatI[n][n]={0};//I矩陣 floatW[n][n]={0};//W矩陣 floatX[n][n]={0}; //X矩陣 floatU[n][n]={0};//U矩陣，即網(wǎng)絡(luò)節(jié)點(diǎn)阻抗矩陣floatUd[n]={0};//各點(diǎn)電壓矩陣floatIL[n][n]={0};//線路電流矩陣3.2程序流程圖3.2.1主程序流程圖3.2.2導(dǎo)納矩陣流程圖3.2.3三角分解法流程圖3.3程序源代碼見(jiàn)附錄1第四章結(jié)果分析測(cè)試系統(tǒng)（例題3-4）的結(jié)果分析測(cè)試系統(tǒng)數(shù)據(jù)與系統(tǒng)圖見(jiàn)附錄3例題計(jì)算結(jié)果節(jié)點(diǎn)導(dǎo)納矩陣節(jié)點(diǎn)3短路的短路電流各節(jié)點(diǎn)電壓任意線路的電流測(cè)試系統(tǒng)的運(yùn)行結(jié)果見(jiàn)附錄4用計(jì)算機(jī)計(jì)算的結(jié)果和例題的結(jié)果一致，說(shuō)明程序在運(yùn)行測(cè)試系統(tǒng)后結(jié)果正確，適用于少節(jié)點(diǎn)系統(tǒng)。第五章收獲與建議兩周的電力系統(tǒng)課程設(shè)計(jì)，使我對(duì)電力系統(tǒng)發(fā)生短路的基礎(chǔ)知識(shí)進(jìn)行了復(fù)習(xí)與鞏固，而且對(duì)于利用計(jì)算機(jī)計(jì)算復(fù)雜系統(tǒng)短路電流的原理計(jì)算有了一定了解。并親自利用其中的用節(jié)點(diǎn)導(dǎo)納矩陣的計(jì)算方法進(jìn)行了編程，并用教材上的例題進(jìn)行了驗(yàn)證。在編程的過(guò)程中遇到了很多問(wèn)題，主要可分為兩方面，理論不清楚的問(wèn)題和編程方面的問(wèn)題。因?yàn)槭巧蠈W(xué)期上的課，所以有些理論不太清楚了，不懂的時(shí)候我就認(rèn)真看書(shū)，看例題，實(shí)在不懂了就請(qǐng)教老師，只有理論弄清楚了，編程時(shí)思路才不會(huì)錯(cuò)。對(duì)于編程方面的問(wèn)題大多也是翻書(shū)，看例題，查資料等，慢慢的就進(jìn)入狀態(tài)了。在編寫(xiě)導(dǎo)納矩陣的時(shí)候，節(jié)點(diǎn)導(dǎo)納矩陣是在電路中學(xué)過(guò)的，我就去翻電路課件，找到了形成節(jié)點(diǎn)導(dǎo)納矩陣的方法：自導(dǎo)納等于與該節(jié)點(diǎn)相關(guān)的所有導(dǎo)納的和，互導(dǎo)納就等于兩節(jié)點(diǎn)之間的負(fù)導(dǎo)納。有了理論編程也就不難了。另外，編寫(xiě)程序時(shí)還應(yīng)該要注意，一個(gè)模塊一個(gè)模塊的編，編好后調(diào)試成功了再編下一個(gè)模塊，這樣就很容易知道錯(cuò)誤出在哪兒，修改起來(lái)也會(huì)簡(jiǎn)單很多。剛開(kāi)始我是采用利用矩陣求逆的方法，由節(jié)點(diǎn)導(dǎo)納矩陣求節(jié)點(diǎn)阻抗矩陣的，但是編程序時(shí)發(fā)現(xiàn)會(huì)用到指針，而我又不會(huì)，編了好長(zhǎng)時(shí)間也沒(méi)編好，最后請(qǐng)教老師，老師說(shuō)可以用三角分解法的方法求節(jié)點(diǎn)阻抗矩陣，求得Y的因子表R和D后，可以采用解方程的方法解出節(jié)點(diǎn)阻抗矩陣。因?yàn)樾纬傻腞、D、L都是特殊矩陣，因此用解方程的方法計(jì)算起來(lái)比用求逆的方法要簡(jiǎn)單得多。雖然中間出現(xiàn)了許多問(wèn)題，但是最后，經(jīng)過(guò)一個(gè)星期的努力，程序終于編好了，經(jīng)過(guò)驗(yàn)證也證明程序正確，但是還是有一些缺憾，我自己對(duì)于編寫(xiě)子函數(shù)還是不在行，簡(jiǎn)單的沒(méi)有返回的函數(shù)可以編寫(xiě)，但是有了返回值，特別返回值不是一個(gè)時(shí)，我就不知道該怎么編了。在編寫(xiě)子函數(shù)和學(xué)會(huì)用指針這兩方面我還須繼續(xù)努力。實(shí)際中的電力系統(tǒng)一般都比較復(fù)雜，人工計(jì)算不僅浪費(fèi)時(shí)間，可能還會(huì)出現(xiàn)錯(cuò)誤，而利用計(jì)算機(jī)進(jìn)行計(jì)算，不僅速度快，而且準(zhǔn)確率也高，因此，掌握計(jì)算機(jī)計(jì)算復(fù)雜系統(tǒng)短路電流的原理對(duì)于以后從事相關(guān)電力專業(yè)的工作有很大幫助。作為為以后工作打基礎(chǔ)來(lái)說(shuō)，我認(rèn)為這兩周的課設(shè)很有意義，不過(guò)，我認(rèn)為還是側(cè)重編程多一些?？傊?，這一次的課程設(shè)計(jì)我復(fù)習(xí)了課本上的知識(shí)，還學(xué)會(huì)了一些編程的技巧，收獲很多。參考文獻(xiàn)[1]陳珩編.電力系統(tǒng)穩(wěn)態(tài)分析（第三版）[M].北京:中國(guó)電力出版社，2008年[2]何仰贊主編.電力系統(tǒng)分析(第三版)[M].武漢:華中科技大學(xué)出版社，2005年.[3]夏道止主編.電力系統(tǒng)分析[M].北京:中國(guó)電力出版社，2004年[4]陳怡主編.電力系統(tǒng)分析[M].北京:中國(guó)電力出版社，2005年[5]譚浩強(qiáng)主編.C語(yǔ)言程序設(shè)計(jì).清華大學(xué)出版社，2000年[6]譚浩強(qiáng)主編.C語(yǔ)言程序設(shè)計(jì)題解與上機(jī)指導(dǎo).清華大學(xué)出版社，2000年附錄附錄1：程序源代碼三相對(duì)稱短路計(jì)算程序源代碼#include<stdio.h>#defineSIZE5#definen3//節(jié)點(diǎn)數(shù)voidshuchu(floatarray[n]);//輸出一維數(shù)組voidshuchu1(floatarray[n][n]);//輸出二維虛數(shù)數(shù)組voidshuchu2(floatarray[n][n]);//輸出二維實(shí)數(shù)數(shù)組structdata_type{intnode1;intnode2;floatx;}data[SIZE];voidmain(){ inti,j,k;intchoose;intf,p,q; floatl,If; floatx[n][n]={0};//原始電抗矩陣 floaty[n][n]={0};//網(wǎng)絡(luò)節(jié)點(diǎn)導(dǎo)納矩陣矩陣 floats[n][n]={0}; floatd[n][n]={0}; floatr[n][n]={0}; floatR[n][n]={0}; floatD[n][n]={0}; floatL[n][n]={0}; floatI[n][n]={0}; floatW[n][n]={0};floatw1[n][n]={0}; floatX[n][n]={0}; floatU[n][n]={0}; floatu1[n][n]={0}; floatUz[n+1]={0};floatUd[n]={0};//各點(diǎn)電壓矩陣floatIL[n][n]={0};//線路電流矩陣floatIL2[n+1][n+1]={0}; FILE*fp; fp=fopen("data.txt","r"); printf("首節(jié)點(diǎn)末節(jié)點(diǎn)電抗值\n"); for(i=0;i<SIZE;i++) { fscanf(fp,"%d,%d,%fi\n",&data[i].node1,&data[i].node2,&data[i].x); printf("%d\t%d\t%fi\n",data[i].node1,data[i].node2,data[i].x);} fclose(fp); for(k=0;k<SIZE;k++) { for(i=0;i<n;i++) for(j=0;j<n;j++) { if(i==(data[k].node1-1)&&j==(data[k].node2-1)) x[i][j]=data[k].x; x[j][i]=x[i][j];} } printf("**********系統(tǒng)原始電抗參數(shù)矩陣x=**********\n");shuchu1(x);//求網(wǎng)絡(luò)節(jié)點(diǎn)導(dǎo)納矩陣 for(i=0;i<n;i++) { for(j=0;j<n;j++) { if(x[i][j]!=0) y[i][i]=y[i][i]-1/x[i][j];} } for(i=0;i<n;i++) { for(j=i+1;j<n;j++) { if(x[i][j]==0) y[i][j]=0; else y[i][j]=1/x[i][j]; y[j][i]=y[i][j];} } printf("**********系統(tǒng)節(jié)點(diǎn)導(dǎo)納矩陣y=**********\n"); shuchu1(y);//求R和D的元素 for(i=0;i<n;i++) {for(k=0;k<i;k++) { s[i][i]=s[i][i]+R[k][i]*R[k][i]*d[k][k]; } d[i][i]=y[i][i]-s[i][i];D[i][i]=-1/d[i][i];//求D矩陣的元素 for(j=i;j<n;j++) {for(k=0;k<i;k++) { r[i][j]=r[i][j]+R[k][i]*R[k][j]*d[k][k]; } R[i][i]=1; R[i][j]=(y[i][j]-r[i][j])/d[i][i];//求R矩陣的元素 }for(j=0;j<n;j++) L[i][j]=R[j][i];//求L矩陣的元素 }/*printf("矩陣D=\n");shuchu2(D);printf("矩陣R=\n");shuchu2(R);printf("矩陣L=\n");shuchu2(L);*///三角分解法求節(jié)點(diǎn)阻抗矩陣 for(i=0;i<n;i++) {for(j=0;j<n;j++) I[i][i]=1; }//printf("矩陣I=\n");//shuchu2(I);//W for(k=0;k<n;k++) for(i=1;i<n;i++) {for(j=0;j<i;j++) {W[0][0]=I[0][0]; w1[i][k]=w1[i][k]+L[i][j]*W[j][k]; } W[i][k]=I[i][k]-w1[i][k]; } //printf("矩陣W=\n");//shuchu2(W);//X for(k=0;k<n;k++) for(i=0;i<n;i++) X[i][k]=D[i][i]*W[i][k];//printf("矩陣X=\n");//shuchu1(X);//U for(k=0;k<n;k++) for(i=n-2;i>=0;i--) {for(j=n-1;j>i;j--) {U[n-1][k]=X[n-1][k]; u1[i][k]=u1[i][k]+R[i][j]*U[j][k]; } U[i][k]=X[i][k]-u1[i][k]; } printf("**********網(wǎng)絡(luò)節(jié)點(diǎn)阻抗矩陣U=**********\n"); shuchu1(U); printf("請(qǐng)輸入短路發(fā)生位置(1-節(jié)點(diǎn)，2-線路):"); scanf("%d",&choose); if(choose==1) { printf("請(qǐng)輸入短路點(diǎn)："); scanf("%d",&f);//求短路點(diǎn)電流 If=-1/U[f-1][f-1]; printf("結(jié)點(diǎn)%d短路電流If=%fi\n",f,If);//求各節(jié)點(diǎn)電壓 for(i=0;i<n;i++) Ud[i]=1+If*U[i][f-1]; printf("**********節(jié)點(diǎn)電壓矩陣Ud=**********\n"); shuchu(Ud);//求線路電流 for(i=0;i<n;i++) { for(j=0;j<n;j++) IL[i][j]=(Ud[j]-Ud[i])*y[i][j]; } printf("**********線路電流矩陣IL=**********\n"); shuchu1(IL); } else { printf("請(qǐng)輸入短路發(fā)生所在線路小節(jié)點(diǎn)號(hào):"); scanf("%d",&p); printf("請(qǐng)輸入短路發(fā)生所在線路大節(jié)點(diǎn)號(hào):"); scanf("%d",&q); printf("請(qǐng)輸入短路發(fā)生位置(用小數(shù)表示):"); scanf("%f",&l);//求短路點(diǎn)相關(guān)的一列阻抗元素 for(i=0;i<n;i++) Uz[i]=(1-l)*U[i][p-1]+l*U[i][q-1]; Uz[n]=Uz[n]+(1-l)*(1-l)*U[p-1][p-1]+l*l*U[q-1][q-1]+2*l*(1-l)*U[p-1][q-1]+l*(1-l)*x[p-1][q-1]; printf("******跟短路點(diǎn)相關(guān)的一列阻抗元素為：****\n"); for(i=0;i<=n;i++) printf("\t%f\n",Uz[i]);//求短路點(diǎn)電流 If=-1/Uz[n]; printf("線路發(fā)生短路的電流If=%fi\n",If);//求各節(jié)點(diǎn)電壓 for(i=0;i<=n;i++) Ud[i]=1+If*Uz[i]; printf("**********節(jié)點(diǎn)電壓矩陣Ud=**********\n"); for(i=0;i<=n;i++) printf("\t%f\n",Ud[i]);//求線路電流 for(i=0;i<n;i++) {for(j=0;j<n;j++) IL2[i][j]=(Ud[j]-Ud[i])*y[i][j]; }IL2[p-1][q-1]=0; IL2[q-1][p-1]=0;IL2[n][p-1]=(Ud[p-1]-Ud[n])/(l*x[p-1][q-1]); IL2[p-1][n]=-IL2[n][p-1]; IL2[n][q-1]=Ud[q-1]/((1-l)*x[p-1][q-1]); IL2[q-1][n]=-IL2[n][q-1]; printf("**********線路電流矩陣IL=**********\n"); for(i=0;i<=n;i++) {for(j=0;j<=n;j++) printf("%fi\t",IL2[i][j]); printf("\n"); } }}/**********輸出一維數(shù)組***********/voidshuchu(floatarray[n]){ inti; for(i=0;i<n;i++) printf("\t%f\n",array[i]);}/********輸出二維虛數(shù)數(shù)組*********/voidshuchu1(floatarray[n][n]){ inti,j; for(i=0;i<n;i++) {for(j=0;j<n;j++) printf("%fi\t",array[i][j]);printf("\n"); }}/********輸出二維實(shí)數(shù)數(shù)組*********/voidshuchu2(floatarray[n][n]){ inti,j; for(i=0;i<n;i++) {for(j=0;j<n;j++) printf("%f\t",array[i][j]);printf("\n"); }}附錄2：測(cè)試系統(tǒng)數(shù)據(jù)與系統(tǒng)圖測(cè)試系統(tǒng)數(shù)據(jù)：線路編號(hào)首節(jié)點(diǎn)號(hào)末節(jié)點(diǎn)號(hào)線路參數(shù)（電抗）111j0.1500212j0.1000313j0.1000422j0.0750523j0.10007節(jié)點(diǎn)網(wǎng)絡(luò)系統(tǒng)數(shù)據(jù)：線路編號(hào)首節(jié)點(diǎn)號(hào)末節(jié)點(diǎn)號(hào)線路參數(shù)(電抗)111j0.2212j0.3323j0.15424j0.2533j0.15634j0.05744j0.4845j0.35946j0.351055j0.51157j0.451266j0.21367j0.0851477j0.2附錄3：測(cè)試系統(tǒng)的運(yùn)行結(jié)果測(cè)試系統(tǒng)結(jié)果輸出網(wǎng)絡(luò)節(jié)點(diǎn)導(dǎo)納矩陣Y=0-26.6667i0+10.0000i0+10.0000i0+10.0000i0-33.3333i0+10.0000i0+10.0000i0+10.0000i0-20.0000i1.請(qǐng)輸入短路發(fā)生位置：1-節(jié)點(diǎn)，2-線路。1請(qǐng)輸入短路點(diǎn)：3If=0-9.8592iV=0.45070.53520.0000I=00+0.8451i0-4.5070i0-0.8451i00-5.3521i0+4.5070i0+5.3521i02.請(qǐng)輸入短路發(fā)生位置：1-節(jié)點(diǎn)，2-線路。2請(qǐng)輸入短路所在線路小節(jié)點(diǎn)號(hào)：2請(qǐng)輸入短路所在線路大節(jié)點(diǎn)號(hào)：3請(qǐng)輸入短路位置（百分?jǐn)?shù)或小數(shù)）：0.3Z1=0+0.0437i0+0.0531i0+0.0634i0+0.0772i短路電流If=0-12.9486i短路電壓V=0.43400.31190.17870線路電流I=00-1.2209i0-2.5527i00+1.2209i000-10.3959i0+2.5527i000-2.5527i00+10.3959i0+2.5527i03.請(qǐng)輸入短路發(fā)生位置：1-節(jié)點(diǎn)，2-線路。2請(qǐng)輸入短路所在線路小節(jié)點(diǎn)號(hào)：1請(qǐng)輸入短路所在線路大節(jié)點(diǎn)號(hào)：3請(qǐng)輸入短路位置（百分?jǐn)?shù)或小數(shù)）：0.2Z1=0+0.