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

lllmos

lllmos

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

玉门关山

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

lllmos

回复 3# 玉门关山

' I) s, f0 ], T3 N. f
$ [. j: C7 D& F4 }* I# K    因为我不知道问题出在哪里了。。
do
        {
7 P+ \* j9 ?; J4 t: i& S1 U                //求解不平衡量
                for(i=0;i<nB;i++)
4 ^! V8 Z3 e' [5 I7 \5 [; b                {
                        if(sB.Type!=2)//假如不是平衡节点
                        {        
/ L, K" x( K( n7 T6 @" i7 s1 F2 l                                DP=sB.GenP-sB.LoadP;
                                DQ=sB.GenQ-sB.LoadQ;
7 s. D2 k2 D8 L; ~! R+ H5 p                               
                                for(j=0;j<nB;j++)
                                {
                                        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节点
$ ~- a/ A% D1 L1 K                                                DQ-=sB.Volt*(sB[j].Volt*(g[j]*sin(A)-b[j]*cos(A)));
                               
                                        else if(sB.Type==1)//PV节点
                                                DQ=0;
                                }
- z9 J; o! u) M8 M( q                        }
                        else if(sB.Type==2)//平衡节点
. s; a& w6 {$ |3 u* a" {                                DP=DQ=0;
                }
                //for(i=0;i<nB;i++)
                //         printf("DP[%d]===%f,DQ[%d]===%f\n",i,DP,i,DQ);
1 [7 ~+ i& m. P* `1 h; }* I
            //求解修正方程
                for(i=0;i<nB-1;i++)
                    AA1=DP[i+1]/sB[i+1].Volt;
/ K. _5 W: ~! G- t, X) u" U; n* J. S                for(i=0;i<nB-1-count_PVnode;i++)
/ A, g8 b2 h' n! O) f                        AA2=DQ[i+1+count_PVnode]/sB[i+1+count_PVnode].Volt;
                calculate_gaosi((double **)b1,BB1,AA1,NBUS-1);//AA是不平衡量，BB是解向量
. N1 j& V! P2 s; f5 N/ y                calculate_gaosi((double **)b2,BB2,AA2,NBUS-1);

. Y1 r2 O/ _7 O7 I% R                max1=fabs(AA1[0]);
& t2 X: n: @7 G1 }! F$ t8 T' s                for(i=1;i<nB-1;i++)
                        if(max1<fabs(AA1))       
7 E" P, c! U+ f                                max1=fabs(AA1);
" L* e/ C1 m4 S9 _8 T                        max2=fabs(AA2[0]);
% Z3 O/ N1 }8 d; V* t                        for(i=1;i<nB-1-count_PVnode;i++)
% v' `( I. B2 y5 K* X                                if(max2<fabs(AA2))       
                                max2=fabs(AA2);
: k; Q% [" U# o8 _0 A( x# P                for(i=0;i<nB-1;i++)
; N* B$ N. {# T/ i. @  T& O                        sB[i+1].Phase+=BB1/sB.Volt;
                for(i=0;i<nB-1-count_PVnode;i++)
                        sB[i+1+count_PVnode].Volt+=BB2;
  r6 B, f2 K1 O* h                for(i=0;i<nB;i++)
# _: b" T, n3 g; P8 |                {
                        printf("sB[%d].Volt=%f,sB[%d].Phase=%f\n",i,sB.Volt,i,sB.Phase*180/PI);
                       
                }
5 b/ a" n& O, R( m" E) U2 D  u# O                printf("\n");
* ?/ T2 M6 h* K* y( ?                        ci++;
        }
, I, ]$ X+ q! w        while(fabs(max1)>0.00001&&fabs(max2)>0.00001&&ci<40);
% G& ~3 U+ J5 q4 z# p4 M: Z这是我求潮流的程序，用的PQ分解法，最后得到的结果是只能精确到小数点后第二位，第三位就不对了。

electri

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

ordinarymengya

回复 2# lllmos


, b) |  ^1 C" ^* z) X2 J   跟你的收敛条件设置有关系啊。

ordinarymengya

回复 4# lllmos
' c1 ^) `& \% B( t
建议楼主换几个软件试试，SPS,PSSE，PowerSystemWorld，特别是PSSE算潮流比较专业。