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

lllmos

lllmos

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

玉门关山

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

lllmos

回复 3# 玉门关山
; `- A9 @& M( ?8 `

    因为我不知道问题出在哪里了。。
! C% G6 Y# y! ?, ~: Odo
        {
                //求解不平衡量
                for(i=0;i<nB;i++)
                {
                        if(sB.Type!=2)//假如不是平衡节点
3 A$ U& e% w) V+ M2 u                        {        
. p9 c; k+ W* ~                                DP=sB.GenP-sB.LoadP;
  g) N8 s% l1 D" U) P                                DQ=sB.GenQ-sB.LoadQ;
                               
0 c' X: i: I; N  |/ g/ k. G, I                                for(j=0;j<nB;j++)
) x! {- I  `4 T1 C3 f1 z% l                                {
/ d. y# o& ?7 H, Q1 a; U                                        A=sB.Phase-sB[j].Phase;
                                        DP-=sB.Volt*(sB[j].Volt*(g[j]*cos(A)+b[j]*sin(A)));
                                       
$ b& S! g; l% r# m: }% @8 p* [( A                                        if(sB.Type==0)//PQ节点
                                                DQ-=sB.Volt*(sB[j].Volt*(g[j]*sin(A)-b[j]*cos(A)));
5 i( G) b, d6 o# Y6 e' F1 f                               
) n# D2 q9 H. f- h                                        else if(sB.Type==1)//PV节点
4 ^6 i0 {' g3 _7 Z9 L                                                DQ=0;
( y9 C+ B  }0 h2 @' K                                }
6 Y5 S1 T' a0 m! e8 s3 m( ~                        }
                        else if(sB.Type==2)//平衡节点
% D+ n# V- {' a0 `                                DP=DQ=0;
! z- ?/ n' a$ r1 X$ s* E8 Z                }
                //for(i=0;i<nB;i++)
                //         printf("DP[%d]===%f,DQ[%d]===%f\n",i,DP,i,DQ);

            //求解修正方程
& L% l  ]- X% u                for(i=0;i<nB-1;i++)
                    AA1=DP[i+1]/sB[i+1].Volt;
                for(i=0;i<nB-1-count_PVnode;i++)
                        AA2=DQ[i+1+count_PVnode]/sB[i+1+count_PVnode].Volt;
, ^1 H! q0 J7 q9 i/ C6 E                calculate_gaosi((double **)b1,BB1,AA1,NBUS-1);//AA是不平衡量，BB是解向量
                calculate_gaosi((double **)b2,BB2,AA2,NBUS-1);
. m0 C2 |/ B* L% X+ ~/ @3 P
                max1=fabs(AA1[0]);
                for(i=1;i<nB-1;i++)
8 u8 H. Z4 B6 x) w/ ~                        if(max1<fabs(AA1))       
                                max1=fabs(AA1);
                        max2=fabs(AA2[0]);
1 x3 B" U& _" H: X                        for(i=1;i<nB-1-count_PVnode;i++)
9 \$ R- A! ^6 O/ d; a! V                                if(max2<fabs(AA2))       
  p+ `! `  M) X' ]3 f' y# L6 ^                                max2=fabs(AA2);
                for(i=0;i<nB-1;i++)
                        sB[i+1].Phase+=BB1/sB.Volt;
5 v, Z# M& @4 O& M                for(i=0;i<nB-1-count_PVnode;i++)
                        sB[i+1+count_PVnode].Volt+=BB2;
                for(i=0;i<nB;i++)
- M+ L1 Z0 P- A! J' r& {& }                {
                        printf("sB[%d].Volt=%f,sB[%d].Phase=%f\n",i,sB.Volt,i,sB.Phase*180/PI);
- e+ U( k  C0 B  ]0 \# h                       
8 N8 _& z& t7 b5 j* ]                }
                printf("\n");
                        ci++;
        }
        while(fabs(max1)>0.00001&&fabs(max2)>0.00001&&ci<40);
7 b1 r# ^: P3 s2 M; y; W1 l3 B4 g这是我求潮流的程序，用的PQ分解法，最后得到的结果是只能精确到小数点后第二位，第三位就不对了。

electri

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

ordinarymengya

回复 2# lllmos


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

ordinarymengya

回复 4# lllmos
& m5 [2 @: F# C) }8 A# x0 n& N# D$ e
5 k4 r- S1 E1 a0 \7 H建议楼主换几个软件试试，SPS,PSSE，PowerSystemWorld，特别是PSSE算潮流比较专业。