【正文】 m++) 。 for(m=0。f)。n,amp。 } fscanf(fp,系統(tǒng)參數(shù)如下所示 :\n\n)。 while((fp=fopen(name,r))==NULL) { printf(輸入有錯誤，請重新輸入 \n)。 printf(請輸入一個已保存的文件名 \n)。 float yc[10]。 float ze[10]。 float z[10][10]。 float b[10]。 int m,n,a,f,i=0,j=0,sd,md,butn,times,timesc。 FILE *fp。 calculate2(name)。 fprintf(fp, %d\ %d\為 %d,n,a,fnew)。 fp=fopen(name,r+)。f)。n,amp。 fscanf(fp,系統(tǒng)參數(shù)如下所示 :\n\n)。fnew)。 printf(請輸入新的短路節(jié)點 )。 } } void change(char name[20]) { FILE *fp。 case 2: quit()。 switch(butn) {case 1: calculate(name)。butn)。 /***********選擇操作 ************/ printf(接下來你希望 \n 短路計算請按 \n 退出請按 \n )。i++) fprintf(fp,Y%d=%f\t,e[i],b[i])。 {for(i=0。b[i])。i++) scanf(%d%f,amp。 for(i=0。times)。 /************節(jié)點對地電容 *********************/ printf(請輸入電容支路數(shù) \n)。itimes。 fprintf(fp,\n)。itimes。 getch()。itimes。 /**************寫入文件 ****************/ fprintf(fp,共 %d臺發(fā)電機進入 系統(tǒng) \n,times)。ae[i],amp。i++) scanf(%d%f%f,amp。 for(i=0。times)。 printf(請輸入發(fā)電機臺數(shù)： \n)。 fprintf(fp,z%d %d=%.4f\t,sd,md,dk)。md,amp。i++) { scanf(%d%d%f,amp。 for(i=0。 fprintf(fp,支路電抗參數(shù)如下 \n)。 fprintf(fp,系統(tǒng)參數(shù)如下所示 :\n\n)。f)。n,amp。 //輸入文件名 /*************輸入系統(tǒng)參數(shù) *************/ printf(請輸入系統(tǒng)的節(jié)點數(shù)，支路數(shù)，故障節(jié)點 \n)。 //發(fā)電機支路電抗 /****************************輸入文件 **************************/ printf(請輸入您想要創(chuàng)建系統(tǒng)名稱 (如 )用以保存系統(tǒng)參數(shù) \n)。 //儲存發(fā)電機接入節(jié)點 float ae[10]。//n,a,f分別為節(jié)點數(shù)、支路數(shù)、短路點 float dk。 FILE *fp。 } } } void create() { char name[20]。 case 4: quit()。 case 3: change(name)。 case 2: load()。 switch(butn) { case 1: create()。butn)。 printf(\n 退出請按 )。 printf(短路計算已經(jīng)完成，請選擇下列操作 \n 新建一個系統(tǒng)請按 )。 break。 fprintf(fp,I%d%d=%.4f\t,count[1][m],count[0][m],I[i])。 fprintf(fp,I%d%d=%.4f\t,count[0][m],count[1][m],I[i])。 j=count[1][m]1。ma。 } fprintf(fp,\n各支路電流 :\n)。i++) if(i+1==e[m]) Ig[e[m]]=(E[m]V[i])/ze[m]。m++) { for(i=0。 for(m=0。 else fprintf(fp,V%d=0\t,i+1,V[i])。in。 fprintf(fp,\n各節(jié)點電壓 :\n)。 } fprintf(fp,\n短路節(jié)點是 :%d,f)。jn。i++) { fprintf(fp,\n)。 for(i=0。 else 信息工程學(xué)院課程設(shè)計（論文） 32 fprintf(fp,% ,Y[i][j])。jn。i++) { fprintf(fp,\n)。 for(i=0。 fp=fopen(outputname,w+)。 printf(你輸入一個文件名用以保存數(shù)據(jù) \n)。g)。 } /***************************儲存文件 **********************************/ printf(是否保存文件，保存請按，不保存按 )。 I[i]=(V[j]V[i])/z[i][j]。m++) { i=count[0][m]1。 for(m=0。 printf(Ig[%d][0]=%.4f\n,e[m],Ig[e[m]])。in。mtimes。 printf(發(fā) 電機支路電流如下 :\n)。 else printf(V%d=\n,i+1,V[i])。in。 printf(短路 電流： \nIf=%.4f\n***********************************\n各節(jié)點電壓為 \n,If)。in。 /***********節(jié)點電壓初值 ***********************/ 信息工程學(xué)院課程設(shè)計（論文） 31 If=V0[f1]/Z[f1][f1]。jn。in。i++) In[i]=((Emid[i])*(ye[i]))。 for(i=0。 system(cls)。j++) printf(%.4f ,Z[i][j])。 for(j=0。in。j++) Z[i][j]=Z[i][j]。i++) for(j=0。 } } for(i=0。kk++) { sum=sumu[i][kk]*Z[kk][j]。 for(kk=i+1。i) if(i==(n1)) Z[i][j]=h[i]。 for(i=n1。i++) if(ij) h[i]=0。 } for(i=0。kk++) sum=sumu[kk][i]*fi[kk]。 for(kk=0。 else if((ij)==0) fi[i]=1。in。jn。 h[j]=0。j10。 float h[10]。 } printf(\n節(jié)點阻抗矩陣如下所示 )。j++) if(u[i][j]==0) printf( )。 for(j=0。in。 } printf(\nU矩陣如下所示 )。j++) if(d[i][j]==0) printf( )。 for(j=0。in。 } } } /************LDU屏幕顯示 *************/ printf(\nD矩陣如下所示 )。kk++) sum=sum+u[kk][i]*u[kk][j]*d[kk][kk]。 信息工程學(xué)院課程設(shè)計（論文） 29 else { for(kk=0。j++) { sum=0。 } for(j=i+1。kk++) sum=sum+u[kk][i]*u[kk][i]*d[kk][kk]。 else { for(kk=0。i++) { sum=0。 } /**************迭代法求解 *************/ for(i=0。i++) { l[i][i]=1。 } for(i=0。 d[i][j]=0。j10。i10。 float sum。 float u[10][10]。 } float l[10][10]。j++) if(Y[i][j]==0) printf( )。 for(j=0。in。 信息工程學(xué)院課程設(shè)計（論文） 28 } printf(\n\n節(jié)點導(dǎo)納矩陣為 )。j++) { Y[i][j]=yjs[i][j]。i++) for(j=i+1。 } for(i=0。j++) Y[i][i]=Y[i][i]yjs[i][j]。i++) { for(j=0。 } /**************導(dǎo)納矩陣 *************/ for(i