如图所示系统,试计算潮流分布,相关数据见《版潮流计算用户手册》P121。#include<>#include<>floatdivRe(floatb1,floatb2,floatb3,floatb4){floata1r;a1r=(b1*b3+b2*b4)/(b3*b3+b4*b4);return(a1r);}floatdivIm(floatb1,floatb2,floatb3,floatb4){floata1i;a1i=(b2*b3-b1*b4)/(b3*b3+b4*b4);return(a1i);}floatmulRe(floatb1,floatb2,floatb3,floatb4){floata2r;a2r=b1*b3-b2*b4;return(a2r);}floatmulIm(floatb1,floatb2,floatb3,floatb4){floata2i;a2i=b2*b3+b1*b4;return(a2i);}floatMax(floata[],intn){inti;floatmax;max=fabs(a[0]);for(i=1;imax)max=fabs(a[i]);}return(max);}voidmain(){inti,j,k,h,km;intT=16;floateps,sumpi1,sumpi2,sumqi1,sumqi2,max,sumir,sumii,I1r,I1i,t,xx,xxx;floatpi0[8],qi0[8],detpi[8],detqi[8],Iir0[8],Iii0[8],J0[16][16],detsi[16],detui[16],delta_p[9][9],delta_q[9][9],a[16][32],ni[16][16],H[8][8],N[8][8],J[8][8],L[8][8],ei1[9],fi1[9],sp[9][9],sq[9][9];staticfloatybr[9][9]={{,,0,0,0,,0,0,0},{,,,0,0,0,0,0,0},{0,,,,0,0,0,0,0},{0,0,,,,0,0,0,0},{0,0,0,,,,0,0,0},{,0,0,0,,,0,0,0},{0},{0},{0}};staticfloatybi[9][9]={{,,0,0,0,,0,0,},{,,,0,0,0,0,0,0},{0,,,,0,0,16,0,0},{0,0,,,,0,0,0,0},{0,0,0,,,,0,,0},{,0,0,0,,,0,0,0},{0,0,16,0,0,0,-16,0,0},{0,0,0,0,,0,0,,0},{,0,0,0,0,0,0,0,}};staticfloatyd[9][9]={{0,,0,0,0,,0,0,0},{,0,,0,0,0,0,0,0},{0,,0,,0,0,0,0,0},{0,0,,0,,0,0,0,0},{0,0,0,,0,,0,0,0},{,0,0,0,,0,0,0,0},{0},{0},{0}};floatei0[9]={,,,,,,,,};floatfi0[9]={};floatpi[9]={0,,0,,0,,,,0};floatqi[9]={0,,0,,0,,0,0,0};h=0;km=15;eps=;do{h+=1;printf("\nNowThe%dthtimes\n",h);for(i=0;i<8;i++){printf("ei0[%d]=%f\t",i,ei0[i]);printf("fi0[%d]=%f\t",i,fi0[i]);}for(i=0;i<8;i++){printf("pi[%d]=%f\t",i,pi[i]);printf("qi[%d]=%f\t",i,qi[i]);}sumpi2=0;sumqi2=0;for(i=0;i<8;i++){for(j=0;j<9;j++){sumpi1=(ei0[i]*(ybr[i][j]*ei0[j]-ybi[i][j]*fi0[j])+fi0[i]*(ybr[i][j]*fi0[j]+ybi[i][j]*ei0[j]));sumpi2+=sumpi1;}pi0[i]=sumpi2;printf("pi0[%d]=%f\t",i,pi0[i]);sumpi2=0;}for(i=0;i<8;i++){for(j=0;j<9;j++){sumqi1=(fi0[i]*(ybr[i][j]*ei0[j]-ybi[i][j]*fi0[j])-ei0[i]*(ybr[i][j]*fi0[j]+ybi[i][j]*ei0[j]));sumqi2+=sumqi1;}qi0[i]=sumqi2;printf("qi0[%d]=%f\t",i,qi0[i]);sumqi2=0;}for(i=0;i<8;i++){detpi[i]=pi[i]-pi0[i];detqi[i]=qi[i]-qi0[i];if(i==6||i==7){qi0[i]=ei0[i]*ei0[i]+fi0[i]*fi0[i];detqi[i]=[i];}printf("detpi[%d]=%f\t",i,detpi[i]);printf("detqi[%d]=%f\t",i,detqi[i]);}//*************************************************************************************************//节点的注入电流表达式for(i=0;i<8;i++){Iii0[i]=0;Iir0[i]=0;}for(i=0;i<8;i++){for(j=0;j<9;j++){Iir0[i]+=ybr[i][j]*ei0[j]-ybi[i][j]*fi0[j];Iii0[i]+=ybr[i][j]*fi0[j]+ybi[i][j]*ei0[j];}}//*********************************************************************************************//求解NHJL矩阵for(i=0;i<8;i++){for(j=0;j<8;j++){if(i==j){if(i==6||i==7){H[i][j]=-ybi[i][j]*ei0[j]+ybr[i][j]*fi0[j]+Iii0[i];N[i][j]=ybr[i][j]*ei0[j]+ybi[i][j]*fi0[j]+Iir0[i];J[i][j]=2*fi0[i];L[i][j]=2*ei0[i];}else{H[i][j]=-ybi[i][j]*ei0[j]+ybr[i][j]*fi0[j]+Iii0[i];N[i][j]=ybr[i][j]*ei0[j]+ybi[i][j]*fi0[j]+Iir0[i];J[i][j]=-ybr[i][j]*ei0[j]-ybi[i][i]*fi0[j]+Iir0[i];L[i][j]=-ybi[i][j]*ei0[j]+ybr[i][j]*fi0[j]-Iii0[i];}}else{if(i==6||i==7){H[i][j]=ybr[i][j]*fi0[j]-ybi[i][j]*ei0[j];N[i][j]=ybr[i][j]*ei0[j]+ybi[i][j]*fi0[j];J[i][j]=0;L[i][j]=0;}else{H[i][j]=ybr[i][j]*fi0[j]-ybi[i][j]*ei0[j];N[i][j]=ybr[i][j]*ei0[j]+ybi[i][j]*fi0[j];J[i][j]=-ybi[i][j]*fi0[j]-ybr[i][j]*ei0[j];L[i][j]=ybr[i][j]*fi0[j]-yb...