【正文】 %第四列存貯支路的對(duì)地導(dǎo)納%第五列存貯變壓器的變比，注意支路為1 %第六列存貯支路的序號(hào)B2=input(‘n請(qǐng)輸入節(jié)點(diǎn)信息：B2=’)； %第一列為電源側(cè)的功率 %第二列為負(fù)荷側(cè)的功率 %第三列為該點(diǎn)的電壓值%第四列為該點(diǎn)的類型：1為PQ，2為PV節(jié)點(diǎn)，3為平衡節(jié)點(diǎn) A=input（‘n請(qǐng)輸入節(jié)點(diǎn)號(hào)及對(duì)地阻抗：A=’）； ip=input（‘n請(qǐng)輸入修正值：ip=’）； %ip為修正值)。For i=1:m S(i)=B2(i1)B2(i2)。[De,Di]=hxf(J,D,F,ph,n,no)。end S(ph)=sum(hh)*v(ph)。end result2=[B1(:,6),B1(:,1),B1(:,2),real(b’),imag(b’),real(c’),imag(c’), real(b’+c’),imag(b’+c’)]。V1=e(i)^2+f(i)^2。end end end Df=Df’。End ElseDX=[DX1,0,0]。%該子程序是用來求取jacci矩陣for i=1:n X3=0。X1=G(i)*f(i)B(i,j)*e(i)。p=2*i1。J(p,q)=X2。X2=C(i)+G(i,j)*e(i)+B(i,j)*f(i)。end J(p,q)=X1。J(:,no)=[]。//i作為整個(gè)程序的循環(huán)變量int N=Bus::ScanfBusNo()。//輸出支路參數(shù)Bus *bus=new Bus[N]。for(i=0。x=new double[NN]。//whether converbence看迭代是否結(jié)束if(icon==1){coutBus::JisuanJacAug(JacAug,X,bus,N)。i{int statemp,endtemp。cc=Complex::productComplex(bb , Complex::getconj(line[i].Y))。}coutBus::JisuanNodeScal(X,bus,N)。coutt(bus[i].)。}include “” include include include include“” include“” include“” using namespace std。//動(dòng)態(tài)分配支路結(jié)構(gòu)體Line::ScanfLineData(line)。//輸出結(jié)點(diǎn)參數(shù)Complex **X。for(i=0。//計(jì)算節(jié)點(diǎn)功率差值 Bus::PrintfNodeScal(X,bus,N)。goto LOOP。aa=Complex::productComplex(Complex::getconj(bus[statemp1].V), B)。dd=Complex::CaddC(aa,cc)。i{coutComplex::PrintfComplex(bus[i].V)。Complex::PrintfComplex(Complex::ComDivRea(line[i].etos,))。//求兩個(gè)復(fù)數(shù)差static Complex productComplex(Complex c1,Complex c2)。//顯示一個(gè)復(fù)數(shù)static void PrintfmultiComplex(Complex C,int N)。=1。=。=(**)/(pow(,2)+pow(,2))。=。}doubleComplex::getCompleximage(Complex c1)//取虛部 {return 。cout”。n=NN。//列n*2矩陣for(i=k。j{al=a[pivrow[k]][j]。i{aik=a[i][k]。j{sum +=a[i][j]*x[j]。====節(jié)點(diǎn)電壓===============發(fā)電機(jī)發(fā)出功率======節(jié)點(diǎn)10。0 0 節(jié)點(diǎn)6。00 0======線路傳輸功率========== 2to103to14to204to35to206to206to407to507to608to606to906to10011to9010to904to12013to12 14to1230to29 15to1228to8 28to6 請(qǐng)按任意鍵繼續(xù)... 20to19 20to10 2。0 00 節(jié)點(diǎn)5。x[pivcol[k][0]]=al。i){sum=0。a[i][k]=al。exit(0)。//行pivcol[k][0]=k。void gauss::gauss_slove(double **a,double *x,int NN){int n,i,j,k,*pivrow,**pivcol。(8)。return Node。} Complex Complex::getconj(Complex c1)//取共軛 {Complex Node。} Complex Complex::divideComplex(Complex c1,Complex c2)//復(fù)數(shù)除法 {Complex Node。}Complex Complex::subComplex(Complex c1,Complex c2)//復(fù)數(shù)減法 {Complex Node。=sqrt(*+*)。//求一個(gè)復(fù)數(shù)實(shí)部static double getCompleximage(Complex c1)。//0表示直角坐標(biāo)，1表示極坐標(biāo)static Complex CaddC(Complex c1,Complex c2)。endtemp=line[i].end。bus[i].V=Complex::Rec2Polar(bus[i].V)。bb=Complex::subComplex(Complex::getconj(bus[endtemp1].V), Complex::getconj(bus[statemp1].V))。statemp=line[i].start。//修正節(jié)點(diǎn)電壓// Bus::PrintfNodeV(bus,N)。//計(jì)算節(jié)點(diǎn)功率Bus::JisuanNodeScal(X,bus,N)。double **JacAug。}Bus::ScanfBusData(bus)。N)==0){return 0。Complex::PrintfComplex(Complex::ComDivRea(line[i].etos,))。i{coutComplex::PrintfComplex(bus[i].V)。dd=Complex::CaddC(aa,cc)。aa=Complex::productComplex(Complex::getconj(bus[statemp1].V), B)。goto LOOP。//計(jì)算節(jié)點(diǎn)功率差值Bus::PrintfNodeScal(X,bus,N)。for(i=0。//輸出結(jié)點(diǎn)參數(shù)Complex **X。//動(dòng)態(tài)分配支路結(jié)構(gòu)體Line::ScanfLineData(line)。謝謝include “” include include include include“” include“” include“” using namespace std。J(no:)=[]。P=2*i1。X1=D(i)+G(i,j)*f(i)B(i,j)*e(i)。m=p+1。X4=X1。ifJ=amp。End End case 2 Function for j=1:n J=jacci(Y,G,B,PQ,e,f,V,C,D,B2,n,ph,no)X1=G(i,j)*f(i)B(i,j)*e(i)。DX(no+1)=0。p=p+1。endP1=C(i)*e(i)+D(i)*f(i)。b(i)=v(B1(i1))*a(i)j*B1(i4)*v(B1(i))^2/2。end v=e+f*j。[C,D,DF]=xxf(G,B,e,f,P,Q,n,B2,ph,V,no)。for i=1：n End e(i)=real(B2(i3))?！尽俊尽縁unction tisco %這是一個(gè)電力系統(tǒng)潮流計(jì)算的程序 n=input（‘n請(qǐng)輸入節(jié)點(diǎn)數(shù)：n=’）； m=input(‘請(qǐng)輸入支路數(shù)：m=’)。Q(i)=U(i)*U(j)*(G(i,j)*sin(a(i)a(j))B(i,j)*cos(a(i)a(j)))+Q(i)。M,L]% J為雅克比矩陣x1=((inv(J))*x2)。L(i,j)=H(i,j)。oQ(i)=0。x1=ones(8,1)。S(3)=。for i=1:5for j=i:5y(j,i)=y(i,j)。% 輸入原始數(shù)據(jù)，求節(jié)點(diǎn)導(dǎo)納矩陣y(1,2)=1/(+)。計(jì)算精度要求各節(jié)點(diǎn)電壓修正量不大于106。outfile()。ee=()。如下“()。char b。量double *A=new double[10*M]。amp。t{tlchildSbegin()。void BinTree::PostOrder(line_part *t){//后序遍歷二叉樹修改潮流if(tget_val()0amp。} } else { if(tget_rchild()get_val()0amp。amp。outfile1”()。double U。/////////////////////////////////////////////////////////////////////////////////////////////////////////////// void BinTree::display(line_part *amp。void display(line_part *amp。void PostOrder(){PostOrder(root)。root = aa。void pass_U2(line_part *a,line_part *b){(*b).U_begin=(*a).U_end。}。}。}。U[1]=U_zong。//U_heng即是△u2 0])/U_begin。//功率傳遞函數(shù)friend void pass_U2(line_part*a,line_part*b)。S0[1]=Q0。X[1]=x。S0[1]=0。X[1]=0。//變壓器勵(lì)磁損耗：S0[0]有功，S0[1]無功double S_org[2]。}。}。//支路首段功率：S_begin[0]有功，S_begin[1]無功double S0[2]。********************************************************************************************************************************************************* 程序中變量定義說明：類中定義的變量class line_part{//定義支路類double U[2]。第一篇：潮流計(jì)算電力系統(tǒng)大作業(yè)(C++編寫)程序設(shè)計(jì)所涉及二叉樹解釋任意一棵樹或一個(gè)森林都能唯一地對(duì)應(yīng)一棵二叉樹，由此而編寫本程序。******************************************************************************************************************************************************* 程序功能說明：本程序可以計(jì)算任意長度線型開始網(wǎng)絡(luò)潮流；支持多電壓等級(jí)下的計(jì)算；可在除供電節(jié)點(diǎn)外任意節(jié)點(diǎn)引出負(fù)載；*******************************************************************************************************************************************************輸入格式說明：以支路為基本單位，按潮流方向輸入數(shù)據(jù)：以下例示意：//支路個(gè)數(shù)1,110,118,0,0,1,0,0//線路等效的支路2,11,110,0,40,30,10,//變壓器等效的支路.//按此方式知道輸入所有的支路.//輸入從上到下的順序?yàn)槌绷髟?//線型開式網(wǎng)絡(luò)中的流動(dòng)方向/*第一行輸入支路個(gè)數(shù)，回車第二行至后輸入各個(gè)支路參數(shù)，回車分隔不同支路；各行輸入的支路參數(shù)順序是：支路編號(hào)，末端電壓，始端電壓，線路等效電阻，線路等效感抗，線路等效容納，末端輸入有功，末端輸入無功，變比，變壓器有功勵(lì)磁損耗，變壓器無功勵(lì)磁損耗對(duì)于線路等效電路：變壓器有功/無功損耗輸入零，變比輸入1；對(duì)于變壓器等效電路：所有的參數(shù)均歸算至高壓側(cè)，Rt,Xt對(duì)應(yīng)輸入線路等效電阻/感抗的位置，線路等效容抗為零；對(duì)于個(gè)節(jié)點(diǎn)的引出負(fù)荷：輸入至以此節(jié)點(diǎn)為末節(jié)點(diǎn)的支路的末端輸入有功/無功部分；*/ *****************************************************************