潮流计算上机报告 华电(3)

课程课程设计报告

{

/* 计算偏移量 */ for(i=1;i

h=1;

for(i=1;i

/* 输出偏移量

if((fp=fopen(\偏移量 .txt\{

printf(\偏移量.txt' \\n\exit(0); }

fprintf(fp,\偏移量 ---\\n\

10

课程课程设计报告

for(i=1;i

fprintf(fp,\if(fp!=NULL) fclose(fp);*/

/* 计算雅克比矩阵 */ for(i=1;i

if(i==j) { H[i][j]=PQ[i][1]-V[i][1]*V[i][1]*YB[i][j]; J[i][j]=-PQ[i][0]-V[i][1]*V[i][1]*YG[i][j]; N[i][j]=-PQ[i][0]+V[i][1]*V[i][1]*YG[i][j]; L[i][j]=-PQ[i][1]-V[i][1]*V[i][1]*YB[i][j]; } else { H[i][j]=V[i][1]*V[j][1]*(YG[i][j]*sin(V[i][0]-V[j][0])-YB[i][j]*cos(V[i][0]-V[j][0])); J[i][j]=-V[i][1]*V[j][1]*(YG[i][j]*cos(V[i][0]-V[j][0])+YB[i][j]*sin(V[i][0]-V[j][0])); N[i][j]=V[i][1]*V[j][1]*(YG[i][j]*cos(V[i][0]-V[j][0])+YB[i][j]*sin(V[i][0]-V[j][0])); L[i][j]=V[i][1]*V[j][1]*(YG[i][j]*sin(V[i][0]-V[j][0])-YB[i][j]*cos(V[i][0]-V[j][0])); } } } h=1;

for(i=1;i

11

课程课程设计报告

}

} } h++; l=1; for(j=1;j

if(sBus[i-1].Type==1) { for(j=1;j

/* 输出雅克比矩阵

if((fp=fopen(\雅克比矩阵 .txt\{

printf(\雅克比矩阵.txt' \\n\

12

课程课程设计报告

exit(0); } fprintf(fp,\雅克比矩阵 ---\\n\ fprintf(fp,\ \ for(i=1;id) { d=fabs(JJ[k][j]); /*在一行中找到一个最大值赋值d,并用JS[K]记住这个最大值所在的列号*/ LH[k]=j; } } if(fabs(d)<0.000001) /*如果d的数值太小，做为被除数将带来很大的误差 */ { l=0; } else if(LH[k]!=k) { for(i=1;i

13

课程课程设计报告

{ t=JJ[i][k]; JJ[i][k]=JJ[i][LH[k]]; /*进行列交换，让最大值始终在对角元上*/ JJ[i][LH[k]]=t; } } if(l==0) { break; } for(j=k+1;j

if(fabs(JJ[h-1][h-1])>0.00001)

{ /*用追赶法，解方程组，求未知数x*/ U[h-1]=dPQ[h-1]; for(i=h-2;i>=0;i--) { t=0.0; for(j=i+1;j

/*输出高斯结果

if((fp=fopen(\高斯.txt\ { printf(\}

for(i=1;i

14