#2

//================================================

//   线性方程组的求解，采用向量运算和行消元方法

//================================================

int System_Analyse::Solute_Y(Matrix_element **Y,_complex *Sii,int dimension)

{

int i,F=0;

double a;

Matrix_element *p1,*p2,*p3,*p4,*p5,*pp;

p1=new Matrix_element;

p2=new Matrix_element;

p3=new Matrix_element;

p4=new Matrix_element;

p5=new Matrix_element;

pp=new Matrix_element;

//矩阵消元开始；

for(i=1;i<=dimension;i++)

{

p1=Y[i];                       //取出要消元的行；

if(p1==NULL)

{

Y[i]=new Matrix_element;

Y[i]->G=0.0;

Y[i]->B=0.0;

Y[i]->J=i;

Y[i]->next=NULL;

continue;

}//end if

if((p1->J==i)&&(p1->G==0.0)&&(p1->B==0.0)) continue;

if(p1==NULL) continue;

while(p1->J

{

p2=Y[p1->J];

while(p2->next!=NULL)            //进行消元；

{

p2=p2->next;

ChengFa(p1,p2,p3);

p4=Search(p1,p2->J,F);

if(F==1)

{

p4->G=p4->G-p3->G;

p4->B=p4->B-p3->B;

}//end if

if(F==2||F==3)

{

p5=new Matrix_element;

p5->next=p4->next;

p4->next=p5;

p5->G=-p3->G;

p5->B=-p3->B;

p5->J=p2->J;

}//end if

}//end while

//计算方程右端的项;

ChengFa(p1,Y[p1->J],p3);

Sii[i].x=Sii[i].x-p3->G;

Sii[i].y=Sii[i].y-p3->B;

p5=p1;

p1=p1->next;

delete(p5);  //释放该单元；

}//end while

//所在行的对角元化为标准的 1+j0 形式；

Y[i]=p1;

a=p1->G*p1->G+p1->B*p1->B;

if(a<=1.0e-15)    //是孤立子块

{

p1->G=1.0;

p1->B=0.0;

a=1.0;

}//end if

// { cout<

//   a=1.0; return 0;  }//if                            //方程出现奇异;

p1->G=p1->G/a;  p1->B=-p1->B/a;  p2=p1->next;

while(p2!=NULL)

{

ChengFa(p1,p2,p2);

p2=p2->next;

}//end while

//计算方程右端的项;

a=p1->G*Sii[i].x-p1->B*Sii[i].y;

Sii[i].y=p1->G*Sii[i].y+p1->B*Sii[i].x;

Sii[i].x=a;

//把方程右边的数放入导纳矩阵的对角元中；

p1->G=Sii[i].x; p1->B=Sii[i].y;

}//end for

//回代消元求解开始；

for(i=dimension-1;i>=1;i--)

{

p1=Y[i];

if(p1==NULL) continue;

// if((p1->G==0.0)&&(p1->B==0.0)) continue;

pp=p1;  p1=p1->next;

while(p1!=NULL)                    //消元,把结果放在对角元中;

{

ChengFa(p1,Y[p1->J],p3);

pp->G-=p3->G; pp->B-=p3->B;

p5=p1;  p1=p1->next;

delete(p5);

}//end while

Y[i]->next=NULL;

}//end for

delete(p3);

return 1;

}//end Solute_Y()

就是对线性方程组进行求解的过程，求解完以后输出Y的值就出现如题所示的问题了