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

lllmos

lllmos

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

玉门关山

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

lllmos

回复 3# 玉门关山
: k& f  u6 T6 l. r0 B; w

+ ?; X. i% z& S2 F    因为我不知道问题出在哪里了。。
2 Y% q% }7 {7 q% d3 u8 P. E9 {% p, tdo
        {
) ~8 N2 I# t, H                //求解不平衡量
                for(i=0;i<nB;i++)
; X& H$ @: O: z$ X" g: G! i7 o- D                {
! |# ~& l/ t9 a: ?4 W                        if(sB.Type!=2)//假如不是平衡节点
3 C2 F, d' x2 i2 \                        {        
                                DP=sB.GenP-sB.LoadP;
4 N. K9 A5 Z% c" s3 d# W6 n                                DQ=sB.GenQ-sB.LoadQ;
                               
4 z( h; v. {- e6 h  M' S5 E' T                                for(j=0;j<nB;j++)
5 g* u7 j8 k8 U" O- F                                {
                                        A=sB.Phase-sB[j].Phase;
                                        DP-=sB.Volt*(sB[j].Volt*(g[j]*cos(A)+b[j]*sin(A)));
                                       
4 T1 F: p# I# z: V" O: C                                        if(sB.Type==0)//PQ节点
/ F& s. g# N+ D                                                DQ-=sB.Volt*(sB[j].Volt*(g[j]*sin(A)-b[j]*cos(A)));
* ~8 q4 C. r. ?2 `5 x! A0 J6 [1 `                               
                                        else if(sB.Type==1)//PV节点
- T% [& J: l) R; D! W                                                DQ=0;
8 q9 e" ~! ?9 d  @                                }
; O# f+ N# L3 D0 f. `                        }
                        else if(sB.Type==2)//平衡节点
                                DP=DQ=0;
                }
                //for(i=0;i<nB;i++)
                //         printf("DP[%d]===%f,DQ[%d]===%f\n",i,DP,i,DQ);

* t) x8 M- S8 t& }% E% q            //求解修正方程
                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;
6 ~0 ]; d3 E2 `* ^3 t                calculate_gaosi((double **)b1,BB1,AA1,NBUS-1);//AA是不平衡量，BB是解向量
) J6 w0 z& x) {0 ~5 ]/ L                calculate_gaosi((double **)b2,BB2,AA2,NBUS-1);
  _% `, w8 Q  A0 ^; W0 P# ?
                max1=fabs(AA1[0]);
. i' K" s" Q5 S                for(i=1;i<nB-1;i++)
5 l1 s$ J0 d' e2 r                        if(max1<fabs(AA1))       
# v# j$ o( |$ l+ N( T' W4 @6 G                                max1=fabs(AA1);
3 A8 D% R$ p4 f$ t2 N$ h  B                        max2=fabs(AA2[0]);
                        for(i=1;i<nB-1-count_PVnode;i++)
8 D* e7 O0 ~; n2 s# ^/ q                                if(max2<fabs(AA2))       
                                max2=fabs(AA2);
/ e, c  R5 w0 M$ H" S                for(i=0;i<nB-1;i++)
                        sB[i+1].Phase+=BB1/sB.Volt;
                for(i=0;i<nB-1-count_PVnode;i++)
; W" a- i& ~* j3 L! l1 T                        sB[i+1+count_PVnode].Volt+=BB2;
! d+ i' I- W$ L                for(i=0;i<nB;i++)
8 Q8 S9 k, F& x- O, Q" C                {
0 C7 A4 C. N3 l                        printf("sB[%d].Volt=%f,sB[%d].Phase=%f\n",i,sB.Volt,i,sB.Phase*180/PI);
! _- s( C# S: `/ G  Z3 w                       
1 b4 |/ m* h: ]$ s+ x) G8 ]8 Z7 s* D                }
                printf("\n");
                        ci++;
% C) Q" K4 A' O        }
/ Y* a1 h% L3 ~- L: g        while(fabs(max1)>0.00001&&fabs(max2)>0.00001&&ci<40);
这是我求潮流的程序，用的PQ分解法，最后得到的结果是只能精确到小数点后第二位，第三位就不对了。

electri

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

ordinarymengya

回复 2# lllmos
1 C0 M- q# ], X. N5 f# j, L9 ^3 v
1 ~, v4 ~" \6 z0 A
   跟你的收敛条件设置有关系啊。

ordinarymengya

回复 4# lllmos
3 p1 ^# w7 p& Q4 f7 p8 @; E* W4 s
建议楼主换几个软件试试，SPS,PSSE，PowerSystemWorld，特别是PSSE算潮流比较专业。