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

lllmos

lllmos

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

玉门关山

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

lllmos

回复 3# 玉门关山
+ `" R* L0 Z" G6 R2 k  q

0 o4 D  Q- E/ I, L0 ~2 s1 q    因为我不知道问题出在哪里了。。
( g# v5 ?& z6 j/ N& Ddo
$ B( a6 f% c2 D% f1 w; O        {
                //求解不平衡量
                for(i=0;i<nB;i++)
8 ]. b0 j& q9 i" b( E5 Z                {
                        if(sB.Type!=2)//假如不是平衡节点
                        {        
                                DP=sB.GenP-sB.LoadP;
) A" @2 F, `" O4 O5 O1 e8 n" T                                DQ=sB.GenQ-sB.LoadQ;
                               
; |6 I3 t* \$ {, z' v8 p                                for(j=0;j<nB;j++)
7 X) I& ]9 Q& O6 S! k4 d! l! P                                {
                                        A=sB.Phase-sB[j].Phase;
                                        DP-=sB.Volt*(sB[j].Volt*(g[j]*cos(A)+b[j]*sin(A)));
7 J- S$ e% x( C) ^, j                                       
                                        if(sB.Type==0)//PQ节点
                                                DQ-=sB.Volt*(sB[j].Volt*(g[j]*sin(A)-b[j]*cos(A)));
# I& C  @; Y, t! S                               
                                        else if(sB.Type==1)//PV节点
1 x8 N3 i3 }( d5 l                                                DQ=0;
0 [! m0 @9 {8 d/ f; e5 [8 a! t                                }
$ l; g. Z, U$ n+ }( B3 ^                        }
" a9 R7 p% w9 r" ^4 k" b                        else if(sB.Type==2)//平衡节点
                                DP=DQ=0;
                }
! e, d; O. A# }- e! T                //for(i=0;i<nB;i++)
- I# q; B( J3 R  L: @% L                //         printf("DP[%d]===%f,DQ[%d]===%f\n",i,DP,i,DQ);

            //求解修正方程
1 i8 g4 ~) r$ A( _, r. Q& P' Z                for(i=0;i<nB-1;i++)
# K+ Q2 M; l+ _6 R                    AA1=DP[i+1]/sB[i+1].Volt;
                for(i=0;i<nB-1-count_PVnode;i++)
0 U2 e9 s! w, {0 n1 a2 i+ h                        AA2=DQ[i+1+count_PVnode]/sB[i+1+count_PVnode].Volt;
$ ~7 e+ Z, J2 C! K& z7 w, l                calculate_gaosi((double **)b1,BB1,AA1,NBUS-1);//AA是不平衡量，BB是解向量
                calculate_gaosi((double **)b2,BB2,AA2,NBUS-1);
" k0 H/ ?% @( i5 _3 k' O& \. w
                max1=fabs(AA1[0]);
                for(i=1;i<nB-1;i++)
                        if(max1<fabs(AA1))       
* U( ]: H, [" \* q                                max1=fabs(AA1);
                        max2=fabs(AA2[0]);
8 t& ?% v6 D3 W6 H3 p                        for(i=1;i<nB-1-count_PVnode;i++)
                                if(max2<fabs(AA2))       
5 M3 c- r: c& W4 y+ \                                max2=fabs(AA2);
                for(i=0;i<nB-1;i++)
                        sB[i+1].Phase+=BB1/sB.Volt;
; _3 y2 U6 v# @                for(i=0;i<nB-1-count_PVnode;i++)
                        sB[i+1+count_PVnode].Volt+=BB2;
8 L/ F& O% |+ p, {2 ]                for(i=0;i<nB;i++)
) }# ?5 ?' U; U% N                {
  m8 I( b. w9 |) L+ n: x: g, V                        printf("sB[%d].Volt=%f,sB[%d].Phase=%f\n",i,sB.Volt,i,sB.Phase*180/PI);
                       
6 ?8 _% {3 b) p9 E, x. e, f# d# l                }
                printf("\n");
2 ?! z9 Z+ ~: Q) [8 r. R                        ci++;
        }
        while(fabs(max1)>0.00001&&fabs(max2)>0.00001&&ci<40);
这是我求潮流的程序，用的PQ分解法，最后得到的结果是只能精确到小数点后第二位，第三位就不对了。

electri

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

ordinarymengya

回复 2# lllmos


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

ordinarymengya

回复 4# lllmos

, B8 k7 b% C& Z( n建议楼主换几个软件试试，SPS,PSSE，PowerSystemWorld，特别是PSSE算潮流比较专业。