電力系統(tǒng)短路故障的計(jì)算機(jī)算法程序設(shè)計(jì)

電力系統(tǒng)短路故障的計(jì)算機(jī)算法程序設(shè)計(jì)DOC課程設(shè)計(jì)報(bào)告——電力系統(tǒng)短路故障的計(jì)算機(jī)算法程序設(shè)計(jì)專業(yè)：電氣工程與其自動(dòng)化班級(jí)：設(shè)計(jì)者：學(xué)號(hào)：華中科技大學(xué)電氣與電子工程學(xué)院電力系統(tǒng)短路故障的計(jì)算機(jī)算法程序設(shè)計(jì)全文共29頁(yè)，當(dāng)前為第1頁(yè)。電力系統(tǒng)短路故障的計(jì)算機(jī)算法程序設(shè)計(jì)全文共29頁(yè)，當(dāng)前為第1頁(yè)。目錄一．課程設(shè)計(jì)說(shuō)明-------------------------------------(3)二．選用所用語(yǔ)言的理由-------------------------------(3)三．程序主框圖子框圖與主要數(shù)據(jù)變量的說(shuō)明-------------(3)四．程序說(shuō)明〔計(jì)算公式與出處〕-----------------------(4)五．設(shè)計(jì)實(shí)際應(yīng)用-------------------------------------(4)六.設(shè)計(jì)體會(huì)-----------------------------------------(10)七．參考資料----------------------------------------(10)八．附錄〔程序源代碼〕-------------------------------(10)電力系統(tǒng)短路故障的計(jì)算機(jī)算法程序設(shè)計(jì)全文共29頁(yè)，當(dāng)前為第2頁(yè)。電力系統(tǒng)短路故障的計(jì)算機(jī)算法程序設(shè)計(jì)全文共29頁(yè)，當(dāng)前為第2頁(yè)。一．課程設(shè)計(jì)說(shuō)明1.課程設(shè)計(jì)題目為：電力系統(tǒng)短路故障的計(jì)算機(jī)算法程序設(shè)計(jì)。2.實(shí)際任務(wù)為：編寫計(jì)算機(jī)程序完成三相對(duì)稱短路電流的實(shí)際計(jì)算。要求形成可執(zhí)行文件，程序合格，輸入輸出數(shù)據(jù)以文件形式。其中，輸入數(shù)據(jù)為：節(jié)點(diǎn)數(shù)，支路數(shù)，故障節(jié)點(diǎn)，支路首段號(hào)，支路末端號(hào)，支路電抗值。輸出數(shù)據(jù)為：節(jié)點(diǎn)導(dǎo)納矩陣，阻抗矩陣，故障點(diǎn)，故障點(diǎn)電流，支路電流，各節(jié)點(diǎn)電壓。并以《電力系統(tǒng)分析》課本例題6-1，習(xí)題6-10,6-11為例進(jìn)展計(jì)算和驗(yàn)算。二．選用所用語(yǔ)言的理由這次課程設(shè)計(jì)我選用的語(yǔ)言是C語(yǔ)言，原因是在大學(xué)階段只接觸過(guò)C語(yǔ)言，C語(yǔ)言有一定的根底，而且C語(yǔ)言完全能勝任這次編程任務(wù)。電力系統(tǒng)短路故障的計(jì)算機(jī)算法程序設(shè)計(jì)全文共29頁(yè)，當(dāng)前為第3頁(yè)。三．程序主框圖子框圖與主要數(shù)據(jù)變量的說(shuō)明電力系統(tǒng)短路故障的計(jì)算機(jī)算法程序設(shè)計(jì)全文共29頁(yè)，當(dāng)前為第3頁(yè)。從duanlu.txt輸入數(shù)據(jù)從duanlu.txt輸入數(shù)據(jù)根據(jù)輸入形成導(dǎo)納矩陣對(duì)導(dǎo)納矩陣求逆，得到阻抗矩陣通過(guò)課本上的公式，通過(guò)阻抗矩陣求故障時(shí)各點(diǎn)電壓，各支路電流用戶輸入短路點(diǎn)和短路前電壓〔設(shè)負(fù)荷電流無(wú)影響〕輸出結(jié)果到result.txt文件主要變量說(shuō)明如下：init();〔數(shù)據(jù)輸入函數(shù)，形成節(jié)點(diǎn)導(dǎo)納矩陣〕brinv(doublea[],intn);〔矩陣求逆函數(shù)〕out_Y();〔形成節(jié)點(diǎn)阻抗矩陣〕short_Select();〔輸入短路點(diǎn)〕Count_If();〔計(jì)算出短路時(shí)各節(jié)點(diǎn)電壓和各支路電流〕floatyy[M][M]〔導(dǎo)納矩陣〕,zz[M][M]〔阻抗矩陣〕;floatIf[M]〔短路節(jié)點(diǎn)電流〕;floatVol[M]〔各節(jié)點(diǎn)電壓〕;floatIII[M]〔各支路電流〕;電力系統(tǒng)短路故障的計(jì)算機(jī)算法程序設(shè)計(jì)全文共29頁(yè)，當(dāng)前為第4頁(yè)。structpowernet電力系統(tǒng)短路故障的計(jì)算機(jī)算法程序設(shè)計(jì)全文共29頁(yè)，當(dāng)前為第4頁(yè)。{BYTEnode;BYTEline;BYTEhanghao;inthead_node;intend_node;floatreact;}net[N]；〔電力網(wǎng)絡(luò)結(jié)構(gòu)體，用來(lái)表示電力網(wǎng)絡(luò)各參數(shù)〕四．程序說(shuō)明〔計(jì)算公式與出處〕短路計(jì)算過(guò)程如下：輸入數(shù)據(jù)－－形成節(jié)點(diǎn)導(dǎo)納矩陣并求逆－－節(jié)點(diǎn)阻抗矩陣－－選擇短路點(diǎn)f，輸入短路前電壓－－通過(guò)公式〔6-10〕計(jì)算短路電流If――通過(guò)公式〔6-11〕計(jì)算節(jié)點(diǎn)電壓－－用公式〔6-9〕計(jì)算短路后支路電流。If＝Vf〔0〕/Zff〔6-10〕Vi＝Vi〔0〕-Zif*If〔6-11〕Ipq=(Vp-Vq)*(-Ypq)(6-9)各公式引自《電力系統(tǒng)分析》上冊(cè)第六章。五．設(shè)計(jì)實(shí)際應(yīng)用電力系統(tǒng)短路故障的計(jì)算機(jī)算法程序設(shè)計(jì)全文共29頁(yè)，當(dāng)前為第5頁(yè)。例6-3題電力系統(tǒng)短路故障的計(jì)算機(jī)算法程序設(shè)計(jì)全文共29頁(yè)，當(dāng)前為第5頁(yè)。輸入數(shù)據(jù)如下：5571,1,0.151,2,0.1052,3,0.0652,4,0.083,4,0.054,5,0.1845,5,0.22.txt文件輸出數(shù)據(jù)：節(jié)點(diǎn)導(dǎo)納矩陣為：-16.19047546*j 9.52380943*j 0.00000000*j 0.00000000*j 0.00000000*j 電力系統(tǒng)短路故障的計(jì)算機(jī)算法程序設(shè)計(jì)全文共29頁(yè)，當(dāng)前為第6頁(yè)。9.52380943*j -37.40842819*j 15.38461590*j 12.50000000*j 0.00000000*j 電力系統(tǒng)短路故障的計(jì)算機(jī)算法程序設(shè)計(jì)全文共29頁(yè)，當(dāng)前為第6頁(yè)。0.00000000*j 15.38461590*j -35.38461685*j 20.00000000*j 0.00000000*j 0.00000000*j 12.50000000*j 20.00000000*j -37.93478394*j 5.43478251*j 0.00000000*j 0.00000000*j 0.00000000*j 5.43478251*j -9.98023701*j 節(jié)點(diǎn)阻抗矩陣為：0.11813837*j 0.09583523*j 0.09017095*j 0.08581381*j 0.04673029*j 0.09583523*j 0.16291989*j 0.15329060*j 0.14588347*j 0.07944150*j 0.09017095*j 0.15329060*j 0.18598767*j 0.16113926*j 0.08774910*j 0.08581381*j 0.14588347*j 0.16113926*j 0.17287450*j 0.09413958*j 0.04673029*j 0.07944150*j 0.08774910*j 0.09413958*j 0.15146215*j 短路節(jié)點(diǎn)為:3電力系統(tǒng)短路故障的計(jì)算機(jī)算法程序設(shè)計(jì)全文共29頁(yè)，當(dāng)前為第7頁(yè)。電力系統(tǒng)短路故障的計(jì)算機(jī)算法程序設(shè)計(jì)全文共29頁(yè)，當(dāng)前為第7頁(yè)。短路前電壓為(忽略負(fù)荷電流的影響)：1.000000短路點(diǎn)的電流為:If[3]=-5.376700*j短路節(jié)點(diǎn)電壓為:Vol[1]=0.515178Vol[2]=0.175802Vol[3]=0.000000Vol[4]=0.133602Vol[5]=0.528199各支路電流為:III[1][1]=-3.232147*jIII[1][2]=-3.232147*j電力系統(tǒng)短路故障的計(jì)算機(jī)算法程序設(shè)計(jì)全文共29頁(yè)，當(dāng)前為第8頁(yè)。III[2][3]=-2.704651*j電力系統(tǒng)短路故障的計(jì)算機(jī)算法程序設(shè)計(jì)全文共29頁(yè)，當(dāng)前為第8頁(yè)。III[2][4]=-0.527499*jIII[3][4]=2.672048*jIII[4][5]=2.144548*jIII[5][5]=-2.144548*j習(xí)題6—10輸入數(shù)據(jù)：5571,1,0.231,3,0.1052,2,0.282,4,0.213,4,0.43554電力系統(tǒng)短路故障的計(jì)算機(jī)算法程序設(shè)計(jì)全文共29頁(yè)，當(dāng)前為第9頁(yè)。3,5,0.2904電力系統(tǒng)短路故障的計(jì)算機(jī)算法程序設(shè)計(jì)全文共29頁(yè)，當(dāng)前為第9頁(yè)。4,5,0.2541輸出數(shù)據(jù)：節(jié)點(diǎn)導(dǎo)納矩陣為：-13.87163544*j0.0000000*j 9.52380943*j 0.00000000*j0.0000000*j 0.00000000*j -8.33333302*j 0.00000000*j 4.76190472*j 0.00000000*j 9.52380943*j 0.000000*j-15.26333618*j2.29600048*j3.44352627*j 0.00000000*j 4.76190472*j 2.29600048*j-10.99336338*j 3.93545866*j 0.00000000*j 0.00000000*j 3.44352627*j 3.93545866*j -7.37898493*j 節(jié)點(diǎn)阻抗矩陣為：0.18042089*j 0.06035720*j 0.15778695*j 0.10562509*j 0.12996729*j 0.06035720*j 0.20652170*j 0.08791157*j 0.15141298*j 0.12177899*j 電力系統(tǒng)短路故障的計(jì)算機(jī)算法程序設(shè)計(jì)全文共29頁(yè)，當(dāng)前為第10頁(yè)。0.15778695*j 0.08791157*j 0.22982012*j 0.15384525*j 0.18930019*j 電力系統(tǒng)短路故障的計(jì)算機(jī)算法程序設(shè)計(jì)全文共29頁(yè)，當(dāng)前為第10頁(yè)。0.10562509*j 0.15141298*j 0.15384525*j 0.26497272*j 0.21311322*j 0.12996729*j 0.12177899*j 0.18930019*j 0.21311322*j 0.33752048*j 短路節(jié)點(diǎn)為:5短路前電壓為(忽略負(fù)荷電流的影響)：1.000000短路點(diǎn)的電流為:If[5]=-2.962783*j短路節(jié)點(diǎn)電壓為:Vol[1]=0.614935Vol[2]=0.639195Vol[3]=0.439145Vol[4]=0.368592電力系統(tǒng)短路故障的計(jì)算機(jī)算法程序設(shè)計(jì)全文共29頁(yè)，當(dāng)前為第11頁(yè)。電力系統(tǒng)短路故障的計(jì)算機(jī)算法程序設(shè)計(jì)全文共29頁(yè)，當(dāng)前為第11頁(yè)。Vol[5]=0.000000各支路電流為:III[1][1]=-1.674195*jIII[1][3]=-1.674195*jIII[2][2]=-1.288588*jIII[2][4]=-1.288588*jIII[3][4]=-0.161989*jIII[3][5]=-1.512206*jIII[4][5]=-1.450577*j習(xí)題6—11電力系統(tǒng)短路故障的計(jì)算機(jī)算法程序設(shè)計(jì)全文共29頁(yè)，當(dāng)前為第12頁(yè)。3電力系統(tǒng)短路故障的計(jì)算機(jī)算法程序設(shè)計(jì)全文共29頁(yè)，當(dāng)前為第12頁(yè)。312輸入數(shù)據(jù)：3341,1,0.21,2,0.22,3,0.13,3,0.15輸出數(shù)據(jù)：節(jié)點(diǎn)導(dǎo)納矩陣為：-10.00000000*j 5.00000000*j 0.00000000*j 5.00000000*j -15.00000000*j 10.00000000*j 0.00000000*j 10.00000000*j -16.66666603*j 節(jié)點(diǎn)阻抗矩陣為：0.13846155*j 0.07692308*j 0.04615385*j 0.07692308*j 0.15384616*j 0.09230770*j 0.04615385*j 0.09230770*j 0.11538462*j 短路節(jié)點(diǎn)為:3電力系統(tǒng)短路故障的計(jì)算機(jī)算法程序設(shè)計(jì)全文共29頁(yè)，當(dāng)前為第13頁(yè)。短路前電壓為(忽略負(fù)荷電流的影響)：1.050000電力系統(tǒng)短路故障的計(jì)算機(jī)算法程序設(shè)計(jì)全文共29頁(yè)，當(dāng)前為第13頁(yè)。短路點(diǎn)的電流為:If[3]=-9.099999*j短路節(jié)點(diǎn)電壓為:Vol[1]=0.630000Vol[2]=0.210000Vol[3]=-0.000000各支路電流為:III[1][1]=-2.100000*jIII[1][2]=-2.100000*jIII[2][3]=-2.100000*jIII[3][3]=-7.000000*j六．設(shè)計(jì)體會(huì)電力系統(tǒng)短路故障的計(jì)算機(jī)算法程序設(shè)計(jì)全文共29頁(yè)，當(dāng)前為第14頁(yè)。這次設(shè)計(jì)體會(huì)良多，首先是通過(guò)短路電流的計(jì)算，熟悉了以前電力系統(tǒng)分析的課程，鞏固了專業(yè)知識(shí)，同時(shí)通過(guò)C語(yǔ)言的使用也使幾乎已經(jīng)忘干凈的C語(yǔ)言技能得以回復(fù)和鞏固。電力系統(tǒng)短路故障的計(jì)算機(jī)算法程序設(shè)計(jì)全文共29頁(yè)，當(dāng)前為第14頁(yè)。七．參考資料《電力系統(tǒng)分析》上冊(cè)，第三版，何仰贊溫增銀華中科技大學(xué)《程序設(shè)計(jì)根底》〔C語(yǔ)言〕王載新曾大亮楊有安X珂梅清華大學(xué)八．附錄〔程序源代碼〕#include<stdio.h>#include<math.h>#include<conio.h>#include<string.h>#include<stdlib.h>#include<malloc.h>#defineN100#defineM100#defineBYTEunsignedcharinit();intbrinv(doublea[],intn);out_Y();電力系統(tǒng)短路故障的計(jì)算機(jī)算法程序設(shè)計(jì)全文共29頁(yè)，當(dāng)前為第15頁(yè)。intshort_Select();電力系統(tǒng)短路故障的計(jì)算機(jī)算法程序設(shè)計(jì)全文共29頁(yè)，當(dāng)前為第15頁(yè)。Count_If();structpowernet{BYTEnode;BYTEline;BYTEhanghao;inthead_node;intend_node;floatreact;}net[N];floatyy[M][M],zz[M][M];floatIf[M];floatVol[M];floatIII[M];voidmain(){init();//輸入原始數(shù)據(jù)電力系統(tǒng)短路故障的計(jì)算機(jī)算法程序設(shè)計(jì)全文共29頁(yè)，當(dāng)前為第16頁(yè)。out_Y();//形成節(jié)點(diǎn)導(dǎo)納、節(jié)點(diǎn)阻抗電力系統(tǒng)短路故障的計(jì)算機(jī)算法程序設(shè)計(jì)全文共29頁(yè)，當(dāng)前為第16頁(yè)。Count_If();//計(jì)算短路電流，節(jié)點(diǎn)電壓，支路電流}init(){inti,j;floattemp=0;FILE*fp;if((fp=fopen("duanlu.txt","r"))==NULL){printf("程序使用方法：\n\n在與程序的同一文件夾下建立duanlu.txt文本文件\nduanlu.txt文本文件第一行輸入的是：節(jié)點(diǎn)數(shù)支路數(shù)輸入數(shù)據(jù)的行數(shù)。\n〔例如：

335〕。表示3個(gè)節(jié)點(diǎn)，3個(gè)支路,txt文件里輸入有5行數(shù)據(jù)〔不包括本行〕。\n回車另起一行為：首節(jié)點(diǎn)，末節(jié)電，支路電抗。〔例如：1，2，0.15〕\n\n");B: gotoB;}電力系統(tǒng)短路故障的計(jì)算機(jī)算法程序設(shè)計(jì)全文共29頁(yè)，當(dāng)前為第17頁(yè)。電力系統(tǒng)短路故障的計(jì)算機(jī)算法程序設(shè)計(jì)全文共29頁(yè)，當(dāng)前為第17頁(yè)。fscanf(fp,"%d%d%d",&net[0].node,&net[0].line,&net[0].hanghao);printf("本系統(tǒng)節(jié)點(diǎn)數(shù)：%d,支路數(shù)：%d\n\n",net[0].node,net[0].line);printf("以下為線路阻抗：\n\n");printf("始節(jié)點(diǎn)末節(jié)電電抗\n");for(i=0;i<net[0].hanghao;i++){fscanf(fp,"%d,%d,%f",&net[i].head_node,&net[i].end_node,&net[i].react);printf("%4d,%8d,%15f*j\n",net[i].head_node,net[i].end_node,net[i].react);}printf("\n");for(i=0;i<=net[0].hanghao-1;i++)//互導(dǎo)納形成{yy[net[i].head_node-1][net[i].end_node-1]=-1/net[i].react;yy[net[i].end_node-1][net[i].head_node-1]=-1/net[i].react;}for(i=0;i<=net[0].hanghao-1;i++)//自導(dǎo)納形成{for(j=0;j<=net[0].hanghao-1;j++)電力系統(tǒng)短路故障的計(jì)算機(jī)算法程序設(shè)計(jì)全文共29頁(yè)，當(dāng)前為第18頁(yè)。{電力系統(tǒng)短路故障的計(jì)算機(jī)算法程序設(shè)計(jì)全文共29頁(yè)，當(dāng)前為第18頁(yè)。if(net[j].head_node==(i+1)||net[j].end_node==(i+1))temp=1/net[j].react+temp;}yy[i][i]=temp;temp=0;}fclose(fp);fp=fopen("result.txt","w");fprintf(fp,"%s","節(jié)點(diǎn)導(dǎo)納矩陣為：\n");fclose(fp);printf("************************節(jié)點(diǎn)導(dǎo)納矩陣為：******************************\n\n");for(i=0;i<net[0].node;i++){for(j=0;j<net[0].node;j++){printf("%.8f*j\t",-yy[i][j]);fp=fopen("result.txt","a"); fprintf(fp,"%.8f*j\t",-yy[i][j]); fclose(fp); } printf("\n");電力系統(tǒng)短路故障的計(jì)算機(jī)算法程序設(shè)計(jì)全文共29頁(yè)，當(dāng)前為第19頁(yè)。 fp=fopen("result.txt","a");電力系統(tǒng)短路故障的計(jì)算機(jī)算法程序設(shè)計(jì)全文共29頁(yè)，當(dāng)前為第19頁(yè)。 fprintf(fp,"\n"); fclose(fp);}printf("\n"); fp=fopen("result.txt","a"); fprintf(fp,"\n"); fclose(fp);}intbrinv(doublea[],intn){int*is,*js,i,j,k,l,u,v;doubled,p;is=(int*)malloc(n*sizeof(int));js=(int*)malloc(n*sizeof(int));for(k=0;k<=n-1;k++){d=0.0;for(i=k;i<=n-1;i++)for(j=k;j<=n-1;j++){l=i*n+j;p=fabs(a[l]);if(p>d){d=p;is[k]=i;js[k]=j;}}電力系統(tǒng)短路故障的計(jì)算機(jī)算法程序設(shè)計(jì)全文共29頁(yè)，當(dāng)前為第20頁(yè)。if(d+1.0==1.0)電力系統(tǒng)短路故障的計(jì)算機(jī)算法程序設(shè)計(jì)全文共29頁(yè)，當(dāng)前為第20頁(yè)。{free(is);free(js);printf("該矩陣為奇異矩陣\n");return(0);}if(is[k]!=k)for(j=0;j<=n-1;j++){u=k*n+j;v=is[k]*n+j;p=a[u];a[u]=a[v];a[v]=p;}if(js[k]!=k)for(i=0;i<=n-1;i++){u=i*n+k;v=i*n+js[k];p=a[u];a[u]=a[v];a[v]=p;}l=k*n+k;a[l]=1.0/a[l];for(j=0;j<=n-1;j++)if(j!=k){u=k*n+j;a[u]=a[u]*a[l];}for(i=0;i<=n-1;i++)if(i!=k)for(j=0;j<=n-1;j++)電力系統(tǒng)短路故障的計(jì)算機(jī)算法程序設(shè)計(jì)全文共29頁(yè)，當(dāng)前為第21頁(yè)。if(j!=k)電力系統(tǒng)短路故障的計(jì)算機(jī)算法程序設(shè)計(jì)全文共29頁(yè)，當(dāng)前為第21頁(yè)。{u=i*n+j;a[u]=a[u]-a[i*n+k]*a[k*n+j];}for(i=0;i<=n-1;i++)if(i!=k){u=i*n+k;a[u]=-a[u]*a[l];}}for(k=n-1;k>=0;k--){if(js[k]!=k)for(j=0;j<=n-1;j++){u=k*n+j;v=js[k]*n+j;p=a[u];a[u]=a[v];a[v]=p;}if(is[k]!=k)for(i=0;i<=n-1;i++){u=i*n+k;v=i*n+is[k];p=a[u];a[u]=a[v];a[v]=p;}}free(is);free(js);return(1);電力系統(tǒng)短路故障的計(jì)算機(jī)算法程序設(shè)計(jì)全文共29頁(yè)，當(dāng)前為第22頁(yè)。}電力系統(tǒng)短路故障的計(jì)算機(jī)算法程序設(shè)計(jì)全文共29頁(yè)，當(dāng)前為第22頁(yè)。out_Y(){inti,j;double**a;a=(float**)malloc(net[0].node*net[0].node*sizeof(float));//為矩陣分配動(dòng)態(tài)存儲(chǔ)空間for(i=0;i<net[0].node;i++) {a[i]=(float*)malloc(net[0].node*sizeof(float)); }FILE*fp;fp=fopen("result.txt","a");fprintf(fp,"%s","節(jié)點(diǎn)阻抗矩陣為：\n");fclose(fp);for(i=0;i<net[0].node;i++)for(j=0;j<net[0].node;j++){a[i][j]=yy[i][j];}i=brinv(*a,net[0].node);電力系統(tǒng)短路故障的計(jì)算機(jī)算法程序設(shè)計(jì)全文共29頁(yè)，當(dāng)前為第23頁(yè)。for(i=0;i<net[0].node;i++)電力系統(tǒng)短路故障的計(jì)算機(jī)算法程序設(shè)計(jì)全文共29頁(yè)，當(dāng)前為第23頁(yè)。for(j=0;j<net[0].node;j++){zz[i][j]=(float)a[i][j];}printf("************************節(jié)點(diǎn)阻抗矩陣為：******************************\n\n");for(i=0;i<net[0].node;i++){for(j=0;j<net[0].node;j++) {printf("%.8f*j\t",zz[i][j]);fp=fopen("result.txt","a"); fprintf(fp,"%.8f*j\t",zz[i][j]); fclose(fp); } printf("\n"); fp=fopen("result.txt","a"); fprintf(fp,"\n"); fclose(fp);}printf("\n");電力系統(tǒng)短路故障的計(jì)算機(jī)算法程序設(shè)計(jì)全文共29頁(yè)，當(dāng)前為第24頁(yè)。 fp=fopen("result.txt","a");電力系統(tǒng)短路故障的計(jì)算機(jī)算法程序設(shè)計(jì)全文共29頁(yè)，當(dāng)前為第24頁(yè)。 fprintf(fp,"\n"); fclose(fp);}intShort_Select(){intss;printf("請(qǐng)選擇短路故障點(diǎn):");scanf("%d",&ss);printf("\n\n");returnss;}//計(jì)算節(jié)點(diǎn)電流,節(jié)電電壓,短路點(diǎn)注入單位電流時(shí)解線性方程得各節(jié)電電壓Count_If(){BYTENum;//網(wǎng)絡(luò)的節(jié)點(diǎn)數(shù)BYTEsss;inti,j;電力系統(tǒng)短路故障的計(jì)算機(jī)算法程序設(shè)計(jì)全文共29頁(yè)，當(dāng)前為第25頁(yè)。 floatE;電力系統(tǒng)短路故障的計(jì)算機(jī)算法程序設(shè)計(jì)全文共29頁(yè)，當(dāng)前為第25頁(yè)。FILE*fp; printf("請(qǐng)輸入短路前電壓，短路前為空載:"); scanf("%f",&E);AAA:sss=Short_Select();fp=fopen("result.txt","a");fprintf(fp,"%s%d\n\n%s%3f\n\n","短路節(jié)點(diǎn)為:",sss,"短路前電壓為(忽略負(fù)荷電流的影響)：",E);fclose(fp);if(sss>net[0].node||sss<1) {printf("短路點(diǎn)選擇錯(cuò)誤,"); gotoAAA; }II[sss]=1;Num=net[0].node;If[sss]=E/zz[sss-1][sss-1]; fp=fopen("result.txt","a"); fprintf(fp,"%sIf[%d]=%f*j\n\n%s\n\n","短路點(diǎn)的電流為:",sss,-If[sss],"短路節(jié)點(diǎn)電壓為:"); fclose(fp);電力系統(tǒng)短路故障的計(jì)算機(jī)算法程序設(shè)計(jì)全文共29頁(yè)，當(dāng)前為第26頁(yè)。電力系統(tǒng)短路故障的計(jì)算機(jī)算法程序設(shè)計(jì)全文共29頁(yè)，當(dāng)前為第26頁(yè)。printf("If[%d]=%f*j\n\n",sss,-If[sss]);printf("****************************短路節(jié)點(diǎn)電壓為:************************\n");for(i=0;i<net[0].node;i++){Vol[i]=(float)(E-zz[sss-1][i]*If[sss]);printf(