基于Matlab的电力系统潮流计算

YOYO玥玥

一部分关于形成雅克比矩阵的程序，看不懂呀，希望各位大侠赐教！
jay = sqrt(-1);
exp_ang = exp(jay*ang);
% Voltage rectangular coordinates
V_rect = V.*exp_ang;
CV_rect=conj(V_rect);
9 ]1 t4 c8 d, B' e' z7 ZY_con = conj(Y);
%vector of conjugate currents
+ P! F1 |* }! t9 e) @i_c=Y_con*CV_rect;
% complex power vector
& X0 b" z1 d; F. Z% G) V- FS=V_rect.*i_c;
S=sparse(diag(S));
Vdia=sparse(diag(V_rect));
+ H! w7 a+ C8 M: ]CVdia=conj(Vdia);
Vmag=sparse(diag(abs(V)));
2 G" T# ^% N1 B$ m+ G6 XS1=Vdia*Y_con*CVdia;
t1=((S+S1)/Vmag)*volt_red';
: {) r7 T; A7 n0 ~, e. i2 st2=(S-S1)*ang_red';
J11=-ang_red*imag(t2);
6 i! @( n$ S' ]: ~( }+ q, BJ12=ang_red*real(t1);
. P* ~3 N3 h4 {$ J7 w2 UJ21=volt_red*real(t2);
J22=volt_red*imag(t1);
if nargout > 3
   Jac11 = J11; clear J11
   Jac12 = J12; clear J12
( W7 R5 [- l4 U, e   Jac21 = J21; clear J21
9 w: A8 `/ E" ~" b9 a8 ^5 B$ g: m1 H   Jac22 = J22; clear J22
9 F$ F$ m+ k8 E0 v  h, Z9 k( Gelse
   Jac11 = [J11 J12;
) G! f, H7 B& W  ^" ?$ O               J21 J22];
end

YOYO玥玥

一部分关于雅克比矩阵的，实在不懂，希望大侠们帮帮忙……
jay = sqrt(-1);
exp_ang = exp(jay*ang);
  U3 D& o9 ?' k* b' U: S, Q* H% Voltage rectangular coordinates
; Q+ S$ X3 ~5 }; i$ l9 ]* oV_rect = V.*exp_ang;
: m$ U2 V1 Y3 g  ECV_rect=conj(V_rect);
* g) d' q" k$ Q- YY_con = conj(Y);
%vector of conjugate currents
i_c=Y_con*CV_rect;
% complex power vector
5 O. w$ M  C# g2 L( oS=V_rect.*i_c;
/ j6 q! W8 Q0 H) WS=sparse(diag(S));
Vdia=sparse(diag(V_rect));
3 ]" q, l' f3 Y7 T% k: e4 SCVdia=conj(Vdia);
Vmag=sparse(diag(abs(V)));
S1=Vdia*Y_con*CVdia;
: C9 O4 f. I5 {) k/ Lt1=((S+S1)/Vmag)*volt_red';
t2=(S-S1)*ang_red';
9 g1 Q- H0 _% }) ]1 V, h: ]J11=-ang_red*imag(t2);
J12=ang_red*real(t1);
J21=volt_red*real(t2);
- a( L- h. Q4 O8 ^) T3 rJ22=volt_red*imag(t1);
if nargout > 3
   Jac11 = J11; clear J11
   Jac12 = J12; clear J12
' b0 X! W+ @2 a! L9 l! q6 [* S   Jac21 = J21; clear J21
   Jac22 = J22; clear J22
: K4 T3 D, K" d4 X- N$ Felse
" B0 Q) o: |0 J1 X) l5 c/ y   Jac11 = [J11 J12;
               J21 J22];
end

wyh4

谢分享的资料,3Q

syh0603010629

好东西啊

百草春生

楼主可以把那个发我吗？下不了，谢谢了wjdf1234@163.com

hz04230807

谢谢楼主，顶，学习下

hotwind_zmo

该文献有硬伤啊。。牛顿拉夫逊法不是Newtown-Raphson 么。。。。

Missing...

貌似不错啊！！！

zibuyu_21

谢谢了

zyy2379

学习了