C语言实现P-Q分解法潮流计算

电力系统课程设计，用C语言实现P-Q分解法潮流计算。其中用到了高斯消去法解矩阵方程，求修正的雅可比矩阵。

#include 
#include 

  #define pi 3.1415926535
      #define S12 S[0][1]
      #define S13 S[0][2]
      #define S14 S[0][3]
      #define S21 S[1][0]
      #define S24 S[1][3]
      #define S31 S[2][0]
      #define S41 S[3][0]
      #define S42 S[3][1]
      #define s12 s[0][1]
      #define s13 s[0][2]
      #define s14 s[0][3]
      #define s21 s[1][0]
      #define s24 s[1][3]
      #define s31 s[2][0]
      #define s41 s[3][0]
      #define s42 s[3][1]
      
      
  
 void xsh2（int mint ndouble a[][n]）  //显示二维矩阵 
   { int ij;
     for（i=0;i     {for（j=0;j     printf（“%8.6f  “a[i][j]）;
     printf（“\n“）;}
     printf（“\n“）;printf（“\n“）; 
    }
 void xsh1（int ndouble a[]）//显示一维数组
        {int i=0;
        while（i        {printf（“%8.6f  “a[i]）;i++;}
        printf（“\n“）;printf（“\n“）;
           }  
           
           
           
           
      #define Y11 Y[0][0]
      #define Y12 Y[0][1]
      #define Y13 Y[0][2]
      #define Y14 Y[0][3]
      #define Y21 Y[1][0]
      #define Y22 Y[1][1]
      #define Y23 Y[1][2]
      #define Y24 Y[1][3]
      #define Y31 Y[2][0]
      #define Y32 Y[2][1]
      #define Y33 Y[2][2]
      #define Y34 Y[2][3]
      #define Y41 Y[3][0]
      #define Y42 Y[3][1]
      #define Y43 Y[3][2]
      #define Y44 Y[3][3]
      #define y11 y[0][0]
      #define y12 y[0][1]
      #define y13 y[0][2]
      #define y14 y[0][3]
      #define y21 y[1][0]
      #define y22 y[1][1]
      #define y23 y[1][2]
      #define y24 y[1][3]
      #define y31 y[2][0]
      #define y32 y[2][1]
      #define y33 y[2][2]
      #define y34 y[2][3]
      #define y41 y[3][0]
      #define y42 y[3][1]
      #define y43 y[3][2]
      #define y44 y[3][3]
      #define y120 yo[0][1]
      #define y130 yo[0][2]
      #define y140 yo[0][3]
      #define y210 yo[1][0]
      #define y240 yo[1][3]
      #define y310 yo[2][0]
      #define y410 yo[3][0]
      #define y420 yo[3][1] 
   double Y[4][4]y[4][4];
   double b[3][3]B[2][2];
   double yo[4][4]={0};  
   double Z12z12Z13z13Z14z14Z24z24 ;
   double k;  //变压器变比k 
  void   juzhen（）
       {  int ij;
   y210=y120;y410=y140;y420=y240;y130=k*（1-k）/z13;y310=（k-1）/z13;
   Y11=Z14/（Z14*Z14+z14*z14）+k*k*Z13/（Z13*Z13+z13*z13）+Z12/（Z12*Z12+z12*z12）;
   y11=y140-z14/（Z14*Z14+z14*z14）-k*k*z13/（Z13*Z13+z13*z13）+y120-z12/（Z12*Z12+z12*z12）;
   Y21=Y12=-Z12/（Z12*Z12+z12*z12）;
   y21=y12=z12/（Z12*Z12+z12*z12）;
   Y13=Y31=-k*Z13/（Z13*Z13+z13*z13）;
   y13=y31=k*z13/（Z13*Z13+z13*z13）;
   Y14=Y41=-Z14/（Z14*Z14+z14*z14）;
   y14=y41=z14/（Z14*Z14+z14*z14）;
   Y22=Z24/（Z24*Z24+z24*z24）+Z12/（Z12*Z12+z12*z12）;
   y22=y240+y210-z24/（Z24*Z24+z24*z24）-z12/（Z12*Z12+z12*z12）;
   Y23=Y32=y23=y32=0;
   Y24=Y42=-Z24/（Z24*Z24+z24*z24）;
   y24=y42=z24/（Z24*Z24+z24*z24）;
   Y33=Z13/（Z13*Z13+z13*z13）;
   y33=-z13/（Z13*Z13+z13*z13）;
   Y34=Y43=y43=y34=0;
   Y44=Z14/（Z14*

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1035  2008-12-22 20:53  电力系统分析课设\P-Q分解法潮流计算C语言实现\Makefile.win

     文件       9598  2008-12-24 10:28  电力系统分析课设\P-Q分解法潮流计算C语言实现\szg.c

     文件        837  2008-12-22 20:49  电力系统分析课设\P-Q分解法潮流计算C语言实现\szg.dev

     文件      23132  2008-12-24 16:33  电力系统分析课设\P-Q分解法潮流计算C语言实现\szg.exe

     文件      12000  2008-12-22 20:52  电力系统分析课设\P-Q分解法潮流计算C语言实现\szg.o

     文件        323  2009-01-06 15:42  电力系统分析课设\readme.txt

     目录          0  2008-12-24 16:33  电力系统分析课设\P-Q分解法潮流计算C语言实现

     目录          0  2009-01-06 15:35  电力系统分析课设

----------- ---------  ---------- -----  ----

                46925                    8