電力系統(tǒng)分析課設 C語言

1、目錄1課程設計說明12編程語言的選擇及理由13程序設計23.1主體流程圖23.2詳細流程圖33.2.1創(chuàng)建系統(tǒng)（create）流程圖33.2.2加載系統(tǒng)函數(shù)（load）流程圖33.2.3計算子函數(shù)（calculate）流程圖43.2.4改變短路節(jié)點（change）流程圖43.3數(shù)據(jù)及變量說明53.4程序代碼及注釋63.5測試算例64設計體會13參考文獻13附錄131 課程設計說明本課程設計根據(jù)電力系統(tǒng)分析課程中所講的基于節(jié)點方程的三相短路計算的原理和方法。采用c語言編程 。實現(xiàn)了三相短路計算的計算機實現(xiàn)。本程序輸入含有網(wǎng)絡拓撲信息和電器元件的參數(shù)，輸出三相短路的短路電流和短路后網(wǎng)絡各節(jié)點的電壓

2、參數(shù)和各支路的電流參數(shù)。并以文件的形式保存，方便用戶的讀取，打印。實現(xiàn)了設計的要求。三相短路在電力系統(tǒng)中雖然發(fā)生的概率很小，但一旦發(fā)生產(chǎn)生的影響非常嚴重，電力系統(tǒng)中設備的動穩(wěn)定和熱穩(wěn)定校驗都要以三相短路電流為依據(jù)，因此我選擇做是電力系統(tǒng)三相短路的計算程序，覺得更有意義。應用計算機進行電力系統(tǒng)計算，首先要建立電力系統(tǒng)相應的數(shù)學模型，如建立電力系統(tǒng)等效網(wǎng)絡和節(jié)點方程；其次是運用合理的計算方法計算結果，如三角分解法，LDU分解法等；第三則是選擇合適的計算機語言編制計算程序，實現(xiàn)程序的準確性，提高精確度。本設計利用節(jié)點方程來實現(xiàn)電力系統(tǒng)三相短路計算。首先，根據(jù)給定的電力系統(tǒng)運行方式制定系統(tǒng)的一相等值網(wǎng)

3、絡，計算出各元器件的參數(shù)，其次根據(jù)網(wǎng)絡方程利用變壓器和線路參數(shù)形成不含發(fā)電機和負荷的節(jié)點導納矩陣Y；然后，利用系統(tǒng)形成的節(jié)點導納矩陣求解系統(tǒng)阻抗矩陣Z；第三則是應用節(jié)點阻抗矩陣計算短路電流。電力系統(tǒng)短路電流的實際工程計算中，許多實際問題的解決并不需要十分精確的結果。因此采用近似計算法，在建立系統(tǒng)節(jié)點方程時，采用網(wǎng)絡模型和參數(shù)做簡化處理，等值電路的制定通常將發(fā)電機作為含源支路，表示為節(jié)點注入電流源和節(jié)點電抗的形式；系統(tǒng)中存在同步發(fā)電機時，按發(fā)電機支路處理；同時，忽略發(fā)電機、變壓器和輸電線路的電阻，不計線路的電容，略去變壓器的勵磁電流，負荷忽略不計；在標幺參數(shù)計算時，選取各級平均額定電壓作為基準電

4、壓，忽略元件的額定電壓和平均電壓的差別，認為變壓器變比的標幺值為1；此外，還假定所有的發(fā)電機的電勢具有相同的相位，所有的元件僅用電抗表示，避免復數(shù)運算。三相對稱短路為金屬性短路。本課程設計的關鍵部分是利用系統(tǒng)節(jié)點導納矩陣形成節(jié)點阻抗矩陣，本程序采用LDU分解法，然后利用節(jié)點導納矩陣和阻抗矩陣的乘積為單位矩陣列寫方程，節(jié)點導納矩陣經(jīng)LDU分解后，回代求解系統(tǒng)阻抗矩陣。利用節(jié)點阻抗矩陣進行短路計算，首先利用故障前的節(jié)點注入電流和系統(tǒng)節(jié)點阻抗矩陣求取故障前系統(tǒng)的各節(jié)點電壓；然后，利用金屬短路時短路點殘壓為0,對地電阻為0，計算短路點電流；最后，利用故障前節(jié)點電壓、節(jié)點阻抗和短路點電流求得短路后各節(jié)點

5、電壓和各支路電流。2 編程語言的選擇及理由能夠?qū)崿F(xiàn)程序編寫的語言有多種，本課程設計中采用C語言進行電力系統(tǒng)對稱短路計算程序的編寫。個人對c語言和metlab相對較熟悉。之所以選擇C語言編程，從主觀上說，是因為我學過C語言課程，相對于其他的編程語言而言，C語言的應用更為熟悉和了解。另外，從客觀上說，C語言是目前世界上最為流行、使用最為廣泛的高級程序設計語言。，C的效率MATLAB比高多了，MATLAB雖然語法簡單，但執(zhí)行效率不高。MATLAB是用于特定的方面的，比如說矩陣運算方面、DSP就很有優(yōu)勢，而C語言不同。很多軟件的底層都是C編寫的。MATLAB提供的是一種基于解釋的語言，雖然也是一種高級

6、語言，但是還是很簡單的，它的目標在于用戶方便，開發(fā)簡單。而實際上這些方便也是有代價的，編譯器承受了這些代價，程序執(zhí)行效率方面也有損失。另一方面MATLAB這么做也是基于其專用于科學計算的目的來考慮的，無法應用于通用的程序設計。c語言這么做保持了其編譯的速度，便于大型程序設計，執(zhí)行效率也很高，因此我選擇用c語言編寫程序 。3 程序設計3.1 主體流程圖開始歡迎使用用戶選擇操作，按鍵123進入創(chuàng)建系統(tǒng)子函數(shù)（create）進入加載系統(tǒng)子函數(shù)（load）進入退出子函數(shù)進入計算子函數(shù)（calculate）改短路點子函數(shù)（change）是否改路點？結束是圖3-1 主體流程圖3.2 詳細流程圖3.2.1

7、創(chuàng)建系統(tǒng)（create）流程圖輸入文件名輸入電抗支路參數(shù)結束輸入節(jié)點數(shù)、支路數(shù)、故障節(jié)點輸入發(fā)電機支路參數(shù)輸入電容支路參數(shù)圖3-2 create函數(shù)流程圖3.2.2 加載系統(tǒng)函數(shù)（load）流程圖輸入系統(tǒng)文件名輸出查看系統(tǒng)參數(shù)選擇進入下一個子函數(shù)圖3-3 節(jié)點導納矩陣形成流程圖3.2.3 計算子函數(shù)（calculate）流程圖 文件讀入電路參數(shù)求取LDU分解的因子表矩陣下三角求取DUZj向量求解完畢？否是結束求取節(jié)點阻抗矩陣列向量根據(jù)支路參數(shù)建立節(jié)點導納矩陣求取短路點電流求取節(jié)點電壓求取支路電流建立文件存儲計算結果結束圖3-4 計算子函數(shù)流程圖3.2.4 改變短路節(jié)點（change）流程圖輸入

8、短路點從新讀取文件中的參數(shù)選擇進入下一個子函數(shù)計算圖3-5 change函數(shù)流程圖3.3 數(shù)據(jù)及變量說明本程序中將電路的支路分為三類：電抗支路、電流支路、發(fā)電機直流，數(shù)據(jù)輸入時分開輸入格式如下：電抗支路：支路數(shù)、節(jié)點首端好、節(jié)點末端好，支路電抗值。發(fā)電機支路：發(fā)電機臺數(shù)、發(fā)電機接入點，發(fā)電機電勢、發(fā)電機之路電抗。電容支路：支路數(shù)、電容節(jié)點號，電容支路電抗值。程序中的變量定義如下： FILE *fp; /文件指針 float If; /短路電流 float V010; /節(jié)點電壓初值 float In10; /注入電流 float V10; /節(jié)點電壓 float I10; /支路電流int c

9、ount210; /數(shù)組用于儲存支路號float z1010; /儲存支路電抗值float yjs1010; /存放計算導納矩陣中的中間值float ye10; /儲存發(fā)電機支路導納float Y1010; /節(jié)點導納矩陣float Z1010; /阻抗矩陣float dk; /電抗值int e110;float yc10; /節(jié)點對地電抗float ycjs10;int e10; /發(fā)電機接入點float E10; float Ig10;float Emid10;float ze10; /儲存發(fā)電機支路導納int m,n,a,f,i=0,j=0,sd,md,butn,times,timesc

10、;/*n是節(jié)點數(shù)，a是支路數(shù)，f是短路點，i為行標，j為列標*/1、文件讀入程序如下： printf(請輸入系統(tǒng)的節(jié)點數(shù)，支路數(shù)，故障節(jié)點n); scanf(%d%d%d,&n,&a,&f); fp=fopen(name,w+); fprintf(fp,系統(tǒng)參數(shù)如下所示:nn); fprintf(fp,1.系統(tǒng)的節(jié)點數(shù)為%dn2.系統(tǒng)的支路數(shù)為%dn3.系統(tǒng)的短路點為%dnn,n,a,f); fprintf(fp,支路電抗參數(shù)如下n);2、計算主程序如下： （1）、 If=V0f-1/Zf-1f-1;（短路電流計算） for(i=0;in;i+) Vi=V0i-Zif-1*If; （2）、fo

11、r(i=0;in;i+)（導納矩陣計算） for(j=0;jn;j+) Yii=Yii-yjsij; Yii=Yii-yei-ycjsi; for(i=0;in;i+) for(j=i+1;jn;j+) Yij=yjsij; Yji=yjsij; 3.4 程序代碼及注釋 具體源程序代碼及其注釋見附錄。3.5 測試算例【例6-3】在如圖3-5所示的電力系統(tǒng)中分別在節(jié)點1和節(jié)點5接入發(fā)電機支路，其標幺值參數(shù)為：。在節(jié)點3發(fā)生三相短路，計算短路電流及網(wǎng)絡中的電流分布。線路的電阻和電容略去不計，變壓器的標幺變比等于1。各元件參數(shù)的標幺值如下：圖3-6 電力系統(tǒng)等值網(wǎng)絡圖圖3-7 三相短路時的等值網(wǎng)絡圖

12、由以上網(wǎng)絡圖可得計算過程如下：（1）、形成輸入文件：系統(tǒng)參數(shù)如下所示:1.系統(tǒng)的節(jié)點數(shù)為52.系統(tǒng)的支路數(shù)為53.系統(tǒng)的短路點為3支路電抗參數(shù)如下z12=0.1050z23=0.0650z24=0.0800z43=0.0500z54=0.1840發(fā)電機接入點及其電抗如下：共2臺發(fā)電機進入系統(tǒng)z1=0.1500z5=0.2200e1=1.0000e5=1.0000共0個節(jié)點有接地電容（2）、計算后得到的輸出文件：計算結果如下：1、節(jié)點導納矩陣:-16.1905 9.5238 9.5238 -37.4084 15.3846 12.5000 15.3846 -35.3846 20.0000 12.5

13、000 20.0000 -37.9348 5.4348 5.4348 -9.9802 2、節(jié)點阻抗矩陣: 0.1181 0.0958 0.0902 0.0858 0.0467 0.0958 0.1629 0.1533 0.1459 0.0794 0.0902 0.1533 0.1860 0.1611 0.0877 0.0858 0.1459 0.1611 0.1729 0.0941 0.0467 0.0794 0.0877 0.0941 0.1515 3、短路節(jié)點是:34、短路電流:If=5.37675、各節(jié)點電壓:V1=0.5152V2=0.1758V3=0V4=0.1336V5=0.528

14、26、發(fā)電機支路電流如下:Ig10=-3.2321Ig50=-2.14457、各支路電流:I12=-3.2321I23=-2.7046I24=-0.5275I43=-2.6720I54=-2.1445圖3-8 對稱短路計算結果圖3-8 改變短路點后計算結果【習題6-10】如圖3-9所示網(wǎng)絡，略去負荷，試用節(jié)點阻抗矩陣法求節(jié)點5發(fā)生三相短路時，短路點的短路電流及線路L-2、L-3的電流。已知各元件參數(shù)如下：發(fā)電機G-1: SN = 120 MVA, x”d = 0.23; G-2: SN = 60 MVA, x”d = 0.14變壓器T-1: SN = 120 MVA, Vs = 10.5%;

15、T-2: SN = 60 MVA, Vs = 10.5%線路參數(shù) x1 = 0.4/km, b1 = 2.810-6 S/km線路長度 L-1: 120 km, L-2: 80 km, L-3: 70 km取SB = 120 MVA, VB = VavGGG-1G-2T-1T-2L-1L-2L-312345圖3-9 電力系統(tǒng)網(wǎng)絡圖圖3-10 三相短路時的等值網(wǎng)絡圖（1）、形成輸入文件：系統(tǒng)參數(shù)如下所示:1.系統(tǒng)的節(jié)點數(shù)為：52.系統(tǒng)的支路數(shù)為：53.系統(tǒng)的短路點為：5支路電抗參數(shù)如下：z3 5=0.29036 z4 5=0.25406 z3 4=0.43554z1 3=0.1050z2 4=0

16、.2100發(fā)電機接入點及其電抗如下：共2臺發(fā)電機進入系統(tǒng)：z1=0.2300z2=0.2800e1=1.0500e2=1.0500共3個節(jié)點有接地電容：Y3=-32.Y4=-34.Y5=-43.（2）、得到的輸出文件：計算結果如下：1、節(jié)點導納矩陣:-13.8716 9.5238 -8.3333 4.7619 9.5238 -15.2329 2.2960 3.4440 4.7619 2.2960 -10.9646 3.9361 3.4440 3.9361 -7.3569 2、節(jié)點阻抗矩陣: 0.1819 0.0616 0.1600 0.1079 0.1326 0.0616 0.2078 0.0

17、898 0.1536 0.1242 0.1600 0.0898 0.2330 0.1571 0.1932 0.1079 0.1536 0.1571 0.2689 0.2174 0.1326 0.1242 0.1932 0.2174 0.3427 3、短路節(jié)點是:54、短路電流：If=3.12635、各節(jié)點電壓:V1=0.6471V2=0.6723V3=0.4632V4=0.3890V5=06、發(fā)電機支路電流如下:Ig10=-1.7516Ig20=-1.34917、各支路電流:I35=-1.5954I45=-1.5310I34=-0.1705I13=-1.7516I24=-1.3491圖3-11

18、 對稱短路計算結果【習題6-11】電力系統(tǒng)等值電路如圖3-13所示，支路阻抗的標幺值已注明圖中。（1）形成節(jié)點導納矩陣（或節(jié)點阻抗矩陣），并用以計算節(jié)點3的三相短路電流。（2）另選一種方法計算短路電流，并用以驗證（1）的計算結果。圖3-12 電力系統(tǒng)等值電路圖圖3-13 經(jīng)Norton變換后的等值電路圖圖3-13 最終計算等值電路圖由以上網(wǎng)絡可得輸入文件如下：系統(tǒng)參數(shù)如下所示:1.系統(tǒng)的節(jié)點數(shù)為：32.系統(tǒng)的支路數(shù)為：23.系統(tǒng)的短路點為：3支路電抗參數(shù)如下：z1 2=0.2000z2 3=0.1000發(fā)電機接入點及其電抗如下：共2臺發(fā)電機進入系統(tǒng)：z1=0.2000z3=0.e1=1.050

19、0e3=0.共1個節(jié)點有接地電容：Y2=-2.（2）、得到的輸出文件：計算結果如下：1、節(jié)點導納矩陣:-10.0000 5.0000 5.0000 -14.5000 10.0000 10.0000 -17.6667 2、節(jié)點阻抗矩陣: 0.1394 0.0789 0.0446 0.0789 0.1577 0.0893 0.0446 0.0893 0.1071 3、短路節(jié)點是:34、短路電流：If=9.18755、各節(jié)點電壓:V1=0.6344V2=0.2188V3=06、發(fā)電機支路電流如下:Ig10=-2.0781Ig30=-7.00007、各支路電流:I12=-2.0781I23=-2.18

20、75圖3-15 對稱短路計算結果通過計算結果與算例答案的對比，計算結果與算例 標準答案一樣，說明了程序計算方法和思路的準確性。說明本次課程設計基本達到了設計的要求。4 設計體會通過本次課程設計讓自己學到了很多東西，包括電力系統(tǒng)分析、和c語言的很知識。通過做課程設計自己又仔細的將兩門課程相關章節(jié)看了一遍，特別是c語言，幾乎把整個課本都看了一遍，并自學了很多相關的知識，重新看書的時候才發(fā)現(xiàn)，雖然以前學過，但幾年不用已經(jīng)忘得差不多啦！當然，重新看書也發(fā)現(xiàn)了許多新的知識，學到不少的東西。本次做課程設計老師給我們留了充足的時間，讓我們耐心去學習知識，而不像以前一樣趕任務。通過這次課程設計，我更好的掌握了

21、電力系統(tǒng)短路計算的方法。我想作為一個學電氣的人這是非常關鍵的。在剛接到題目時，我感覺很迷茫，覺得很難，當我仔細準備的時候發(fā)現(xiàn)這也不是想象中的難，最終通過自己的努力和通過學的幫助完成了課程設計。這次設計給我最大的收獲并不僅僅是相關的專業(yè)知識，這是我第一次做這種要幾門課結合的課程設計，我覺得這使自己應用知識的能力得到了鍛煉，本次課程設計用c語言編程序?qū)崿F(xiàn)了電力系統(tǒng)的三相短路計算。雖然程序能夠計算，但我也知道還有很多的問題。還需要努力。第一次不自己學的知識結合在一起，覺得很有用，自己沒白學。在答辯的過程老師對程序提出了修改意見，是自己的程序更加實用，在這里感謝老師的指導，讓我受益匪淺。同時在做的過程

22、中我也和很多的同學交流、討論，大家交流自己的想法，覺得是很開心的一件事。參考文獻1 何仰贊，溫增銀電力系統(tǒng)分析(上冊)(第三版)武漢：華中科技大學出版社，20022 何仰贊，溫增銀電力系統(tǒng)分析題解(上、下)(第三版)武漢：華中科技大學出版社，20023 鄭平安，曾大亮編程序設計基礎(C語言)第二版北京：清華大學出版社，2006附錄主程序及注釋（這里只給出主要計算程序）：#include #include #include 聲明外部函數(shù) extern void create(); extern void load(); extern void quit();extern void change(

23、char name20); extern void calculate(char name20) FILE *fp; /文件指針 float If; /短路電流 float V010; /節(jié)點電壓初值 float In10; /注入電流 float V10; /節(jié)點電壓 float I10; /支路電流int count210; /數(shù)組用于儲存支路號float z1010; /儲存支路電抗值float yjs1010; /存放計算導納矩陣中的中間值float ye10; /儲存發(fā)電機支路導納float Y1010; /節(jié)點導納矩陣float Z1010; /阻抗矩陣float dk; /電抗值

24、int e110;float yc10; /節(jié)點對地電抗float ycjs10;int e10; /發(fā)電機接入點float E10; float Ig10;float Emid10;float ze10; /儲存發(fā)電機支路導納int m,n,a,f,g,i=0,j=0,sd,md,butn,times,timesc;/*n是節(jié)點數(shù)，a是支路數(shù)，f是短路點，i為行標，j為列標for(i=0;i10;i+)yci=0; ycjsi=0; Ei=0; Emidi=0;for(i=0;i10;i+) /初始賦 V0i=0; for(j=0;j10;j+) zij=0; Yij=0; /*讀取系統(tǒng)參數(shù)

25、*/ fp=fopen(name,r); fscanf(fp,系統(tǒng)參數(shù)如下所示:nn); fscanf(fp,1.系統(tǒng)的節(jié)點數(shù)為%dn2.系統(tǒng)的支路數(shù)為%dn3.系統(tǒng)的短路點為%dnn,&n,&a,&f); /*讀電抗*/ fscanf(fp,支路電抗參數(shù)如下n); for(m=0;ma;m+) fscanf(fp,z%d %d=%ft,&sd,&md,&dk); i=sd-1; j=md-1; count0m=sd; count1m=md; zij=dk; zji=dk; printf(系統(tǒng)支路電抗如下：n); for(m=0;ma;m+) i=count0m-1; j=count1m-1;

26、 printf(z%d%d = %.4fn,count0m,count1m,zij); fscanf(fp,nn發(fā)電機接入點及其電抗如下：); fscanf(fp,共%d臺發(fā)電機進入系統(tǒng)n,×); for(m=0;mtimes;m+) fscanf(fp,z%d=%ft,&em,&zem); fscanf(fp,n); for(m=0;mtimes;m+) fscanf(fp,e%d=%ft,&em,&Em); printf(輸入的發(fā)電機支路電抗為n); for(m=0;mtimes;m+) printf(z%d=%ft,em,zem); for(m=0;mtimes;m+) for

27、(m=0;mtimes;m+) printf(e%d=%ft,em,Em); for(m=0;m0) for(m=0;mtimesc;m+) fscanf(fp,Y%d=%ft,&e1m,&ycm); printf(n系統(tǒng)有%d個節(jié)點電容n,timesc); for(m=0;mtimesc;m+) printf(z%d0=%f ,e1m,ycm); fclose(fp);/*計算節(jié)點阻抗矩陣*/ for(m=0;m10;m+) yem=0; for(m=0;mtimes;m+) for(j=0;jn;j+) if(j=(em-1) yej=1/zem; Emidj=Em; for(m=0;mt

28、imesc;m+) for(j=0;jn;j+) if(j=(e1m-1) ycjsj=1/ycm; /*系統(tǒng)支路導納*/ for(i=0;in;i+) for(j=0;jn;j+) if(zij=0) yjsij=0; else yjsij=1/zij; /*導納矩陣*/ for(i=0;in;i+) for(j=0;jn;j+) Yii=Yii-yjsij; Yii=Yii-yei-ycjsi; for(i=0;in;i+) for(j=i+1;jn;j+) Yij=yjsij; Yji=yjsij; printf(nn節(jié)點導納矩陣為); for(i=0;in;i+) printf(n);

29、 for(j=0;jn;j+) if(Yij=0) printf( ); else printf(%8.4f ,Yij); float l1010; float d1010; float u1010; int kk; float sum; for(i=0;i10;i+) for(j=0;j10;j+) lij=0; dij=0; uij=0; for(i=0;i10;i+) lii=1; uii=1; /*迭代法求解*/ for(i=0;in;i+) sum=0; if(i1) dii=Yii; else for(kk=0;kki;kk+) sum=sum+ukki*ukki*dkkkk; d

30、ii=Yii-sum; for(j=i+1;jn;j+) sum=0; if(i1) uij=Yij/dii; else for(kk=0;kki;kk+) sum=sum+ukki*ukkj*dkkkk; uij=(Yij-sum)/dii; /*LDU屏幕顯示*/ printf(nD矩陣如下所示); for(i=0;in;i+) printf(n); for(j=0;jn;j+) if(dij=0) printf( ); else printf(%8.4f ,dij); printf(nU矩陣如下所示); for(i=0;in;i+) printf(n); for(j=0;jn;j+) i

31、f(uij=0) printf( ); else printf(%8.4f,uij); printf(n節(jié)點阻抗矩陣如下所示); float fi10; float h10; for(j=0;j10;j+) fij=0; hj=0; for(j=0;jn;j+) for(i=0;in;i+) if(i-j)0) fii=0; else if(i-j)=0) fii=1; else sum=0; for(kk=0;kki;kk+) sum=sum-ukki*fikk; fii=sum; for(i=0;in;i+) if(i=0;i-) if(i=(n-1) Zij=hi; else sum=0

32、; for(kk=i+1;kkn;kk+) sum=sum-uikk*Zkkj; Zij=hi+sum ; for(i=0;in;i+) for(j=0;jn;j+) Zij=-Zij; for(i=0;in;i+) printf(n); for(j=0;jn;j+) printf(%.4f ,Zij); getch(); system(cls);/*計算短路電壓電流網(wǎng)絡中電流分布*/ /*注入電流計算*/ printf(tttt短路計算結果如下nn); for(i=0;in;i+) Ini=(Emidi)*(yei);/*注入節(jié)點電壓初值*/ for(i=0;in;i+) for(j=0;j

33、n;j+) V0i=V0i+Inj*Zji; /*節(jié)點電壓初值*/ If=V0f-1/Zf-1f-1; for(i=0;in;i+) Vi=V0i-Zif-1*If; printf(1、短路電流：nIf=%.4fn*n2、各節(jié)點電壓為n,If); for(i=0;in;i+) if(if-1) printf(V%d=%.4fn,i+1,Vi); else printf(V%d=0.0000n,i+1,Vi); printf(*n); printf(3、發(fā)電機支路電流如下:n); for(m=0;mtimes;m+) for(i=0;in;i+) if(i+1=em)Igem=-(Em-Vi)/

34、zem;printf(Ig%d0=%.4fn,em,Igem); printf(*n4、各支路電流為n); for(m=0;ma;m+) i=count0m-1; j=count1m-1; Ii=(Vj-Vi)/zij; printf(I%d%d=%.4fn,count1m,count0m,Ii); /*儲存文件*/ printf(是否保存文件，保存請按，不保存按);scanf(%d,&g);switch(g) case 1: char outputname20; printf(5、你輸入一個文件名用以保存數(shù)據(jù)n); scanf(%s,outputname); fp=fopen(outputn

35、ame,w+); fprintf(fp,計算結果如下：n1、節(jié)點導納矩陣:); for(i=0;in;i+) fprintf(fp,n); for(j=0;jn;j+) if(Yij=0) fprintf(fp, ); else fprintf(fp,%8.4f ,Yij); fprintf(fp,n2、節(jié)點阻抗矩陣:);for(i=0;in;i+) fprintf(fp,n); for(j=0;jn;j+) fprintf(fp,%7.4f ,Zij); fprintf(fp,n3、短路節(jié)點是:%d,f); fprintf(fp,n4、短路電流:If=%.4f,If);fprintf(fp,n5、各節(jié)點電壓:n); for(i=0;in;i+) if(if-1) fprintf(fp,V%d=%.4ft,i+1,Vi); else fprintf(fp,V%d=0t,i+1,Vi); fprintf(fp,n6、發(fā)電機支路電流如下:n); for(m=0;mtimes;m+) for(i=0;in;i+) if(i+1=em)Igem=-(Em-Vi)/zem;fprintf(fp,Ig%d0=%.4fn,em,Igem); fprintf(fp,n6、各支路電流:n); for(m=0;mVj) Ii=(Vj-Vi)/zij; fprintf(fp,I%d%d=%.4ft,