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

lllmos

lllmos

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

玉门关山

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

lllmos

回复 3# 玉门关山


$ H6 z" Q. a! f! a: r3 d3 K    因为我不知道问题出在哪里了。。
4 A( l& a/ g$ a+ y7 |do
        {
                //求解不平衡量
6 H, j6 j1 w% ]* W                for(i=0;i<nB;i++)
                {
# T! l7 x/ N# D7 {1 [  ~( b, [                        if(sB.Type!=2)//假如不是平衡节点
                        {        
                                DP=sB.GenP-sB.LoadP;
                                DQ=sB.GenQ-sB.LoadQ;
- v" {! O0 S; V% ]& U7 i% I- T                               
/ f4 x8 S4 g9 i) s% o- l                                for(j=0;j<nB;j++)
" G: T0 u- H. ^                                {
/ V# T- }8 ?* O, K- k                                        A=sB.Phase-sB[j].Phase;
                                        DP-=sB.Volt*(sB[j].Volt*(g[j]*cos(A)+b[j]*sin(A)));
                                       
" C/ n( [) `; D( k: I$ g, ^' Z, @& ]                                        if(sB.Type==0)//PQ节点
                                                DQ-=sB.Volt*(sB[j].Volt*(g[j]*sin(A)-b[j]*cos(A)));
) Y; O, X  [- _! j* [' V, X8 z: n                               
8 Z, l* U1 O) I/ V                                        else if(sB.Type==1)//PV节点
                                                DQ=0;
                                }
+ `" B" X: F( g& V! i                        }
/ s5 _, U. I$ K( F6 i  Z0 b; V                        else if(sB.Type==2)//平衡节点
                                DP=DQ=0;
8 Q6 U9 q" m, c  w0 c  [                }
                //for(i=0;i<nB;i++)
                //         printf("DP[%d]===%f,DQ[%d]===%f\n",i,DP,i,DQ);

/ |& h* [# B6 u+ s$ s            //求解修正方程
                for(i=0;i<nB-1;i++)
4 [  m2 p7 p' p; d% @( z! R  E# }                    AA1=DP[i+1]/sB[i+1].Volt;
/ ?2 B( j) @& O' A7 |% q/ F                for(i=0;i<nB-1-count_PVnode;i++)
                        AA2=DQ[i+1+count_PVnode]/sB[i+1+count_PVnode].Volt;
& M# K: j+ q) }% E                calculate_gaosi((double **)b1,BB1,AA1,NBUS-1);//AA是不平衡量，BB是解向量
                calculate_gaosi((double **)b2,BB2,AA2,NBUS-1);

, ~: k; `( L* @8 V                max1=fabs(AA1[0]);
2 ?  t; O* L4 @4 U3 |                for(i=1;i<nB-1;i++)
! y" g4 {7 Q* p                        if(max1<fabs(AA1))       
& t) U* C$ @6 A6 k( @& M6 N( K& i                                max1=fabs(AA1);
2 V9 `* P: _9 U9 D* ?* @3 l/ N                        max2=fabs(AA2[0]);
                        for(i=1;i<nB-1-count_PVnode;i++)
                                if(max2<fabs(AA2))       
+ H$ Q5 P; }, z                                max2=fabs(AA2);
                for(i=0;i<nB-1;i++)
9 U8 R) Z' R" b$ S                        sB[i+1].Phase+=BB1/sB.Volt;
                for(i=0;i<nB-1-count_PVnode;i++)
                        sB[i+1+count_PVnode].Volt+=BB2;
$ V8 y- s' a  z+ [- \! l                for(i=0;i<nB;i++)
- L- E2 a2 K, Y8 i7 U. D                {
                        printf("sB[%d].Volt=%f,sB[%d].Phase=%f\n",i,sB.Volt,i,sB.Phase*180/PI);
                       
6 ]6 r5 A+ V' \# s  q                }
  Y# m, |- Q/ \. }                printf("\n");
  i6 w7 q  x2 x. e3 P/ H7 U                        ci++;
) ^" y# x/ l- `        }
! H. N. A4 u  E  @1 A$ H' n        while(fabs(max1)>0.00001&&fabs(max2)>0.00001&&ci<40);
这是我求潮流的程序，用的PQ分解法，最后得到的结果是只能精确到小数点后第二位，第三位就不对了。

electri

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

ordinarymengya

回复 2# lllmos


8 L1 F$ C: d" q+ F4 X' |   跟你的收敛条件设置有关系啊。

ordinarymengya

回复 4# lllmos
2 D8 q, m0 H4 B1 k$ y0 h! Q; H
建议楼主换几个软件试试，SPS,PSSE，PowerSystemWorld，特别是PSSE算潮流比较专业。