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

lllmos

lllmos

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

玉门关山

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

lllmos

回复 3# 玉门关山
8 @. x+ s3 N/ C3 |% \( A- _' @
- a; O. h) R! i) P+ q, l( U
& g2 }# x9 t" U: s/ F5 T& s    因为我不知道问题出在哪里了。。
do
* Q" x1 Q  X, b$ c& w        {
6 j  @; w' G; g$ r, M6 g: F1 J1 k                //求解不平衡量
, Z7 ^9 y- H. Z                for(i=0;i<nB;i++)
                {
4 X1 T; ~3 Q- ^5 S/ |9 |+ @; X3 |                        if(sB.Type!=2)//假如不是平衡节点
                        {        
                                DP=sB.GenP-sB.LoadP;
0 b  x" g7 F0 q7 V) [) G) C                                DQ=sB.GenQ-sB.LoadQ;
                               
                                for(j=0;j<nB;j++)
                                {
8 G8 z/ s% \: e4 r7 y6 w7 Y                                        A=sB.Phase-sB[j].Phase;
                                        DP-=sB.Volt*(sB[j].Volt*(g[j]*cos(A)+b[j]*sin(A)));
& ?8 }7 ^9 c; R( j; {2 b# P0 I+ a                                       
6 K' T7 u5 T2 \( A& F6 ~8 K9 o                                        if(sB.Type==0)//PQ节点
                                                DQ-=sB.Volt*(sB[j].Volt*(g[j]*sin(A)-b[j]*cos(A)));
0 y4 x: C* v% q) L! R. j0 b                               
                                        else if(sB.Type==1)//PV节点
                                                DQ=0;
2 J% f4 v/ x4 W; `$ y                                }
                        }
                        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);
! s4 C& t: U& z
            //求解修正方程
                for(i=0;i<nB-1;i++)
& I6 C9 v9 r+ }; k- t$ G                    AA1=DP[i+1]/sB[i+1].Volt;
$ b: r0 e; w( h6 U* T                for(i=0;i<nB-1-count_PVnode;i++)
                        AA2=DQ[i+1+count_PVnode]/sB[i+1+count_PVnode].Volt;
0 @8 T% s6 p2 t7 E$ s1 h9 r                calculate_gaosi((double **)b1,BB1,AA1,NBUS-1);//AA是不平衡量，BB是解向量
2 b1 u: _, X+ e, w" c  _" U2 V                calculate_gaosi((double **)b2,BB2,AA2,NBUS-1);
! h  ~4 d" @6 k
, C; \0 p: L4 F3 c+ x                max1=fabs(AA1[0]);
                for(i=1;i<nB-1;i++)
( |  q) Y( k; Z; g4 ?" D, P                        if(max1<fabs(AA1))       
; O: v) C) m# ?) f! y+ ?% z% O" {                                max1=fabs(AA1);
                        max2=fabs(AA2[0]);
$ s7 ^% _1 ]: |3 N                        for(i=1;i<nB-1-count_PVnode;i++)
                                if(max2<fabs(AA2))       
                                max2=fabs(AA2);
& K* p# b& x+ G& O( ^& Y# }- ~  N                for(i=0;i<nB-1;i++)
1 j2 d/ @% @, x2 r. l6 F" k                        sB[i+1].Phase+=BB1/sB.Volt;
                for(i=0;i<nB-1-count_PVnode;i++)
                        sB[i+1+count_PVnode].Volt+=BB2;
, z! i4 d/ l. \5 v1 q                for(i=0;i<nB;i++)
                {
                        printf("sB[%d].Volt=%f,sB[%d].Phase=%f\n",i,sB.Volt,i,sB.Phase*180/PI);
. C7 y# _% Y) @* A2 {                       
                }
                printf("\n");
                        ci++;
        }
        while(fabs(max1)>0.00001&&fabs(max2)>0.00001&&ci<40);
这是我求潮流的程序，用的PQ分解法，最后得到的结果是只能精确到小数点后第二位，第三位就不对了。

electri

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

ordinarymengya

回复 2# lllmos


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

ordinarymengya

回复 4# lllmos
; a, f* e' S6 K9 Q3 \
建议楼主换几个软件试试，SPS,PSSE，PowerSystemWorld，特别是PSSE算潮流比较专业。