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

lllmos

lllmos

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

玉门关山

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

lllmos

回复 3# 玉门关山

7 ^/ t* Q5 B8 S# _0 Z: S" K( c
% j# |* J0 W( y: w/ h$ R2 u    因为我不知道问题出在哪里了。。
$ }6 L$ ^6 G% x  P% v- K/ \7 m0 Mdo
        {
                //求解不平衡量
$ v& W1 G( e' @5 D; l                for(i=0;i<nB;i++)
6 v  F  h# m5 w6 q0 G, Y                {
                        if(sB.Type!=2)//假如不是平衡节点
                        {        
, p2 b- O% c/ y+ ^. n+ w& G8 m                                DP=sB.GenP-sB.LoadP;
, [$ P9 U2 A$ t. I1 Q) K. z                                DQ=sB.GenQ-sB.LoadQ;
. n( L$ _- i( V                               
) _9 @- E3 h# f! R! \                                for(j=0;j<nB;j++)
; q8 k4 d# ~5 N* q) g' K# E5 r& h                                {
$ _3 |! x5 z. y' y2 @. n                                        A=sB.Phase-sB[j].Phase;
                                        DP-=sB.Volt*(sB[j].Volt*(g[j]*cos(A)+b[j]*sin(A)));
' E# V8 x* j$ n. ], Q$ S: j* q                                       
" i0 y9 p. J, H& @. o3 w                                        if(sB.Type==0)//PQ节点
                                                DQ-=sB.Volt*(sB[j].Volt*(g[j]*sin(A)-b[j]*cos(A)));
                               
                                        else if(sB.Type==1)//PV节点
                                                DQ=0;
( Y2 J1 m6 X- C5 T+ ?                                }
1 j& K8 g' l2 C- x, y                        }
% F) Z7 N6 O) |# H% f                        else if(sB.Type==2)//平衡节点
# E" `. I1 X7 v2 {  e! B                                DP=DQ=0;
                }
3 U! }6 t. K: l: W                //for(i=0;i<nB;i++)
) }8 r8 b+ O1 k! G  ^( @  _. O                //         printf("DP[%d]===%f,DQ[%d]===%f\n",i,DP,i,DQ);

            //求解修正方程
1 }! @* h6 a2 c) X- Q. I6 o. B                for(i=0;i<nB-1;i++)
' {  _& ]2 {9 @. H2 t                    AA1=DP[i+1]/sB[i+1].Volt;
  V$ f8 G# h" j+ ^; B1 f# A                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是解向量
                calculate_gaosi((double **)b2,BB2,AA2,NBUS-1);
. ~! E( o- p# A3 e1 O
                max1=fabs(AA1[0]);
                for(i=1;i<nB-1;i++)
- {, L: W7 `: z; W9 _! c                        if(max1<fabs(AA1))       
1 x3 N+ V8 \: V) E- X                                max1=fabs(AA1);
                        max2=fabs(AA2[0]);
                        for(i=1;i<nB-1-count_PVnode;i++)
. @. Q; r5 C2 N2 f                                if(max2<fabs(AA2))       
* D8 h/ g6 Q( j; N                                max2=fabs(AA2);
# F, b' P( u% }; i% \; l                for(i=0;i<nB-1;i++)
                        sB[i+1].Phase+=BB1/sB.Volt;
                for(i=0;i<nB-1-count_PVnode;i++)
- F0 b' ^4 w) C8 e' q, @& ^                        sB[i+1+count_PVnode].Volt+=BB2;
                for(i=0;i<nB;i++)
- L# O; Q7 R/ H- o) r" l0 j+ _" L                {
! J# l: i  Q% l6 a) n                        printf("sB[%d].Volt=%f,sB[%d].Phase=%f\n",i,sB.Volt,i,sB.Phase*180/PI);
6 }% H- F3 B5 F; q                       
                }
                printf("\n");
                        ci++;
, v/ _" y3 t# K  t* g3 ^        }
$ e: D( @9 N; ]2 p" Q        while(fabs(max1)>0.00001&&fabs(max2)>0.00001&&ci<40);
这是我求潮流的程序，用的PQ分解法，最后得到的结果是只能精确到小数点后第二位，第三位就不对了。

electri

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

ordinarymengya

回复 2# lllmos
6 Z" Y' z0 N, Q0 o; C, B3 G

; u1 ]6 P6 q. R! v' L2 G5 n1 b9 X   跟你的收敛条件设置有关系啊。

ordinarymengya

回复 4# lllmos

建议楼主换几个软件试试，SPS,PSSE，PowerSystemWorld，特别是PSSE算潮流比较专业。