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

 找回密码
 立即加入
搜索
查看: 1685|回复: 1

夏道止老师《电力系统分析》N-R迭代过程

[复制链接]

该用户从未签到

尚未签到

发表于 2009-11-14 11:38:18 | 显示全部楼层 |阅读模式
文档文案
文档文件名称: 夏道止老师《电力系统分析》N-R迭代过程
内容简介: 夏道止老师《电力系统分析》N-R迭代过程
所属专业方向: 电力
文件类型: txt
来源: 原创

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

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

×
=========================导纳矩阵输出开始:=========================
# x- m4 q# I2 d9 F  J2.9049 - j11.5030 0.0000 + j5.3182 -1.6600 + j3.1618 -1.2450 + j2.3714 5 O; H% h) u# s4 R
0.0000 + j5.3182 0.0000 - j4.6638 0.0000 + j0.0000 0.0000 + j0.0000
9 T2 i: {9 C7 T-1.6600 + j3.1618 0.0000 + j0.0000 2.4900 - j4.7039 -0.8300 + j1.5809 # L( |+ J! u# `# O
-1.2450 + j2.3714 0.0000 + j0.0000 -0.8300 + j1.5809 2.0750 - j3.9092 , H" }5 L  M; S9 K4 E/ C6 b$ F: t
#########################导纳矩阵输出完毕!#########################, t4 S; {- O- @5 g& l, f9 D
8 m& ]/ J$ o$ c1 `
=========================初始功率向量输出开始:=========================, Y& _& Z, {5 r
0.0000 + j0.0000 -0.5000 + j-0.3000 0.2000 + j0.0000 0.0000 + j0.0000 ( X- ~! e+ D- c
#########################初始功率向量输出完毕!#########################
  p5 O/ q' h% i2 d: @6 Z* V5 t5 ?$ B& r- N# c( q% r7 l
=========================初始电压向量输出开始:=========================
1 w3 i0 O3 s) T0 D% S% Y1.0000 + j0.0000 1.0000 + j0.0000 1.0500 + j0.0000 1.0500 + j0.0000
0 Y9 c% p- C$ }5 B0 b) s' s#########################初始电压向量输出完毕!#########################; E6 [5 A- Y9 V5 V7 m, q' W0 I
9 F/ \) h1 F$ k& P: Q
===================Iteration  k = 0 Begin=========================! g% o! \7 D8 p& R5 y2 i9 G
P[0] = -0.14525# b3 g5 w6 o$ M  M& C# E
Q[0] = 0.37497! X; W& W% M2 e! l: n- Y" u( w# g% C
P[1] = 0.00000
. \" h. e" `( Y) I0 Y; c- ?+ m+ GQ[1] = -0.65434# ?1 o7 Y7 G# a1 t$ V8 g( r6 a
P[2] = 0.08715
7 x0 Q; G: Z7 J0 b8 XQ[2] = 1.10250
4 u) Y5 u5 u* a/ e; ?7 P' ^误差向量 p、q:, B- e; Q( W% T( S1 k& f  o3 v+ N' y
0.14524743, 0.14524743, -0.37496691, -0.50000000, 0.35433884, 0.11285154, 0.00000000,
" ~1 y" H) d3 f% H5 T0 ?- @& \电流 aii、bii:: ?1 R0 G' p" L; G
aii[0] = -0.1452474277   bii[0] = -0.3749669147 aii[1] = 0.0000000000   bii[1] = 0.6543388430 aii[2] = 0.0829985301   bii[2] = 0.0000000000
0 V* S% L: ]+ ~0 `! t( e! {========================雅克比矩阵========================
5 S+ V* v' e" a9 Q& y: m4 _-2.7597000    -11.128100    0.00000000    5.31820000    1.66000000    3.16180000    ! v* W2 k* W) f) y0 E' q
-11.878000    3.05020000    5.31820000    0.00000000    3.16180000    -1.6600000    : {9 G+ u1 Y( N/ G
0.00000000    5.31820000    0.00000000    -5.3182000    0.00000000    0.00000000   
1 F9 \0 J8 m0 _/ j5.31820000    0.00000000    -4.0095000    0.00000000    0.00000000    0.00000000   
& u* j8 b7 ?, [( ?1 Y% f# x' l2 F1.74300000    3.31990000    0.00000000    0.00000000    -2.6975000    -4.9391000      H& E. V/ @( j7 b" W7 L& G
0.00000000    0.00000000    0.00000000    0.00000000    -2.1000000    0.00000000    7 [# t& T" U. r% p
  T  J2 @  s" O  {/ Q; d3 M
所求矩阵的逆矩阵为:  u. e/ ^8 x! y! e3 z& j% J' k
-0.08369     -0.15942     -0.08369     -0.21145     0.00000     -0.30618     9 |3 Q/ b( P, a$ e9 i5 G
-0.23406     0.07113     -0.23406     0.09434     -0.17374     0.14524     1 `  \7 p/ I' T6 G# T. i
-0.11101     -0.21145     -0.11101     -0.52987     0.00000     -0.40612     
  y- J/ x0 z9 ?/ B5 P-0.23406     0.07113     -0.42209     0.09434     -0.17374     0.14524     , U. r9 A/ ?3 y* k6 q
0.00000     -0.00000     0.00000     0.00000     0.00000     -0.47619     5 f, A( O0 E' g/ k7 [" C3 J4 E
-0.18686     -0.00845     -0.18686     -0.01121     -0.31925     0.24965     " z, I3 p2 U; l) a" S* |
方程的解:+ [) `. G- L! F. s- n& d+ {
0.014542  0.070185  -0.069086  0.164202  0.000000  0.029460  
. O( R2 ?% R1 R1 \$ O$ x; U新的电压值(e+jf):
1 D, v; ]+ ~( e, Z: I0.98546 - j0.07019, 1.06909 - j0.16420, 1.05000 - j0.02946, ( _- f9 V# t6 w. O
===================Iteration  k = 0 End  =========================
! s6 Y& ?* v* v( \: O- L! N+ L7 J# }===================Iteration  k = 1 Begin=========================
2 n/ g. m; o  D$ |9 N' H% oP[0] = -0.02835
4 l. f/ d# N' R5 o9 _  u+ MQ[0] = -0.00268
  ^/ H. U$ F" ?9 Q; F" |% \* IP[1] = -0.46151
2 _& S' y% A" U+ B9 gQ[1] = -0.20794& ?: a- `  X, N+ f- F
P[2] = 0.20354
# z3 p2 Q8 y0 p( NQ[2] = 1.10337
8 t7 g) L+ ?  y! L' `误差向量 p、q:
, V7 @- |5 n& i& h4 A0.02834863, 0.02834863, 0.00267722, -0.03848582, -0.09206171, -0.00353918, -0.00086787,
1 V0 K, c: N; g6 x9 {电流 aii、bii:- w7 o9 I; a# M1 j) d2 M, P
aii[0] = -0.0284292626   bii[0] = 0.0047414768 aii[1] = -0.3925561267   bii[1] = 0.2547940506 aii[2] = 0.1642347318   bii[2] = 0.0000000000
! Q0 |; h. m! S$ [  B" i========================雅克比矩阵========================) j/ ~; f% T! J7 G8 m6 Q
-3.6416000    -11.136600    0.37330000    5.24080000    1.85770000    2.99940000    ( _/ H' }1 g0 j% U7 t8 T) b1 u
-11.127100    3.69850000    5.24080000    -0.3733000    2.99940000    -1.8577000   
9 U7 w- _4 H! C; f0.87330000    5.68560000    -0.3733000    -5.2408000    0.00000000    0.00000000   
9 b1 q3 @. F8 F" v9 j5.68560000    -0.8733000    -4.7313000    1.15840000    0.00000000    0.00000000    ' r$ h$ @# E! _, X  E0 {$ J6 H
1.83610000    3.27100000    0.00000000    0.00000000    -2.9173000    -4.8658000   
! u, o6 ~$ S2 U; u/ d+ G# C9 p$ P0.00000000    0.00000000    0.00000000    0.00000000    -2.1000000    0.05890000   
  {% A" A- _& c* [6 |
0 K. ^; Q9 J$ u3 [4 }所求矩阵的逆矩阵为:4 K) ^+ Q$ j* {1 R+ W/ k
-0.11168     -0.15094     -0.13743     -0.16516     -0.01477     -0.29385     
" B7 Z6 P& S- [  T+ b7 D6 N-0.23803     0.07191     -0.22576     0.07869     -0.17259     0.13189     ( h5 @: ]# K* ?
-0.15534     -0.17860     -0.23173     -0.40317     -0.03179     -0.34835     0 O+ g4 ]+ o1 C5 h8 f" P
-0.26578     0.06559     -0.44212     0.08656     -0.18743     0.11893     " e5 }2 R. N( r" ^, P
-0.00558     -0.00024     -0.00562     -0.00026     -0.00903     -0.46893     
+ r  n3 t& c7 y-0.19882     -0.00847     -0.20026     -0.00927     -0.32170     0.25892     * Z  l& u8 ?- @! v5 {3 i
方程的解:- {$ ~0 z7 ?% K" E
0.017232  -0.004615  0.041568  0.002248  0.000520  0.003815  1 ]* O# V7 Y" Z9 n) s
新的电压值(e+jf):
* i4 L! b3 u2 Y/ J0.96823 - j0.06557, 1.02752 - j0.16645, 1.04948 - j0.03328,
) {' I# Y3 C- v% d===================Iteration  k = 1 End  =========================
" F& S0 }, g" K" {===================Iteration  k = 2 Begin=========================
5 I+ ^9 o/ F" }! IP[0] = -0.00050
5 l7 r$ z$ H5 v# |1 y4 I* \. KQ[0] = 0.00005, q- \- y  w  K; o/ z" b# a! |
P[1] = -0.49877
$ R4 h, x0 h; aQ[1] = -0.29567+ V3 F4 x) ?  ~5 G0 |
P[2] = 0.20071
. F3 s- y0 }% D2 k7 c& KQ[2] = 1.10251/ v  Z: g! w; n$ q4 ^2 `
误差向量 p、q:% }' o7 p9 F  f0 c
0.00050286, 0.00050286, -0.00004686, -0.00122615, -0.00432794, -0.00071146, -0.00001483,
) k+ N' u5 I( ?9 _7 b电流 aii、bii:
- r- w% o: k4 P$ c9 V% j+ F2 q/ Yaii[0] = -0.0005202507   bii[0] = -0.0000131695 aii[1] = -0.4275816134   bii[1] = 0.3570182626 aii[2] = 0.1577813702   bii[2] = 0.0000000000 0 Q1 Q! R3 C3 u5 `) X
========================雅克比矩阵========================& R9 j9 ~- @* w( {% a! [
-3.5664000    -10.947100    0.34870000    5.14920000    1.81460000    2.95250000    . F' g) @+ S) |( `3 ~1 d- m1 D3 A
-10.947100    3.56740000    5.14920000    -0.3487000    2.95250000    -1.8146000   
! w* W( Q; }' a" g6 L, z0.88520000    5.46450000    -0.3487000    -5.1492000    0.00000000    0.00000000    / [% G, @. r* `* M" Z& ^! P
5.46450000    -0.8852000    -4.4352000    1.20390000    0.00000000    0.00000000    " F" M1 {. T$ t% o4 F' d5 n
1.84730000    3.26310000    0.00000000    0.00000000    -2.9275000    -4.8538000    . c2 j; }9 S& r3 ?8 c; a+ ?
0.00000000    0.00000000    0.00000000    0.00000000    -2.0990000    0.06660000    * ~& o% w/ p: w
8 X/ k1 H* f9 i% A# g8 z2 y$ i: B
所求矩阵的逆矩阵为:
) Y; |6 T4 D+ n- g2 [-0.11450     -0.15942     -0.14639     -0.18258     -0.01421     -0.30342     
- q1 [" j. l1 ]) T) P-0.23498     0.07099     -0.22078     0.08130     -0.16768     0.13058     
3 |0 r6 x/ x5 N5 @& Y-0.16419     -0.19402     -0.25476     -0.44360     -0.03241     -0.36965     
. B* A5 R4 A) p, d7 l, e% r, ^-0.25794     0.06106     -0.43642     0.08493     -0.17820     0.11145     
/ Q3 ~+ F- z* p  g4 D-0.00627     -0.00040     -0.00635     -0.00046     -0.01009     -0.46835     # k% z8 v8 d/ e2 @4 R; C2 D
-0.19777     -0.01271     -0.20031     -0.01456     -0.31808     0.25478     ; w! c, F! L; y: M$ C) w- \5 [
方程的解:  l+ B1 D2 U$ u' E5 w8 [$ r
0.000934  -0.000085  0.002187  0.000160  0.000021  0.000432  
7 g  i) T; `0 ]6 T" G2 c( V新的电压值(e+jf):
1 ?; i9 X- p$ d1 h' I2 P0.96729 - j0.06549, 1.02533 - j0.16661, 1.04946 - j0.03371,
  R& t3 X" d+ i===================Iteration  k = 2 End  =========================
! M2 s, a' m1 L7 @+ a; x===================Iteration  k = 3 Begin=========================
5 X) k, w$ M5 LP[0] = -0.00000
- ?: G0 z' P4 \3 z9 `0 }# vQ[0] = 0.00000: @2 k: c* i: J2 v1 ?( x
P[1] = -0.50000
# }1 v5 M9 {" A6 e6 b- NQ[1] = -0.29999
# D* \( X5 s8 \2 \2 ?P[2] = 0.20008
/ V7 C. M3 i& p* cQ[2] = 1.10250) q% [* w2 T/ v/ b8 D* {
误差向量 p、q:
7 ?. ~; m2 @" z4 g0.00000048, 0.00000048, -0.00000006, -0.00000179, -0.00001164, -0.00008305, -0.00000019,
. {/ x: }( w  x! u电流 aii、bii:
1 d1 h: ?/ \. p8 f6 Z0 Taii[0] = -0.0000005025   bii[0] = -0.0000000245 aii[1] = -0.4287817842   bii[1] = 0.3622513393 aii[2] = 0.1567496678   bii[2] = 0.0000000000 . D  i; h6 R/ b- T% E' r1 _- R8 K
========================雅克比矩阵========================# Q# w2 w8 y3 k) D. t2 d1 t1 S* h
-3.5632000    -10.936600    0.34830000    5.14420000    1.81270000    2.94970000    * c) z& F% k8 {  O
-10.936600    3.56320000    5.14420000    -0.3483000    2.94970000    -1.8127000   
0 b9 x5 w( `$ k; B6 l- J9 Z+ p0.88610000    5.45290000    -0.3483000    -5.1442000    0.00000000    0.00000000    4 ]$ F5 p. C8 n, e2 v2 X& ]9 Z
5.45290000    -0.8861000    -4.4197000    1.20580000    0.00000000    0.00000000   
7 N$ @0 Z  E+ w2 g$ o, o4 s1.84860000    3.26230000    0.00000000    0.00000000    -2.9284000    -4.8527000   
# ?: M1 l, L* d# x( _* E0.00000000    0.00000000    0.00000000    0.00000000    -2.0989000    0.06740000    + K; ]6 q, O: u- U1 |& a

* P4 y0 P7 ^' M% X0 t5 N所求矩阵的逆矩阵为:
2 i. C( ^+ P, I* a; m5 h2 j3 E( M-0.11472     -0.15989     -0.14693     -0.18356     -0.01423     -0.30393     4 G& W2 ?9 R7 Y% Z
-0.23481     0.07095     -0.22052     0.08145     -0.16742     0.13050     . v8 S0 w$ W, p; A6 C
-0.16472     -0.19489     -0.25605     -0.44590     -0.03248     -0.37084     
- x0 S+ k( A* e* w8 b-0.25751     0.06086     -0.43612     0.08491     -0.17772     0.11109     
! T. T1 a! g. `- S-0.00635     -0.00042     -0.00643     -0.00048     -0.01021     -0.46826     
, ~8 \7 }6 @+ G6 t% k, q) B-0.19773     -0.01296     -0.20034     -0.01488     -0.31788     0.25453     
3 p: Z, x$ A- W! c$ b5 W; @) ?4 Q方程的解:
7 }3 j! y1 y  R* j, }7 Z- w  P- m! o0.000004  0.000013  0.000008  0.000014  0.000001  0.000027  5 h! b% @+ Y# _# _3 d0 T& ~
新的电压值(e+jf):
! B$ |% }7 y4 c& {+ [0.96729 - j0.06550, 1.02532 - j0.16662, 1.04946 - j0.03373,
6 ~: l' v  h# I: B===================Iteration  k = 3 End  =========================& e8 x5 }/ g  F
===================Iteration  k = 4 Begin=========================
% J% L8 s' O: hP[0] = -0.00000( p& N: {8 n4 L" e9 `
Q[0] = -0.00000* G! [3 X; [& F: H6 Y* p, N  o* Y
P[1] = -0.500006 P, ?0 q; _/ ~5 X
Q[1] = -0.30000
' L" G/ k/ v& j  T+ o9 OP[2] = 0.20001
4 n( E9 M: b1 @- C! `( h) y$ kQ[2] = 1.10250
% M) V" l3 t+ v误差向量 p、q:
7 w8 e$ V* m& u# {# i0.00000000, 0.00000000, 0.00000000, 0.00000000, -0.00000000, -0.00000516, -0.00000000, 7 T$ Y, S  x# g4 m" k
电流 aii、bii:9 {2 B' u' ~# n: T% d3 ?  n# ^" s/ Z
aii[0] = -0.0000000000   bii[0] = 0.0000000000 aii[1] = -0.4287787078   bii[1] = 0.3622711617 aii[2] = 0.1566485901   bii[2] = 0.0000000000
/ g+ V3 M# ~/ ?3 L========================雅克比矩阵========================
! M5 P& _( ]/ j1 c-3.5634000    -10.936500    0.34830000    5.14420000    1.81280000    2.94970000      I& b( b9 u2 T" Y
-10.936500    3.56340000    5.14420000    -0.3483000    2.94970000    -1.8128000    , f/ h$ {$ |% y( j
0.88610000    5.45290000    -0.3483000    -5.1442000    0.00000000    0.00000000   
1 g6 l) [( |1 `3 }, C5.45290000    -0.8861000    -4.4197000    1.20590000    0.00000000    0.00000000    3 \5 [8 O& \  D& G# i6 T* \
1.84870000    3.26220000    0.00000000    0.00000000    -2.9284000    -4.8526000   
% u) D5 V- _7 S2 S+ D5 N8 @0.00000000    0.00000000    0.00000000    0.00000000    -2.0989000    0.06750000   
: Z3 [5 {2 w! E  H. H% Z. F" c  ^) T- f! U- T% l
所求矩阵的逆矩阵为:' c8 {' e8 Y# C0 X3 `& X
-0.11473     -0.15989     -0.14693     -0.18357     -0.01423     -0.30393     7 K1 J9 g& `% n
-0.23481     0.07095     -0.22052     0.08145     -0.16742     0.13050     
4 g% Q0 U2 }7 a" r3 a3 p-0.16473     -0.19489     -0.25606     -0.44590     -0.03248     -0.37084     - Y* V6 u; j- i3 D7 Z9 y- w
-0.25751     0.06086     -0.43611     0.08491     -0.17772     0.11108     2 d( h' u5 v: t
-0.00636     -0.00042     -0.00644     -0.00048     -0.01022     -0.46825     . R1 j$ q( A4 ?% a/ c
-0.19773     -0.01297     -0.20034     -0.01489     -0.31788     0.25452     ) u* o7 ]- _* d
方程的解:9 {/ C8 F/ r/ Q, K* U
0.000000  0.000001  0.000000  0.000001  0.000000  0.000002  
# w+ R9 a1 P6 p4 \& x! ?# R3 h$ ~( @新的电压值(e+jf):+ x  }  Y4 B$ I! p
0.96729 - j0.06550, 1.02532 - j0.16663, 1.04946 - j0.03374, 5 V* }. }' @" b1 n3 t! P8 ^
===================Iteration  k = 4 End  =========================2 R- H) d$ x, C7 \2 N
===================Iteration  k = 5 Begin=========================" q- U7 A! m  l0 U0 F6 v
P[0] = 0.00000
2 K4 Z) b, o2 C' E5 x; |3 I, MQ[0] = -0.00000
7 H5 V) I* t7 C- v% e+ kP[1] = -0.500003 E1 T: v3 D  `6 K1 X
Q[1] = -0.30000
# E& I1 r" _, s3 b9 YP[2] = 0.20000
# b/ t2 p: T2 x; @Q[2] = 1.102501 J, ~: a! H( ?6 B
误差向量 p、q:
5 N# ~! q5 H5 K; c-0.00000000, -0.00000000, 0.00000000, 0.00000000, 0.00000000, -0.00000032, -0.00000000,
+ w/ @2 C8 n( ~. U* f5 W=================k=5时,结果收敛=======================
7 d6 m$ y0 P# D8 {0 x================牛顿法计算共耗时:0秒:609毫秒=========

评分

参与人数 1学分 +1 收起 理由
真龍虎王 + 1 实际应用..GOOD

查看全部评分

"真诚赞赏,手留余香"
还没有人打赏,支持一下
楼主热帖
帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】

该用户从未签到

尚未签到

发表于 2009-11-14 12:16:05 | 显示全部楼层
很详细  呵呵
"真诚赞赏,手留余香"
还没有人打赏,支持一下
帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】
您需要登录后才可以回帖 登录 | 立即加入

本版积分规则

招聘斑竹

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

GMT+8, 2025-4-6 07:55

Powered by Discuz! X3.5 Licensed

© 2001-2025 Discuz! Team.

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