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

lllmos

lllmos

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

玉门关山

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

lllmos

回复 3# 玉门关山

& @% P) q0 s* m/ `9 b
  }! M$ j& ^, t: w7 k    因为我不知道问题出在哪里了。。
9 b' L% E7 u, o0 W9 A: J0 ~! sdo
        {
                //求解不平衡量
                for(i=0;i<nB;i++)
                {
) l) v4 p/ {1 Q# B                        if(sB.Type!=2)//假如不是平衡节点
                        {        
                                DP=sB.GenP-sB.LoadP;
/ h5 x3 s' [4 {" I% M                                DQ=sB.GenQ-sB.LoadQ;
0 l$ y: f+ ]' _" R* U: N  g$ l                               
8 R) i: n. ?1 f! q0 H                                for(j=0;j<nB;j++)
                                {
) [/ t( j: W7 O6 C2 T3 o                                        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)));
                               
                                        else if(sB.Type==1)//PV节点
3 J" D/ q% Q9 Z) w- g' W, E                                                DQ=0;
                                }
7 ~2 {+ a$ N( N9 [, H                        }
                        else if(sB.Type==2)//平衡节点
' l5 X7 f2 Y  u' U) i! l. G$ o                                DP=DQ=0;
                }
' v( i  h( ]! F0 g                //for(i=0;i<nB;i++)
                //         printf("DP[%d]===%f,DQ[%d]===%f\n",i,DP,i,DQ);

            //求解修正方程
- \0 k' D' h# [' P& j7 k! T& E                for(i=0;i<nB-1;i++)
% W9 K  X' X% C7 Y; B& k- E                    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是解向量
                calculate_gaosi((double **)b2,BB2,AA2,NBUS-1);
/ v5 ]' ?7 Y1 b
                max1=fabs(AA1[0]);
                for(i=1;i<nB-1;i++)
1 \& ]: u  W. ^$ t- _                        if(max1<fabs(AA1))       
                                max1=fabs(AA1);
3 v; \4 h' I9 {                        max2=fabs(AA2[0]);
4 i4 S$ V0 i" K/ Z" z                        for(i=1;i<nB-1-count_PVnode;i++)
                                if(max2<fabs(AA2))       
                                max2=fabs(AA2);
                for(i=0;i<nB-1;i++)
9 o, b' h  \! D                        sB[i+1].Phase+=BB1/sB.Volt;
" f$ U7 S# {' P/ U% b! T                for(i=0;i<nB-1-count_PVnode;i++)
! F( r, u& Q& [. x$ s! s9 H! W                        sB[i+1+count_PVnode].Volt+=BB2;
9 K% b7 Y' D3 `0 U' s1 j/ {                for(i=0;i<nB;i++)
                {
% }% }% q! @1 r# h2 K                        printf("sB[%d].Volt=%f,sB[%d].Phase=%f\n",i,sB.Volt,i,sB.Phase*180/PI);
                       
9 L) l$ y0 Y, U                }
                printf("\n");
7 y. L! b8 v3 v. c8 D                        ci++;
. @0 g! `# J! `, T        }
        while(fabs(max1)>0.00001&&fabs(max2)>0.00001&&ci<40);
这是我求潮流的程序，用的PQ分解法，最后得到的结果是只能精确到小数点后第二位，第三位就不对了。

electri

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

ordinarymengya

回复 2# lllmos


+ K$ Z+ M; C1 z1 e   跟你的收敛条件设置有关系啊。

ordinarymengya

回复 4# lllmos
% {5 Q$ H' `8 L3 @
/ e4 g7 R$ o  D  L建议楼主换几个软件试试，SPS,PSSE，PowerSystemWorld，特别是PSSE算潮流比较专业。