牛頓-拉夫遜法潮流計(jì)算課程設(shè)計(jì)說(shuō)明書-資料下載頁(yè)
2025-08-02 21:06本頁(yè)面
【正文】 rintf(有功功率增量P[2]=%f,P[2])。 printf( ,)。printf(\n)。printf(無(wú)功功率增量Q[1]=%f,Q[1])。 printf( ,)。printf(\n)。printf(電壓增量V22=%f,V22)。printf(\n)。//************篩選出最大值***********************max=fabs(P[1])fabs(P[2])?fabs(P[1]):fabs(P[2])。max=maxfabs(Q[1])?max:fabs(Q[1])。max=maxfabs(V22)?max:fabs(V22)。printf(max=%f\n,max)。//********************************************while (k1=4){//****形成雅克比矩陣**********************for(j=1。j3。j++){if(1==j){float c=0,d=0。int m。for(m=1。m4。m++){c+=G[1][m]*e[m]B[1][m]*f[m]。 d+=G[1][m]*f[m]+B[1][m]*e[m]。}J[1*N1][j*N1]=cG[1][j]*e[1]B[1][j]*f[1]。J[1*N1][j*N]=d+B[1][j]*e[1]G[1][j]*f[1]。J[1*N][j*N1]=d+B[1][j]*e[1]G[1][j]*f[1]。J[1*N][j*N]=c+G[1][j]*e[1]+B[1][j]*f[1]。}else{J[1*N1][j*N1]=G[1][j]*e[1]B[1][j]*f[1]。 J[1*N][j*N]=G[1][j]*e[1]+B[1][j]*f[1]。 J[1*N1][j*N]=B[1][j]*e[1]G[1][j]*f[1]。 J[1*N][j*N1]=B[1][j]*e[1]G[1][j]*f[1]。}}for(j=1。j3。j++){if(2==j){float c=0,d=0。int m。for(m=1。m4。m++){c+=G[2][m]*e[m]B[2][m]*f[m]。 d+=G[2][m]*f[m]+B[2][m]*e[m]。}J[2*N1][j*N1]=cG[2][j]*e[2]B[2][j]*f[2]。J[2*N1][j*N]=d+B[2][j]*e[2]G[2][j]*f[2]。J[2*N][j*N1]=2*e[2]。J[2*N][j*N]=2*f[2]。}else{J[2*N1][j*N1]=G[2][j]*e[2]B[2][j]*f[2]。J[2*N][j*N]=0。J[2*N1][j*N]=B[2][j]*e[2]G[2][j]*f[2]。J[2*N][j*N1]=0。}}printf(雅克比矩陣是:\n)。for(i=1。i5。i++){for(j=1。j5。j++){printf(%f,J[i][j])。printf( )。}printf(\n)。}//********計(jì)算修正方程*************for(i=1。iM。i++){L[i][i]=1。}for(i=1。iM。i++){U[1][i]=J[1][i]。L[i][1]=J[i][1]/U[1][1]。}for(n=2。nM。n++){for(j=n。jM。j++){sigma1=0。for(s=0。s=n1。s++)sigma1+=L[n][s]*U[s][j]。U[n][j]=J[n][j]sigma1。}for(i=n。iM。i++){sigma2=0。for(s=0。s=n1。s++)sigma2+=L[i][s]*U[s][n]。L[i][n]=(J[i][n]sigma2)/U[n][n]。}}b[1]=P[1]。b[2]=Q[1]。b[3]=P[2]。b[4]=V22。for(i=1。iM。i++){sigma1=0。for(n=1。n=i1。n++)sigma1+=L[i][n]*y[n]。y[i]=b[i]sigma1。}for(i=M1。i=1。i){sigma2=0。for(n=i+1。nM。n++)sigma2+=U[i][n]*x[n]。x[i]=(y[i]sigma2)/U[i][i]。}xe[1]=x[1]。xe[2]=x[3]。 xf[1]=x[2]。xf[2]=x[4]。 printf(節(jié)點(diǎn)電壓:\n)。for(i=1。i3。i++){e[i]+=xe[i]。 f[i]+=xf[i]。}for(i=1。i3。i++){printf(e[%d]=,i)。 printf(%f,e[i])。 printf( ,)。}for(i=1。i3。i++){printf(f[%d]=,i)。 printf(%f,f[i])。 printf( ,)。}printf(\n)。k1=k1+1。goto loop1。}for(j=1。j4。j++){a1+=G[3][j]*e[j]B[3][j]*f[j]。 b1+=G[3][j]*f[j]+B[3][j]*e[j]。}P3=e[3]*a1+f[3]*b1。Q3=f[3]*a1e[3]*b1。printf(P3+jQ3=%f+j%f,P3,Q3)。printf(\n)。 33
法律信息相關(guān)推薦
鄂ICP備17016276號(hào)-1