【文章內(nèi)容簡(jiǎn)介】 。//變壓器勵(lì)磁損耗：S0[0]有功，S0[1]無功double S_org[2]。//支路末端負(fù)載：S_org[0]有功，S_org[1]無功line_part *lchild,*rchild。//支路的后繼兩個(gè)支路 }class BinTree{//定義樹類void PreOrder(){PreOrder(root)。}。//樹的先序遍歷修改電壓void PostOrder(){PostOrder(root)。}。//樹的后序遍歷修改潮流void display(){display(root)。}。//樹的先序遍歷顯示數(shù)據(jù)line_part *root。//樹的根，是一個(gè)支路類}。主函數(shù)中定義的數(shù)據(jù)ofstream outfile。//輸出數(shù)據(jù)流定ifstream infile。//輸入數(shù)據(jù)流定const int M(a)。//支路個(gè)數(shù)常量主程序（復(fù)制粘貼到C++就能用）**************************************************************************************************************************************************************include include includedouble p_cost_all=0。double sqr(double x){//平方計(jì)算函數(shù) return x*x。}。class line_part{//定義支路類 private:double val。//支路排序double U_end。//支路末端電壓double U_begin。//支路首段電壓double X[3]。//支路等效阻抗: X[0]電阻，X[1]感抗，X[2]容納double S_end[2]。//支路末端功率：S_end[0]有功，S_end[1]無功double S_begin[2]。//支路首段功率：S_begin[0]有功，S_begin[1]無功double S0[2]。//變壓器勵(lì)磁損耗：S0[0]有功，S0[1]無功double S_org[2]。//支路末端負(fù)載：S_org[0]有功，S_org[1]無功line_part *lchild,*rchild。public:double U[2]。//支路電壓降落：U[0]電壓降落橫分量，U[1]電壓降落縱分量double k。//變壓器變比public: line_part(){ val=0。U_end=0。U_begin=0。X[0]=0。X[1]=0。X[2]=0。S_end[0]=S_org[0]=0。S_end[1]=S_org[1]=0。S_begin[0]=0。S_begin[1]=0。k=1。S0[0]=0。S0[1]=0。U[0]=0。U[1]=0。lchild = rchild= NULL。}line_part(double vall,double u_end=0,double u_begin=0,double r=0，double x=0,double b=0,double Pe=0,double Xe=0,double K=0,double P0=0,double Q0=0){val=vall。U_end=u_end。U_begin=u_begin。X[0]=r。X[1]=x。X[2]=b。S_end[0]=S_org[0]=Pe。S_end[1]=S_org[1]=Xe。S_begin[0]=0。S_begin[1]=0。k=K。S0[0]=P0。S0[1]=Q0。U[0]=0。U[1]=0。lchild = rchild= NULL。}。friend class BinTree。friend void pass_U(line_part*a,line_part*b,line_part*c)。//電壓傳遞函數(shù)friend void pass_w(line_part*a,line_part*b,line_part*c)。//功率傳遞函數(shù)friend void pass_U2(line_part*a,line_part*b)。//電壓傳遞函數(shù)friend void pass_w2(line_part*a,line_part*b)。//功率傳遞函數(shù)void Sbegin(){(U_end)*X[2]/2))/sqr(Uend)。sqr(U_begin)*X[2]/2sqr(U_end)*X[2]/2。}。void Uend(){double U_heng(0),U_zong(0)。X[2]/2S0[1]。1])/U_begin。//U_heng即是△u2 0])/U_begin。//U_zong即是δu2 sqr(U_zong))/k。double get_val(){}。//支路首段功率計(jì)算函數(shù) double Uend,I2。Uend=k*U_end。I2=(sqr(S_end[0])+sqr(S_end[1]sqrS_begin[0]=S_end[0]+I2*X[0]+S0[0]S_begin[1]=S_end[1]+I2*X[1]+S0[1]//支路末端電壓計(jì)算函數(shù) double p_begin,q_begin。p_begin=S_begin[0]S0[0]。q_begin=S_begin[1]+sqr(U_begin)*U_heng=(p_begin*X[0]+q_begin*X[U_zong=(p_begin*X[1]q_begin*X[U_end=sqrt(sqr(U_beginU_heng)+U[0]=U_heng。U[1]=U_zong。}。//返回支路編號(hào) if(this==0){return1。}else{ if(val0amp。amp。valdouble get_Uend(){//返回支路末端電壓return U_end。}。double get_Ubegin(){//返回支路首段電壓return U_begin。}。double get_Pbegin(){//返回支路首段有功return S_begin[0]。}。double get_Pend(){//返回支路末端有功return S_end[0]。}。double get_Qbegin(){//返回支路首段無功return S_begin[1]。}。double get_Qend(){//返回支路末端無功return S_end[1]。}。double get_Pcost(){//返回支路有功損耗return S_begin[0]S_end[0]。}。double get_Qcost(){//返回支路無功損耗return S_begin[1]S_end[1]。}。line_part *get_lchild(){//返回支路無功損耗return lchild。}。line_part *get_rchild(){//返回支路無功損耗return rchild。}。}。void pass_U(line_part *a,line_part *b,line_part *c){(*c).U_begin=(*a).U_end。(*b).U_begin=(*a).U_end。}。void pass_w(line_part *a,line_part *b,line_part *c){(*a).S_end[0]=(*b).S_begin[0]+(*c).S_begin[0]+(*a).S_org[0]。(*a).S_end[1]=(*b).S_begin[1]+(*c).S_begin[1]+(*a).S_org[1]。}。void pass_U2(line_part *a,line_part *b){(*b).U_begin=(*a).U_end。}。void pass_w2(line_part *a,line_part *b){(*a).S_end[0]=(*b).S_begin[0]+(*a).S_org[0]。(*a).S_end[1]=(*b).S_begin[1]+(*a).S_org[1]。}。class BinTree{ public:friend void pass_U(line_part*,line_part*,line_part*)。//電壓傳遞函數(shù)pass_w(line_part*,line_part*,line_part*)。line_part(0,0,0,0,0,0,0,0,0,0,0)。root = aa。}vall,double u_end=0,double u_begin=0,double r=0, Pe=0,double Xe=0, P0=0,double Q0=0){ egin,r,x,b,Pe,Xe,K,P0,Q0)。private:*amp。t,double vall,double u_end,double u_begin,double r, Pe,double Xe, Q0)。}。friend void //功率傳遞函數(shù) BinTree(){line_part *aa=new line_part *Getroot(){return root。} void insertline_part(doubledouble x=0,double b=0,doubledouble K=0,double insertline_part(root,vall,u_end,u_b} void PreOrder(){PreOrder(root)。}。void PostOrder(){PostOrder(root)。}。void display(){display(root)。}。line_part *root。void insertline_part(line_partdouble x,double b,doubledouble K,double P0,double void PreOrder(line_part *amp。t)。void PostOrder(line_part *t)。void display(line_part *amp。t)。void BinTree::insertline_part(line_part *amp。t, double vall,double u_end=0,double u_begin=0,double r=0，double x=0,double b=0,double Pe=0,double Xe=0，double K=0,double P0=0,double Q0=0){//插入節(jié)點(diǎn)double ass=tget_val()。if(t==0||tget_val()t=new line_part(vall,u_end,u_begin,r,x,b,Pe,Xe,K,P0,Q0)。} else if(vallget_val()){insertline_part(tlchild, vall,u_end,u_begin,r,x,b,Pe,Xe,K,P0,Q0)。