|
|
文档文案
| 文档文件名称: |
夏道止老师《电力系统分析》N-R迭代过程 |
| 内容简介: |
夏道止老师《电力系统分析》N-R迭代过程 |
| 所属专业方向: |
电力 |
| 文件类型: |
txt |
| 来源: |
原创 |
马上加入，结交更多好友，共享更多资料，让你轻松玩转电力研学社区！
您需要 登录 才可以下载或查看，没有账号？立即加入
×
=========================导纳矩阵输出开始:=========================
! E) o" U# n' J2.9049 - j11.5030 0.0000 + j5.3182 -1.6600 + j3.1618 -1.2450 + j2.3714 2 {! O3 K2 a' T) z
0.0000 + j5.3182 0.0000 - j4.6638 0.0000 + j0.0000 0.0000 + j0.0000
; k, v! m' m4 C: `6 a7 U d: R$ d1 {5 x& r-1.6600 + j3.1618 0.0000 + j0.0000 2.4900 - j4.7039 -0.8300 + j1.5809
, U2 h$ {, }- U9 n; x-1.2450 + j2.3714 0.0000 + j0.0000 -0.8300 + j1.5809 2.0750 - j3.9092 : G X% c( o0 b
#########################导纳矩阵输出完毕!#########################+ t1 J* `: E+ B2 ~ {
* h$ r; |8 X0 J3 A ]+ N+ x2 ~- d=========================初始功率向量输出开始:=========================
; x+ m8 i& \* m% ]: E" K6 B& r E0.0000 + j0.0000 -0.5000 + j-0.3000 0.2000 + j0.0000 0.0000 + j0.0000
7 j2 h$ Q/ B. D: f# x#########################初始功率向量输出完毕!#########################. R- u* [. \+ _# x
, n I) W) O3 {) z. M$ c=========================初始电压向量输出开始:=========================% _5 I2 O( b5 l, C
1.0000 + j0.0000 1.0000 + j0.0000 1.0500 + j0.0000 1.0500 + j0.0000 ( y/ j# E2 w0 M
#########################初始电压向量输出完毕!#########################
% c1 y/ @ q+ P8 `( o- d( D4 N- b( P) ^7 W( B
===================Iteration k = 0 Begin=========================
F4 b4 G: p6 n) V4 y! |P[0] = -0.145259 S- T( q2 ?5 z" W) Z$ X: `3 v
Q[0] = 0.37497: z5 e$ e2 Y) T" p: M5 ~( P
P[1] = 0.00000
8 ?6 W9 p: g7 x/ \Q[1] = -0.65434
) ?4 n$ D5 s( }8 r5 UP[2] = 0.08715
% u5 {: N3 T; W3 v v& T* tQ[2] = 1.10250& S5 f" y( L* M% Q
误差向量 p、q:
% q) Y8 T. _/ N- C5 _0.14524743, 0.14524743, -0.37496691, -0.50000000, 0.35433884, 0.11285154, 0.00000000,
$ a$ s0 ~3 t$ A8 S电流 aii、bii:& {4 n" c S/ Z0 ~7 R% c
aii[0] = -0.1452474277 bii[0] = -0.3749669147 aii[1] = 0.0000000000 bii[1] = 0.6543388430 aii[2] = 0.0829985301 bii[2] = 0.0000000000 - X9 L, R- Z* O9 Q4 P# ^) z
========================雅克比矩阵========================
+ |) B% T8 N7 ~; @! T-2.7597000 -11.128100 0.00000000 5.31820000 1.66000000 3.16180000 9 @ P& P C5 M* c$ b+ ^) ?3 j, Y0 H
-11.878000 3.05020000 5.31820000 0.00000000 3.16180000 -1.6600000
/ }0 ?! H8 @! V: r) W0.00000000 5.31820000 0.00000000 -5.3182000 0.00000000 0.00000000
0 \2 N) p$ J. B A4 R- M# b5.31820000 0.00000000 -4.0095000 0.00000000 0.00000000 0.00000000 $ S; t6 [1 ?7 E
1.74300000 3.31990000 0.00000000 0.00000000 -2.6975000 -4.9391000
! y: c4 ~- C& [9 M0.00000000 0.00000000 0.00000000 0.00000000 -2.1000000 0.00000000
% M- P8 ^2 x! a
" b! h; [' [& M2 S+ C7 m0 ?所求矩阵的逆矩阵为：
, {! a! F" C' @) G; ^-0.08369 -0.15942 -0.08369 -0.21145 0.00000 -0.30618 ! M, N8 s8 z1 t6 r- @" I4 w( t6 D
-0.23406 0.07113 -0.23406 0.09434 -0.17374 0.14524
6 Y3 w1 w4 Y) p' k. p5 k-0.11101 -0.21145 -0.11101 -0.52987 0.00000 -0.40612
& O6 f+ U/ ?1 @" G; s9 K-0.23406 0.07113 -0.42209 0.09434 -0.17374 0.14524 6 r$ W- @! i6 ^* h E( \5 r
0.00000 -0.00000 0.00000 0.00000 0.00000 -0.47619 % @! a$ e; W; s3 [% |& M6 G5 ^
-0.18686 -0.00845 -0.18686 -0.01121 -0.31925 0.24965 5 x* p. z! N' f- v4 v8 e3 h
方程的解:
7 g9 T& f1 m, d* N- z# B z0.014542 0.070185 -0.069086 0.164202 0.000000 0.029460 * D5 @# N# i6 E- _* k; f! x
新的电压值(e+jf):
; L6 M [& h8 h- i5 i8 q0.98546 - j0.07019, 1.06909 - j0.16420, 1.05000 - j0.02946, 4 o2 n/ c) j! W+ Y8 w
===================Iteration k = 0 End =========================! v" l, X$ n" t0 c: p4 ~1 @
===================Iteration k = 1 Begin=========================) B E0 ]9 a8 ^( o5 w) `
P[0] = -0.02835
7 @2 t( ^* j! K! [' B" XQ[0] = -0.00268+ E& q2 u! \6 g- J9 s
P[1] = -0.46151
; u6 j* i! `: W1 B; tQ[1] = -0.207948 _* @: `: f0 @8 A. v. `$ V7 Y* `! I/ ~
P[2] = 0.20354+ \2 p0 g) r5 K: v
Q[2] = 1.10337
: u: w3 _) L$ Q误差向量 p、q:4 y1 L. ?4 I4 G9 h( |
0.02834863, 0.02834863, 0.00267722, -0.03848582, -0.09206171, -0.00353918, -0.00086787, ' i% }/ {0 I. L1 o
电流 aii、bii:
/ [" n2 {% O- saii[0] = -0.0284292626 bii[0] = 0.0047414768 aii[1] = -0.3925561267 bii[1] = 0.2547940506 aii[2] = 0.1642347318 bii[2] = 0.0000000000
5 b0 E) K* C" o6 a9 u========================雅克比矩阵========================. e% i( @2 I( \, [! y) d4 g) z
-3.6416000 -11.136600 0.37330000 5.24080000 1.85770000 2.99940000
W; o$ t) p; R$ [-11.127100 3.69850000 5.24080000 -0.3733000 2.99940000 -1.8577000 # H( [3 X3 b) O% V5 R$ ?/ Z
0.87330000 5.68560000 -0.3733000 -5.2408000 0.00000000 0.00000000 6 ]* u; I- B# P) l5 C
5.68560000 -0.8733000 -4.7313000 1.15840000 0.00000000 0.00000000 " Y! G3 [9 a/ }( q6 H( V
1.83610000 3.27100000 0.00000000 0.00000000 -2.9173000 -4.8658000 . F: d9 d7 d" B
0.00000000 0.00000000 0.00000000 0.00000000 -2.1000000 0.05890000
! b6 U4 A" ~1 W5 F9 X: [8 H% z0 y
8 u" a4 G5 g4 L; r! @7 t所求矩阵的逆矩阵为：
9 C- t. q8 ^4 v6 L-0.11168 -0.15094 -0.13743 -0.16516 -0.01477 -0.29385
5 g2 R* F7 F9 Z" G L-0.23803 0.07191 -0.22576 0.07869 -0.17259 0.13189 * W: ]- b/ ] d1 {, s6 s
-0.15534 -0.17860 -0.23173 -0.40317 -0.03179 -0.34835
/ e3 A( L* F% E. j( p; [% _4 q-0.26578 0.06559 -0.44212 0.08656 -0.18743 0.11893
8 Z. H4 O0 i1 \$ e-0.00558 -0.00024 -0.00562 -0.00026 -0.00903 -0.46893 9 x5 P5 C! e/ f+ U
-0.19882 -0.00847 -0.20026 -0.00927 -0.32170 0.25892 ; K6 m8 i' k( C1 x* q* u
方程的解:
( b, Z6 Y( e) ~# { b0.017232 -0.004615 0.041568 0.002248 0.000520 0.003815
: [4 ^* j5 A1 g3 _5 s- Y& h新的电压值(e+jf):
6 I5 [" J' V* g2 _( a# [/ Y, F0.96823 - j0.06557, 1.02752 - j0.16645, 1.04948 - j0.03328, , q" u4 Q7 n" E# w" V( v1 i
===================Iteration k = 1 End =========================
1 f+ j# |( T: Y1 A7 u% C& V0 W===================Iteration k = 2 Begin=========================$ `9 @, s& t% i r7 M4 k3 _
P[0] = -0.00050
( E8 s0 F7 h4 d& T8 E% d* OQ[0] = 0.00005* j3 n+ U) J- t+ ~
P[1] = -0.49877
$ K1 K+ I0 h4 o! Z7 r4 xQ[1] = -0.29567
% S ~) v/ M9 n4 ]0 FP[2] = 0.20071
: E3 f* |" o2 q1 S& a: S( ?Q[2] = 1.10251" E4 X7 [% D- ` ` D/ }* [
误差向量 p、q:( [; Y& `8 Z- V4 T
0.00050286, 0.00050286, -0.00004686, -0.00122615, -0.00432794, -0.00071146, -0.00001483,
0 `5 |! H( e2 X6 Q: h电流 aii、bii:5 \- `7 V* {" e6 G1 D
aii[0] = -0.0005202507 bii[0] = -0.0000131695 aii[1] = -0.4275816134 bii[1] = 0.3570182626 aii[2] = 0.1577813702 bii[2] = 0.0000000000
* V, a. H+ D* f) h4 M========================雅克比矩阵========================$ c7 E( K+ J# w |: t7 x! x
-3.5664000 -10.947100 0.34870000 5.14920000 1.81460000 2.95250000 ( X+ ?! E, p Z
-10.947100 3.56740000 5.14920000 -0.3487000 2.95250000 -1.8146000
0 b- ?0 y; _9 B0.88520000 5.46450000 -0.3487000 -5.1492000 0.00000000 0.00000000 0 ?: }. w: T( Q" ^ ?$ u1 z0 |4 Q
5.46450000 -0.8852000 -4.4352000 1.20390000 0.00000000 0.00000000 % g. H6 m. [7 S: }
1.84730000 3.26310000 0.00000000 0.00000000 -2.9275000 -4.8538000 # b D' h4 B8 ^2 e
0.00000000 0.00000000 0.00000000 0.00000000 -2.0990000 0.06660000
8 s; Q- r. p4 W+ D# k0 b; ?# B' z
- t' x$ {5 v/ p* k0 G* l: [所求矩阵的逆矩阵为：& f! |- O D7 k
-0.11450 -0.15942 -0.14639 -0.18258 -0.01421 -0.30342
! ^2 m* t1 \3 f- l/ R-0.23498 0.07099 -0.22078 0.08130 -0.16768 0.13058 - V, _( P/ k! L/ ^2 K$ x8 Z5 C5 C
-0.16419 -0.19402 -0.25476 -0.44360 -0.03241 -0.36965 4 E G2 [1 j6 J; Y1 _5 q1 s
-0.25794 0.06106 -0.43642 0.08493 -0.17820 0.11145 3 ]) H+ r/ @9 y) M7 ^
-0.00627 -0.00040 -0.00635 -0.00046 -0.01009 -0.46835
/ I, V: m& g$ ]/ ^: D5 S-0.19777 -0.01271 -0.20031 -0.01456 -0.31808 0.25478 ) ^, D) Q( n' Y
方程的解:& H2 Q0 P. ~- M: R8 u: u7 E
0.000934 -0.000085 0.002187 0.000160 0.000021 0.000432 ! P7 o5 @' }6 |' G6 V5 ]& `/ h
新的电压值(e+jf):
& B: W ^2 j( [( b7 J' l/ i2 }5 ]0.96729 - j0.06549, 1.02533 - j0.16661, 1.04946 - j0.03371,
& f7 d3 G4 W! y- u9 ]4 {. N" j===================Iteration k = 2 End =========================; \ S6 F" L6 o2 \* E) @
===================Iteration k = 3 Begin=========================
% x3 x4 z" F) s; I% nP[0] = -0.00000
. a6 F) r- l# h. e2 e) Q @Q[0] = 0.00000
' J! H: x& O; F( H- B- k5 ]; fP[1] = -0.50000
$ D" E' Y3 D% f6 eQ[1] = -0.29999
) k0 S' T& D- n& E4 [P[2] = 0.20008
% Q6 V+ W; O6 B2 pQ[2] = 1.10250
& ?; E- j0 d3 x% t& F) c# Y5 J误差向量 p、q:, t2 F; r9 j5 q* A& T% a K' b
0.00000048, 0.00000048, -0.00000006, -0.00000179, -0.00001164, -0.00008305, -0.00000019,
+ s; M$ a4 ^0 H- M, }电流 aii、bii:$ T8 \& {0 [0 ?' `' A" R
aii[0] = -0.0000005025 bii[0] = -0.0000000245 aii[1] = -0.4287817842 bii[1] = 0.3622513393 aii[2] = 0.1567496678 bii[2] = 0.0000000000 - l h. @; [1 J0 C; ~, |" }
========================雅克比矩阵========================
( [5 a5 }1 ]) V-3.5632000 -10.936600 0.34830000 5.14420000 1.81270000 2.94970000
1 D3 h4 `, a& b) E-10.936600 3.56320000 5.14420000 -0.3483000 2.94970000 -1.8127000 : b o @: n2 F8 l' _
0.88610000 5.45290000 -0.3483000 -5.1442000 0.00000000 0.00000000 2 d8 R) ~0 [- w+ G8 R# F
5.45290000 -0.8861000 -4.4197000 1.20580000 0.00000000 0.00000000
6 W5 m( F: r6 k8 j4 ?4 @1.84860000 3.26230000 0.00000000 0.00000000 -2.9284000 -4.8527000 3 y) g( b, A( n, Z+ `+ Y
0.00000000 0.00000000 0.00000000 0.00000000 -2.0989000 0.06740000 ! n3 [/ @" \( }6 w4 T
; d6 C: v5 B' S3 S4 j/ e& I# q, F; |
所求矩阵的逆矩阵为：
$ `2 o; l( x2 ~9 H3 h* ]-0.11472 -0.15989 -0.14693 -0.18356 -0.01423 -0.30393 9 h2 p0 \4 u( E; \1 x! x
-0.23481 0.07095 -0.22052 0.08145 -0.16742 0.13050
! `, @' V, I' G-0.16472 -0.19489 -0.25605 -0.44590 -0.03248 -0.37084 ; k' n7 H, w0 |& d0 f: z$ q5 h) w
-0.25751 0.06086 -0.43612 0.08491 -0.17772 0.11109 - [% M- H B# z7 s
-0.00635 -0.00042 -0.00643 -0.00048 -0.01021 -0.46826
, S/ H8 U* n7 u-0.19773 -0.01296 -0.20034 -0.01488 -0.31788 0.25453 % W7 z3 _) P- B% A8 i/ ^, v
方程的解:+ `. ~! J. |8 y9 e( S; e, k, _* D! z
0.000004 0.000013 0.000008 0.000014 0.000001 0.000027 # {" c. S; N, t$ e; R" _8 k/ ^# i
新的电压值(e+jf):
- z; r6 M" }4 k, z7 ]; Q0.96729 - j0.06550, 1.02532 - j0.16662, 1.04946 - j0.03373,
3 z( L4 }' s( o' K===================Iteration k = 3 End =========================# Y' b0 c9 L5 C$ E/ e. I8 D1 J" N1 b
===================Iteration k = 4 Begin=========================% m3 S+ N# |" ]8 E
P[0] = -0.00000, E2 k8 ?% {1 R' C0 x! f' O
Q[0] = -0.00000
% l' r+ `! x+ V @P[1] = -0.500006 S( P4 _. O7 |6 p( n2 }" s5 I
Q[1] = -0.300006 ?, j; w4 ?, \( L: x% I5 \
P[2] = 0.20001
8 P! r6 `0 X+ J) k) o" qQ[2] = 1.10250
7 n9 |6 E9 w5 j9 z, [' ~误差向量 p、q:
% \5 q" L. g& ?* d5 r+ w4 b" D0.00000000, 0.00000000, 0.00000000, 0.00000000, -0.00000000, -0.00000516, -0.00000000, ; C- ?' _9 n' L) l1 E0 Q
电流 aii、bii:
1 ^( g3 r- W4 @0 q7 u l* Q) Q- vaii[0] = -0.0000000000 bii[0] = 0.0000000000 aii[1] = -0.4287787078 bii[1] = 0.3622711617 aii[2] = 0.1566485901 bii[2] = 0.0000000000 4 `9 S# \0 E- _) k
========================雅克比矩阵========================
9 V! ^& O3 v8 u5 a-3.5634000 -10.936500 0.34830000 5.14420000 1.81280000 2.94970000
D% l: o M! h& I" d-10.936500 3.56340000 5.14420000 -0.3483000 2.94970000 -1.8128000
5 ` `5 P* a" D( c* x' b) i% O0.88610000 5.45290000 -0.3483000 -5.1442000 0.00000000 0.00000000
; i# S! ~. ]$ c5 \$ n5.45290000 -0.8861000 -4.4197000 1.20590000 0.00000000 0.00000000 ; A! S* d2 w T; D$ J# _+ ~
1.84870000 3.26220000 0.00000000 0.00000000 -2.9284000 -4.8526000
" l0 t$ q _2 D& d0.00000000 0.00000000 0.00000000 0.00000000 -2.0989000 0.06750000
3 z9 v9 l% Y$ u0 _- H% I' F3 C* \# K7 P0 k; b) _, t' k u$ |) P1 g
所求矩阵的逆矩阵为：
4 O, A! T4 B( I0 R-0.11473 -0.15989 -0.14693 -0.18357 -0.01423 -0.30393
; }8 A' \, s5 Z- M( g-0.23481 0.07095 -0.22052 0.08145 -0.16742 0.13050
9 p3 i+ @! G- _3 @$ ^" L* W-0.16473 -0.19489 -0.25606 -0.44590 -0.03248 -0.37084
1 I; N9 {4 o3 r3 Y% x, T9 A, |-0.25751 0.06086 -0.43611 0.08491 -0.17772 0.11108
( v( I7 R5 l8 `( b: E7 \- D$ t-0.00636 -0.00042 -0.00644 -0.00048 -0.01022 -0.46825 / ]( n2 h& {0 Q: o N! E
-0.19773 -0.01297 -0.20034 -0.01489 -0.31788 0.25452
* |. T: A# A/ d方程的解:" C! _* K: P# Z' d
0.000000 0.000001 0.000000 0.000001 0.000000 0.000002 5 D! d. F/ X$ F/ H4 S) @/ h' Z
新的电压值(e+jf):
, v7 m) a3 c1 r$ R3 ` K0.96729 - j0.06550, 1.02532 - j0.16663, 1.04946 - j0.03374,
# {6 D) y/ L% f# e===================Iteration k = 4 End =========================
/ v& _3 o3 {$ H. l: Q8 N# s: h===================Iteration k = 5 Begin=========================
$ f4 h/ k. ~0 q9 v! D* u* Z2 TP[0] = 0.000007 u% n! a& e, x* {
Q[0] = -0.000002 V& }2 z0 G' v' ]3 t) c/ S% S2 i
P[1] = -0.50000
" F3 I! x1 Y3 v0 [& [; CQ[1] = -0.30000
9 \) p, @0 z9 g1 x; [P[2] = 0.20000
: M# F7 Q+ [2 o& H" T; Q7 ?- \Q[2] = 1.10250( i: g: I6 E& E" h1 | M* n
误差向量 p、q:0 A# K) c' U D8 u1 ]5 K4 f; w, G
-0.00000000, -0.00000000, 0.00000000, 0.00000000, 0.00000000, -0.00000032, -0.00000000,
, b' a+ J5 A }' V( ?" |) e=================k=5时，结果收敛=======================
+ N6 @/ N7 L( ^ F7 K8 F& p3 R% g================牛顿法计算共耗时:0秒:609毫秒========= |
评分
-
查看全部评分
|
赣公网安备36040302000210号 )|网站地图
狗仔卡
发表于 2009-11-14 11:38:18
提升卡
置顶卡
沉默卡
喧嚣卡
变色卡
显身卡