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

lllmos · 发表于 2015-6-16 10:44:58

马上加入，结交更多好友，共享更多资料，让你轻松玩转电力研学社区！

您需要登录才可以下载或查看，没有账号？立即加入

×

lllmos · 发表于 2015-6-16 11:17:56

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

玉门关山 · 发表于 2015-6-16 11:33:18

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

lllmos · 发表于 2015-6-16 15:02:40

回复 3# 玉门关山

因为我不知道问题出在哪里了。。
do
{
//求解不平衡量
for(i=0;i<nB;i++)
{
if(sB.Type!=2)//假如不是平衡节点2 V9 H" t, a2 L# ^# Z# H
{
3 l: q( P  R% r# s$ q DP=sB.GenP-sB.LoadP;. S; c  v# ]  U' u
DQ=sB.GenQ-sB.LoadQ;
- g+ D; s" N# f3 h# s- k$ m- \
1 ^9 D3 A; Y  x  k# |: @$ V for(j=0;j<nB;j++)
* R. \. b- V# b0 a& c {
0 R6 Z$ E: n+ m2 F! j0 J* Y A=sB.Phase-sB[j].Phase;
3 ~% |, m# T/ a. ?/ @ DP-=sB.Volt*(sB[j].Volt*(g[j]*cos(A)+b[j]*sin(A)));
5 C5 N* r  P- L, o$ y1 N) m* X$ h
6 N  K3 P" M+ o1 n8 |1 d# O if(sB.Type==0)//PQ节点
4 k. V# N) I+ z; k7 ~ DQ-=sB.Volt*(sB[j].Volt*(g[j]*sin(A)-b[j]*cos(A)));
( T! [" j6 g  M0 {
, U4 l& q/ t6 q7 P else if(sB.Type==1)//PV节点: [% s$ I1 P8 g9 ^( ?
DQ=0;% Q% C# x: W. K4 F% I: V9 x2 E2 N
}4 X/ j8 d) N3 C8 r% ^# v
}
1 O$ k4 \' J5 Z* u7 ^ else if(sB.Type==2)//平衡节点
" t- u" L" K: j3 Y4 } DP=DQ=0;2 T2 H' w" M" h2 i9 B- @) H
}
  E( ?6 y' i7 ~9 a" \ //for(i=0;i<nB;i++)
3 g; G( S* Q! N6 ]0 F4 p //    printf("DP[%d]===%f,DQ[%d]===%f\n",i,DP,i,DQ);7 \. L, n) ?: A. A7 l, ]1 g
  |3 B* H6 ^/ D7 M) a7 Q6 S
      //求解修正方程
0 `- s/ b6 u7 a for(i=0;i<nB-1;i++)
' n# `8 K/ E& h* x* {) ^       AA1=DP[i+1]/sB[i+1].Volt;
  ^% m3 E, F& z/ e" M for(i=0;i<nB-1-count_PVnode;i++)& Y1 v- F$ _! h5 y: H' }
AA2=DQ[i+1+count_PVnode]/sB[i+1+count_PVnode].Volt;
$ Z$ U+ y+ f, }: I: a, Q4 a( w$ k3 g- ] calculate_gaosi((double **)b1,BB1,AA1,NBUS-1);//AA是不平衡量，BB是解向量- I# |* q4 g# B; o% t6 T# E
calculate_gaosi((double **)b2,BB2,AA2,NBUS-1);
; U% {. P$ l2 Y( {4 A$ Z( z. W3 ?  ?, I) t5 T; i$ v4 c( R. O
max1=fabs(AA1[0]);" r+ Z3 ^( n5 h# r) [! ~1 x# X
for(i=1;i<nB-1;i++)
' E  D$ N, B/ b+ E" r if(max1<fabs(AA1))
: [" }/ L9 l) r. }' w9 S max1=fabs(AA1);
3 C0 Y$ h: ?7 x/ u max2=fabs(AA2[0]);
( z+ E4 b" p+ H  h$ u for(i=1;i<nB-1-count_PVnode;i++)6 U; M  x2 z. W/ T( y
if(max2<fabs(AA2)) ( M" t6 R! M# I  ^3 B; w3 k3 _" m
max2=fabs(AA2);& U& {4 Y  q& p' T+ m6 G( _
for(i=0;i<nB-1;i++)5 R. ^7 d! t9 I; u0 O: t* g
sB[i+1].Phase+=BB1/sB.Volt;
1 u6 P: N. e1 H* X' W for(i=0;i<nB-1-count_PVnode;i++)
9 z+ S3 Q" q& E2 E sB[i+1+count_PVnode].Volt+=BB2;/ ]  d: K2 p4 n/ f" ^: C
for(i=0;i<nB;i++)
8 u) s& i2 `" d) Z8 Z9 n3 n {
% b7 c2 U* [' N( E9 L* x, X. S& } printf("sB[%d].Volt=%f,sB[%d].Phase=%f\n",i,sB.Volt,i,sB.Phase*180/PI);
0 P& Y6 U3 ?& f7 y0 b5 H : n9 k0 E% \* \, N6 x+ n) ?
}. ]- h6 D, Y% ~" G+ Z' E' ?' q
printf("\n");
; m! R0 c2 w& B' y4 n" R' M* K ci++;/ h0 L) M" b# s
}+ o  r: d6 T3 Y$ u4 q7 @
while(fabs(max1)>0.00001&&fabs(max2)>0.00001&&ci<40);
8 B/ g% O# j  @, X4 @这是我求潮流的程序，用的PQ分解法，最后得到的结果是只能精确到小数点后第二位，第三位就不对了。

electri · 发表于 2015-6-16 19:38:59

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

ordinarymengya · 发表于 2015-10-4 18:59:11

回复 2# lllmos

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

ordinarymengya · 发表于 2015-10-4 19:00:33

回复 4# lllmos

建议楼主换几个软件试试，SPS,PSSE，PowerSystemWorld，特别是PSSE算潮流比较专业。

账号		自动登录	找回密码
密码			立即加入

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

马上加入，结交更多好友，共享更多资料，让你轻松玩转电力研学社区！

相关帖子