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

YOYO玥玥

一部分关于形成雅克比矩阵的程序，看不懂呀，希望各位大侠赐教！
jay = sqrt(-1);
exp_ang = exp(jay*ang);
% Voltage rectangular coordinates
V_rect = V.*exp_ang;
  s. h4 ^! O: v" TCV_rect=conj(V_rect);
Y_con = conj(Y);
%vector of conjugate currents
i_c=Y_con*CV_rect;
& d1 k  r6 g* d2 [* @& [& S% complex power vector
  {& S$ k6 D( z/ j7 JS=V_rect.*i_c;
: E  p8 c. ]$ M9 \S=sparse(diag(S));
) c, Q) c- R* m6 R3 BVdia=sparse(diag(V_rect));
CVdia=conj(Vdia);
$ D% n$ R' f1 @, w" c: E) oVmag=sparse(diag(abs(V)));
2 w/ I" h4 @6 O* L2 ^S1=Vdia*Y_con*CVdia;
t1=((S+S1)/Vmag)*volt_red';
t2=(S-S1)*ang_red';
J11=-ang_red*imag(t2);
J12=ang_red*real(t1);
4 \# U  O, ?8 x) yJ21=volt_red*real(t2);
J22=volt_red*imag(t1);
" H* r. C. S1 G' r  Gif nargout > 3
   Jac11 = J11; clear J11
( x- R1 ~1 a% {1 G7 v8 z$ }   Jac12 = J12; clear J12
; M  O  t$ y1 `8 v# e3 r$ ?; q- n   Jac21 = J21; clear J21
   Jac22 = J22; clear J22
  H9 }1 z) D6 J5 belse
   Jac11 = [J11 J12;
) C* F# [3 i6 l. G; T! F& V               J21 J22];
end

YOYO玥玥

一部分关于雅克比矩阵的，实在不懂，希望大侠们帮帮忙……
jay = sqrt(-1);
exp_ang = exp(jay*ang);
% Voltage rectangular coordinates
V_rect = V.*exp_ang;
CV_rect=conj(V_rect);
1 o4 C& K' K2 ]/ N, I1 _, l. uY_con = conj(Y);
%vector of conjugate currents
1 r- F4 W" A% j7 u8 U1 a% di_c=Y_con*CV_rect;
% complex power vector
6 ~% u% s, k- I' f/ n. Z8 cS=V_rect.*i_c;
S=sparse(diag(S));
Vdia=sparse(diag(V_rect));
$ o/ m* [. T) `$ v- c' r: u4 |5 B4 PCVdia=conj(Vdia);
Vmag=sparse(diag(abs(V)));
* L# d* R/ x9 iS1=Vdia*Y_con*CVdia;
" N- T! L$ h# N  b$ o- ft1=((S+S1)/Vmag)*volt_red';
t2=(S-S1)*ang_red';
J11=-ang_red*imag(t2);
: `; \  \! l" ?" z  @J12=ang_red*real(t1);
' m  M6 G9 a+ n, w8 R# ]7 kJ21=volt_red*real(t2);
8 X& W1 h9 I& h; J! jJ22=volt_red*imag(t1);
: F1 o2 I) ~$ Iif nargout > 3
   Jac11 = J11; clear J11
   Jac12 = J12; clear J12
   Jac21 = J21; clear J21
) _7 M- M" Q  L8 ~6 `' f   Jac22 = J22; clear J22
else
   Jac11 = [J11 J12;
               J21 J22];
end

wyh4

谢分享的资料,3Q

syh0603010629

好东西啊

百草春生

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

hz04230807

谢谢楼主，顶，学习下

hotwind_zmo

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

Missing...

貌似不错啊！！！

zibuyu_21

谢谢了

zyy2379

学习了