设为首页收藏本站|繁體中文 快速切换版块
切换到宽版

 找回密码
 立即加入
搜索

电力研学论坛»交流社区 › 学术与仿真软件优秀版区 › MATLAB › 配电系统三相潮流计算 matlab仿真 结果不正确 求解答 ...
返回列表 发新帖
查看: 3446|回复: 5

配电系统三相潮流计算 matlab仿真 结果不正确 求解答

 荐 [复制链接]
咻咻永远

该用户从未签到

尚未签到
发表于 2013-5-17 11:12:13 | 显示全部楼层 |阅读模式

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

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

×
各位~~我根据下图的一个系统用Matlab编写了一段程序，是电流型前推回代法的三相配电系统潮流计算，但是现在结果不是我想要的，不收敛。请各位帮我看一下，提提建议。图片见附件了。还有个问题：已知的线路都是三相对称线路的阻抗值，那么线路的阻抗矩阵中的互阻抗应该如何计算，是什么样的形式，应用怎样的公式？这里我不是很清楚。
  1. %程序名：qiantuihuitui_I_3.m . j5 b! s8 [) C7 E8 @" K/ }
  2. %功能：支路电流前推回推法求解潮流 * M( F, v6 ?3 j+ x
  3. clc
    ; K9 p' M% G& Q9 Y' g6 S; d1 ^8 @) f
  4. clear all;
    3 @* ?! l0 N* d
  5. %--------------输入网络参数--------------
    & r( ^/ I9 p& H0 A6 d
  6. %1-支路编号，2-首节点，3-尾节点，4-自阻抗，5-尾节点复功率，6-支路性质（1-馈线段支路，2-变压器支路）,7-尾节点是否带负荷
    4 W! s% U2 O  j& ~/ @/ M
  7. DB=[1 1 2 0.000167+j*0.000208 0.42+j*0.31 1 1
    4 I" Z  d" j- m* \
  8.     2 2 3 0.000151+j*0.000188 6.15 1 0
    : p% n* x) ~5 U
  9.     3 2 4 0.000066+j*0.000082 0.38+j*0.29 1 1
    3 O8 d  M" F5 V$ v
  10.     4 2 5 0.000249+j*0.000310 0 1 0
    6 Q- L/ {3 e# e8 R
  11.     5 2 6 0.000172+j*0.000215 0 1 00 {" D" v! ]& C8 G( s# N0 M4 c3 n
  12.     6 4 7 0.000156+j*0.000195 6.06 1 06 |: _; A9 W+ H3 l* E7 a/ e
  13.     7 4 8 0.000162+j*0.000202 6.04 1 0
    ! J1 @8 H4 P2 P+ c: c
  14.     8 4 9 0.000345+j*0.000430 0 1 0 & j0 k2 T2 f2 e
  15.     9 4 10 0.000287+j*0.000358 0 1 0
    # ?/ o9 Q# Q8 }+ c4 m+ w" J
  16.     10  5 11 0.020563+j*0.321594 0 2 0
    ) t7 ?% ], v1 \: s* L2 f5 P
  17.     11  6 12 0.020563+j*0.321594 0 2 0
    1 t+ G- a" J( l3 E
  18.     12  9 13 0.020563+j*0.321594 0 2 0 5 T4 {9 z5 V: i2 p4 b" m  S
  19.     13  10 14 0.020563+j*0.321594 0 2 0  4 r  Q3 ^# v! ~. N
  20.     14  11 15 0.000237+j*0.000408 5.72+j*0.12 1 1 ( `7 i: A8 [- n
  21.     15  12 16 0.000237+j*0.000408 5.76+j*0.09 1 1
    6 t6 e' K' {+ U& e7 e' y, d
  22.     16  13 17 0.000292+j*0.000502 5.86+j*0.11 1 1+ f  T# s+ j& f* c/ d' }/ |. j
  23.     17  14 18 0.000274+j*0.000470 5.81+j*0.14 1 1];
    1 f8 a0 g0 ^9 H$ q* l
  24. [n,m]=size(DB);  * C! e5 G- @& ]7 b/ {6 a5 W6 @3 C
  25. B=[1 sin(2*pi/3)+j*cos(2*pi/3) sin(4*pi/3)+j*cos(4*pi/3)];" C' \' J1 f. W# b& Y# [
  26. C=[1.02 1 1.02 1 1 1 1.02 1.02 1 1 1 1 1 1 1 1 1 1];: M& q7 \  k5 I
  27. U(:,1)=B(1,1)*C';7 Q' b8 F" K* B7 V, X: g
  28. U(:,2)=B(1,2)*C';
    2 N; ?0 p) L1 M0 k3 o. @4 i* d) }
  29. U(:,3)=B(1,3)*C';
    + v3 P5 J6 H! G% l) o
  30. %-------------------求解潮流----------------- 2 l. d6 m: P+ f, `9 X* q+ \
  31. for k=1:15 ; x# \( c# B& {2 G$ I# f
  32. %   I(:,k)=((DB(:,5).')*(diag(1./U(:,k))))';
    & I9 v/ V) v/ p+ _: r2 B3 ^! }
  33. for i=n:-1:1
    " \& @1 j( c6 A) N( |2 D
  34.   %如果尾节点带恒功率负荷，需计算节点上负荷注入的电流4 s6 p+ V; s" n# |
  35. if DB(i,7)==1
    8 \% x0 h: Z! P( q
  36. c=DB(i,5)/3;
    ) g. Z0 h4 u. Q
  37. d=c/U(i+1,(3*k-2));. t' Q9 w% w( E8 n& D/ V3 G, Y9 e
  38. IL(i,3*k-2)=conj(d);( ]7 j7 [8 T. I- G5 }& E
  39. d=c/U(i+1,3*k-1);  H, v. m4 \4 `/ @5 D
  40. IL(i,3*k-1)=conj(d);
    . e" r0 D% q% Y" ^0 c! h8 C% q
  41. d=c/U(i+1,3*k);8 F+ t5 D+ J; H" s8 ]
  42. IL(i,3*k)=conj(d);1 B$ v+ ?, \2 b* Z: D, _
  43. else
    + w" J& S8 I; A  c2 H
  44. IL(i,3*k-2)=0;4 a# }4 E, f1 L4 o: ]
  45. IL(i,3*k-1)=0;5 R7 q  P8 O: b1 L- U5 ^
  46. IL(i,3*k)=0;. Q6 |% N/ @  D
  47. end
    & v6 h: V. t) o3 ]- y: X" L
  48. %找出所有与尾节点相连的支路,计算进支电流（末端电流），存IKj
    & c& w$ e3 @( s
  49. A=(find(DB(:,2)==(i+1)));
    % m4 r: t3 Z, B1 N" l
  50. if isempty(A)~=1
    : o3 a6 }" E! |( i9 J8 ^! \. D
  51. IKj(i,3*k-2)=IL(i,3*k-2)-sum(IKj(A,3*k-2));
    5 h  Z% A: ?. I6 W, ^( ~2 h* F$ ~; W
  52. IKj(i,3*k-1)=IL(i,3*k-1)-sum(IKj(A,3*k-1));! y8 [) O4 n1 l4 _$ ]! k( `
  53. IKj(i,3*k)=IL(i,3*k)-sum(IKj(A,3*k));
      [! q; @, D6 _& v7 ]* O5 `
  54. else. j$ }. _6 j0 o, U" W, I
  55. IKj(i,3*k-2)=IL(i,3*k-2);
    ! Q* V& M7 K* F, }  M" Y2 v
  56. IKj(i,3*k-1)=IL(i,3*k-1);
    % M% H2 U' `* u7 h' s  P$ a
  57. IKj(i,3*k)=IL(i,3*k);5 M8 b; Z  r! _4 M% P
  58. end# I" Z! P2 I9 M2 J4 X6 Y  ]
  59. %计算出支电流（始端电流），存IKi0 a( }0 n/ b( W6 P) Y+ {- M7 D2 g
  60. a=DB(i,2);
    ! K2 s: R# `2 w) b0 T3 c  ^
  61. b=DB(i,3);
    ' ^" ?) [6 D7 U' r# c* r
  62. Yi=1/DB(i,4)*eye(3);
    0 d. |2 Q& F3 k, ~! o
  63. if DB(A,6)==1
    # E) {4 Z  T7 O; X1 k
  64. F=0.5*Yi*[(U(a,3*k-2)+U(b,3*k-2)) (U(a,3*k-1)+U(b,3*k-1)) (U(a,3*k)+U(b,3*k))]'+[IKj(i,3*k-2) IKj(i,3*k-1) IKj(i,3*k)]';! A" o: y( I* _. z6 q! a% l
  65. IKi(i,3*k-2)=F(1,1);
    " T1 X4 g9 v' ^' B7 E
  66. IKi(i,3*k-1)=F(2,1);
    ! _. Z* ]. |, W
  67. IKi(i,3*k)=F(3,1);
    , h/ d7 n% ?5 k# m  U* {1 ~
  68. else( @- G+ i* J4 `/ @& j+ v* p/ c! ]
  69. YT=1/real(DB(i,4))+j*(1/imag(DB(i,4)));. f) T: ^! Y* T: F! N
  70. YI=YT*eye(3);
    - C* w' L3 c: T8 W
  71. YII=YT*eye(3);
    6 X5 @3 u& f5 v( r
  72. YIII=-YT*eye(3);1 {9 l. ^/ H2 K+ b& B
  73. D=inv(YIII)*(-[IKj(i,3*k-2) IKj(i,3*k-1) IKj(i,3*k)]'-YII*([U(b,3*k-2) U(b,3*k-1) U(b,3*k)]'));. w, j# [9 b4 i2 S. C
  74. U(a,3*k+1)=D(1,1);4 Z) E; @$ U# l% {3 y  j" F
  75. U(a,3*k+2)=D(2,1);& @; f$ f/ I0 s$ P0 e8 v  _6 Y
  76. U(a,3*k+3)=D(3,1);4 ^5 b8 L& N6 \. c% b- F5 L. j/ G; f
  77. E=YI*[U(a,3*k+1) U(a,3*k+2) U(a,3*k+3)]'-YIII*[U(b,3*k-2) U(b,3*k-1) U(b,3*k)]';
    8 Z5 a6 F9 k5 K& c: e
  78. IKi(i,3*k-2)=E(1,1);! \$ M/ c( R9 H
  79. IKi(i,3*k-1)=E(2,1);4 |2 R; y4 U) p" d& x7 y
  80. IKi(i,3*k)=E(3,1);
    ( _) c2 c2 x/ U6 ?. q& e( g
  81. end
    * L  o: a3 Q) R1 C1 K1 x
  82. end
    1 e. l0 ~0 r7 Z1 o& q
  83. %前推电压 " o- D( V6 d4 f$ V
  84. for j=2:n
    , D. {4 a. [6 \$ L$ [5 {+ ?5 r: @
  85. U(1,3*k-2)=1.02;
    3 E+ |' T- [2 m" w! I2 m" [
  86. U(1,3*k-1)=1.02*(sin(2*pi/3)+j*cos(2*pi/3));, _# P6 _$ }& o/ ?( ~
  87. U(1,3*k)=1.02*(sin(4*pi/3)+j*cos(4*pi/3));' x3 M1 d& l* Z2 O, p  M
  88. a=DB(j-1,2);4 A1 e0 }7 y4 W
  89. if DB(j-1,6)==1
    - h  Y. [& v4 f' L# q# s
  90. Yi=1/DB(j-1,4)*eye(3);
    * {: p$ @8 j; |! Z1 L; j4 j5 ?
  91. G1=[IKi(a,3*k-2) IKi(a,3*k-1) IKi(a,3*k)]';: I1 |# X0 I; ^8 t" ^0 `. j+ S
  92. G=[U(a,3*k-2) U(a,3*k-1) U(a,3*k)]'-DB(j-1,4)*eye(3)*(G1-0.5*Yi*[U(a,3*k-2) U(a,3*k-1) U(a,3*k)]');
    8 f; g  e& `, q0 b6 ^: `* `
  93. U(j,3*k+1)=G(1,1); 9 G, [) E% X! \" p; d
  94. U(j,3*k+2)=G(2,1);/ w1 ]: {" ~% R
  95. U(j,3*k+3)=G(3,1); : p# ^3 K2 E: P9 |/ q! q
  96. else
    . z. l3 }6 B- H$ \' n4 t
  97. YT=1/real(DB(i,4))+j*(1/imag(DB(i,4)));
      c7 `/ Q0 t: v7 @/ q% I. q
  98. YI=YT*eye(3);
    " u( I, h  G1 ?8 d
  99. YII=YT*eye(3);
    ! D6 ~5 {; s' |
  100. YIII=-YT*eye(3);  @" N. b2 }* }1 Z3 q3 Z
  101. H=inv(YIII)*([IKi(a,3*k-2) IKi(a,3*k-1) IKi(a,3*k)]'-YI*[U(a,3*k+1) U(a,3*k+2) U(a,3*k+3)]');
    0 ?2 q3 J1 R% F9 j& N, g, X
  102. U(j,3*k+1)=H(1,1);" P% p4 K$ K5 \. |  A; h3 G) @
  103. U(j,3*k+2)=H(2,1);* d# I. d/ G2 a" g4 D0 c! m* i( x
  104. U(j,3*k+3)=H(3,1);
    * }" l; P5 R; A) o" y6 [5 W  [* C
  105. end $ z! c, M  Z0 S% @4 y6 J: S- {9 B) x
  106. end5 T& h8 i  j8 V1 R: O: E
  107. dU1=U(:,3*k+3)-U(:,3*k); ) b7 z/ w% Q/ c" ?6 E: e5 P& q
  108. dU2=U(:,3*k+2)-U(:,3*k-1); 3 L/ A$ S0 i3 H4 o1 V0 ^6 q
  109. dU3=U(:,3*k+1)-U(:,3*k-2);
    ' {6 ^1 D$ r0 O- w1 y% q: o
  110. if (max(abs(dU1))<0.001&&max(abs(dU2))<0.001&&max(abs(dU2))<0.001) & g/ Q1 j: k& M: |! B2 F
  111. break;
    2 Z# u6 ]# b% x  c, x+ s5 ~
  112. else             * V, W0 p- r% F" I) r
  113. k=k+1;
    $ C/ o7 U& N+ e/ |
  114. end + G% _) F" }; i' v0 F
  115. end
    1 G8 [& r% G# ^+ o5 w% h# D( ~
  116. U4 }. @7 i% m! k/ v" C
  117. IL
    ; v- T/ d' ?$ E! Y8 p7 M
  118. IKj
    5 @2 }% a$ |0 H0 h( {& e+ ]- q4 P# _
  119. IKi2 a, a- b4 P5 T4 b6 R
  120. k
    % {7 t8 K( Z- u! m1 ^! I2 t" X& E9 k
复制代码
matlab潮流计算仿真, matlab潮流计算, matlab潮流计算程序, 基于matlab的潮流计算
QQ截图20130517110314.jpg
"真诚赞赏，手留余香"
打赏
还没有人打赏，支持一下
楼主热帖

相关帖子

帖文化：【文明发帖 和谐互动】 社区精神：【创新、交流、互助、共享】
参与“宣传中心”推广活动！让更多的人参与进来，让更多的人可以帮助您！+学分+金币
回复

使用道具 举报

redplum
  • TA的每日心情
    愤怒
    2021-6-12 00:00

    • 签到天数: 1657 天

    连续签到: 28 天

    [LV.Master]伴坛终老

    累计签到：3132 天
    连续签到：3 天
    发表于 2013-5-17 14:45:43 | 显示全部楼层
    我运行了楼主的这段代码，发现最后结果矩阵中很多元素都是Nan或Inf，是不是潮流计算阻抗矩阵中多数情况是这样的，还是不应该很多元素是NaN或Inf
    "真诚赞赏，手留余香"
    打赏
    还没有人打赏，支持一下
    帖文化：【文明发帖 和谐互动】 社区精神：【创新、交流、互助、共享】
    电力研学论坛招聘版主、专家--欢迎您加入中国电力研学网管理团队！
    回复 推荐 踩下

    使用道具 举报

    咻咻永远

    该用户从未签到

    尚未签到
     楼主| 发表于 2013-5-18 08:05:22 | 显示全部楼层
    回复 2# redplum
    " x' ]& o& q2 P7 I" M: w
    * w5 d) ?  j: V* i7 e% y. r8 e$ ^1 x" P% r2 a; `
        不应该是这样的，这是一个迭代计算，结果不收敛，所以还是哪个环节有问题，但是我还没有找到。
    "真诚赞赏，手留余香"
    打赏
    还没有人打赏，支持一下
    帖文化：【文明发帖 和谐互动】 社区精神：【创新、交流、互助、共享】
    小贴士：资料共享版块上传附件“专家分”值+5，“学分”+1
    回复 推荐 踩下

    使用道具 举报

    redplum
  • TA的每日心情
    愤怒
    2021-6-12 00:00

    • 签到天数: 1657 天

    连续签到: 28 天

    [LV.Master]伴坛终老

    累计签到：3132 天
    连续签到：3 天
    发表于 2013-5-18 13:43:28 | 显示全部楼层
    这段代码的原始数据有17行7列，K=1:15计算了15次循环。建议将数据弄少一些，这样才能更好判断到底是哪里不收敛
    "真诚赞赏，手留余香"
    打赏
    还没有人打赏，支持一下
    帖文化：【文明发帖 和谐互动】 社区精神：【创新、交流、互助、共享】
    参与“宣传中心”推广活动！让更多的人参与进来，让更多的人可以帮助您！+学分+金币
    回复 推荐 踩下

    使用道具 举报

    咻咻永远

    该用户从未签到

    尚未签到
     楼主| 发表于 2013-5-21 16:38:26 | 显示全部楼层
    回复 4# redplum
    ; I' O5 d: |% n1 f  C2 `3 n就是即使我弄成迭代次数为2，结果也不啊~~
    "真诚赞赏，手留余香"
    打赏
    还没有人打赏，支持一下
    帖文化：【文明发帖 和谐互动】 社区精神：【创新、交流、互助、共享】
    小贴士：资料共享版块上传附件“专家分”值+5，“学分”+1
    回复 推荐 踩下

    使用道具 举报

    魏丹

    该用户从未签到

    尚未签到
    发表于 2014-6-17 20:04:48 | 显示全部楼层
    这个问题解决了吗 求答案~
    帖文化：【文明发帖 和谐互动】 社区精神：【创新、交流、互助、共享】
    参与“宣传中心”推广活动！让更多的人参与进来，让更多的人可以帮助您！+学分+金币
    回复 推荐 踩下

    使用道具 举报

    返回列表 发新帖
    高级模式
    B Color Image Link Quote Code Smilies
    您需要登录后才可以回帖 登录 | 立即加入

    本版积分规则

    招聘斑竹

    小黑屋|手机版|APP下载(beta)|Archiver|电力研学网 ( 赣ICP备12000811号-1|赣公网安备36040302000210号 )|网站地图

    GMT+8, 2026-8-12 04:48

    Powered by Discuz! X3.5 Licensed

    © 2001-2026 Discuz! Team.

    快速回复 返回顶部 返回列表