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

lllmos

lllmos

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

玉门关山

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

lllmos

回复 3# 玉门关山
8 D. @, B1 u: u  I6 o
- Y  ]" F0 T+ C! s
, F! D: \# n- q0 }0 }    因为我不知道问题出在哪里了。。
do
        {
                //求解不平衡量
                for(i=0;i<nB;i++)
0 F9 y' U0 b8 ]. Z* W: R                {
                        if(sB.Type!=2)//假如不是平衡节点
                        {        
+ U" \2 C4 ~# S8 Q$ {                                DP=sB.GenP-sB.LoadP;
1 G8 ~3 D* t1 [7 g6 b                                DQ=sB.GenQ-sB.LoadQ;
                               
. M0 L8 i& l3 ~( }. Z                                for(j=0;j<nB;j++)
                                {
( R6 i0 k0 O8 x8 S9 o" L1 |% _9 n! v, K                                        A=sB.Phase-sB[j].Phase;
9 _0 I7 f  ~, _0 G! ?2 S, {1 ?                                        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)));
2 x; V) C9 p0 w. S. c; D                               
                                        else if(sB.Type==1)//PV节点
7 L+ B3 \  @; z                                                DQ=0;
                                }
- T; K  |4 G# k7 H$ D7 W                        }
0 j- d/ x7 C; R$ m# _4 S                        else if(sB.Type==2)//平衡节点
                                DP=DQ=0;
                }
$ f* r0 v! Z! t8 a% W; i) u; V                //for(i=0;i<nB;i++)
                //         printf("DP[%d]===%f,DQ[%d]===%f\n",i,DP,i,DQ);
$ E% c1 Y& @) D+ f/ }8 d& M
            //求解修正方程
. s" f' V" n% ^7 Q* X" J0 Q                for(i=0;i<nB-1;i++)
                    AA1=DP[i+1]/sB[i+1].Volt;
                for(i=0;i<nB-1-count_PVnode;i++)
7 ]0 Z* o  a5 n) Z                        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);
* Y2 T* m  e7 e8 o6 G' F, B- g, U
# b; l) D) i' e& L/ B* z( l2 Y                max1=fabs(AA1[0]);
                for(i=1;i<nB-1;i++)
                        if(max1<fabs(AA1))       
- E1 g- `  j* y* m2 u( s1 w% m7 {                                max1=fabs(AA1);
& ^' r9 M1 F: f9 q7 g5 q                        max2=fabs(AA2[0]);
2 M6 N* S4 I0 \5 J7 K* q                        for(i=1;i<nB-1-count_PVnode;i++)
                                if(max2<fabs(AA2))       
3 f" W1 k% d7 \3 h; m7 G7 q$ }; U! o% m                                max2=fabs(AA2);
                for(i=0;i<nB-1;i++)
                        sB[i+1].Phase+=BB1/sB.Volt;
                for(i=0;i<nB-1-count_PVnode;i++)
; b1 T) q" f. r8 M0 Z' E$ g" B! U                        sB[i+1+count_PVnode].Volt+=BB2;
9 T$ s3 `# A' Z+ X( b                for(i=0;i<nB;i++)
                {
                        printf("sB[%d].Volt=%f,sB[%d].Phase=%f\n",i,sB.Volt,i,sB.Phase*180/PI);
0 n6 n: e6 [1 b, F# }4 z% N3 V                       
                }
                printf("\n");
; k7 ?$ g8 i+ v) O! H9 W                        ci++;
2 y$ l+ O7 J8 [7 r3 q, k9 E* m0 }        }
        while(fabs(max1)>0.00001&&fabs(max2)>0.00001&&ci<40);
这是我求潮流的程序，用的PQ分解法，最后得到的结果是只能精确到小数点后第二位，第三位就不对了。

electri

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

ordinarymengya

回复 2# lllmos


8 _( R+ n/ ^0 n# R. p4 v: e   跟你的收敛条件设置有关系啊。

ordinarymengya

回复 4# lllmos
, {$ _1 A' {% R$ f9 E0 R
1 j3 i/ R4 q* _* q  m2 ~3 S: E0 ^2 J建议楼主换几个软件试试，SPS,PSSE，PowerSystemWorld，特别是PSSE算潮流比较专业。