故障上機(jī)計(jì)算電力系統(tǒng)故障分析課程設(shè)計(jì)-資料下載頁
2025-06-07 08:52本頁面
【正文】 for(i=0。i3。i++) fprintf(wf,%f ,IF[i])。 for(i=0。i3。i++) fprintf(wf,%f ,IFX[i])。 fprintf(wf,\n)。 fprintf(wf,各支路的序 (1,2,0)，相 (A,B,C)電流為： \n)。 prog5(IL,I,m1)。 fprintf(wf,各結(jié)點(diǎn)序 (1,2,0),相電壓 (A,B,C)為 :\n)。 prog5(UN,U,n1)。 } void prog1(float Y[N1][N1],float z[M1],int n1[M1],int n2[M1],int n,int m) /*形成結(jié)點(diǎn)導(dǎo)納矩陣 */ { float y[M1]。 int i,j,k。 for(i=0。in。i++) for(j=0。jn。j++) Y[i][j]=0。 for(i=0。im。i++) y[i]=1/z[i]。 for(i=0。in。i++) for(j=0。jm。j++) if(n1[j]==i+1||n2[j]==i+1) Y[i][i]+=y[j]。 for(i=0。in。i++) for(j=0。jn。j++) for(k=0。km。k++) if(i!=jamp。amp。((n1[k]==i+1amp。amp。n2[k]==j+1)||(n1[k]==j+1amp。amp。n2[k]==i+1))) Y[i][j]=y[k]。 } void prog2(float Y[][N1],float Z[][N1],int n) /*形成結(jié)點(diǎn)阻抗矩陣 */ { /*三角分解法 */ float L[N1][N1],D[N1],I[N1],U[N1]。 int i,j,k。 for(i=0。in。i++) { D[i]=Y[i][i]。 for(j=0。ji。j++) L[i][j]=Y[i][j]。 } for(k=0。kn。k++) { L[k][k]=1。 for(j=0。jk。j++) { for(i=0。ij。i++) L[k][j]=L[k][i]*D[i]*L[j][i]。 L[k][j]/=D[j]。 } for(i=0。ik。i++) D[k]=L[k][i]*L[k][i]*D[i]。 } for(i=0。in。i++) { for(j=0。jn。j++) I[j]=U[j]=0。 I[i]=1。 for(k=0。kn。k++) { for(j=0。jk。j++) U[k]+=L[k][j]*U[j]。 U[k]=I[k]U[k]。 } for(k=0。kn。k++) U[k]/=D[k]。 for(k=n1。k=0。k) for(j=k+1。jn。j++) U[k]=L[j][k]*U[j]。 for(k=0。kn。k++) Z[k][i]=U[k]。 } } void prog3(float A[][3],float B[][3],int n) /*由序分量形成相分量 */ { float R[M1][3],I[M1][3]。 int i,j,k。 for(i=0。in。i++) for(j=0。j3。j++) R[i][j]=I[i][j]=0。 for(i=0。in。i++) for(j=0。j3。j++) { for(k=0。k3。k++) { R[i][j]+=SR[j][k]*A[i][k]。 I[i][j]+=SI[j][k]*A[i][k]。 } B[i][j]=R[i][j]*R[i][j]+I[i][j]*I[i][j]。 B[i][j]=sqrt(B[i][j])。 } } void prog4(float Y[][N1],int n) /*輸出矩陣 */ { int i,j。 for(i=0。in。i++) { for(j=0。jn。j++) fprintf(wf,%f ,Y[i][j])。 fprintf(wf,\n)。 } } void prog5(float A[][3],float B[][3],int n) /*輸出矩陣 */ { int i,j。 for(i=0。in。i++) { fprintf(wf,%d ,i+1)。 for(j=0。j3。j++) fprintf(wf,%f ,A[i][j])。 for(j=0。j3。j++) fprintf(wf,%f ,B[i][j])。 fprintf(wf,\n)。 } }
畢業(yè)設(shè)計(jì)相關(guān)推薦
鄂ICP備17016276號-1