fm_errv.m

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

function fm_errv(Vold,msg,busidx)
% FM_ERRV compute absolute voltage rate errors
%
% FM_ERRV(VRATE,MSG,BUSIDX)
%       VRATE  -> component voltage rate
%       MSG    -> string for displaying message
%       BUSIDX -> component indexes
%
%Author:    Federico Milano
%Date:      25-Dic-2005
%Version:   1.0.0
%
%E-mail:    [e-mail]
%Web-site:  http://thunderbox.uwaterloo.ca/~fmilano
%
% Copyright (C) 2002-2006 Federico Milano

global Varname Bus

Vbus = Bus.con(busidx,2);
Verr = abs((Vbus-Vold)./Vbus);
idx = find(Verr > 0.1);
for iii = 1:length(idx)
  kkk = idx(iii);
  fm_disp(['Voltage rate of ',msg,' #', num2str(kkk), ...
           ' at Bus ',Varname.bus{busidx(kkk)}, ...
           ' differs more than 10% from Bus voltage rate'],2)
end