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

lllmos

lllmos

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

玉门关山

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

lllmos

回复 3# 玉门关山


  ^& b6 K# y5 O3 g: A7 H9 |8 d    因为我不知道问题出在哪里了。。
( X, I3 M% l6 V+ l1 }do
        {
$ L3 m! `8 n7 u                //求解不平衡量
% l0 C+ q, A' `- \2 c                for(i=0;i<nB;i++)
! J$ R3 y9 [5 E! h' m! [" p7 l                {
                        if(sB.Type!=2)//假如不是平衡节点
                        {        
3 L* C3 L' ?3 P1 R' c                                DP=sB.GenP-sB.LoadP;
                                DQ=sB.GenQ-sB.LoadQ;
                               
                                for(j=0;j<nB;j++)
; o2 T+ p. s* i& G                                {
# u  J% B2 Z/ s" {# i9 [                                        A=sB.Phase-sB[j].Phase;
                                        DP-=sB.Volt*(sB[j].Volt*(g[j]*cos(A)+b[j]*sin(A)));
  a& r/ ~7 I; Y( F% W                                       
" o& h+ X4 c* @; Y3 h                                        if(sB.Type==0)//PQ节点
  g+ g* v* S9 u  P0 n- p4 N: K" y4 S                                                DQ-=sB.Volt*(sB[j].Volt*(g[j]*sin(A)-b[j]*cos(A)));
7 t/ L* Y1 B( N, c                               
                                        else if(sB.Type==1)//PV节点
6 ]+ W7 ]: A! ~: a! L; o! @                                                DQ=0;
                                }
                        }
                        else if(sB.Type==2)//平衡节点
                                DP=DQ=0;
                }
1 i. c: G, r% Z4 Y                //for(i=0;i<nB;i++)
) b, Z7 W/ k. ~# z( j                //         printf("DP[%d]===%f,DQ[%d]===%f\n",i,DP,i,DQ);

            //求解修正方程
3 I* B& S) g& p  F                for(i=0;i<nB-1;i++)
, T# o- C2 F' X! e8 O3 a9 i, M                    AA1=DP[i+1]/sB[i+1].Volt;
1 `& R; a. }5 U* [( g  t                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是解向量
% X' H& g" `4 A: B& P# f  @                calculate_gaosi((double **)b2,BB2,AA2,NBUS-1);

8 a! C% }, n' A& f, C, t: X                max1=fabs(AA1[0]);
: A- Q+ O5 c. a0 |' n, F. d/ v! G& i                for(i=1;i<nB-1;i++)
                        if(max1<fabs(AA1))       
                                max1=fabs(AA1);
                        max2=fabs(AA2[0]);
& @/ V( F1 W8 D$ u; X# [8 B                        for(i=1;i<nB-1-count_PVnode;i++)
3 ^+ z) R( D" p) n0 x& }                                if(max2<fabs(AA2))       
                                max2=fabs(AA2);
                for(i=0;i<nB-1;i++)
4 c% g# _- N; o; |9 }                        sB[i+1].Phase+=BB1/sB.Volt;
! n& @7 o+ @8 m3 k1 j' ?# y" {" l                for(i=0;i<nB-1-count_PVnode;i++)
0 E+ G8 J: k- G& ^+ S) @                        sB[i+1+count_PVnode].Volt+=BB2;
                for(i=0;i<nB;i++)
                {
                        printf("sB[%d].Volt=%f,sB[%d].Phase=%f\n",i,sB.Volt,i,sB.Phase*180/PI);
1 `0 K8 Q7 x$ A2 H: E6 {2 Q                       
1 ?9 @+ G9 j1 b% V9 L+ ?                }
" e7 }, D! w7 n, x+ _                printf("\n");
                        ci++;
        }
        while(fabs(max1)>0.00001&&fabs(max2)>0.00001&&ci<40);
: X: q% N( |1 f; I" @2 j$ L这是我求潮流的程序，用的PQ分解法，最后得到的结果是只能精确到小数点后第二位，第三位就不对了。

electri

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

ordinarymengya

回复 2# lllmos


. ?$ Z5 p# E9 e, T' \" x3 P+ a, y   跟你的收敛条件设置有关系啊。

ordinarymengya

回复 4# lllmos
, N/ k. J; D1 Q# W" M! d
4 E/ S9 K! @/ |: i* D& E! @6 b3 U建议楼主换几个软件试试，SPS,PSSE，PowerSystemWorld，特别是PSSE算潮流比较专业。