交流社区 › 电力系统与电网技术 › 基于Matlab的电力系统潮流计算

YOYO玥玥 发表于 2010-4-19 20:39:41

一部分关于形成雅克比矩阵的程序，看不懂呀，希望各位大侠赐教！
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);
%vector of conjugate currents
i_c=Y_con*CV_rect;
% 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;
t1=((S+S1)/Vmag)*volt_red';
t2=(S-S1)*ang_red';
J11=-ang_red*imag(t2);
J12=ang_red*real(t1);
J21=volt_red*real(t2);
J22=volt_red*imag(t1);
if nargout > 3
   Jac11 = J11; clear J11
   Jac12 = J12; clear J12
   Jac21 = J21; clear J21
   Jac22 = J22; clear J22
else
   Jac11 = [J11 J12;
             J21 J22];
end

YOYO玥玥 发表于 2010-4-19 20:44:51

一部分关于雅克比矩阵的，实在不懂，希望大侠们帮帮忙……
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);
%vector of conjugate currents
i_c=Y_con*CV_rect;
% 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;
t1=((S+S1)/Vmag)*volt_red';
t2=(S-S1)*ang_red';
J11=-ang_red*imag(t2);
J12=ang_red*real(t1);
J21=volt_red*real(t2);
J22=volt_red*imag(t1);
if nargout > 3
   Jac11 = J11; clear J11
   Jac12 = J12; clear J12
   Jac21 = J21; clear J21
   Jac22 = J22; clear J22
else
   Jac11 = [J11 J12;
             J21 J22];
end

wyh4 发表于 2010-4-21 22:30:39

谢分享的资料,3Q

syh0603010629 发表于 2010-4-25 15:20:02

:shutup:好东西啊

百草春生 发表于 2010-5-17 19:12:08

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

hz04230807 发表于 2010-5-22 10:37:23

谢谢楼主，顶，学习下

hotwind_zmo 发表于 2010-5-23 16:12:23

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

Missing... 发表于 2010-5-24 16:50:03

貌似不错啊！！！

zibuyu_21 发表于 2016-4-11 16:45:49

谢谢了

zyy2379 发表于 2016-4-19 15:09:14

学习了

查看完整版本: 基于Matlab的电力系统潮流计算