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

20031374

网上找到计算潮流的版本之一如下：
看到一个帖子：http://www.ilovematlab.cn/viewth ... =page%3D1#pid230344
$ u; V. Q: a; j: _4 H5 ?有谁可以解释一下B2各列参数的意义吗？谢谢!
请输入节点数:n=12
: u2 e. @. \4 I  w# W1 t  _请输入支路数:nl=11
请输入平衡母线节点号:isb=1
请输入误差精度: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]
请输入各节点参数形成的矩阵: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         
-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         
3 G( A6 {2 E/ W# b& V-9.6108 +19.7021i        0                  0             9.7768 -25.1143i        0                  0                  0         
, H' T3 V- `2 x" @4 T-6.7873 +14.3288i        0                  0                  0             6.9534 -19.6955i        0                  0         
0 g1 `# ?, d/ I1 V( s# A-4.2052 +13.8772i        0                  0                  0                  0            10.0733 -37.9479i -5.7121 +18.6596i
        0                  0                  0                  0                  0            -5.7121 +18.6596i   5.7121 -18.5071i
        0                  0            -0.1641 + 5.8760i        0                  0                  0                  0         
; u1 D6 ]1 u" J4 d# y* l8 ?: u        0                  0                  0            -0.1578 + 5.2584i        0                  0                  0         
3 r$ C  F0 L! D+ ^1 d4 F5 G+ r        0                  0                  0                  0            -0.1578 + 5.2584i        0                  0         
0 e6 Q( w, \+ W7 V) F+ k        0                  0                  0                  0                  0            -0.1559 + 5.5822i        0         
) N! u- \7 M" r-0.1686 + 6.0349i        0                  0                  0                  0                  0                  0         
! j" T1 ]  Y7 A- p, fColumns 8 through 12
        0                  0                  0                  0            -0.1686 + 6.0349i
        0                  0                  0                  0                  0         
-0.1641 + 5.8760i        0                  0                  0                  0         
        0            -0.1578 + 5.2584i        0                  0                  0         
  ]$ Y- O1 ?' H, S+ G1 `3 S8 Z        0                  0            -0.1578 + 5.2584i        0                  0         
        0                  0                  0            -0.1559 + 5.5822i        0         
1 B2 q- s( Q( T/ I$ Y% O        0                  0                  0                  0                  0         
& u# N. g$ o6 @+ x2 `   0.1559 - 5.5822i        0                  0                  0                  0         
        0             0.1499 - 4.9955i        0                  0                  0         
. i5 N4 Y* n6 V, v! s        0                  0             0.1499 - 4.9955i        0                  0         
1 }4 E7 q+ y; Y0 t; z% w- Z        0                  0                  0             0.1559 - 5.5822i        0         
        0                  0                  0                  0             0.1822 - 6.5242i
迭代次数
1 m# F, N$ I* r7 C3 L     4
0 ~% d  _& I, }$ \& l5 Q# R没有达到精度要求的个数
    21    21    14     0
各节点的实际电压标么值E为(节点号从小到大排列):
Columns 1 through 7
1 e7 |9 C  d: z1 Y   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
Columns 8 through 12
; W& H) @( d' ]9 k   1.0225 - 0.0954i   1.0044 - 0.0833i   0.9983 - 0.0890i   1.0110 - 0.0197i   1.0211 + 0.0930i
各节点的电压大小V为(节点号从小到大排列):
( \6 @1 |- v9 a: 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
+ Z9 Q: ^9 x" v) S0 E# i各节点的电压相角时shita为(节点号从小到大排列):
! s( u* K% y, L+ ~: m: i6 h* Q" i         0   -5.6958   -1.0221   -0.5827   -0.8925    3.3184    6.6965   -5.3305   -4.7429   -5.0961   -1.1182    5.2027
各节点的功率S为(节点号从小到大排列):
0 s$ [. Z* f9 c* `# uColumns 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
) G) e2 q; {) T( N2 k6 O+ FColumns 8 through 12
-0.4753 - 0.2947i -0.4074 - 0.3608i -0.4074 - 0.3608i -0.4753 - 0.2947i   0.6000 + 0.3720i
" U: y# Q$ n& N2 g各条支路的首端功率Si为(顺序同您输入B1时一样):
   1.6349 + 0.2888i
   0.4827 + 0.0380i
   0.4139 + 0.1626i
   0.4150 + 0.0691i
5 A+ ~( Y+ z$ d9 h$ p-0.9770 + 0.0037i
-0.4753 - 0.2947i
! p* m& w! s9 Y# V: Q-0.4074 - 0.3608i
-0.4074 - 0.3608i
-0.4753 - 0.2947i
-1.4711 - 0.3611i
2 {5 v, S. B  j* f; [-0.5980 - 0.2994i
各条支路的末端功率Sj为(顺序同您输入B1时一样):
- R  z/ H# U% s! n. `-1.6000 - 0.5772i
+ `$ J7 y9 E' i$ Q-0.4768 - 0.3478i
9 W, [" E6 m/ l4 f/ X-0.4091 - 0.4192i
-0.4092 - 0.4198i
3 S5 w% U& M/ @( @! P# n& }! }   0.9943 + 0.0116i
; h! `$ [& y! M  C   0.4768 + 0.3478i
; ~4 F- Y' v+ x/ p+ X, |   0.4091 + 0.4192i
$ B7 X& f' d: s. e4 K1 `7 q, k   0.4092 + 0.4198i
   0.4768 + 0.3495i
   1.5000 + 0.0966i
   0.6000 + 0.3720i
各条支路的功率损耗DS为(顺序同您输入B1时一样):
   0.0349 - 0.2884i
- F- G( d* M& ^3 c5 @) S! `4 Q$ z& Q   0.0060 - 0.3098i
   0.0047 - 0.2566i
   0.0058 - 0.3508i
2 \) d7 y: G( q+ E7 `   0.0173 + 0.0153i
( L) H* L& {7 b   0.0015 + 0.0531i
2 ~3 U1 G( M( d9 }  F8 p   0.0017 + 0.0583i
/ V, |: S" @- ^7 H- @+ a   0.0018 + 0.0590i
   0.0015 + 0.0548i
- C- A/ N& L$ e. c- e8 a# d   0.0289 - 0.2644i
2 I1 b8 h1 }% |* \3 T! t' X8 l; [   0.0020 + 0.0726i

20031374

clc
3 e' Y1 S' f' O! ^, G; uclear
2 @! Q' ?) N* B5 r  vn=input('请输入节点数:n=');
3 j/ ~0 K9 y$ |nl=input('请输入支路数:nl=');
isb=input('请输入平衡母线节点号:isb=');
* f/ N, M% N, B) ~pr=input('请输入误差精度:pr=');
B1=input('请输入由各支路参数形成的矩阵:B1=');
B2=input('请输入各节点参数形成的矩阵:B2=');
Y=zeros(n); e=zeros(1,n);f=zeros(1,n);V=zeros(1,n);
5 I' L2 l' {  oO=zeros(1,n);S1=zeros(nl);
% W5 S3 N4 e2 |: T%求导纳矩阵
for i=1:nl
    if B1(i,6)==0 %NO.6
        p=B1(i,1);q=B1(i,2);
; ^( ]8 O+ Q5 C; [+ V    else p=B1(i,2);q=B1(i,1);
    end
    Y(p,q)=Y(p,q)-1./(B1(i,3)*B1(i,5));%非对角元素
% M' P# u7 Z4 g) D! M    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变比
  G4 a% ^. o8 b" w* o& Q    Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2;
end
: l5 [* _4 [$ i& H: S" fdisp('导纳矩阵Y=');
4 Q" S( c7 T- w1 b* q2 |disp(Y);
G=real(Y);B=imag(Y);
8 N9 {( n1 T  K) D2 ]! P; Rfor i=1:n%n=节点数
# P4 B4 U/ W& `9 X8 K" ^' ?9 L    e(i)=real(B2(i,3));%NO.3节点电压初值
, X6 W7 I  ]" i. m9 j% B    f(i)=imag(B2(i,3));
    V(i)=B2(i,4);%？
4 z1 u) @6 b. G% @8 ^end
for i=1:n
    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？
' p$ l8 [. A& o" l( U, nend
P=real(S);Q=imag(S);

20031374

B2的第4，5列代表什么含义？
请输入各节点参数形成的矩阵: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 w2 H/ y& j9 b4 D5 M" C0        0                          1           0 0 2;
4 t4 B  f. M  o5 j6 T' O0        0                         1           0 0 2;
0.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;
) f, ]7 d+ k& s0       0.66+0.4092i     1       0 0 2;
0 o0 N. |  c; F( y& M0 q1 v0       0.77+0.4774i     1        0 0 2;
8 b1 L) K& h+ n0 C: w- }1+0.62i  0.4+0.248i       1         0 0 2
0 u7 V/ n  K" R) i3 k]

wsz310

谢谢楼主啊,很详细!