電力系統(tǒng)短路故障的計(jì)算機(jī)算法程序設(shè)計(jì)(編輯修改稿)
2025-08-03 15:47 本頁(yè)面
【文章內(nèi)容簡(jiǎn)介】 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。 inet[0].hanghao。 i++ ) { fscanf (fp,%d,%d,%f,amp。net[i].head_node,amp。net[i].end_node,amp。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].hanghao1。i++ ) // 互導(dǎo)納形成 { yy[net[i].head_node1][net[i].end_node1]=1/net[i].react。 yy[net[i].end_node1][net[i].head_node1]=1/net[i].react。 } for (i=0。i=net[0].hanghao1。i++ ) //自導(dǎo)納形成 { for (j=0。j=net[0].hanghao1。j++ ) { 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 (,w)。 fprintf(fp,%s,節(jié)點(diǎn)導(dǎo)納矩陣為：\n)。 fclose(fp)。 printf (************************節(jié)點(diǎn)導(dǎo)納矩陣為：******************************\n\n)。 for (i=0。 inet[0].node。 i++ ) { for (j=0。 jnet[0].node。 j++ ) {printf (%.8f*j\t,yy[i][j])。 fp = fopen (,a)。 fprintf(fp,%.8f*j\t,yy[i][j])。 fclose(fp)。 } printf (\n)。 fp = fopen (,a)。 fprintf(fp,\n)。 fclose(fp)。 } printf (\n)。 fp = fopen (,a)。 fprintf(fp,\n)。 fclose(fp)。}int brinv(double a[],int n) { int *is,*js,i,j,k,l,u,v。 double d,p。 is=(int*)malloc(n*sizeof(int))。 js=(int*)malloc(n*sizeof(int))。 for (k=0。 k=n1。 k++) { d=。 for (i=k。 i=n1。 i++) for (j=k。 j=n1。 j++) { l=i*n+j。 p=fabs(a[l])。 if (pd) { d=p。 is[k]=i。 js[k]=j。} } if (d+==) { free(is)。 free(js)。 printf(該矩陣為奇異矩陣\n)。 return(0)。 } if (is[k]!=k) for (j=0。 j=n1。 j++) { u=k*n+j。 v=is[k]*n+j。 p=a[u]。 a[u]=a[v]。 a[v]=p。 } if (js[k]!=k)
規(guī)章制度相關(guān)推薦
鄂ICP備17016276號(hào)-1