基于MATLAB牛顿拉夫逊法进行夏季最小潮流计算

20031374

网上找到计算潮流的版本之一如下：
! ?7 I3 B$ P; b2 p5 L5 @看到一个帖子：http://www.ilovematlab.cn/viewth ... =page%3D1#pid230344
有谁可以解释一下B2各列参数的意义吗？谢谢!
! d# d/ `4 Z: j0 F% Y: Z请输入节点数:n=12
请输入支路数:nl=11
请输入平衡母线节点号:isb=1
请输入误差精度:pr=0.00001
% B5 Z$ m/ L  [+ x  J& R请输入由各支路参数形成的矩阵:B1=[1 2 0.013+0.068i 0.448i 1 0;1 3 0.024+0.050i 0.298i 1 0;1 4 0.020+0.041i 0.246i 1 0;1 5 0.027+0.057i 0.337i 1 0;1 6 0.020+0.066i 0.0373i 1 0;3 8 0.005+0.179i 0 0.95 1;4 9 0.006+0.200i 0 0.95 1;5 10 0.006+0.200i 0 0.95 1;6 11 0.005+0.179i 0 1.0 1;6 7 0.015+0.049i 0.305i 1 0;1 12 0.005+0.179i 0 0.925 0]
请输入各节点参数形成的矩阵:B2=[0 0 1.05 0 0 1;0 1.6+0.992i 1.05 1 0 3;0 0 1 0 0 2;0 0 1 0 0 2;0 0 1 0 0 2;0 0 1 0 0 2;1.5 0 1.1 1.1 0 3;0 0.4753+0.2947i 1 0 0 2;0 0.4074+0.36084i 1 0 0 2;0 0.4074+0.36084i 1 0 0 2;0 0.4753+0.2947i 1 0 0 2;1+0.62i 0.4+0.248i 1 0 0 2]
导纳矩阵Y=
Columns 1 through 7
31.2739 -83.2494i -2.7123 +14.1874i -7.8023 +16.2549i -9.6108 +19.7021i -6.7873 +14.3288i -4.2052 +13.8772i        0         
( ~+ T5 ]( o0 N1 r1 g( O; r-2.7123 +14.1874i   2.7123 -13.9634i        0                  0                  0                  0                  0         
. U6 u0 g- T4 A8 W' `: Z) r-7.8023 +16.2549i        0             7.9751 -22.2912i        0                  0                  0                  0         
-9.6108 +19.7021i        0                  0             9.7768 -25.1143i        0                  0                  0         
-6.7873 +14.3288i        0                  0                  0             6.9534 -19.6955i        0                  0         
8 _% f9 G+ i# v+ Y2 I-4.2052 +13.8772i        0                  0                  0                  0            10.0733 -37.9479i -5.7121 +18.6596i
* G: p+ ]6 l1 {* V4 }        0                  0                  0                  0                  0            -5.7121 +18.6596i   5.7121 -18.5071i
' W$ Y* H5 ~& S/ m& L: b" V7 f        0                  0            -0.1641 + 5.8760i        0                  0                  0                  0         
/ }3 G+ z0 P/ d# Z! b+ Y) Y" a        0                  0                  0            -0.1578 + 5.2584i        0                  0                  0         
7 k" G  f. g7 L  D* u+ d        0                  0                  0                  0            -0.1578 + 5.2584i        0                  0         
' Q" Y7 h6 I$ ^1 V" w( t        0                  0                  0                  0                  0            -0.1559 + 5.5822i        0         
-0.1686 + 6.0349i        0                  0                  0                  0                  0                  0         
Columns 8 through 12
        0                  0                  0                  0            -0.1686 + 6.0349i
/ D$ i2 W2 Q- t& B        0                  0                  0                  0                  0         
-0.1641 + 5.8760i        0                  0                  0                  0         
3 i8 [  J* O. F; z* x        0            -0.1578 + 5.2584i        0                  0                  0         
" Q, L9 W' A- w  `5 {/ X  b) F' e        0                  0            -0.1578 + 5.2584i        0                  0         
1 T  e& v8 |) e6 u, g! C        0                  0                  0            -0.1559 + 5.5822i        0         
$ ^  ]+ y! U6 M. k5 o6 W. o& l2 V        0                  0                  0                  0                  0         
: h2 ^# [% |4 D   0.1559 - 5.5822i        0                  0                  0                  0         
- k0 c* H6 N! I  V4 d1 k        0             0.1499 - 4.9955i        0                  0                  0         
2 |, G& o: p! v8 a        0                  0             0.1499 - 4.9955i        0                  0         
        0                  0                  0             0.1559 - 5.5822i        0         
6 o7 \$ \, W) E* v# n        0                  0                  0                  0             0.1822 - 6.5242i
迭代次数
     4
0 n8 f! E  P$ D4 M# Y没有达到精度要求的个数
& ^! ^* g1 G3 {8 U, B  c- M    21    21    14     0
各节点的实际电压标么值E为(节点号从小到大排列):
% h9 x8 p( V0 s' PColumns 1 through 7
   1.0500             0.9951 - 0.0992i   1.0293 - 0.0184i   1.0305 - 0.0105i   1.0255 - 0.0160i   1.0671 + 0.0619i   1.0925 + 0.1283i
0 w$ J% G8 H2 c0 b- _; F0 B$ F: {& `% [Columns 8 through 12
   1.0225 - 0.0954i   1.0044 - 0.0833i   0.9983 - 0.0890i   1.0110 - 0.0197i   1.0211 + 0.0930i
! v' O5 l1 O3 e8 `8 v# Y7 d( Z各节点的电压大小V为(节点号从小到大排列):
1 o" H/ M/ f: O    1.0500    1.0000    1.0295    1.0305    1.0256    1.0689    1.1000    1.0269    1.0079    1.0023    1.0112    1.0253
- ?, f6 p; P2 Y/ B2 M0 u9 c各节点的电压相角时shita为(节点号从小到大排列):
4 A$ E; d$ b' v# e! Q         0   -5.6958   -1.0221   -0.5827   -0.8925    3.3184    6.6965   -5.3305   -4.7429   -5.0961   -1.1182    5.2027
各节点的功率S为(节点号从小到大排列):
3 c, W! J3 `/ F  L- w* f  YColumns 1 through 7
- d0 n8 ~7 b9 j6 }$ a( \   1.3716 + 0.2627i -1.6000 - 0.5772i -0.0000 + 0.0000i   0.0000 - 0.0000i   0.0000 - 0.0000i -0.0000 - 0.0000i   1.5000 + 0.0966i
8 [! h3 P/ y# FColumns 8 through 12
3 S7 c# M+ u8 y4 _! o" S3 H-0.4753 - 0.2947i -0.4074 - 0.3608i -0.4074 - 0.3608i -0.4753 - 0.2947i   0.6000 + 0.3720i
各条支路的首端功率Si为(顺序同您输入B1时一样):
. _/ h5 z& d* Z3 \& X$ {   1.6349 + 0.2888i
5 s* \; ^% D& l7 J& D* C/ N0 G   0.4827 + 0.0380i
   0.4139 + 0.1626i
& Z9 P+ D: o. b4 ^   0.4150 + 0.0691i
-0.9770 + 0.0037i
-0.4753 - 0.2947i
- k+ h0 z* A/ g2 F4 [0 b-0.4074 - 0.3608i
3 ^# _  J0 d' ^/ j1 p, j9 T, R0 n-0.4074 - 0.3608i
-0.4753 - 0.2947i
-1.4711 - 0.3611i
-0.5980 - 0.2994i
各条支路的末端功率Sj为(顺序同您输入B1时一样):
  c3 r  t: @7 L. v-1.6000 - 0.5772i
-0.4768 - 0.3478i
-0.4091 - 0.4192i
-0.4092 - 0.4198i
3 j2 m7 ~' q$ u, ?   0.9943 + 0.0116i
8 s& r) E4 L- M$ s   0.4768 + 0.3478i
7 N6 s. D$ g0 V2 G" ^   0.4091 + 0.4192i
. j$ E5 A; K$ K: S# K8 c! l   0.4092 + 0.4198i
( X6 W: c8 c; B: D4 b- o   0.4768 + 0.3495i
   1.5000 + 0.0966i
2 e1 L8 g! o9 ^& E* ~   0.6000 + 0.3720i
各条支路的功率损耗DS为(顺序同您输入B1时一样):
' `2 B* r4 b3 L$ M  |  l   0.0349 - 0.2884i
* f$ ~2 z* K$ w( @4 F/ {   0.0060 - 0.3098i
   0.0047 - 0.2566i
, _$ E/ a; d) _" t3 {; Q   0.0058 - 0.3508i
% }+ u/ A$ u3 Y# V8 L   0.0173 + 0.0153i
   0.0015 + 0.0531i
: q% W$ G) R2 L7 w( n/ E1 a   0.0017 + 0.0583i
   0.0018 + 0.0590i
# Z  m- H: i6 z3 q$ W7 K   0.0015 + 0.0548i
   0.0289 - 0.2644i
( o; U2 M2 w+ W   0.0020 + 0.0726i

20031374

clc
clear
n=input('请输入节点数:n=');
nl=input('请输入支路数:nl=');
+ O, G8 N7 X5 K6 P5 c1 u* I3 k8 Zisb=input('请输入平衡母线节点号:isb=');
9 s# r, b- T8 t8 W1 [- v0 Q# C( @3 a  Opr=input('请输入误差精度:pr=');
B1=input('请输入由各支路参数形成的矩阵:B1=');
: G/ j# J9 v$ N& }1 hB2=input('请输入各节点参数形成的矩阵:B2=');
1 a6 @* F/ T, M& |Y=zeros(n); e=zeros(1,n);f=zeros(1,n);V=zeros(1,n);
O=zeros(1,n);S1=zeros(nl);
& \9 ?* r4 v; t, |%求导纳矩阵
for i=1:nl
2 {% ~. J) I- r8 @9 `1 d' G    if B1(i,6)==0 %NO.6
3 H+ I# o% L* K% u, q        p=B1(i,1);q=B1(i,2);
    else p=B1(i,2);q=B1(i,1);
7 ]& N0 i/ P0 v    end
    Y(p,q)=Y(p,q)-1./(B1(i,3)*B1(i,5));%非对角元素
+ b* t1 O+ t! o8 S+ b# S    Y(q,p)=Y(p,q);%非对角元素
' H! _8 |& b+ X" H% R  R    Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)^2)+B1(i,4)./2;%变压器支路 NO.5变比
    Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2;
4 C+ B% @1 W8 o" e1 K# s) nend
# n- M, S6 m: A: Y* u; ddisp('导纳矩阵Y=');
2 @# V' _; k7 Zdisp(Y);
G=real(Y);B=imag(Y);
, c( C/ V7 C& gfor i=1:n%n=节点数
7 e' y2 O9 }7 B; \# `& K9 v; x2 D    e(i)=real(B2(i,3));%NO.3节点电压初值
    f(i)=imag(B2(i,3));
    V(i)=B2(i,4);%？
/ Z# g9 s1 N5 d9 Eend
for i=1:n
, E, d. U& v7 {# ~/ e" y    S(i)=B2(i,1)-B2(i,2);% NO.1节点发出的功率，NO.2节点负荷大小，S（i）负荷数据
4 }; a  v; b8 k) M$ t8 z7 u    B(i,i)=B(i,i)+B2(i,5);% B=imag(Y)，NO.5？
: G5 Y; \& w' send
, `3 H: f) [) T8 c0 CP=real(S);Q=imag(S);

20031374

B2的第4，5列代表什么含义？
0 `* Z8 ?1 @5 D# }9 n请输入各节点参数形成的矩阵:B2=[
0        0                          1.05       0 0 1;
2 J' `4 c9 b* D" Z2 N: j0        1.6+0.992i           1.05       1 0 3;  
4 ~7 A& W" S' o/ ]3 w* c6 a$ i6 s0        0                         1             0 0 2;
3 Y* |) u2 U8 f; u& A0        0                         1            0 0 2;
0        0                          1           0 0 2;
' r" E. F# h! a0        0                         1           0 0 2;
; g0 ?; \) ^' J7 ~. x0.7708   0                     1.1        1.1 0 3;
* ~' S# r" P( \& K4 R1 A0       0.539+0.33418i   1        0 0 2;   
0       0.66+0.4092i     1       0 0 2;
2 e2 X( W$ ]- j5 x8 Z0       0.66+0.4092i     1       0 0 2;
1 R5 q5 m6 s6 L/ s/ e0       0.77+0.4774i     1        0 0 2;
! j2 V3 O: d7 u1+0.62i  0.4+0.248i       1         0 0 2
]

wsz310

谢谢楼主啊,很详细!