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

 找回密码
 立即加入
搜索

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

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

 荐 [复制链接]
咻咻永远

该用户从未签到

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

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

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

×
各位~~我根据下图的一个系统用MATLAB编写了一段程序，是电流型前推回代法的三相配电系统潮流计算，但是现在结果不是我想要的，不收敛。请各位帮我看一下，提提建议。图片见附件了。还有个问题：已知的线路都是三相对称线路的阻抗值，那么线路的阻抗矩阵中的互阻抗应该如何计算，是什么样的形式，应用怎样的公式？这里我不是很清楚。
  1. %程序名：qiantuihuitui_I_3.m $ ]9 [8 f7 L/ U
  2. %功能：支路电流前推回推法求解潮流
    ! s" Q1 \4 \, n7 w
  3. clc ; j% b6 e1 J6 |# }  e! g
  4. clear all; $ h4 B9 C3 T' d5 I% |9 r  r# g9 m
  5. %--------------输入网络参数--------------
    ) t, g; V( U. z5 p( X
  6. %1-支路编号，2-首节点，3-尾节点，4-自阻抗，5-尾节点复功率，6-支路性质（1-馈线段支路，2-变压器支路）,7-尾节点是否带负荷
    ' `! l4 p" ]: W2 B8 ^, Z( u: K
  7. DB=[1 1 2 0.000167+j*0.000208 0.42+j*0.31 1 1
    0 |  p$ C- q5 K- a; h
  8.     2 2 3 0.000151+j*0.000188 6.15 1 09 @0 e2 l. e- h3 \7 x2 r
  9.     3 2 4 0.000066+j*0.000082 0.38+j*0.29 1 1
    8 E+ o+ D- S9 Q7 K% t# s6 q
  10.     4 2 5 0.000249+j*0.000310 0 1 0
    ( F9 L0 {# t# u1 ~: a
  11.     5 2 6 0.000172+j*0.000215 0 1 0: V* {; A# C5 ~
  12.     6 4 7 0.000156+j*0.000195 6.06 1 0" J' @" J( A4 P: ~( u) |6 C
  13.     7 4 8 0.000162+j*0.000202 6.04 1 0
    : r; b/ v0 b2 y$ Z& n
  14.     8 4 9 0.000345+j*0.000430 0 1 0
    - T3 v" V* \7 r- J
  15.     9 4 10 0.000287+j*0.000358 0 1 0# S0 R4 R& _3 G  e
  16.     10  5 11 0.020563+j*0.321594 0 2 0 / d7 ~, @. ^: O, K) m& E! c
  17.     11  6 12 0.020563+j*0.321594 0 2 0" v& {% G2 @- V' M7 y4 }& n  [
  18.     12  9 13 0.020563+j*0.321594 0 2 0 % [2 G! J% r" b0 `& t; d* u: M
  19.     13  10 14 0.020563+j*0.321594 0 2 0  
    ' V7 d8 K& g1 a, {
  20.     14  11 15 0.000237+j*0.000408 5.72+j*0.12 1 1 / o( l$ a& n1 a- L' H, @
  21.     15  12 16 0.000237+j*0.000408 5.76+j*0.09 1 15 }: a# z3 n2 W, N. t
  22.     16  13 17 0.000292+j*0.000502 5.86+j*0.11 1 1: u8 ~0 Q# Q0 X
  23.     17  14 18 0.000274+j*0.000470 5.81+j*0.14 1 1]; ; e4 d  P! J0 e9 f5 ^1 k$ o- o
  24. [n,m]=size(DB);  4 [" X# H0 ]  M# b6 y2 ^
  25. B=[1 sin(2*pi/3)+j*cos(2*pi/3) sin(4*pi/3)+j*cos(4*pi/3)];' j, q. K9 P8 W$ o% O8 H+ H6 `
  26. C=[1.02 1 1.02 1 1 1 1.02 1.02 1 1 1 1 1 1 1 1 1 1];- n) n. X, `8 k) N/ y
  27. U(:,1)=B(1,1)*C';6 x) C4 Q, g/ _" T3 b. A: z* \$ H
  28. U(:,2)=B(1,2)*C';
    6 t$ f3 ^, }5 M* K; k* V1 k7 m
  29. U(:,3)=B(1,3)*C';0 Q+ w# p& C% i0 P
  30. %-------------------求解潮流-----------------
    1 F7 u* }6 J: @
  31. for k=1:15 # |' Q1 N! G7 r5 ~& }- \
  32. %   I(:,k)=((DB(:,5).')*(diag(1./U(:,k))))'; ) n% B, x# D' \+ h- ~
  33. for i=n:-1:1 5 j6 L2 g2 ~  Z4 B& T2 }8 C# T
  34.   %如果尾节点带恒功率负荷，需计算节点上负荷注入的电流) H9 [6 }; f+ K3 E! M" N; W
  35. if DB(i,7)==1( {. w2 O8 |/ C" ~
  36. c=DB(i,5)/3;
    2 Z- v2 X* i: ]$ p1 R0 |1 z
  37. d=c/U(i+1,(3*k-2));
    , ^+ I9 R# t8 z( _- @. Y
  38. IL(i,3*k-2)=conj(d);! g5 x! P9 P" M5 B; r; R% I" t. ]
  39. d=c/U(i+1,3*k-1);
    5 K& M5 c0 G( n$ ~8 _
  40. IL(i,3*k-1)=conj(d);6 v1 p, C  M8 Q& i' m$ x8 W' K
  41. d=c/U(i+1,3*k);
    7 @+ d5 P6 ]1 v& Z' G$ a
  42. IL(i,3*k)=conj(d);( V' z( U- f& d' m6 f% A8 ~. w
  43. else
    4 a" ]2 C. I3 C9 n
  44. IL(i,3*k-2)=0;' Z; _. ?1 r" O- V& k# Q7 m2 \
  45. IL(i,3*k-1)=0;0 K2 E7 j& D$ q7 y; ]) J7 F
  46. IL(i,3*k)=0;! e  m2 V% m: ]3 |0 J
  47. end- g. j0 Q: N% _# M3 ?0 }+ i
  48. %找出所有与尾节点相连的支路,计算进支电流（末端电流），存IKj
    0 D9 G$ }9 P0 [! l; ?8 p; h8 a, Z. w
  49. A=(find(DB(:,2)==(i+1)));7 j4 V4 Q$ a  Z  C# \
  50. if isempty(A)~=1" T) `+ S& ]( A
  51. IKj(i,3*k-2)=IL(i,3*k-2)-sum(IKj(A,3*k-2));
    & q& v/ t' e8 n; ^+ ^
  52. IKj(i,3*k-1)=IL(i,3*k-1)-sum(IKj(A,3*k-1));* q4 d* z5 T4 V* j8 Y5 G7 N
  53. IKj(i,3*k)=IL(i,3*k)-sum(IKj(A,3*k));- Y+ {4 T+ ~* ~' }, x! x
  54. else
    7 }9 X9 c% l) v  l
  55. IKj(i,3*k-2)=IL(i,3*k-2);
    , e7 J' r, ~2 b1 V" F
  56. IKj(i,3*k-1)=IL(i,3*k-1);  z: [; A/ J! _+ ~
  57. IKj(i,3*k)=IL(i,3*k);# \; H) P1 H+ C" ?* e% E
  58. end# n% [$ u% z  j4 s1 l
  59. %计算出支电流（始端电流），存IKi, h0 r/ @+ m3 }
  60. a=DB(i,2);
    ( V* d' S" f. t- E6 W( q
  61. b=DB(i,3);
    3 C; q( T' W8 h
  62. Yi=1/DB(i,4)*eye(3);1 p" I9 g9 X( S7 Y5 O+ u; U: f8 D
  63. if DB(A,6)==11 V. m2 ]$ q' p& z6 e
  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)]';
    8 w4 o, Y; n. F/ Y
  65. IKi(i,3*k-2)=F(1,1);
    % _" U  M9 @: s
  66. IKi(i,3*k-1)=F(2,1);
    # e6 J  ?! A+ }8 O) n7 F
  67. IKi(i,3*k)=F(3,1);) x! Q8 g8 H3 Y) O& r
  68. else: M+ \: S1 P: M  L9 m
  69. YT=1/real(DB(i,4))+j*(1/imag(DB(i,4)));
    " V! H3 Y# y+ W- J! a. J
  70. YI=YT*eye(3);
    7 c- W/ U" x5 G4 b% G4 [3 B
  71. YII=YT*eye(3);
    4 @% Z* \7 a- x+ O* [
  72. YIII=-YT*eye(3);
    , S' U# Q1 @: J3 M
  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)]'));
    8 n" f" i7 l8 [& @- c  n  B
  74. U(a,3*k+1)=D(1,1);
    ( h: }7 z0 v, b$ X8 }% H$ i
  75. U(a,3*k+2)=D(2,1);
      w1 F& j/ E1 G5 R2 Y* V
  76. U(a,3*k+3)=D(3,1);
    ( t8 Y; D! z- j) S7 ?# K* C
  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)]';
    # S9 Q: E( ^2 |1 z, Z# o
  78. IKi(i,3*k-2)=E(1,1);, ^  Z; [2 ^9 f5 A5 b6 h! Q
  79. IKi(i,3*k-1)=E(2,1);0 v* m8 D4 F! D' U. m
  80. IKi(i,3*k)=E(3,1);$ ~& G2 M& @7 ~# M/ l* q: v
  81. end 5 Q$ X$ u8 P9 R6 {
  82. end
    . X( J. p7 G4 ]# Z
  83. %前推电压
    # z/ G+ t, g4 _: |0 `1 x) k8 T
  84. for j=2:n
    % ^7 j6 s9 O3 Q3 \: a. R3 d+ q4 o0 W/ P
  85. U(1,3*k-2)=1.02;/ Z3 t& E- u: ^) a$ N
  86. U(1,3*k-1)=1.02*(sin(2*pi/3)+j*cos(2*pi/3));$ g' Q% u: d9 h7 y; o' _& y& v6 B
  87. U(1,3*k)=1.02*(sin(4*pi/3)+j*cos(4*pi/3));
    ! P* A  P! a7 E: i. Q5 p
  88. a=DB(j-1,2);
    # a8 P/ ^6 T5 N$ e9 s2 B# ~/ D
  89. if DB(j-1,6)==1
    ) C+ ~, p+ X) ^
  90. Yi=1/DB(j-1,4)*eye(3);
    : y' R  E! c4 p/ ^
  91. G1=[IKi(a,3*k-2) IKi(a,3*k-1) IKi(a,3*k)]';. H7 Q. ~  s) s# r( Y$ _
  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)]');/ ^0 ?3 m; _1 d9 L
  93. U(j,3*k+1)=G(1,1);
    ) `! M# f& _% I, h: S- k, A% v
  94. U(j,3*k+2)=G(2,1);
    0 Q+ G' ~5 v' h& X1 l
  95. U(j,3*k+3)=G(3,1);
    ( O+ d+ x& t: s, K8 J+ G% y
  96. else: o/ L1 l$ h% ]# h$ U* L, [1 W9 d) p
  97. YT=1/real(DB(i,4))+j*(1/imag(DB(i,4)));
    + g* M, `4 m7 {, u) E& W4 G
  98. YI=YT*eye(3);# k2 p8 e8 g* u0 {' _  I; J
  99. YII=YT*eye(3);* d' U6 c, w& C" Z9 S
  100. YIII=-YT*eye(3);  M; d  K* y. q# ^
  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)]');
    9 O( d) v  D6 A' m; j
  102. U(j,3*k+1)=H(1,1);* ~. w7 i; V* y4 e6 v
  103. U(j,3*k+2)=H(2,1);& f  v) }) G/ K- c; o1 u0 \3 E
  104. U(j,3*k+3)=H(3,1);% Q. ^' V% q. i6 z' E' t
  105. end 2 Z4 ^' Z# M8 ~2 O6 M5 b
  106. end
    & c' A! b3 _( D: v" F& Z
  107. dU1=U(:,3*k+3)-U(:,3*k);
    4 T4 B9 s. m9 v! S  V
  108. dU2=U(:,3*k+2)-U(:,3*k-1);
    $ k# W- u7 J, @& M0 j8 W- f- c
  109. dU3=U(:,3*k+1)-U(:,3*k-2); # }" m8 q6 _! D! g* R8 H
  110. if (max(abs(dU1))<0.001&&max(abs(dU2))<0.001&&max(abs(dU2))<0.001)
    # J/ C4 a  A8 w# L
  111. break; & T; l, J7 E- W3 ^' M  t9 E) h: k& P: O
  112. else            
    $ b& G# O2 ~4 \: o
  113. k=k+1;
    & n8 Q( G" g( U! j% M: j8 A$ y
  114. end
    4 ?- `( \6 ^" \# L+ P9 o; E
  115. end
    2 G$ D0 P* ?. D7 H
  116. U
    1 U, Z# R/ J" ]: z, `7 V2 Y# U4 h
  117. IL
    - P* k; K$ d5 ~1 z
  118. IKj
    " v: |1 H3 O) h8 J8 ]
  119. IKi1 d2 a+ R5 X7 R* U: g1 X4 W
  120. k
    5 F$ g+ m! x; K5 r
复制代码
matlab潮流计算仿真, matlab潮流计算, matlab潮流计算程序, 基于matlab的潮流计算
QQ截图20130517110314.jpg
"真诚赞赏，手留余香"
打赏
还没有人打赏，支持一下
楼主热帖

相关帖子

帖文化：【文明发帖 和谐互动】 社区精神：【创新、交流、互助、共享】
电力研学论坛招聘版主、专家--欢迎您加入中国电力研学网管理团队！
回复

使用道具 举报

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

    • 签到天数: 1657 天

    连续签到: 28 天

    [LV.Master]伴坛终老

    累计签到：3028 天
    连续签到：6 天
    发表于 2013-5-17 14:45:43 | 显示全部楼层
    我运行了楼主的这段代码，发现最后结果矩阵中很多元素都是Nan或Inf，是不是潮流计算阻抗矩阵中多数情况是这样的，还是不应该很多元素是NaN或Inf
    "真诚赞赏，手留余香"
    打赏
    还没有人打赏，支持一下
    帖文化：【文明发帖 和谐互动】 社区精神：【创新、交流、互助、共享】
    小贴士：资料共享版块上传附件“专家分”值+5，“学分”+1
    回复 推荐 踩下

    使用道具 举报

    咻咻永远

    该用户从未签到

    尚未签到
     楼主| 发表于 2013-5-18 08:05:22 | 显示全部楼层
    回复 2# redplum
    0 n/ T) [0 n# j+ c+ O& d, U) p8 G4 D0 n4 f1 U  [4 w4 A) ^
    ) @* J( M/ b, E* C/ L
        不应该是这样的，这是一个迭代计算，结果不收敛，所以还是哪个环节有问题，但是我还没有找到。
    "真诚赞赏，手留余香"
    打赏
    还没有人打赏，支持一下
    帖文化：【文明发帖 和谐互动】 社区精神：【创新、交流、互助、共享】
    电力研学论坛招聘版主、专家--欢迎您加入中国电力研学网管理团队！
    回复 推荐 踩下

    使用道具 举报

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

    • 签到天数: 1657 天

    连续签到: 28 天

    [LV.Master]伴坛终老

    累计签到：3028 天
    连续签到：6 天
    发表于 2013-5-18 13:43:28 | 显示全部楼层
    这段代码的原始数据有17行7列，K=1:15计算了15次循环。建议将数据弄少一些，这样才能更好判断到底是哪里不收敛
    "真诚赞赏，手留余香"
    打赏
    还没有人打赏，支持一下
    帖文化：【文明发帖 和谐互动】 社区精神：【创新、交流、互助、共享】
    电力研学论坛招聘版主、专家--欢迎您加入中国电力研学网管理团队！
    回复 推荐 踩下

    使用道具 举报

    咻咻永远

    该用户从未签到

    尚未签到
     楼主| 发表于 2013-5-21 16:38:26 | 显示全部楼层
    回复 4# redplum ) a; ^# O0 R2 n- Y+ C2 ]
    就是即使我弄成迭代次数为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-3-17 01:38

    Powered by Discuz! X3.5 Licensed

    © 2001-2025 Discuz! Team.

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