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

lllmos

lllmos

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

玉门关山

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

lllmos

回复 3# 玉门关山

# i8 T) j+ S' R% d7 w( K0 W5 I& F
    因为我不知道问题出在哪里了。。
$ w/ ], j. q, \7 y% bdo
        {
( v8 R4 \( b2 r                //求解不平衡量
                for(i=0;i<nB;i++)
                {
                        if(sB.Type!=2)//假如不是平衡节点
                        {        
                                DP=sB.GenP-sB.LoadP;
                                DQ=sB.GenQ-sB.LoadQ;
" [) j# J+ Y: g+ m! Z5 E                               
% J* u% Z1 F8 A" y                                for(j=0;j<nB;j++)
                                {
' F# m5 }: P; ?* u! R& }                                        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节点
# n3 f- K7 `0 N                                                DQ-=sB.Volt*(sB[j].Volt*(g[j]*sin(A)-b[j]*cos(A)));
                               
                                        else if(sB.Type==1)//PV节点
& L& a+ w% n5 [7 a+ N% I0 O* G) N                                                DQ=0;
) V8 A  }; O* w! C3 o: @6 I                                }
/ ]1 D8 e. K# x1 H) r                        }
                        else if(sB.Type==2)//平衡节点
                                DP=DQ=0;
                }
                //for(i=0;i<nB;i++)
- z. h* z. {- d! [; q                //         printf("DP[%d]===%f,DQ[%d]===%f\n",i,DP,i,DQ);
+ X3 M& _( ~7 R* j
            //求解修正方程
                for(i=0;i<nB-1;i++)
                    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);
; [  a/ `3 g' d; q. b
9 r; a* v  b$ o; _, \% [, k                max1=fabs(AA1[0]);
! r0 G( P: a" g% {& w5 T                for(i=1;i<nB-1;i++)
                        if(max1<fabs(AA1))       
                                max1=fabs(AA1);
* H1 H1 Z8 e4 u9 x* X7 C. `4 v  ]                        max2=fabs(AA2[0]);
" M6 H( Y# L  y7 u: H" I. q                        for(i=1;i<nB-1-count_PVnode;i++)
                                if(max2<fabs(AA2))       
                                max2=fabs(AA2);
                for(i=0;i<nB-1;i++)
                        sB[i+1].Phase+=BB1/sB.Volt;
: o1 p3 m" C' f' k                for(i=0;i<nB-1-count_PVnode;i++)
  h' k5 E+ ?" U                        sB[i+1+count_PVnode].Volt+=BB2;
- U5 ]/ t& Y6 D. U2 W; G4 ~: R                for(i=0;i<nB;i++)
                {
9 E) h8 H, y9 P6 V                        printf("sB[%d].Volt=%f,sB[%d].Phase=%f\n",i,sB.Volt,i,sB.Phase*180/PI);
                       
/ f- {- n3 r5 E; q6 d  D/ ~2 r7 D                }
+ O1 T/ c8 H" U# T2 ]9 Q  Y                printf("\n");
                        ci++;
        }
        while(fabs(max1)>0.00001&&fabs(max2)>0.00001&&ci<40);
这是我求潮流的程序，用的PQ分解法，最后得到的结果是只能精确到小数点后第二位，第三位就不对了。

electri

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

ordinarymengya

回复 2# lllmos


   跟你的收敛条件设置有关系啊。

ordinarymengya

回复 4# lllmos
' H( ], s2 y% u; S( l- B2 B" O
建议楼主换几个软件试试，SPS,PSSE，PowerSystemWorld，特别是PSSE算潮流比较专业。