%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%% Modelos Minto Remifentanil %%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%                                                                         %
% Modelos de Minto para o Remifentanil                                    %
%                                                                         %
% Miguel Teixeira, Novembro, 2012                                         %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%% Variáveis
V1 = 5.1-0.0201*(Age-40)+0.072*(LBM-55);                                   % [L]
V2 = 9.82-0.0811*(Age-40)+0.108*(LBM-55);                                  % [L]
V3 = 5.42;                                                                 % [L]

k10 = (2.6-0.0162*(Age-40)+0.0191*(LBM-55))/V1;                            % [min^-1]
k12 = (2.05-0.0301*(Age-40))/V1;                                           % [min^-1]
k13 = (0.076-0.00113*(Age-40))/V1;                                         % [min^-1]
k21 = k12*V1/V2;                                                           % [min^-1]
k31 = k13*V1/V2;                                                           % [min^-1]
ke0 = 0.595-0.007*(Age-40);                                                % [min^-1]

%% Converter variavéis
V1 = V1*1000;                                                              % [mL]
V2 = V2*1000;                                                              % [mL]
V3 = V3*1000;                                                              % [mL]

%% Matrizes Compartimentais
Matriz_A = [-(k10+k12+k13) (V2/V1)*k21 (V3/V1)*k31 0;
           (V1/V2)*k12 -k21 0 0;
           (V1/V3)*k13 0 -k31 0;
           ke0 o o -ke0];
Matriz_B = [1/V1 0 0 0]';
Matriz_C = [0 0 0 1];
Matriz_D = 0;

% Discretização das Matrizes
[Matriz_A, Matriz_B] = c2d(Matriz_A, Matriz_B, Ts);                        % As matrizes C e D discretizadas ficam iguais às contínuas, logo não é necessária a sua discretização.
                                                                           % Ts é o tempo de amostragem em minutos (Ex: 5s é 1/12 de minuto).
%% Sistema de Espaço de Estados
Estado = Matriz_A*Estado+Matriz_B*U;                                       % Estado são as concentrações nos diferentes compartimentos [ug/mL]. U é o Input de Remifentanil [ug/min].
Ce = Matriz_C*Estado+Matriz_D;                                             % Obtenção da Concentração de Efeito [ug/mL].