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

20031374

网上找到计算潮流的版本之一如下：
/ `; k* N4 t/ s3 U6 U看到一个帖子：http://www.ilovematlab.cn/viewth ... =page%3D1#pid230344
有谁可以解释一下B2各列参数的意义吗？谢谢!
7 @/ W1 y- C9 n3 r7 y9 y请输入节点数:n=12
0 N  u3 j- O5 X7 V. j请输入支路数:nl=11
5 s$ u4 Q# M2 w请输入平衡母线节点号:isb=1
! U  Z$ H. q' i0 I7 B% ~2 G% H请输入误差精度:pr=0.00001
请输入由各支路参数形成的矩阵: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]
# G) R5 P( u+ s7 ~请输入各节点参数形成的矩阵: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]
+ @- M' {/ m' v% z  G+ q导纳矩阵Y=
Columns 1 through 7
& h% X  Y6 y- ~5 K31.2739 -83.2494i -2.7123 +14.1874i -7.8023 +16.2549i -9.6108 +19.7021i -6.7873 +14.3288i -4.2052 +13.8772i        0         
9 n" N' g' E% f& g8 ]+ M-2.7123 +14.1874i   2.7123 -13.9634i        0                  0                  0                  0                  0         
-7.8023 +16.2549i        0             7.9751 -22.2912i        0                  0                  0                  0         
1 x( \+ L, P' ?, Z/ p! g-9.6108 +19.7021i        0                  0             9.7768 -25.1143i        0                  0                  0         
0 @3 j. O9 L; g6 s-6.7873 +14.3288i        0                  0                  0             6.9534 -19.6955i        0                  0         
6 @! Z# q1 a  }-4.2052 +13.8772i        0                  0                  0                  0            10.0733 -37.9479i -5.7121 +18.6596i
3 G0 P' F  o* y3 F) Z- H        0                  0                  0                  0                  0            -5.7121 +18.6596i   5.7121 -18.5071i
8 [/ F; }9 ]" C" g3 E        0                  0            -0.1641 + 5.8760i        0                  0                  0                  0         
% V0 S* c  r4 g7 z2 `        0                  0                  0            -0.1578 + 5.2584i        0                  0                  0         
3 r. {- v- G% N4 c        0                  0                  0                  0            -0.1578 + 5.2584i        0                  0         
6 U* V! K: M1 \* a' E5 \# W( }3 C9 n        0                  0                  0                  0                  0            -0.1559 + 5.5822i        0         
-0.1686 + 6.0349i        0                  0                  0                  0                  0                  0         
; y; ?8 ^, }6 j+ B% q6 A) s  h3 q& GColumns 8 through 12
        0                  0                  0                  0            -0.1686 + 6.0349i
3 b7 t$ Z; c, T+ c" b& C- I        0                  0                  0                  0                  0         
8 G  ^6 F+ [8 p. ?6 c5 G-0.1641 + 5.8760i        0                  0                  0                  0         
5 Q/ X  C$ A/ N7 M; B. L5 G$ X2 [        0            -0.1578 + 5.2584i        0                  0                  0         
, N, A, q, L% {        0                  0            -0.1578 + 5.2584i        0                  0         
7 T* |/ y9 C* Y8 {1 e: L        0                  0                  0            -0.1559 + 5.5822i        0         
        0                  0                  0                  0                  0         
   0.1559 - 5.5822i        0                  0                  0                  0         
        0             0.1499 - 4.9955i        0                  0                  0         
        0                  0             0.1499 - 4.9955i        0                  0         
3 r- ]0 E; w; D( X# ^: [, c7 o        0                  0                  0             0.1559 - 5.5822i        0         
4 P# O* M  m( i8 _# \1 k3 Z+ X8 g# ]        0                  0                  0                  0             0.1822 - 6.5242i
% H) y# F. M! u' X# W迭代次数
& I$ M" E8 J7 s6 F) k1 X$ O     4
没有达到精度要求的个数
) I$ V3 k7 v& K5 {    21    21    14     0
: h. Q5 `$ ?! y6 d各节点的实际电压标么值E为(节点号从小到大排列):
% w* Q( a9 P3 G, o2 R* JColumns 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
1 D. e3 n) ]# qColumns 8 through 12
   1.0225 - 0.0954i   1.0044 - 0.0833i   0.9983 - 0.0890i   1.0110 - 0.0197i   1.0211 + 0.0930i
各节点的电压大小V为(节点号从小到大排列):
    1.0500    1.0000    1.0295    1.0305    1.0256    1.0689    1.1000    1.0269    1.0079    1.0023    1.0112    1.0253
' L4 }4 k/ W; i( q7 k2 T" p各节点的电压相角时shita为(节点号从小到大排列):
         0   -5.6958   -1.0221   -0.5827   -0.8925    3.3184    6.6965   -5.3305   -4.7429   -5.0961   -1.1182    5.2027
8 [) y" [4 s- Z各节点的功率S为(节点号从小到大排列):
, @' v& v+ D5 m% C; B. O- QColumns 1 through 7
   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
: j, D& O2 B  DColumns 8 through 12
. z6 O# Q& |4 E% L. k. g; p5 \-0.4753 - 0.2947i -0.4074 - 0.3608i -0.4074 - 0.3608i -0.4753 - 0.2947i   0.6000 + 0.3720i
各条支路的首端功率Si为(顺序同您输入B1时一样):
. \2 x; F) B0 i2 M+ s3 S  [   1.6349 + 0.2888i
0 {. {% O8 @; S+ K7 e! j" t   0.4827 + 0.0380i
   0.4139 + 0.1626i
   0.4150 + 0.0691i
-0.9770 + 0.0037i
: Q4 s# K: I( ]/ D& j+ B( b; L-0.4753 - 0.2947i
; b; O1 @/ r$ E1 j-0.4074 - 0.3608i
-0.4074 - 0.3608i
/ b0 k( u2 F' q( T- Q-0.4753 - 0.2947i
-1.4711 - 0.3611i
0 J+ H7 T# }( B-0.5980 - 0.2994i
! v* d, L8 _/ l5 }! L; M% E+ I% |( q各条支路的末端功率Sj为(顺序同您输入B1时一样):
' E/ M9 X+ N; G5 w7 S2 g-1.6000 - 0.5772i
-0.4768 - 0.3478i
7 ~5 `. R9 M) V2 W. O; g  l-0.4091 - 0.4192i
-0.4092 - 0.4198i
$ T1 W# V5 ^/ ]/ I5 Q7 W, G) H   0.9943 + 0.0116i
6 V$ A  }$ I9 Y: ]# a+ U   0.4768 + 0.3478i
  y4 f7 }" v( d   0.4091 + 0.4192i
, H" Z5 C! j, S1 R: K   0.4092 + 0.4198i
5 Q* a8 ~: D/ z9 X- R' o8 H' d   0.4768 + 0.3495i
9 V* P* L1 w8 @  S) {! k: o2 W2 m   1.5000 + 0.0966i
+ h; |- c1 ~$ q2 k' C/ R   0.6000 + 0.3720i
各条支路的功率损耗DS为(顺序同您输入B1时一样):
   0.0349 - 0.2884i
   0.0060 - 0.3098i
   0.0047 - 0.2566i
   0.0058 - 0.3508i
   0.0173 + 0.0153i
   0.0015 + 0.0531i
; @5 k  y( D3 _) \9 n, R   0.0017 + 0.0583i
" A8 j8 R) a' o* \3 N   0.0018 + 0.0590i
   0.0015 + 0.0548i
$ I% K  q; z$ V! }+ g5 O+ Y   0.0289 - 0.2644i
: ^8 V* Z. Y/ B, @   0.0020 + 0.0726i

20031374

clc
: r( }/ Y# o& b+ \9 }2 D1 gclear
: C% T( M, ~5 z& `8 ^3 U5 d5 yn=input('请输入节点数:n=');
nl=input('请输入支路数:nl=');
isb=input('请输入平衡母线节点号:isb=');
pr=input('请输入误差精度:pr=');
B1=input('请输入由各支路参数形成的矩阵:B1=');
6 T: P- h, K) p/ ~" x- L! SB2=input('请输入各节点参数形成的矩阵:B2=');
0 ]" l% I; [' g. FY=zeros(n); e=zeros(1,n);f=zeros(1,n);V=zeros(1,n);
O=zeros(1,n);S1=zeros(nl);
%求导纳矩阵
2 n, ~5 x  \# ^9 f5 v. ~for i=1:nl
    if B1(i,6)==0 %NO.6
  C0 b2 t% Y6 C9 a9 F- d        p=B1(i,1);q=B1(i,2);
8 b. g! K5 e( X3 Z1 V2 }    else p=B1(i,2);q=B1(i,1);
    end
    Y(p,q)=Y(p,q)-1./(B1(i,3)*B1(i,5));%非对角元素
, _& N6 l" m. T) y' Q, z  h; @    Y(q,p)=Y(p,q);%非对角元素
    Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)^2)+B1(i,4)./2;%变压器支路 NO.5变比
+ [, I( b; C3 z+ t' A% K# O" l! @    Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2;
end
disp('导纳矩阵Y=');
9 ^+ z' r' N6 |$ l3 h/ Xdisp(Y);
G=real(Y);B=imag(Y);
for i=1:n%n=节点数
8 }/ P0 @+ Q1 ^- v    e(i)=real(B2(i,3));%NO.3节点电压初值
    f(i)=imag(B2(i,3));
    V(i)=B2(i,4);%？
end
for i=1:n
! v4 }) I0 @9 o" U    S(i)=B2(i,1)-B2(i,2);% NO.1节点发出的功率，NO.2节点负荷大小，S（i）负荷数据
    B(i,i)=B(i,i)+B2(i,5);% B=imag(Y)，NO.5？
5 u! c1 {) I  e3 U# E- L9 ^end
P=real(S);Q=imag(S);

20031374

B2的第4，5列代表什么含义？
( M5 T$ h7 k( w$ |请输入各节点参数形成的矩阵:B2=[
0        0                          1.05       0 0 1;
0        1.6+0.992i           1.05       1 0 3;  
4 U0 {: ]4 V/ E) I0        0                         1             0 0 2;
0        0                         1            0 0 2;
/ l9 `! T: Y! r9 k$ {" |0        0                          1           0 0 2;
0        0                         1           0 0 2;
1 `  _  a( M% }& U4 A0.7708   0                     1.1        1.1 0 3;
0       0.539+0.33418i   1        0 0 2;   
0       0.66+0.4092i     1       0 0 2;
; A+ h+ X, H1 _' l# z0 r8 g* i0       0.66+0.4092i     1       0 0 2;
0 e# b" {) D+ _' A5 z9 T0       0.77+0.4774i     1        0 0 2;
1+0.62i  0.4+0.248i       1         0 0 2
]

wsz310

谢谢楼主啊,很详细!