電力系統(tǒng)潮流計(jì)算程序

1、電力系統(tǒng)潮流計(jì)算c語(yǔ)言程序，兩行，大家可以看看，仔細(xì)研究，然后在這個(gè)基礎(chǔ)上修改。謝謝#include stdafx.h#include #include #include #includeComplex.h#includewanjing.h#includegauss.husing namespace std;int _tmain(int argc, _TCHAR* argv)int i;/i作為整個(gè)程序的循環(huán)變量int N=Bus:ScanfBusNo();/輸入節(jié)點(diǎn)個(gè)數(shù)int L=Line:ScanflineNo();/輸入支路個(gè)數(shù)if(L&N)=0)return 0;/如果找不到兩個(gè)文件中

2、的任意一個(gè)，退出Line *line=new LineL;/動(dòng)態(tài)分配支路結(jié)構(gòu)體Line:ScanfLineData(line);/輸入支路參數(shù)Line:PrintfLineData(line,L);/輸出支路參數(shù)Bus *bus=new BusN;/動(dòng)態(tài)分配結(jié)點(diǎn)結(jié)構(gòu)體for(int i=0;iN;i+)busi.Sdelta.real=0;busi.Sdelta.image=0;Bus:ScanfBusData(bus);/輸入節(jié)點(diǎn)參數(shù)Bus:PrintfBusData(bus,N);/輸出結(jié)點(diǎn)參數(shù)Complex *X; X=new Complex *N; for(i=0;iN;i+) Xi=

3、new ComplexN;/動(dòng)態(tài)分配二維數(shù)組作節(jié)點(diǎn)導(dǎo)納陣Bus:JisuanNodeDnz(X,line,bus,L,N);/計(jì)算節(jié)點(diǎn)導(dǎo)納矩陣Bus:PrintfNodeDnz(X,N);/輸出節(jié)點(diǎn)導(dǎo)納矩陣int NN=(N-1)*2;double *JacAug; JacAug=new double *NN; for(i=0;iNN;i+) JacAugi=new doubleNN+1;/為雅可比增廣動(dòng)態(tài)分配一個(gè)二維數(shù)組double *x;x=new doubleNN;int count=1;LOOP:Bus:JisuanNodeI(X,bus,N);/計(jì)算節(jié)點(diǎn)注入電流Bus:JisuanN

4、odeScal(X,bus,N);/計(jì)算節(jié)點(diǎn)功率Bus:JisuanNodeScal(X,bus,N);/計(jì)算節(jié)點(diǎn)功率Bus:JisuanNodeSdelta(bus,N);/計(jì)算節(jié)點(diǎn)功率差值Bus:PrintfNodeScal(X,bus,N);/輸出節(jié)點(diǎn)功率差值int icon=wehcon1(bus,N);/whether converbence看迭代是否結(jié)束if(icon=1)couticon=icon進(jìn)行第count次迭代endl;Bus:JisuanJacAug(JacAug,X,bus,N);/計(jì)算雅可比增廣矩陣/Bus:PrintfJacAug(JacAug,N);gauss:

5、gauss_slove(JacAug,x,NN);/解方程組求出電壓差值Bus:ReviseNodeV(bus,x,N);/修正節(jié)點(diǎn)電壓/Bus:PrintfNodeV(bus,N);count+;goto LOOP;else for(i=0;iL;i+)/計(jì)算輸電線路上的傳輸功率int statemp,endtemp;Complex aa,bb,cc,dd,B;B.real=0;B.image=-linei.B;statemp=linei.start;endtemp=linei.end;aa=Complex:productComplex (Complex:getconj(busstatemp

6、-1.V), B);bb=Complex:subComplex (Complex:getconj(busstatemp-1.V), Complex:getconj(busendtemp-1.V);cc=Complex:productComplex (bb , Complex:getconj(linei.Y);dd=Complex:CaddC(aa,cc);linei.stoe=Complex:productComplex(busstatemp-1.V,dd);aa=Complex:productComplex (Complex:getconj(busendtemp-1.V), B);bb=Co

7、mplex:subComplex (Complex:getconj(busendtemp-1.V), Complex:getconj(busstatemp-1.V);cc=Complex:productComplex (bb , Complex:getconj(linei.Y);dd=Complex:CaddC(aa,cc);linei.etos=Complex:productComplex(busendtemp-1.V,dd);couticon=icon,迭代結(jié)束。endl;Bus:JisuanNodeScal(X,bus,N);/計(jì)算節(jié)點(diǎn)功率for(i=0;iN;i+)busi.Scal.

8、real = busi.Scal.real + busi.Load.real;/發(fā)電機(jī)功率=注入功率+負(fù)荷功率busi.Scal.image= busi.Scal.image+ busi.Load.image;busi.V=Complex:Rec2Polar(busi.V);cout=節(jié)點(diǎn)電壓=發(fā)電機(jī)發(fā)出功率=endl;for(i=0;iN;i+)cout節(jié)點(diǎn)(i+1)t;Complex:PrintfComplex(busi.V);coutt(busi.Scal.real);coutt(busi.Scal.image);coutendl;cout=線路傳輸功率=endl;for(i=0;iL;

9、i+)/輸出線路上的傳輸功率int statemp,endtemp;statemp=linei.start;endtemp=linei.end;coutstatemptoendtempt;Complex:PrintfComplex(Complex:ComDivRea(linei.stoe,0.01);Complex:PrintfComplex(Complex:ComDivRea(linei.etos,0.01);coutendl;return 0;#include stdafx.h#include #include #include #includeComplex.h#includewanji

10、ng.h#includegauss.husing namespace std;int _tmain(int argc, _TCHAR* argv)int i;/i作為整個(gè)程序的循環(huán)變量int N=Bus:ScanfBusNo();/輸入節(jié)點(diǎn)個(gè)數(shù)int L=Line:ScanflineNo();/輸入支路個(gè)數(shù)if(L&N)=0)return 0;/如果找不到兩個(gè)文件中的任意一個(gè)，退出Line *line=new LineL;/動(dòng)態(tài)分配支路結(jié)構(gòu)體Line:ScanfLineData(line);/輸入支路參數(shù)Line:PrintfLineData(line,L);/輸出支路參數(shù)Bus *bus=n

11、ew BusN;/動(dòng)態(tài)分配結(jié)點(diǎn)結(jié)構(gòu)體for(int i=0;iN;i+)busi.Sdelta.real=0;busi.Sdelta.image=0;Bus:ScanfBusData(bus);/輸入節(jié)點(diǎn)參數(shù)Bus:PrintfBusData(bus,N);/輸出結(jié)點(diǎn)參數(shù)Complex *X; X=new Complex *N; for(i=0;iN;i+) Xi=new ComplexN;/動(dòng)態(tài)分配二維數(shù)組作節(jié)點(diǎn)導(dǎo)納陣Bus:JisuanNodeDnz(X,line,bus,L,N);/計(jì)算節(jié)點(diǎn)導(dǎo)納矩陣Bus:PrintfNodeDnz(X,N);/輸出節(jié)點(diǎn)導(dǎo)納矩陣int NN=(N-1)*

12、2;double *JacAug; JacAug=new double *NN; for(i=0;iNN;i+) JacAugi=new doubleNN+1;/為雅可比增廣動(dòng)態(tài)分配一個(gè)二維數(shù)組double *x;x=new doubleNN;int count=1;LOOP:Bus:JisuanNodeI(X,bus,N);/計(jì)算節(jié)點(diǎn)注入電流Bus:JisuanNodeScal(X,bus,N);/計(jì)算節(jié)點(diǎn)功率Bus:JisuanNodeScal(X,bus,N);/計(jì)算節(jié)點(diǎn)功率Bus:JisuanNodeSdelta(bus,N);/計(jì)算節(jié)點(diǎn)功率差值Bus:PrintfNodeScal(X

13、,bus,N);/輸出節(jié)點(diǎn)功率差值int icon=wehcon1(bus,N);/whether converbence看迭代是否結(jié)束if(icon=1)couticon=icon進(jìn)行第count次迭代endl;Bus:JisuanJacAug(JacAug,X,bus,N);/計(jì)算雅可比增廣矩陣/Bus:PrintfJacAug(JacAug,N);gauss:gauss_slove(JacAug,x,NN);/解方程組求出電壓差值Bus:ReviseNodeV(bus,x,N);/修正節(jié)點(diǎn)電壓/Bus:PrintfNodeV(bus,N);count+;goto LOOP;else fo

14、r(i=0;iL;i+)/計(jì)算輸電線路上的傳輸功率int statemp,endtemp;Complex aa,bb,cc,dd,B;B.real=0;B.image=-linei.B;statemp=linei.start;endtemp=linei.end;aa=Complex:productComplex (Complex:getconj(busstatemp-1.V), B);bb=Complex:subComplex (Complex:getconj(busstatemp-1.V), Complex:getconj(busendtemp-1.V);cc=Complex:product

15、Complex (bb , Complex:getconj(linei.Y);dd=Complex:CaddC(aa,cc);linei.stoe=Complex:productComplex(busstatemp-1.V,dd);aa=Complex:productComplex (Complex:getconj(busendtemp-1.V), B);bb=Complex:subComplex (Complex:getconj(busendtemp-1.V), Complex:getconj(busstatemp-1.V);cc=Complex:productComplex (bb , C

16、omplex:getconj(linei.Y);dd=Complex:CaddC(aa,cc);linei.etos=Complex:productComplex(busendtemp-1.V,dd);couticon=icon,迭代結(jié)束。endl;Bus:JisuanNodeScal(X,bus,N);/計(jì)算節(jié)點(diǎn)功率for(i=0;iN;i+)busi.Scal.real = busi.Scal.real + busi.Load.real;/發(fā)電機(jī)功率=注入功率+負(fù)荷功率busi.Scal.image= busi.Scal.image+ busi.Load.image;busi.V=Comp

17、lex:Rec2Polar(busi.V);cout=節(jié)點(diǎn)電壓=發(fā)電機(jī)發(fā)出功率=endl;for(i=0;iN;i+)cout節(jié)點(diǎn)(i+1)t;Complex:PrintfComplex(busi.V);coutt(busi.Scal.real);coutt(busi.Scal.image);coutendl;cout=線路傳輸功率=endl;for(i=0;iL;i+)/輸出線路上的傳輸功率int statemp,endtemp;statemp=linei.start;endtemp=linei.end;coutstatemptoendtempt;Complex:PrintfComplex(

18、Complex:ComDivRea(linei.stoe,0.01);Complex:PrintfComplex(Complex:ComDivRea(linei.etos,0.01);coutendl;return 0;#include using namespace std;class Complex/定義復(fù)數(shù)類(lèi)public:double real;double image;int RecPolar;/0表示直角坐標(biāo)，1表示極坐標(biāo)static Complex CaddC(Complex c1,Complex c2); /求兩個(gè)復(fù)數(shù)和static Complex subComplex(Comp

19、lex c1,Complex c2);/求兩個(gè)復(fù)數(shù)差static Complex productComplex(Complex c1,Complex c2); /求兩個(gè)復(fù)數(shù)積 static Complex divideComplex(Complex c1,Complex c2); /求兩個(gè)復(fù)數(shù)商static Complex ComDivRea(Complex c1,double r2);/除數(shù)static Complex getconj(Complex c1); /求一個(gè)復(fù)數(shù)共軛static Complex getinverse(Complex c1);/取倒數(shù)static double ge

20、tComplexReal(Complex c1); /求一個(gè)復(fù)數(shù)實(shí)部 static double getCompleximage(Complex c1);/求一個(gè)復(fù)數(shù)虛部static void PrintfComplex(Complex c1);/顯示一個(gè)復(fù)數(shù)static void PrintfmultiComplex(Complex C,int N);/顯示多個(gè)復(fù)數(shù)static void zeroComplex(Complex c1);/將復(fù)數(shù)復(fù)零static Complex Rec2Polar(Complex c1);/取極坐標(biāo)Complex()RecPolar=0; Complex Co

21、mplex:Rec2Polar(Complex c1)/極坐標(biāo)表示Complex Node;Node.real=sqrt(c1.real*c1.real+c1.image*c1.image);Node.image=atan2(c1.image,c1.real)*180/3.1415926;Node.RecPolar=1;return Node;Complex Complex:CaddC(Complex c1,Complex c2) /復(fù)數(shù)加法 Complex Node; Node.real=c1.real+c2.real; Node.image=c1.image+c2.image; retur

22、n Node; Complex Complex:subComplex(Complex c1,Complex c2) /復(fù)數(shù)減法 Complex Node; Node.real=c1.real-c2.real; Node.image=c1.image-c2.image; return Node; Complex Complex:productComplex(Complex c1,Complex c2) /復(fù)數(shù)乘法 Complex Node; Node.real=c1.real*c2.real-c1.image*c2.image; Node.image=c1.image*c2.real+c2.im

23、age*c1.real; return Node; Complex Complex:divideComplex(Complex c1,Complex c2) /復(fù)數(shù)除法 Complex Node; Node.real=(c1.real*c2.real+c1.image*c2.image)/(pow(c2.real,2)+pow(c2.image,2);Node.image=(c1.image*c2.real-c1.real*c2.image)/(pow(c2.real,2)+pow(c2.image,2);return Node; Complex Complex:ComDivRea(Compl

24、ex c1,double r1)/復(fù)數(shù)除數(shù)Complex Node;Node.real=c1.real/(r1);Node.image=c1.image/(r1);return Node;Complex Complex:getconj(Complex c1) /取共軛 Complex Node; Node.real=c1.real;Node.image=-c1.image; return Node; Complex Complex:getinverse(Complex c1)/取倒數(shù)Complex Node;Node.real=1;Node.image=0;Node=(Complex:divi

25、deComplex(Node,c1);return Node;double Complex:getComplexReal(Complex c1) /取實(shí)部 return c1.real; double Complex:getCompleximage(Complex c1) /取虛部 return c1.image; void Complex:PrintfComplex(Complex c1)/按直角坐標(biāo)輸出if(c1.RecPolar=0)cout.precision(6);cout.width(8);cout.setf(ios:right);coutc1.real ;cout.precisi

26、on(6);cout.width(8);cout.setf(ios:left);coutc1.imagei ;elsecoutc1.realtc1.image。t;/按極坐標(biāo)輸出void Complex:zeroComplex(Complex c1)/清零c1.real=0;c1.image=0;class gausspublic:static void gauss_slove(double *a,double *x,int NN);static void gauss_output();void gauss:gauss_slove(double *a,double *x,int NN)int

27、n,i,j,k,*pivrow,*pivcol;double eps,pivot,sum,aik,al;n=NN;pivrow=new intn;pivcol=new int *n;for(i=0;in;i+)pivcoli=new int2;eps=0.00001;for(k=0;k(n-1);k+)pivot= fabs(akk);pivrowk=k;/行pivcolk0=k;pivcolk1=k;/列n*2矩陣for(i=k;in;i+)/找最大值放到主對(duì)角線for(j=k;jn;j+)if (pivotfabs(aij)pivot=fabs(aij);pivrowk=i;/行pivco

28、lk1=j;/列if(pivoteps)couterrorendl;getchar();exit(0);if(pivrowk!=k)/行變換for(j=k;j(n+1);j+)al=apivrowkj;apivrowkj=akj;akj=al;if(pivcolk1!=k)/列變換for(i=0;in;i+)al=aipivcolk1;aipivcolk1=aik;aik=al;if (k!=(n-1)/將矩陣化為上三角形式for(i=(k+1);in;i+)aik=aik;for(j=k;j=0;i-)sum=0;for(j=(i+1);j=0;k-)al=xpivcolk1;xpivcol

29、k1=xpivcolk0;xpivcolk0=al;cout節(jié)點(diǎn)電壓修正量endl=endl;for(i=0;in;i+)coutxiendl;coutendlendl;=節(jié)點(diǎn)功率計(jì)算值=0.935261 -0.1590480.573909 0.0789973-0.00289889 -0.00796623-0.0791247 -0.0168362-0.436255 -0.05803920.0359139 -0.0106592-0.229118 -0.0885419-0.136179 -0.1482070.0446243 0.0111298-0.0223764 -0.006957750.1827

30、09 0.016894-0.0310701 -0.04020510.156702 -0.0355909-0.0668055 -0.00703229-0.0886481 -0.0129814-0.0390805 -0.0135062-0.1023 -0.0460568-0.0342827 -0.00382402-0.102896 -0.0184062-0.0237482 -0.00347744-0.198318 -0.0869347-5.24266e-015 -9.48909e-015-0.0354071 -0.0110778-0.0925078 -0.0538236-1.05629e-015

31、-7.53784e-016-0.0391348 -0.01680970.014529 7.049e-0050.00158644 -0.00146487-0.0258771 -0.00458276-0.109514 0.00251645icon=1進(jìn)行第2次迭代節(jié)點(diǎn)電壓修正量=-0.00164889-0.000540034-0.00261067-0.00532027-0.00235315-0.00600971-0.00189677-0.00643874-0.0023631-0.00650659-0.00170949-0.0074907-0.00164545-0.00485415-0.004939

32、77-0.0119042-0.00331285-0.0175611-0.00207908-0.00336375-0.00530645-0.0147816-0.000326161-0.00640487-0.00251701-0.0169829-0.00175286-0.0174333-0.00307326-0.0163809-0.00232773-0.01758068.74293e-005-0.01920180.000558996-0.0197776-0.000247851-0.0193784-0.00115346-0.0185848-0.00127275-0.0186244-0.0001010

33、8-0.01889660.000553585-0.0200901-3.76315e-005-0.02083030.00308341-0.0219386-0.00195916-0.0205356-0.00184757-0.00764010.00197593-0.02455340.00434657-0.027534=節(jié)點(diǎn)功率計(jì)算值=0.98623 -0.1341630.583136 0.166278-0.0239063 -0.0119192-0.076014 -0.0160104-0.441997 -0.07502850.000250012 3.72542e-005-0.228052 -0.108

34、844-0.100078 -0.1056340.000410707 0.000378067-0.057497 -0.01958790.200039 0.0582563-0.111173 -0.07398460.199792 0.0227868-0.0621041 -0.0158709-0.0821379 -0.0248173-0.0350785 -0.0179447-0.0902383 -0.0578368-0.0320461 -0.00890719-0.0951562 -0.0337091-0.0220362 -0.00693706-0.175458 -0.1116014.72557e-01

35、5 1.21429e-014-0.0320661 -0.0159145-0.0871134 -0.0667319-7.03489e-017 9.24355e-016-0.0350769 -0.02285920.000273455 7.10354e-0051.51804e-005 -6.6188e-006-0.0240417 -0.00889343-0.10604 -0.0184098icon=1進(jìn)行第3次迭代節(jié)點(diǎn)電壓修正量=-2.67079e-005-2.30128e-006-2.20543e-005-6.00686e-005-2.33043e-005-6.85601e-005-3.22294

36、e-005-2.61107e-005-2.80198e-005-6.6167e-005-2.34528e-005-5.66132e-005-4.4646e-005-1.74668e-005-4.50947e-005-0.000181763-3.81763e-006-0.000286581-6.68993e-005-1.28441e-005-5.17172e-005-0.000223284-4.54717e-005-2.47586e-0054.32335e-007-0.0002584941.82635e-005-0.000272051-6.95195e-006-0.0002519691.1131

37、8e-005-0.0002794185.74737e-005-0.0003073686.86998e-005-0.0003202745.38112e-005-0.000314473.59531e-005-0.000304943.37607e-005-0.0003074495.26532e-005-0.0003107216.92761e-005-0.0003503735.60942e-005-0.000409770.000123641-0.0004402591.36149e-005-0.000426973-1.70227e-005-9.37794e-0050.000113675-0.000544

38、0110.000176034-0.000636202=節(jié)點(diǎn)功率計(jì)算值=0.986878 -0.1339790.583 0.167193-0.024 -0.012-0.076 -0.016-0.442 -0.07486061.43501e-008 1.07366e-008-0.228 -0.109-0.0999999 -0.1040494.51318e-008 8.98835e-008-0.0579999 -0.01999990.2 0.0591018-0.112 -0.07499970.2 0.0242519-0.062 -0.016-0.082 -0.025-0.035 -0.018-0.0

39、900001 -0.058-0.032 -0.00899997-0.095 -0.0339999-0.022 -0.00699998-0.175 -0.112-6.07156e-015 -1.19217e-014-0.032 -0.016-0.087 -0.06699997.03078e-017 -9.23979e-016-0.035 -0.02299991.09492e-007 4.45699e-0081.54958e-009 -2.01531e-010-0.024 -0.00899994-0.106 -0.0189996icon=0,迭代結(jié)束。=節(jié)點(diǎn)電壓=發(fā)電機(jī)發(fā)出功率=節(jié)點(diǎn)1 1.05

40、0。 98.6878 -13.3979節(jié)點(diǎn)2 1.045 -1.846。 80 29.4193節(jié)點(diǎn)3 1.02384 -3.83352。 0 0節(jié)點(diǎn)4 1.01637 -4.55698。 0 0節(jié)點(diǎn)5 1.01 -6.48617。 50 11.5139節(jié)點(diǎn)6 1.01332 -5.38073。 0 0節(jié)點(diǎn)7 1.00489 -6.38368。 0 0節(jié)點(diǎn)8 1.01 -5.62974。 20 19.5951節(jié)點(diǎn)9 1.03905 -6.78143。 0 0節(jié)點(diǎn)10 1.03595 -8.69362。 0 0節(jié)點(diǎn)11 1.05 -4.5962。 20 5.91018節(jié)點(diǎn)12 1.04711 -7

41、.80323。 0 0節(jié)點(diǎn)13 1.05 -6.34392。 20 2.42519節(jié)點(diǎn)14 1.03242 -8.7401。 0 0節(jié)點(diǎn)15 1.02788 -8.86784。 0 0節(jié)點(diǎn)16 1.03458 -8.45044。 0 0節(jié)點(diǎn)17 1.03051 -8.83678。 0 0節(jié)點(diǎn)18 1.01845 -9.5141。 0 0節(jié)點(diǎn)19 1.01604 -9.70326。 0 0節(jié)點(diǎn)20 1.02022 -9.50938。 0 0節(jié)點(diǎn)21 1.0237 -9.17478。 0 0節(jié)點(diǎn)22 1.02432 -9.17024。 0 0節(jié)點(diǎn)23 1.01802 -9.36719。 0 0節(jié)點(diǎn)

42、24 1.01339 -9.68362。 0 0節(jié)點(diǎn)25 1.01216 -9.68486。 0 0節(jié)點(diǎn)26 0.994393 -10.1089。 0 0節(jié)點(diǎn)27 1.02012 -9.42025。 0 0節(jié)點(diǎn)28 1.00992 -5.86244。 0 0節(jié)點(diǎn)29 1.00022 -10.6579。 0 0節(jié)點(diǎn)30 0.988705 -11.5464。 0 0=線路傳輸功率=2to1 -57.7373 5.41674i 58.3454 -15.1827i3to1 -39.659 -7.75964i 40.3424 1.78481i4to2 -30.87 -9.74186i 31.4153 3.58352i4to3 -37.0772 -7.78596i 37.259 6.55964i5to2 -44.3717 -9.78456i 45.2968 4.84242i6to2 -38.4766 -8.22625i 39.3252 2.87667i6to4 -34.946 1.92384i 35.0885 -3.28202i7to5 -0.16304 -6.41767i 0.1