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

lllmos

lllmos

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

玉门关山

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

lllmos

回复 3# 玉门关山
& n0 c& a0 g. I! i* z$ `  r! A
% o2 E; \+ b3 h# ]; M) \3 y
    因为我不知道问题出在哪里了。。
- J0 l2 C% Z: V) qdo
        {
                //求解不平衡量
- B# Q) i6 q. K! g1 ~6 r                for(i=0;i<nB;i++)
                {
& x5 L3 I; ~+ k1 X: Q! w/ ?                        if(sB.Type!=2)//假如不是平衡节点
( t9 j) G# C3 Z% n                        {        
                                DP=sB.GenP-sB.LoadP;
                                DQ=sB.GenQ-sB.LoadQ;
& C- ]3 v# e9 C                               
                                for(j=0;j<nB;j++)
                                {
                                        A=sB.Phase-sB[j].Phase;
% x0 u4 ^* G' c6 D! p9 M4 G                                        DP-=sB.Volt*(sB[j].Volt*(g[j]*cos(A)+b[j]*sin(A)));
                                       
                                        if(sB.Type==0)//PQ节点
2 H% G2 [# Q" m4 |; i                                                DQ-=sB.Volt*(sB[j].Volt*(g[j]*sin(A)-b[j]*cos(A)));
                               
4 W( ^3 n: Y* U! R                                        else if(sB.Type==1)//PV节点
3 t" W$ [; Z* Q3 ?! U1 A: d9 s                                                DQ=0;
8 `7 L& `' U4 z2 E5 Q" A3 o                                }
                        }
# g1 T" f, G0 f                        else if(sB.Type==2)//平衡节点
7 g& y3 u, C8 A                                DP=DQ=0;
                }
& J& h6 Z$ R9 e! x6 A                //for(i=0;i<nB;i++)
7 ^' t/ F; H9 j- J/ v6 \  L                //         printf("DP[%d]===%f,DQ[%d]===%f\n",i,DP,i,DQ);
3 s8 o5 V6 e1 Z
            //求解修正方程
. t7 z- }3 w; c9 v                for(i=0;i<nB-1;i++)
  {0 S0 a' A0 M                    AA1=DP[i+1]/sB[i+1].Volt;
$ Z5 A8 U, |- a+ K  m* c                for(i=0;i<nB-1-count_PVnode;i++)
+ [3 e2 x( y0 W1 W* J3 G                        AA2=DQ[i+1+count_PVnode]/sB[i+1+count_PVnode].Volt;
                calculate_gaosi((double **)b1,BB1,AA1,NBUS-1);//AA是不平衡量，BB是解向量
0 Z) m9 }* N. x4 k: m                calculate_gaosi((double **)b2,BB2,AA2,NBUS-1);
) Z3 q" c$ |, d' q5 d
' a7 y$ A6 v& s                max1=fabs(AA1[0]);
7 w5 r& c! j5 F  v& p                for(i=1;i<nB-1;i++)
                        if(max1<fabs(AA1))       
                                max1=fabs(AA1);
' u! H3 g+ u/ w8 u2 O+ B2 N                        max2=fabs(AA2[0]);
                        for(i=1;i<nB-1-count_PVnode;i++)
                                if(max2<fabs(AA2))       
) V( y  A0 g) k  F& ^                                max2=fabs(AA2);
                for(i=0;i<nB-1;i++)
                        sB[i+1].Phase+=BB1/sB.Volt;
- [) \" s& F. x7 S4 n6 s                for(i=0;i<nB-1-count_PVnode;i++)
                        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);
' l; v; O8 D) ]( Z2 z                       
                }
                printf("\n");
                        ci++;
        }
+ Q" ]( G. h. V& v* Y6 N        while(fabs(max1)>0.00001&&fabs(max2)>0.00001&&ci<40);
这是我求潮流的程序，用的PQ分解法，最后得到的结果是只能精确到小数点后第二位，第三位就不对了。

electri

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

ordinarymengya

回复 2# lllmos
/ f( l# l6 |0 y/ @" {1 n( T# D
! h, ~) d. n8 N
/ _6 U8 ?1 z% v   跟你的收敛条件设置有关系啊。

ordinarymengya

回复 4# lllmos

建议楼主换几个软件试试，SPS,PSSE，PowerSystemWorld，特别是PSSE算潮流比较专业。