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

lllmos

lllmos

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

玉门关山

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

lllmos

回复 3# 玉门关山

& g  z7 Q. p* A! G; a1 v( A
    因为我不知道问题出在哪里了。。
- \* p4 }3 M; ?1 M: Ydo
. J+ ?" G/ H; m6 O, V% d1 Q% a' w        {
. a7 ?3 U$ }0 G0 j; w                //求解不平衡量
; ]; w8 ]/ l# P  N; Z                for(i=0;i<nB;i++)
                {
                        if(sB.Type!=2)//假如不是平衡节点
                        {        
/ x! y) f( x+ A: @# z                                DP=sB.GenP-sB.LoadP;
                                DQ=sB.GenQ-sB.LoadQ;
9 S& d) u, |; D, n7 x, j5 D                               
" U+ U6 b; Q8 @. l5 q. V" {                                for(j=0;j<nB;j++)
                                {
                                        A=sB.Phase-sB[j].Phase;
9 g1 I6 ^, L& A. v5 n                                        DP-=sB.Volt*(sB[j].Volt*(g[j]*cos(A)+b[j]*sin(A)));
                                       
+ K& p# B; w+ c                                        if(sB.Type==0)//PQ节点
3 w$ I& ?6 V1 O+ u                                                DQ-=sB.Volt*(sB[j].Volt*(g[j]*sin(A)-b[j]*cos(A)));
                               
                                        else if(sB.Type==1)//PV节点
/ c0 X  h8 ~1 j# r                                                DQ=0;
0 R& t  R9 n" ]. Q                                }
                        }
                        else if(sB.Type==2)//平衡节点
1 `7 W3 F9 k  c8 b: E* z                                DP=DQ=0;
                }
5 Q7 J' b! n* Z5 I2 a                //for(i=0;i<nB;i++)
) I+ Y+ Y' p* _" f  ^! {                //         printf("DP[%d]===%f,DQ[%d]===%f\n",i,DP,i,DQ);

9 H! x) a) [; _" H            //求解修正方程
' e  ~8 o* ]/ B+ K. i" e2 x! T+ H                for(i=0;i<nB-1;i++)
3 I9 R: F- L3 `8 s  A2 Z. @/ a                    AA1=DP[i+1]/sB[i+1].Volt;
, X) I; _! j/ _+ ?: n                for(i=0;i<nB-1-count_PVnode;i++)
                        AA2=DQ[i+1+count_PVnode]/sB[i+1+count_PVnode].Volt;
( f  o  e4 M" @3 B" ^7 K3 w+ ]                calculate_gaosi((double **)b1,BB1,AA1,NBUS-1);//AA是不平衡量，BB是解向量
  [( n% y1 d6 O2 _6 N                calculate_gaosi((double **)b2,BB2,AA2,NBUS-1);

2 l& |: p' u* p  o                max1=fabs(AA1[0]);
                for(i=1;i<nB-1;i++)
                        if(max1<fabs(AA1))       
4 x) r* N6 @& Y) u                                max1=fabs(AA1);
9 V, F) z" G( x4 l$ P. k                        max2=fabs(AA2[0]);
                        for(i=1;i<nB-1-count_PVnode;i++)
0 E; q2 m. [1 `* X# \5 o) |3 Z                                if(max2<fabs(AA2))       
) _; t: ]" a5 ~! `4 u7 Y4 S0 `                                max2=fabs(AA2);
                for(i=0;i<nB-1;i++)
                        sB[i+1].Phase+=BB1/sB.Volt;
! s, K+ h5 s) F3 ?                for(i=0;i<nB-1-count_PVnode;i++)
3 T+ j* B' [% ^4 g                        sB[i+1+count_PVnode].Volt+=BB2;
$ d4 J/ G  v2 c4 N, s, u/ s# V                for(i=0;i<nB;i++)
                {
                        printf("sB[%d].Volt=%f,sB[%d].Phase=%f\n",i,sB.Volt,i,sB.Phase*180/PI);
                       
                }
                printf("\n");
- Y; u( I4 C# b" i$ ^                        ci++;
        }
        while(fabs(max1)>0.00001&&fabs(max2)>0.00001&&ci<40);
这是我求潮流的程序，用的PQ分解法，最后得到的结果是只能精确到小数点后第二位，第三位就不对了。

electri

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

ordinarymengya

回复 2# lllmos

, D& [! G8 i$ F2 A1 z
5 U$ ~  @( o' q5 g1 l& ]. ?   跟你的收敛条件设置有关系啊。

ordinarymengya

回复 4# lllmos

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