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

YOYO玥玥

一部分关于形成雅克比矩阵的程序，看不懂呀，希望各位大侠赐教！
jay = sqrt(-1);
exp_ang = exp(jay*ang);
% Voltage rectangular coordinates
V_rect = V.*exp_ang;
CV_rect=conj(V_rect);
Y_con = conj(Y);
$ m+ X5 @, a* @& P4 B%vector of conjugate currents
9 v- ^5 s3 U) C4 V) r3 Q; a: @4 O: t" {i_c=Y_con*CV_rect;
, S  X/ m8 b, v# h+ \9 d) r0 X% complex power vector
S=V_rect.*i_c;
S=sparse(diag(S));
Vdia=sparse(diag(V_rect));
CVdia=conj(Vdia);
Vmag=sparse(diag(abs(V)));
S1=Vdia*Y_con*CVdia;
% b  n! Y9 X: I/ z* R4 qt1=((S+S1)/Vmag)*volt_red';
t2=(S-S1)*ang_red';
J11=-ang_red*imag(t2);
: W$ O4 D4 y6 X! Z- V, CJ12=ang_red*real(t1);
J21=volt_red*real(t2);
J22=volt_red*imag(t1);
if nargout > 3
/ P, [$ ^3 V5 A+ {# \   Jac11 = J11; clear J11
# z1 S! I, `6 N) Q) W   Jac12 = J12; clear J12
   Jac21 = J21; clear J21
   Jac22 = J22; clear J22
8 B3 X' g0 V1 L0 `! Yelse
   Jac11 = [J11 J12;
               J21 J22];
# K  h8 }/ a9 Q+ S! g. Send

YOYO玥玥

一部分关于雅克比矩阵的，实在不懂，希望大侠们帮帮忙……
jay = sqrt(-1);
exp_ang = exp(jay*ang);
% Voltage rectangular coordinates
V_rect = V.*exp_ang;
( r- C; I- q  ~# E* u5 v/ HCV_rect=conj(V_rect);
2 X0 {  v3 I) b! G* j, |7 q* HY_con = conj(Y);
. K, }7 J  X9 u/ G- ]%vector of conjugate currents
$ I. i, |; |* g. |3 [i_c=Y_con*CV_rect;
. l. N! o# q  ^- \3 ?. H8 V. i% complex power vector
S=V_rect.*i_c;
S=sparse(diag(S));
Vdia=sparse(diag(V_rect));
+ |0 v- g2 I% P$ n6 [* O8 z7 e8 NCVdia=conj(Vdia);
  C4 T1 f$ y* q4 w& V7 g* d1 zVmag=sparse(diag(abs(V)));
  T& A( r; {: D- US1=Vdia*Y_con*CVdia;
) ]. v# d: y3 {3 Z1 c  Pt1=((S+S1)/Vmag)*volt_red';
9 c" N% r* D8 Z, ]( xt2=(S-S1)*ang_red';
- [- R2 v+ F( R0 @J11=-ang_red*imag(t2);
J12=ang_red*real(t1);
J21=volt_red*real(t2);
J22=volt_red*imag(t1);
1 t$ ]0 Z2 U. V$ ?6 ^if nargout > 3
   Jac11 = J11; clear J11
5 {* N1 Q4 T9 d' R1 s' F   Jac12 = J12; clear J12
: m! J% Y/ `$ v   Jac21 = J21; clear J21
" ~2 B: s7 C4 I% T   Jac22 = J22; clear J22
9 E: b/ b) a3 M: u" J$ Oelse
3 u$ R7 w9 k1 i: ~; `  {   Jac11 = [J11 J12;
' |; D5 o: K* a8 q) O0 N4 @6 T               J21 J22];
end

wyh4

谢分享的资料,3Q

syh0603010629

好东西啊

百草春生

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

hz04230807

谢谢楼主，顶，学习下

hotwind_zmo

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

Missing...

貌似不错啊！！！

zibuyu_21

谢谢了

zyy2379

学习了