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

lllmos

lllmos

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

玉门关山

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

lllmos

回复 3# 玉门关山

& ?$ O4 c/ T) a: \& E
    因为我不知道问题出在哪里了。。
do
        {
                //求解不平衡量
/ V" a5 s- H5 L' K, I0 H                for(i=0;i<nB;i++)
0 {! \2 P& Q* J5 q; K+ s7 P. r                {
                        if(sB.Type!=2)//假如不是平衡节点
: `* H4 s) N" m- Q& a- W                        {        
                                DP=sB.GenP-sB.LoadP;
! j7 g" m: b" p0 b/ ~. s3 ]                                DQ=sB.GenQ-sB.LoadQ;
8 `* D+ |6 I8 m                               
6 h. ]+ R) ]- ^6 H                                for(j=0;j<nB;j++)
                                {
                                        A=sB.Phase-sB[j].Phase;
                                        DP-=sB.Volt*(sB[j].Volt*(g[j]*cos(A)+b[j]*sin(A)));
+ j  P% K" i5 v! x9 c; o4 i9 O                                       
                                        if(sB.Type==0)//PQ节点
                                                DQ-=sB.Volt*(sB[j].Volt*(g[j]*sin(A)-b[j]*cos(A)));
, `# J. i/ a" c2 D% v/ A" \# p4 A$ Y                               
                                        else if(sB.Type==1)//PV节点
                                                DQ=0;
                                }
                        }
! ]1 P' I0 p% k                        else if(sB.Type==2)//平衡节点
5 U% w) W# j( r: w, R' Y. D                                DP=DQ=0;
                }
, D6 N7 J; {1 Y# W                //for(i=0;i<nB;i++)
                //         printf("DP[%d]===%f,DQ[%d]===%f\n",i,DP,i,DQ);

5 z% q3 w8 s/ f  M. x            //求解修正方程
3 y& K/ v7 Q0 x& I! ?5 k) g( s" J                for(i=0;i<nB-1;i++)
                    AA1=DP[i+1]/sB[i+1].Volt;
3 A: o  K* [3 Q+ B2 X                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是解向量
9 u. K' @/ U8 q( S4 C7 O                calculate_gaosi((double **)b2,BB2,AA2,NBUS-1);

                max1=fabs(AA1[0]);
8 f0 g4 {& \! ^. z" o7 ?                for(i=1;i<nB-1;i++)
                        if(max1<fabs(AA1))       
0 L4 O" i/ M7 M/ ~, T' [, `                                max1=fabs(AA1);
                        max2=fabs(AA2[0]);
) r. W( A; J- ?" K) n. c; z                        for(i=1;i<nB-1-count_PVnode;i++)
                                if(max2<fabs(AA2))       
+ O/ ]  N( O( b+ o) m  Y/ U, E/ S                                max2=fabs(AA2);
& z5 j- ~- s! ~% c                for(i=0;i<nB-1;i++)
' `4 S$ }; H1 K; B! ?* Z                        sB[i+1].Phase+=BB1/sB.Volt;
                for(i=0;i<nB-1-count_PVnode;i++)
5 d& G3 t) S( ]                        sB[i+1+count_PVnode].Volt+=BB2;
1 b% _, M) E* N$ S/ t                for(i=0;i<nB;i++)
. v8 V/ l, R& d: g: m* c                {
. @% q/ p1 j' o1 }& p( T                        printf("sB[%d].Volt=%f,sB[%d].Phase=%f\n",i,sB.Volt,i,sB.Phase*180/PI);
                       
& b3 ?, n, r( c8 [" [                }
) P9 h* l3 c! k) ^  n4 u                printf("\n");
                        ci++;
        }
        while(fabs(max1)>0.00001&&fabs(max2)>0.00001&&ci<40);
这是我求潮流的程序，用的PQ分解法，最后得到的结果是只能精确到小数点后第二位，第三位就不对了。

electri

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

ordinarymengya

回复 2# lllmos
( S4 w0 a/ ~0 A

* P" m5 f4 V3 P7 {/ T6 }# d- v0 ~   跟你的收敛条件设置有关系啊。

ordinarymengya

回复 4# lllmos
8 F/ C' {3 S, ?; O" W
建议楼主换几个软件试试，SPS,PSSE，PowerSystemWorld，特别是PSSE算潮流比较专业。