【正文】 %第四列存貯支路的對地導納%第五列存貯變壓器的變比，注意支路為1 %第六列存貯支路的序號B2=input(‘n請輸入節(jié)點信息：B2=’)； %第一列為電源側的功率 %第二列為負荷側的功率 %第三列為該點的電壓值%第四列為該點的類型：1為PQ，2為PV節(jié)點，3為平衡節(jié)點 A=input（‘n請輸入節(jié)點號及對地阻抗：A=’）； ip=input（‘n請輸入修正值：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作為整個程序的循環(huán)變量int N=Bus::ScanfBusNo()。//輸出支路參數(shù)Bus *bus=new Bus[N]。for(i=0。x=new double[NN]。//whether converbence看迭代是否結束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。//動態(tài)分配支路結構體Line::ScanfLineData(line)。//輸出結點參數(shù)Complex **X。for(i=0。//計算節(jié)點功率差值 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,))。//求兩個復數(shù)差static Complex productComplex(Complex c1,Complex c2)。//顯示一個復數(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é)點電壓===============發(fā)電機發(fā)出功率======節(jié)點10。0 0 節(jié)點6。00 0======線路傳輸功率========== 2to103to14to204to35to206to206to407to507to608to606to906to10011to9010to904to12013to12 14to1230to29 15to1228to8 28to6 請按任意鍵繼續(xù)... 20to19 20to10 2。0 00 節(jié)點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)//復數(shù)除法 {Complex Node。}Complex Complex::subComplex(Complex c1,Complex c2)//復數(shù)減法 {Complex Node。=sqrt(*+*)。//求一個復數(shù)實部static double getCompleximage(Complex c1)。//0表示直角坐標，1表示極坐標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é)點電壓// Bus::PrintfNodeV(bus,N)。//計算節(jié)點功率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。//計算節(jié)點功率差值Bus::PrintfNodeScal(X,bus,N)。for(i=0。//輸出結點參數(shù)Complex **X。//動態(tài)分配支路結構體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 %這是一個電力系統(tǒng)潮流計算的程序 n=input（‘n請輸入節(jié)點數(shù)：n=’）； m=input(‘請輸入支路數(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é)點導納矩陣y(1,2)=1/(+)。計算精度要求各節(jié)點電壓修正量不大于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。//變壓器勵磁損耗：S0[0]有功，S0[1]無功double S_org[2]。}。}。//支路首段功率：S_begin[0]有功，S_begin[1]無功double S0[2]。********************************************************************************************************************************************************* 程序中變量定義說明：類中定義的變量class line_part{//定義支路類double U[2]。第一篇：潮流計算電力系統(tǒng)大作業(yè)(C++編寫)程序設計所涉及二叉樹解釋任意一棵樹或一個森林都能唯一地對應一棵二叉樹，由此而編寫本程序。******************************************************************************************************************************************************* 程序功能說明：本程序可以計算任意長度線型開始網(wǎng)絡潮流；支持多電壓等級下的計算；可在除供電節(jié)點外任意節(jié)點引出負載；*******************************************************************************************************************************************************輸入格式說明：以支路為基本單位，按潮流方向輸入數(shù)據(jù)：以下例示意：//支路個數(shù)1,110,118,0,0,1,0,0//線路等效的支路2,11,110,0,40,30,10,//變壓器等效的支路.//按此方式知道輸入所有的支路.//輸入從上到下的順序為潮流在.//線型開式網(wǎng)絡中的流動方向/*第一行輸入支路個數(shù)，回車第二行至后輸入各個支路參數(shù)，回車分隔不同支路；各行輸入的支路參數(shù)順序是：支路編號，末端電壓，始端電壓，線路等效電阻，線路等效感抗，線路等效容納，末端輸入有功，末端輸入無功，變比，變壓器有功勵磁損耗，變壓器無功勵磁損耗對于線路等效電路：變壓器有功/無功損耗輸入零，變比輸入1；對于變壓器等效電路：所有的參數(shù)均歸算至高壓側，Rt,Xt對應輸入線路等效電阻/感抗的位置，線路等效容抗為零；對于個節(jié)點的引出負荷：輸入至以此節(jié)點為末節(jié)點的支路的末端輸入有功/無功部分；*/ *****************************************************************