numdiff.m

这是一个很适合研究和学习用的电力系统仿真软件

function numdiff(a,ffn,xa,Va,anga)

global NLA Settings DAE 

if ~a.n, return, end

% B matrix using numerical differentiation

fm_nrlf(Settings.lfmit,Settings.lftol,0);

incY = 0; 

for j = 1:a.n

  deltaU = max(NLA.tol,abs(NLA.tol*a.Vref(j)));
  old = a.Vref(j);
  a.Vref(j) = a.Vref(j)+deltaU;
  % update algebraic variables
  fm_nrlf(Settings.lfmit,Settings.lftol,0);
  a = fcall(a);
  ff = DAE.f;
  NLA.b_svc(:,j)=(ff-ffn)./deltaU;
  DAE.x = xa;
  DAE.a = anga;
  DAE.V = Va;   
  a.Vref(j) = old;
  
end