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

lllmos

lllmos

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

玉门关山

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

lllmos

回复 3# 玉门关山
3 U3 y  t3 ]# W+ g8 G1 L

- U9 z' R; a* d3 P+ g8 P2 N    因为我不知道问题出在哪里了。。
6 w2 S! B$ s0 edo
) G7 J; h/ K7 X% w        {
; f9 y; O3 J7 T2 U$ d                //求解不平衡量
. w/ S; I  |# x4 h: E                for(i=0;i<nB;i++)
                {
/ S3 f. g( y7 L  F2 J. G! |. X7 ^$ ?                        if(sB.Type!=2)//假如不是平衡节点
) @; I, x7 V9 ~7 N                        {        
                                DP=sB.GenP-sB.LoadP;
3 O1 s+ m8 J' d" X                                DQ=sB.GenQ-sB.LoadQ;
                               
+ T# G+ ^8 Y6 L0 Z                                for(j=0;j<nB;j++)
                                {
: h! v) T9 g) W# a  m  O3 H$ y' I                                        A=sB.Phase-sB[j].Phase;
                                        DP-=sB.Volt*(sB[j].Volt*(g[j]*cos(A)+b[j]*sin(A)));
                                       
                                        if(sB.Type==0)//PQ节点
                                                DQ-=sB.Volt*(sB[j].Volt*(g[j]*sin(A)-b[j]*cos(A)));
% s; s) }' ~7 W7 F* Q- l- \                               
                                        else if(sB.Type==1)//PV节点
                                                DQ=0;
/ B, N6 q! c  b5 W: D1 k                                }
8 Q& c& \4 s2 |: B- z+ r0 G                        }
                        else if(sB.Type==2)//平衡节点
; R4 V! X& @$ d                                DP=DQ=0;
                }
                //for(i=0;i<nB;i++)
                //         printf("DP[%d]===%f,DQ[%d]===%f\n",i,DP,i,DQ);

            //求解修正方程
                for(i=0;i<nB-1;i++)
: B; [% M  Y% H& Z                    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;
                calculate_gaosi((double **)b1,BB1,AA1,NBUS-1);//AA是不平衡量，BB是解向量
+ R; x. N5 T. _; i$ s3 X0 M                calculate_gaosi((double **)b2,BB2,AA2,NBUS-1);
6 w1 ^2 v% [$ M) o% Z
                max1=fabs(AA1[0]);
                for(i=1;i<nB-1;i++)
) {% j: Q: x& p$ T) z: w                        if(max1<fabs(AA1))       
( \* N; K7 o0 u                                max1=fabs(AA1);
                        max2=fabs(AA2[0]);
                        for(i=1;i<nB-1-count_PVnode;i++)
1 `+ h1 Y+ c5 z8 \! c+ @4 q                                if(max2<fabs(AA2))       
                                max2=fabs(AA2);
, ?9 v7 J1 P6 t8 F                for(i=0;i<nB-1;i++)
                        sB[i+1].Phase+=BB1/sB.Volt;
# E. I$ D! ]1 x3 ~0 ~                for(i=0;i<nB-1-count_PVnode;i++)
( R: I" a7 `- V5 r1 C' e                        sB[i+1+count_PVnode].Volt+=BB2;
                for(i=0;i<nB;i++)
                {
                        printf("sB[%d].Volt=%f,sB[%d].Phase=%f\n",i,sB.Volt,i,sB.Phase*180/PI);
                       
. y5 O: |, D6 x" T( ?                }
                printf("\n");
                        ci++;
6 L) }  e8 Q  Y. d9 b( p        }
( S4 D! Q- l7 t8 _- V+ o& ^        while(fabs(max1)>0.00001&&fabs(max2)>0.00001&&ci<40);
2 F+ Y) v& h8 g! e7 i: t, q4 T这是我求潮流的程序，用的PQ分解法，最后得到的结果是只能精确到小数点后第二位，第三位就不对了。

electri

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

ordinarymengya

回复 2# lllmos
! T4 }* E, I- Y+ G: q

/ \8 H0 k' E: F1 b' R   跟你的收敛条件设置有关系啊。

ordinarymengya

回复 4# lllmos
5 |0 l/ z9 u) b
" P4 x. O7 M5 Y5 \- E4 S建议楼主换几个软件试试，SPS,PSSE，PowerSystemWorld，特别是PSSE算潮流比较专业。