潮流计算程序有误差，不知道是哪里出了问题。求大神解救

lllmos

lllmos

潮流程序计算结果又误差是什么问题造成的啊？

玉门关山

你给这么些信息量，，真得靠大神来答了。。

lllmos

回复 3# 玉门关山

+ ^  X: k) O# g2 c$ Q; L  Y
0 w! H( x2 m% {; i7 Q5 Q    因为我不知道问题出在哪里了。。
do
        {
- P/ C% E( ~# \4 P" R                //求解不平衡量
* F9 b7 u5 I- b3 l; W% |                for(i=0;i<nB;i++)
( W, \- M% Q. s/ V- V4 x& z3 ~                {
                        if(sB.Type!=2)//假如不是平衡节点
                        {        
                                DP=sB.GenP-sB.LoadP;
                                DQ=sB.GenQ-sB.LoadQ;
8 V+ N* f4 X/ f                               
1 M3 S$ u  T) C$ o( L                                for(j=0;j<nB;j++)
                                {
3 b) f8 X- v$ E# p5 G; d" V9 C6 D                                        A=sB.Phase-sB[j].Phase;
2 e$ I- }- L& {( a, M                                        DP-=sB.Volt*(sB[j].Volt*(g[j]*cos(A)+b[j]*sin(A)));
                                       
. K' V# t, n: c% q; n- B4 C                                        if(sB.Type==0)//PQ节点
+ l0 y+ C! y6 C; f( {& q# }. ~  q                                                DQ-=sB.Volt*(sB[j].Volt*(g[j]*sin(A)-b[j]*cos(A)));
                               
                                        else if(sB.Type==1)//PV节点
                                                DQ=0;
                                }
; @* Y- Q( b# ]+ b& _/ w                        }
$ ?6 u" E  ~# U2 f0 H/ p5 R                        else if(sB.Type==2)//平衡节点
1 E% ]. o6 q$ [+ w$ b9 A% q4 X- O                                DP=DQ=0;
0 L& T: \& N0 C8 D9 J/ |                }
# d; L! u5 T: d! ^+ w                //for(i=0;i<nB;i++)
, y! I: F" |- S+ I                //         printf("DP[%d]===%f,DQ[%d]===%f\n",i,DP,i,DQ);

9 b5 B$ @( K  @5 s- F            //求解修正方程
0 K8 Z3 [  t& r- n* l                for(i=0;i<nB-1;i++)
                    AA1=DP[i+1]/sB[i+1].Volt;
7 b3 _( H4 p( j8 V2 A                for(i=0;i<nB-1-count_PVnode;i++)
                        AA2=DQ[i+1+count_PVnode]/sB[i+1+count_PVnode].Volt;
5 ^7 j* _4 M: V                calculate_gaosi((double **)b1,BB1,AA1,NBUS-1);//AA是不平衡量，BB是解向量
                calculate_gaosi((double **)b2,BB2,AA2,NBUS-1);

                max1=fabs(AA1[0]);
$ ^: I) o% K% q; c! }# R2 t0 c3 K                for(i=1;i<nB-1;i++)
                        if(max1<fabs(AA1))       
                                max1=fabs(AA1);
$ s" B/ h( s0 x5 \4 Q* m3 \                        max2=fabs(AA2[0]);
                        for(i=1;i<nB-1-count_PVnode;i++)
                                if(max2<fabs(AA2))       
                                max2=fabs(AA2);
; N2 `! m2 f2 ~- x; Q6 p                for(i=0;i<nB-1;i++)
                        sB[i+1].Phase+=BB1/sB.Volt;
7 ~" M; X0 S9 N9 G                for(i=0;i<nB-1-count_PVnode;i++)
. G7 U1 N. v" L/ q                        sB[i+1+count_PVnode].Volt+=BB2;
# ^* V% `$ f# _$ s                for(i=0;i<nB;i++)
                {
' A+ E3 T! p+ r                        printf("sB[%d].Volt=%f,sB[%d].Phase=%f\n",i,sB.Volt,i,sB.Phase*180/PI);
- }4 r" z  m  i4 v                       
                }
" q9 X, }& K; Z( g! ~+ n3 [$ t; B                printf("\n");
9 N7 _! v( a8 v. C) L- d                        ci++;
        }
! n0 i$ l* i1 ^: u        while(fabs(max1)>0.00001&&fabs(max2)>0.00001&&ci<40);
: }9 D* V# N! m# C这是我求潮流的程序，用的PQ分解法，最后得到的结果是只能精确到小数点后第二位，第三位就不对了。

electri

别急躁，好好找找一定能好到原因。

ordinarymengya

回复 2# lllmos


   跟你的收敛条件设置有关系啊。

ordinarymengya

回复 4# lllmos

9 T, i4 h  O1 A9 C7 e! H4 w建议楼主换几个软件试试，SPS,PSSE，PowerSystemWorld，特别是PSSE算潮流比较专业。