故障上機計算電力系統(tǒng)故障分析課程設計(參考版)
2025-01-21 17:14本頁面
【正文】 } }。j++) fprintf(wf,%f ,B[i][j])。 for(j=0。j3。i++) { fprintf(wf,%d ,i+1)。 for(i=0。 fprintf(wf,\n)。jn。in。 } } void prog4(float Y[][N1],int n) /*輸出矩陣*/ { int i,j。 } B[i][j]=R[i][j]*R[i][j]+I[i][j]*I[i][j]。k++) { R[i][j]+=SR[j][k]*A[i][k]。j++) { for(k=0。i++) for(j=0。 for(i=0。j3。in。 int i,j,k。k++) Z[k][i]=U[k]。 for(k=0。jn。k=0。k++) U[k]/=D[k]。 } for(k=0。j++) U[k]+=L[k][j]*U[j]。k++) { for(j=0。 for(k=0。j++) I[j]=U[j]=0。i++) { for(j=0。 } for(i=0。ik。 L[k][j]/=D[j]。ij。jk。k++) { L[k][k]=1。 } for(k=0。ji。i++) { D[i]=Y[i][i]。 for(i=0。 } void prog2(float Y[][N1],float Z[][N1],int n) /*形成結(jié)點阻抗矩陣*/ { /*三角分解法*/ float L[N1][N1],D[N1],I[N1],U[N1]。amp。amp。amp。km。jn。in。j++) if(n1[j]==i+1||n2[j]==i+1) Y[i][i]+=y[j]。i++) for(j=0。 for(i=0。im。j++) Y[i][j]=0。i++) for(j=0。 for(i=0。 } void prog1(float Y[N1][N1],float z[M1],int n1[M1],int n2[M1],int n,int m) /*形成結(jié)點導納矩陣*/ { float y[M1]。 fprintf(wf,各結(jié)點序(1,2,0),相電壓(A,B,C)為:\n)。 fprintf(wf,各支路的序(1,2,0)，相(A,B,C)電流為：\n)。i++) fprintf(wf,%f ,IFX[i])。 for(i=0。i3。 fprintf(wf,故障點的序(1,2,0)，相(A,B,C)電流為：\n)。 else if(kp==4) fprintf(wf,k(1,1)\n)。 else if(kp==2) fprintf(wf,k(2)\n)。 fprintf(wf,故障類型為: )。 fprintf(wf,零序結(jié)點阻抗矩陣為：\n)。 fprintf(wf,零序結(jié)點導納矩陣為：\n)。 fprintf(wf,正序結(jié)點阻抗矩陣為：\n)。 fprintf(wf,正序結(jié)點導納矩陣為：\n)。 /*輸出計算結(jié)果*/ fprintf(wf,零序網(wǎng)絡的結(jié)點數(shù)，支路數(shù)分別為：%d %d\n,n0,m0)。 } prog3(IL,I,m1)。 else if(p==0) IL[i][2]=UN[t1][2]/z0[i]。 p=n02[i]。im0。 IL[i][1]=(UN[t1][1]UN[p1][1])/z1[i]。 IL[i][1]=UN[t1][1]/z1[i]。 IL[i][1]=UN[p1][1]/z1[i]。 p=n12[i]。im1。 } prog3(UN,U,n1)。 UN[i][1]=Z1[i][sign1]*IF[1]。j++) UN[i][0]+=Z1[i][j]*I0[j]。i++) { for(j=0。 for(i=0。j3。in1。 } prog3(IF,IFX,1)。 IF[1]=(IF[0]*ZP0)/(ZP1+ZP0)。 } else if(kp==3) IF[0]=IF[1]=IF[2]=UP1/(ZP2+ZP1+ZP0)。 IF[1]=IF[0]。 IF[1]=IF[2]=0。 ZP0=Z0[sign1][sign1]。i++) UP1+=Z1[sign1][i]*I0
研究報告相關推薦
鄂ICP備17016276號-1