fm_vs.m

用于电力系统的一个很好的分析软件

function fm_vs(k,idx0)
% FM_VS computes determinants and eigenvalues during
%       time domain integrations
%
% FM_VS(K)
%      K    integration step number
%      IDX0 starting index of the output array
%
%Author:    Federico Milano
%Date:      11-Nov-2002
%Version:   1.0.0
%
%E-mail:    [e-mail]
%Web-site:  http://www.uclm.es/area/gsee/Web/Federico
%
% Copyright (C) 2002-2008 Federico Milano

global Settings Bus DAE Varout Line

% As
As = DAE.Fx - DAE.Fy*(DAE.Gy\DAE.Gx);
AutoState = round(eig(full(As))*1e5)/1e5]';

% Jlf
Jlf = build_gy(Line);
detJlf = det(Jlf);
Jlfr = (Jlf(Bus.v,Bus.v)-Jlf(Bus.v,Bus.a)* ...
        (Jlf(Bus.a,Bus.a)\Jlf(Bus.a,Bus.v)));
AutoJlfr = eig(full(Jlfr))';

% Jlfv
detGy = det(DAE.Gy);
Gyr = (DAE.Gy(Bus.v,Bus.v)-DAE.Gy(Bus.v,Bus.a)*...
       (DAE.Gy(Bus.a,Bus.a)\DAE.Gy(Bus.a,Bus.v)));
AutoJlfv = eig(full(Jlfr))';

% Jlfd
Jlfd = (DAE.Gy - DAE.Gx*(DAE.Fx\DAE.Fy));
detJlfd = det(Jlfd);
Jlfdr = (Jlfd(Bus.v,Bus.v)-Jlfd(Bus.v,Bus.a)* ...
         (Jlfd(Bus.a,Bus.a)\Jlfd(Bus.a,Bus.v)));
AutoJlfd = eig(full(Jlfdr))';

vars = [detJlf,detGy,detJlfd,AutoState,AutoJlfr,AutoJlfv,AutoJlfd];
idxs = Varout.idx(find(Varout.idx>idx0))-idx0;

Varout.vars(k,[idx0+1:end]) = vars(idxs);