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

lllmos

lllmos

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

玉门关山

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

lllmos

回复 3# 玉门关山


+ A0 P5 M/ [4 y3 z- D6 i    因为我不知道问题出在哪里了。。
0 [; s. s7 \3 Q* t. Qdo
7 d4 N  @$ Z# L; d. Q7 @9 H        {
( E: W% V& ?% W8 v8 E0 @3 l7 h                //求解不平衡量
                for(i=0;i<nB;i++)
: u7 h+ v( l' H. b6 l: Q                {
                        if(sB.Type!=2)//假如不是平衡节点
, w5 `  v5 Z- Z/ ]                        {        
' p# H2 K9 Z' q2 i" u                                DP=sB.GenP-sB.LoadP;
6 |- }+ P: I2 ^                                DQ=sB.GenQ-sB.LoadQ;
                               
                                for(j=0;j<nB;j++)
) X  O7 d2 f5 K/ v                                {
                                        A=sB.Phase-sB[j].Phase;
                                        DP-=sB.Volt*(sB[j].Volt*(g[j]*cos(A)+b[j]*sin(A)));
  H" ?7 h# e. K) M5 y                                       
                                        if(sB.Type==0)//PQ节点
6 R4 c) F/ ]. K! \( L4 T+ F) W9 \                                                DQ-=sB.Volt*(sB[j].Volt*(g[j]*sin(A)-b[j]*cos(A)));
& k- ~* G8 ~0 Y' t5 i                               
" p( d3 Y- u! C                                        else if(sB.Type==1)//PV节点
% F2 v$ Z/ I8 g% ^$ y  o# a' r                                                DQ=0;
  _) [2 U( c6 N7 x9 a6 [+ c                                }
( N7 J3 E" O& l; A3 t/ p$ c* y                        }
: ?0 O, @' {1 H' x5 b& |7 b6 }& B( f                        else if(sB.Type==2)//平衡节点
. e# f$ ]8 h" m( U! F, [: p) m* Q6 N                                DP=DQ=0;
8 n5 Q7 k9 v0 S, W, I2 e# ]- O, W; ^                }
% s8 {8 z# n0 r, }" J                //for(i=0;i<nB;i++)
                //         printf("DP[%d]===%f,DQ[%d]===%f\n",i,DP,i,DQ);

  }9 \8 h+ w4 s! n* M/ [+ C            //求解修正方程
                for(i=0;i<nB-1;i++)
                    AA1=DP[i+1]/sB[i+1].Volt;
) n+ A: D8 {# ~# K# r                for(i=0;i<nB-1-count_PVnode;i++)
                        AA2=DQ[i+1+count_PVnode]/sB[i+1+count_PVnode].Volt;
                calculate_gaosi((double **)b1,BB1,AA1,NBUS-1);//AA是不平衡量，BB是解向量
: K5 v0 z; o9 i2 l                calculate_gaosi((double **)b2,BB2,AA2,NBUS-1);

! }  @, C0 B. a3 ]. G) k" R6 w, b                max1=fabs(AA1[0]);
1 n4 E- C: |2 Z! A4 S: V                for(i=1;i<nB-1;i++)
                        if(max1<fabs(AA1))       
/ w- g$ ?$ l6 ]1 _                                max1=fabs(AA1);
                        max2=fabs(AA2[0]);
                        for(i=1;i<nB-1-count_PVnode;i++)
: t7 ?; c; s6 E" B+ q% i: K                                if(max2<fabs(AA2))       
1 b5 y# {' h1 h# o6 \                                max2=fabs(AA2);
                for(i=0;i<nB-1;i++)
4 `: |  W. S* r" S+ T4 t, ?                        sB[i+1].Phase+=BB1/sB.Volt;
3 U1 q" K4 q  S1 S6 B4 T. A                for(i=0;i<nB-1-count_PVnode;i++)
                        sB[i+1+count_PVnode].Volt+=BB2;
                for(i=0;i<nB;i++)
1 {( k# O& T* H, B/ Z+ @                {
& h- Q% u! K! E0 s' A4 `+ V                        printf("sB[%d].Volt=%f,sB[%d].Phase=%f\n",i,sB.Volt,i,sB.Phase*180/PI);
9 X  q" o" \1 |8 G4 Q                       
# r7 A0 N, V2 I& D                }
! q2 k% F1 h5 s                printf("\n");
                        ci++;
        }
        while(fabs(max1)>0.00001&&fabs(max2)>0.00001&&ci<40);
7 P2 B. K" B& V. @- L8 W这是我求潮流的程序，用的PQ分解法，最后得到的结果是只能精确到小数点后第二位，第三位就不对了。

electri

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

ordinarymengya

回复 2# lllmos
2 q2 Z9 e* m0 E+ |
0 A9 ^; r! |( @) z! U
   跟你的收敛条件设置有关系啊。

ordinarymengya

回复 4# lllmos

" ~' G# d" j- B1 q' z8 H, B建议楼主换几个软件试试，SPS,PSSE，PowerSystemWorld，特别是PSSE算潮流比较专业。