write.m

电力系统分析计算程序

function write(a,fid,buslist)
% write voltage dependent loads loads

if ~a.n, return, end

% filter loads using bus list
idx = [];
for i = 1:a.n
  jdx = find(buslist == a.bus(i)*a.u(i));
  if ~isempty(jdx), idx = [idx; i]; end
end

if isempty(idx), return, end

% write Mn data
fprintf(fid,'Mn.con = [ Mn.con; ...\n');
fprintf(fid,['   ',a.format,';\n'],a.con(idx,:)');
fprintf(fid,'   ];\n\n');