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

 找回密码
 立即加入
搜索
查看: 3576|回复: 8

[讨论] IEEE14节点潮流程序

[复制链接]

该用户从未签到

尚未签到

发表于 2010-10-28 11:34:36 | 显示全部楼层 |阅读模式

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

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

×
clear
: ^% |6 H+ [8 s/ `/ G1 Mr(14,10)=0.01938;x(14,10)=0.05917;r(10,11)=0.04699;x(10,11)=0.19797;r(10,1)=0.05811;x(10,1)=0.17632;r(14,2)=0.05403;x(14,2)=0.22304;r(10,2)=0.05695;x(10,2)=0.17388;
7 A- Y/ B: o- R8 zr(10,14)=0.01938;x(10,14)=0.05917;r(11,10)=0.04699;x(11,10)=0.19797;r(1,10)=0.05811;x(1,10)=0.17632;r(2,14)=0.05403;x(2,14)=0.22304;r(2,10)=0.05695;x(2,10)=0.17388;: [- j$ x! b/ [5 s
r(11,1)=0.06701;x(11,1)=0.17103;r(1,2)=0.01335;x(1,2)=0.04211;r(2,12)=0.00000;x(2,12)=0.23488;r(1,3)=0.00000;x(1,3)=0.20452;r(3,13)=0.00000;x(3,13)=0.17615;1 j2 n6 h0 N, s# f  b
r(1,11)=0.06701;x(1,11)=0.17103;r(2,1)=0.01335;x(2,1)=0.04211;r(12,2)=0.00000;x(12,2)=0.23488;r(3,1)=0.00000;x(3,1)=0.20452;r(13,3)=0.00000;x(13,3)=0.17615;
6 f$ w% d5 z5 k- z+ y+ A3 mr(1,4)=0.00000;x(1,4)=0.53894;r(3,4)=0.00000;x(3,4)=0.11001;r(4,5)=0.03181;x(4,5)=0.08450;r(12,6)=0.09498;x(12,6)=0.19890;r(12,7)=0.12291;x(12,7)=0.25581;. n% }/ l/ i: [6 Z  H* K
r(4,1)=0.00000;x(4,1)=0.53894;r(4,3)=0.00000;x(4,3)=0.11001;r(5,4)=0.03181;x(5,4)=0.08450;r(6,12)=0.09498;x(6,12)=0.19890;r(7,12)=0.12291;x(7,12)=0.25581;" H( ^9 R; ]: O0 x* X
r(12,8)=0.06615;x(12,8)=0.13027;r(4,9)=0.12711;x(4,9)=0.27038;r(5,6)=0.08205;x(5,6)=0.19207;r(7,8)=0.22092;x(7,8)=0.19988;r(8,9)=0.17093;x(8,9)=0.34802;
1 W' V5 V+ \* br(8,12)=0.06615;x(8,12)=0.13027;r(9,4)=0.12711;x(9,4)=0.27038;r(6,5)=0.08205;x(6,5)=0.19207;r(8,7)=0.22092;x(8,7)=0.19988;r(9,8)=0.17093;x(9,8)=0.34802;8 o: T2 q1 }, n( n; G" \) S0 P! I
/ g5 T! Y8 G; s* w: D* i
y(14,10)=0.02640;y(10,11)=0.02190;y(10,1)=0.01870;y(14,2)=0.02460;y(10,2)=0.01700;y(11,1)=0.01730;y(1,2)=0.00640;y(2,12)=-0.31063;y(1,3)=-0.10999;y(1,4)=-0.05936;* `# |- O4 {) A  t0 }0 ]
y(10,14)=0.02640;y(11,10)=0.02190;y(1,10)=0.01870;y(2,14)=0.02460;y(2,10)=0.01700;y(1,11)=0.01730;y(2,1)=0.00640;y(12,2)=0.28951;y(3,1)=0.10757;y(4,1)=0.05752;! V( Q9 w5 r/ O# }$ s, }

5 s/ Y) |7 t" @for m=1:14
2 }$ g7 D( ~5 F2 _* C; q    t(m)=0;
' R4 Z4 t' ]+ m    for n=1:14
' z  J. c8 c5 |/ ^5 T4 a1 v7 {        t(m)=t(m)+y(m,n);6 U$ s( d' |$ g! U) h; w
    end
* v1 o8 R' V1 M& {6 e$ g    if m==40 y* x( q) ^! X9 w
        t(m)=t(m)+0.190;
  z* J; A$ O8 A* m! z; ]$ F    end
: o/ Q8 L- d, [end
) h1 ]2 _' k* x" J- _for m=1:142 q) C- h) ~: a. v4 f6 ~
    d(m)=0;( `) S6 h" X& \+ j- s5 |# Q9 q5 D7 Y6 M
    for n=1:14, ?* @. Y- s" m( z6 z: u  F
        if m==n  z5 o  \* I+ C2 F- v% w  s
            
+ O: i$ F% m; I$ e1 k        elseif (r(m,n)==0)&(x(m,n)==0)
4 r" l- t; ^  \5 K4 C, w$ ], R           d(m)=d(m)+0;
) |% b  E$ r7 W1 G        else
; Z$ A" f. K/ j. d7 w           d(m)=d(m)+1/(r(m,n)+j*x(m,n));   
& X! b+ L) r2 u) \8 ]# r        end- A) U5 D( G0 ]2 B, i9 z8 ?
    end
9 E4 W4 k) x6 c, o! l1 A" C7 u/ f; ?end
8 K( Z7 R* L$ X4 X* lfor m=1:14
* h. q, a' u: D$ j    for n=1:14
  f* e, G* j8 s- Y$ F3 |$ p& K0 d        if m==n
4 `% c( Y# i: |8 O           Y(m,n)=j*t(m)+d(m);$ ~; {% C2 o/ M8 d2 f2 M+ c* ~  r
        elseif (r(m,n)==0)&(x(m,n)==0)
5 L- X0 |% c8 Q1 Z) B: E( p$ M           Y(m,n)=0;
9 C5 E* h7 P8 `1 o7 E  l. u        else, e- }5 X: ]9 @; D$ j4 }4 h
           Y(m,n)=-1/(r(m,n)+j*x(m,n));
- }4 f7 k( o: g5 k; F  S        end6 P& B  g0 m) W9 f; r& w. T
    end, c" n# e2 w6 I
end
6 n9 y1 w7 i2 w2 f9 V5 p3 {8 h( wY;
! R# S5 n$ S5 K. f: m0 ~$ z* K: vG=real(Y);B=imag(Y);
: a% o7 [! i, _%给定节点电压初值和循环次数k=06 d7 \5 S- h8 Y$ z
delt(10)=0;delt(11)=0;u(1)=1.0;delt(1)=0;u(2)=1.0;delt(2)=0;delt(12)=0;u(3)=1.0;delt(3)=0;delt(13)=0;u(4)=1.0;delt(4)=0;u(5)=1.0;delt(5)=0;u(6)=1.0;delt(6)=0;u(7)=1.0;delt(7)=0;u(8)=1.0;delt(8)=0;u(9)=1.0;delt(9)=0;7 B' X% G2 v/ x6 w! S
p(10)=0.183;p(11)=-0.942;p(1)=-0.478;q(1)=0.039;p(2)=-0.076;q(2)=-0.016;p(12)=-0.112;p(3)=0;q(3)=0;p(13)=0;p(4)=-0.295;p(5)=-0.09;q(5)=-0.058;p(6)=-0.035;q(6)=-0.018;p(7)=-0.061;q(7)=-0.016;p(8)=-0.135;q(8)=-0.058;p(9)=-0.149;q(9)=-0.05;9 L' }9 T. u% u4 F
k=0;precision=1;N1=13;, Q; X) u2 l% M* y+ l. R
%求节点功率的不平衡量�6�2P和�6�2Q
. T, u' z8 V. k3 B4 c) p% Nwhile precision>0.00001: N0 {- c: S: ]  K& u2 z7 t( I+ j
u(10)=1.0450;u(11)=1.01;u(12)=1.07;u(13)=1.09;u(14)=1.06;delt(14)=0;q(4)=-0.166+(u(4))^2*0.190;7 J+ w: Z# n9 Z: k2 t0 G& M; A8 ]
     for m=1:N1     
7 }& b+ [3 b( Z5 _$ k+ U! X         if m<=9
" L: F$ o* A0 ?2 @             for n=1:N1+1* D) \1 o2 B" t. `" I% _
pt(n)=u(m)*u(n)*(G(m,n)*cos(delt(m)-delt(n))+B(m,n)*sin(delt(m)-delt(n)));
; n/ R' g* b- Qqt(n)=u(m)*u(n)*(G(m,n)*sin(delt(m)-delt(n))-B(m,n)*cos(delt(m)-delt(n)));     - Z1 J$ s, Y  b/ O; n
             end
" B7 y5 i/ Z; ~2 a% g             pp(m)=p(m)-sum(pt);qq(m)=q(m)-sum(qt);
6 V. p5 x  z* u8 a& w9 Q% B         else
( b, `5 S8 e9 l4 Y4 U4 N+ ~             for n=1:N1+1: C5 x* `/ c$ P: F& j8 c
f(n)=u(m)*u(n)*(G(m,n)*cos(delt(m)-delt(n))+B(m,n)*sin(delt(m)-delt(n)));
+ k% o3 f5 \3 K9 }3 L5 q" ~& Z             end2 h* ~6 z: E3 M" M; {8 W/ Z
             pp(m)=p(m)-sum(f);qq(m)=0;+ V0 Q$ R' {* x, C: b2 V
         end      
. w) C2 w- M1 h     end
6 P- G  H7 Z5 ~: ^, t- C# ^%计算雅克比矩阵各元素
: {! B3 S! T& l  B: Q  Q, Z! {1 }for m=1:N1
/ K! \) X( E; [3 ]; o3 u    for n=1:N1+1& F4 m+ U) u  h6 S. V
         h0(n)=u(m)*u(n)*(G(m,n)*sin(delt(m)-delt(n))-B(m,n)*cos(delt(m)-delt(n)));( j% E7 V0 z* w" V
         n0(n)=-u(m)*u(n)*(G(m,n)*cos(delt(m)-delt(n))+B(m,n)*sin(delt(m)-delt(n)));
. K  k- ]& h% A* w  [6 O% p( C  I         j0(n)=-u(m)*u(n)*(G(m,n)*cos(delt(m)-delt(n))+B(m,n)*sin(delt(m)-delt(n)));: x, j$ A, _. O9 U
         L0(n)=-u(m)*u(n)*(G(m,n)*sin(delt(m)-delt(n))-B(m,n)*cos(delt(m)-delt(n)));
) R6 Z( {; ~4 v' p    end9 q  C( ]" m: q/ s& g. m
    if m<=9: a% S( K& O3 n# t
H(m,m)=sum(h0)- u(m)*u(m)*(G(m,m)*sin(delt(m)-delt(m))-B(m,m)*cos(delt(m)-delt(m)));
0 T/ A1 R, N6 C( bN(m,m)=sum(n0)-2*u(m)^2*G(m,m)+u(m)*u(m)*(G(m,m)*cos(delt(m)-delt(m))+B(m,m)*sin(delt(m)-delt(m)));; F+ g0 G9 `2 M: r: I0 M" l
J(m,m)=sum(j0)+u(m)*u(m)*(G(m,m)*cos(delt(m)-delt(m))+B(m,m)*sin(delt(m)-delt(m)));
9 W+ ?' Y" X& n2 uL(m,m)=sum(L0)+2*u(m)^2*B(m,m)+u(m)*u(m)*(G(m,m)*sin(delt(m)-delt(m))-B(m,m)*cos(delt(m)-delt(m)));5 R1 X) ^7 `( T/ J, ^! }
    else8 c9 n; L$ M3 {8 _* B6 F
H(m,m)=sum(h0)- u(m)*u(m)*(G(m,m)*sin(delt(m)-delt(m))-B(m,m)*cos(delt(m)-delt(m)));
* W' w$ h# i: V, W6 x: Y. GN(m,m)=0;
& y. M9 C4 g5 y3 h, a' o. N% [# tJ(m,m)=0;
/ h) Q" ^: x$ O( p2 `L(m,m)=0;9 h2 U/ v6 P9 m) |6 p
    end7 I# l" Q, n* i3 U4 V/ y( l
end
# u9 [0 _1 ^$ v+ a0 }for m=1:N10 @, n2 m- v5 L% }" C
    JJ(2*m-1,2*m-1)=H(m,m);JJ(2*m-1,2*m)=N(m,m);& V1 T) R; }8 d
    JJ(2*m,2*m-1)=J(m,m);JJ(2*m,2*m)=L(m,m);8 q0 C# l9 p' m. B! M
end
  b2 h( z8 A$ s# k ' Z2 j$ W& C2 O+ b6 m
for m=1:N1  y; s+ }, A8 g& V4 l
    for n=1:N1   ; D! j; Q0 m2 J6 Q
        if     (m<=9)&(n<=9)&(m~=n)                       6 s% ]& \" p$ q+ X% H' m
H(m,n)= -u(m)*u(n)*(G(m,n)*sin(delt(m)-delt(n))-B(m,n)*cos(delt(m)-delt(n)));
' m  Y$ o3 W) _J(m,n)=u(m)*u(n)*(G(m,n)*cos(delt(m)-delt(n))+B(m,n)*sin(delt(m)-delt(n)));
& g, P' J3 Q. ]$ ^* z! IN(m,n)=-J(m,n);L(m,n)=H(m,n);2 k& V+ z! B! b
JJ(2*m-1,2*n-1)=H(m,n);JJ(2*m-1,2*n)=N(m,n);
* c, E( B$ j; v5 y  Z# t7 GJJ(2*m,2*n-1)=J(m,n);JJ(2*m,2*n)=L(m,n);. @, h, ?  A6 ?: l3 T$ u4 @
        elseif ((m==10)&(n<10))|((m==11)&(n<10))|((m==12)&(n<10))|((m==13)&(n<10))
# W# K( Z  l! d! c! zH(m,n)= -u(m)*u(n)*(G(m,n)*sin(delt(m)-delt(n))-B(m,n)*cos(delt(m)-delt(n)));
2 {' c5 I& A- A; E+ yJ(m,n)=u(m)*u(n)*(G(m,n)*cos(delt(m)-delt(n))+B(m,n)*sin(delt(m)-delt(n)));
( o% d6 T# s% {$ Q; [/ P7 ^: AN(m,n)=-J(m,n);
% m  m; ~/ L! y1 q% O6 \9 P$ @JJ(2*m-1,2*n-1)=H(m,n);JJ(2*m-1,2*n)=N(m,n);& M, E+ D$ S5 K0 `
JJ(2*m,2*n-1)=0;JJ(2*m,2*n)=0;
0 h+ b  x- ]! [0 k5 q9 _        elseif ((n==10)&(m<10))|((n==11)&(m<10))|((n==12)&(m<10))|((n==13)&(m<10))
" v$ h0 H4 H/ U/ l# YH(m,n)= -u(m)*u(n)*(G(m,n)*sin(delt(m)-delt(n))-B(m,n)*cos(delt(m)-delt(n)));
9 c/ i0 y0 z! P! s  \2 I8 J. EJ(m,n)=u(m)*u(n)*(G(m,n)*cos(delt(m)-delt(n))+B(m,n)*sin(delt(m)-delt(n)));. m9 ^) M% n! s8 X; I$ @1 B0 J
JJ(2*m-1,2*n-1)=H(m,n);JJ(2*m-1,2*n)=0;
4 y7 G' Z. h1 W/ x* A+ l: yJJ(2*m,2*n-1)=J(m,n);JJ(2*m,2*n)=0;' r6 U' ?( R, d7 t
        elseif (m>=10)&(n>=10)&(m~=n), V4 [$ i& J3 z# a
H(m,n)= -u(m)*u(n)*(G(m,n)*sin(delt(m)-delt(n))-B(m,n)*cos(delt(m)-delt(n)));2 ?0 f$ V# ], |
JJ(2*m-1,2*n-1)=H(m,n);JJ(2*m-1,2*n)=0;6 Y0 B7 p4 X# y2 J) h7 z
JJ(2*m,2*n-1)=0;JJ(2*m,2*n)=0;" Q" v/ q2 r1 Z9 R+ e4 Z
        end
* a! B/ s8 I3 n7 h6 ~/ o& n: A    end. z# A! I  U! d8 B+ R0 }
end- ~' N! D" Y' L
% 对雅克比矩阵进行修正,使之变成非奇异阵,以便后续求逆计算
! Z6 y% }5 ]! X/ a" o2 db=0;4 J2 G8 C8 _7 j7 j5 h
for m=1:22
. @) K" Z! X* ~    for n=1:222 ~" ?# r+ i) {0 U' A) ?
        if (m<=18)&(n<=18)
- O4 U. Y0 T, b. D8 u. {        A(m,n)=JJ(m,n);
3 w7 s- O4 T4 ?/ E! N' p        elseif (m>18)&(n<=18). N" a: G4 m# F
        A(m,n)=JJ(m+b,n);! [+ M9 S: d/ U, r' k" `' {
        end   
+ Y7 U: u/ F# [6 {/ T) s    end! F) G6 [- _! p1 z* J2 u  v' X
    if m>18
! ?3 }2 g+ a5 H$ R        b=b+1;: a6 e5 ?$ Y" c7 n! R
    end  s! [# |1 q! X  Z' B: }
end
5 m( j+ {: o4 k  N
  U9 O3 k" g4 I8 [b=0;
& d& c. m* s' g6 O) c( z8 w. Bfor n=1:22
! w$ [4 T. k- A( ?* Q' s% n* Q& ~: b7 {    for m=1:22, o2 o9 e2 I9 M- W2 i6 H
        if n<=18
9 v) E0 k$ b/ F' u6 N        A(m,n)=A(m,n);' V4 [$ t5 n' F
        elseif (m<=18)&(n>18)
2 w1 z' n9 L! }) n$ d- w3 ?        A(m,n)=JJ(m,n+b);
9 q4 [6 U5 H( Z8 q/ x        end6 h3 Y) o9 I, x3 @# D
    end! c% X* k( j9 j
    if n>18
. i( ], K. P9 ^1 C! I/ N. c& A" D        b=b+1;& @8 I4 T- f, k) _3 ?
    end
" I9 }; O' L, Y6 O8 hend9 l, k- \  ]4 F" e
5 ~! f4 O% F# J6 X
b=0;
- \4 Y) v, q  \/ Y. F0 C% z7 `for m=1:22
# \3 q+ @& E( Q4 y0 f    for n=1:26
8 u2 b! N' E5 e! [        if (m>=19)&(n>=19)
# g, B  X8 W. l7 t. p/ b        D(m,n)=JJ(m+b,n);/ ~- x, z; y& T) n) u. M- Z
        end
$ a, l- N- o$ k) B    end
' E2 ^8 x& ?; w! X    if m>18
6 @& N2 \- ?& Z        b=b+1;
/ i  E9 o, u" j0 Q) l    end3 }6 s8 q6 u5 z8 E8 X
end
" {) j" j" J4 g. A! \1 T& a+ i! ~ 5 P6 T& X5 U' z
b=0;
! I  \8 I# |# Vfor n=1:22
( p0 k1 g, z# @- A# Y  r8 W1 ^    for m=1:22
" b9 x3 Q  U) W" S        if (m>=19)&(n>=19)
" J. p3 c# G8 k8 y+ @3 ^4 n! F2 n        D(m,n)=D(m,n+b);
$ V& T' c  f4 A8 O" ^' Z7 T; E        end; x- _# O1 r9 z3 c  F
    end
% I' A* G% `8 S" B: r' A/ L/ u    if n>182 c# u& _! Q* Y
        b=b+1;
& X! j$ [- A% i$ ?) {1 w9 P. {    end
7 w' v$ S4 y: mend
' t" v7 r3 R; B- K; b! E+ h
( r7 _8 c# {, f* Y# Yfor m=1:22  m" R4 d- M; l+ }
    for n=1:22+ p  s# K3 R& u
        if (m>=19)&(n>=19) ! v$ z* @1 f: V0 x! u" f+ B! m8 Q( Z
        A(m,n)=D(m,n);/ w2 }/ d! `! |( B6 Q, d" I
        else$ Q- Q0 E' E" A0 H2 {4 f5 I" q1 Y
        A(m,n)=A(m,n);9 o2 @( L! b4 U! \# V
        end
  V- L; l4 B* l8 o% @    end
5 x3 Y7 h3 p$ g) A& I( R/ w4 C  lend, ]! @# K7 h& ~
% 形成不平衡量的列矩阵$ |. \2 K6 S7 `1 V, N5 z' L6 B  q
for m=1:N1
+ Q% E) g  C* B# O3 Y, l    PP(2*m-1)=pp(m);PP(2*m)=qq(m);
% L0 w/ {5 X! U/ C0 A5 t; Eend
! }' o* T: Y9 a: I# d" T. XC=PP(1:1:18);5 |8 o$ G9 C6 N9 o3 Y/ D" _9 J
b=0;
& G1 ?7 Y: L4 c& k1 }6 c+ Q1 h1 O  xfor m=1:22# G# E& E& f6 A0 T$ D- W
    if m<=18
9 ]" v  n$ e$ E% s    C=C;
2 b1 P) G3 ?( \    else+ q( G  I! L5 Z6 i- C8 y7 X
    C(m)=PP(m+b);
4 Y+ i" j+ s: M    b=b+1;8 i/ \9 z% H9 n0 d! l
    end! t; D# I, ~4 X' F! U( B0 }! Z
end; B, q$ O7 b) v" j  L
% 解修正方程式,得到修正量
9 ^% {. S! N6 T' B* Xuu=-inv(A)*C';, U$ w* K. y! S: e
precision=max(abs(uu));
# V; p* T9 I( H$ A2 W% 计算个节点电压新值,即修正后值3 B& h' ^& O( D# X6 s- z
for n=1:N1
4 Z3 A- a( E: y0 J  v' `! U& N/ D' I+ u9 q    if n<=9
4 W8 I6 f) B- L+ |" a" f( D    delt(n)=delt(n)+uu(2*n-1);
2 r5 P3 x' ^& r. N    u(n)=u(n)+uu(2*n);8 Z: z0 a4 h, a5 Y* w4 e3 S8 J
    else5 N. x, @- X# X: R! i
    delt(n)=delt(n)+uu(n+9);
  j3 U- t: m% S2 I* }    u(n)=u(n);( q8 s. S0 b, \6 V
    end/ C/ Q6 \6 q8 \5 \' D9 A. _$ Z* l
end
9 }; m+ \) W" W, y) \# p, L  J+ V5 v% 循环次数k加一,带入新值进入下一次迭代# a' d: s7 ~5 ]# E2 `! p' N
k=k+1;2 y7 Q6 \8 ?& `% c( T  O  E
end
1 Q! {, Z5 J$ m4 z; Y4 b! k% 迭代完成后,求迭代次数、节点母线电压幅值与相角- i4 D: l+ C; T0 L1 p0 ?' ?2 l
k-1,(delt*180/pi)',u'$ d2 A' W: k% h; F0 r
% 计算各节点的功率和线路功率8 V+ V0 i& F6 `: T/ g
for n=1:N1+1! |  ?8 X+ M! F% L
    U(n)=u(n)*(cos(delt(n))+j*sin(delt(n)));# |0 Y8 G9 I" S* m: t
end
- A! k9 T5 a9 i9 u5 F' bfor m=1:N1+1
/ v+ O8 y7 ?$ p& _0 q0 H    for n=1:N1+1
! z% _1 u+ j6 y3 S/ x    I(n)=Y(m,n)*U(n);
* y! }7 Y* K, B  k/ V5 h  Y3 w    end
& ^9 i( A3 i8 }) y    S(m)=U(m)*sum(conj(I));
9 [& c# f' Z$ s1 X4 f2 X# T1 yend
( Y% T2 d$ D8 Z& |7 |- TS
' |# G+ t0 u  Y" J3 Zfor m=1:N1+1/ j% ^2 D% H" S5 J" ^
    for n=1:N1+1; l! W; x& Q" N/ ?( C- q  p
        T(m,n)=U(m)*conj((U(m)-U(n))*(-Y(m,n)))+U(m)*conj(U(m)*j*y(m,n));
9 p& Y$ k) G, g" a3 U    end" s5 i5 ^( U: `7 i  @
end% |# ]0 D: Q$ s; y/ m( p' l1 I) i/ {: b
T
"真诚赞赏,手留余香"
还没有人打赏,支持一下
楼主热帖
帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】

该用户从未签到

尚未签到

发表于 2010-10-28 14:31:14 | 显示全部楼层
多谢楼主分享。。。。。。。。。
"真诚赞赏,手留余香"
还没有人打赏,支持一下
帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】

该用户从未签到

尚未签到

 楼主| 发表于 2010-10-28 20:40:12 | 显示全部楼层
这个好像算的结果不是很准确,但能学习下思想
"真诚赞赏,手留余香"
还没有人打赏,支持一下
帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】

该用户从未签到

尚未签到

发表于 2010-11-12 16:42:56 | 显示全部楼层
楼主的分享精神不错。
"真诚赞赏,手留余香"
还没有人打赏,支持一下
帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】

该用户从未签到

尚未签到

发表于 2010-12-24 11:13:09 | 显示全部楼层
现在我们编程的缺点是不能用矩阵的思想全局的看待系统问题,会用很多for循环去实现功能
帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】
  • TA的每日心情
    郁闷
    2017-3-30 20:24
  • 签到天数: 2 天

    连续签到: 1 天

    [LV.1]初来乍到

    累计签到:2 天
    连续签到:1 天
    发表于 2011-11-17 19:26:20 | 显示全部楼层
    太复杂了,看不懂、有14节点编的算节点导纳矩阵的程序么
    "真诚赞赏,手留余香"
    还没有人打赏,支持一下
    帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】

    该用户从未签到

    尚未签到

    发表于 2012-4-13 14:04:15 | 显示全部楼层
    牛人啊,佩服佩服
    "真诚赞赏,手留余香"
    还没有人打赏,支持一下
    帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】

    该用户从未签到

    尚未签到

    发表于 2012-4-22 19:12:09 | 显示全部楼层
    这个程序matlab论坛里就有
    "真诚赞赏,手留余香"
    还没有人打赏,支持一下
    帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】
  • TA的每日心情
    奋斗
    2017-3-13 16:44
  • 签到天数: 10 天

    连续签到: 1 天

    [LV.3]偶尔看看II

    累计签到:10 天
    连续签到:1 天
    发表于 2017-1-12 20:05:26 | 显示全部楼层
    多谢楼主,可以借鉴来学习啦
    "真诚赞赏,手留余香"
    还没有人打赏,支持一下
    帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】
    回复 推荐 踩下

    使用道具 举报

    您需要登录后才可以回帖 登录 | 立即加入

    本版积分规则

    招聘斑竹

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

    GMT+8, 2025-4-4 02:49

    Powered by Discuz! X3.5 Licensed

    © 2001-2025 Discuz! Team.

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