C++实现的IEEE33节点潮流计算程序详解

#include <iostream>
#include <math.h>
#include <iomanip>
using namespace std;
int main()
{ int i,j,k=0,n=0,X1[32],X2[32],sort[32];
float V=1,V1,V2,I,R,X,Pload,Qload,Ploss,Qloss,Vbus[33],Vbus1[33],Pj[33],Qj[33],P[32],Q[32]; 
for(i=0;i<33;i++)  {Vbus[i]=12.66; Vbus1[i]=12.66;}
float branch[32][4] = {   { 0,  1,  0.0922, 0.0407 },{ 1,  2,  0.4930, 0.2511 },{ 2,  3,  0.3660, 0.1864 },{ 3,  4,  0.3811, 0.1941 },
  { 4,  5,  0.8190, 0.7070 },{ 5,  6,  0.1872, 0.6188 },{ 6,  7,  0.7144, 0.2351 },{ 7,  8,  1.0300, 0.7400 },
  { 8,  9,  1.0440, 0.7400 },{ 9,  10, 0.1966, 0.0650 },{ 10, 11, 0.3744, 0.1238 },{ 11, 12, 1.4680, 1.1550 },
  { 12, 13, 0.5416, 0.7129 },{ 13, 14, 0.5910, 0.5260 },{ 14, 15, 0.7463, 0.5450 },{ 15, 16, 1.2890, 1.7210 },
  { 16, 17, 0.7320, 0.5740 },{ 1,  18, 0.1460, 0.1565 },{ 18, 19, 1.5042, 1.3554 },{ 19, 20, 0.4095, 0.4784 },
  { 20, 21, 0.7089, 0.9373 },{ 2,  22, 0.4512, 0.3083 },{ 22, 23, 0.8980, 0.7091 },{ 23, 24, 0.8960, 0.7011 },   
  { 5,  25, 0.2030, 0.1034 },{ 25, 26, 0.2842, 0.1447 },{ 26, 27, 1.0590, 0.9337 },{ 27, 28, 0.8042, 0.7006 },
  { 28, 29, 0.5075, 0.2585 },{ 29, 30, 0.9744, 0.9630 },{ 30, 31, 0.3105, 0.3619 },{ 31, 32, 0.3410, 0.5302 }};
float bus[33][2] = {  { 0,   0   },{ 100, 60  },{ 90,  40  },{ 120, 80  },{ 60,  30  },{ 60,  20  },
  { 200, 100 },{ 200, 100 },{ 60,  20  },{ 60,  20  },{ 45,  30  },{ 60,  35  },
  { 60,  35  },{ 120, 80  },{ 60,  10  },{ 60,  20  },{ 60,  20  },{ 90,  40  },
  { 90,  40  },{ 90,  40  },{ 90,  40  },{ 90,  40  },{ 90,  50  },{ 420, 200 },
  { 420, 200 },{ 60,  25  },{ 60,  25  },{ 60,  20  },{ 120, 70  },{ 200, 600 },
  { 150, 70  },{ 210, 100 },{ 60,  40  }};
for(i=0;i<33;i++)  
{bus[i][0]=0.001*bus[i][0]; bus[i][1]=0.001*bus[i][1];}
for (i=0;i<32;i++)
{
X1[i]= int(branch[i][0]);
X2[i]= int(branch[i][1]);
}
while (X2[0]==1)