function [ mu_korr, i_korr, nu_korr ] = get_korr( rez, t )
% Возвращает скорректированные кэффициенты по номеру режима и времени
% регулирования

% Корректированный = Опорный - макс. отклонение

A_dop = 0.05;
Fi_dop= 0.1; % -0.1<Fi<0

%dA_tek = 0;

% из графика АФЧХ опорного режима при Fi=-80
w_min = 0.6;
w_max = 1.5;
w_dop = w_max-w_min;


% Передаточная функция для опорного режима
[ a1, a2, a3, a4, a5, a6, a7, b1, b2, b3, b4, b5, b6, b7 ]=make_koefs(2);
mu_o = (18-b1*t)/(b3*t);
i_o=108/(b3*t^2);
nu_o=216/(b3*t^3);
G1 = get_tf_kren_zamkn(mu_o, i_o, nu_o, b1, b3);   % t min

dW_tek = 0.08; % просто участок частоты
% значение производной по амплитуде dA/dw в точке w_min
%awo = get_Awo(G1, w_min, dW_tek);
% значение производной по фазе dFi/dw в точке w_min
%fwo = get_Fwo(G1, w_min, dW_tek);
[ awo, fwo ] = getMaxAFwo( G1, w_min,w_max );

% графики dA/dmu, dA/di, dA/dnu на опорном режиме
[X1,Y1,X2,Y2,X3,Y3] = get_dAdz_arrays(2, t);

% максимальные значения по амплитуде с учетом знака
[x, a_mu] = get_max_val(X1, Y1, w_min, w_max);
[x, a_i] = get_max_val(X2, Y2, w_min, w_max);
[x, a_nu] = get_max_val(X3, Y3, w_min, w_max);


[X1,Y1,X2,Y2,X3,Y3] = get_dFdz_arrays(2, t);
[x, f_mu] = get_max_val(X1, Y1, w_min, w_max);
[x, f_i] = get_max_val(X2, Y2, w_min, w_max);
[x, f_nu] = get_max_val(X3, Y3, w_min, w_max);


%fprintf('dA/dmu max = %f\n', a_mu);
%fprintf('dA/dnu max = %f\n', a_nu);




[ a1, a2, a3, a4, a5, a6, a7, b1, b2, b3, b4, b5, b6, b7 ]=make_koefs(rez);

% номинальные параметры на режиме
mu_nom = (18-b1*t)/(b3*t);
if (mu_nom<0)
    mu_nom = 0.01;
end
i_nom=108/(b3*t^2);
nu_nom=216/(b3*t^3);



% максимальные отклонения, больше которых быть не может
max_dmui = mu_nom - mu_o; % расчетное значение минус значение оп. режима
max_dii = i_nom - i_o;
max_dnui = nu_nom - nu_o;
%if (max_dmui<0 || max_dii<0)
    
%end

% Записываем векторы и матрицы для функции linprog
% f - 
% нужно чтобы сумма dmu+di+dnu была максимальной при опред. условиях
% или же

%    dmu di dnu dw_тек
f = [1;  1; 1;  0];
f = -f; % так как linprog ищет минимум, а не максимум

A =  [a_mu  a_i  a_nu   awo
      -a_mu -a_i -a_nu  -awo
      f_mu  f_i  f_nu   fwo
      -f_mu  -f_i  -f_nu   -fwo];
%A = [a_mu  a_i  a_nu   awo
%     f_mu  f_i  f_nu   fwo];

b = [ A_dop; -A_dop+0.01; -Fi_dop; Fi_dop-0.01];
%b = [ A_dop; Fi_dop];

% векторы нижних и верхних границ решений
%lb = [0.01; 0.01; 0.01; 0];
%ub = [max_dmui; max_dii; max_dnui; w_dop];
lb = [0; 0; 0; 0];
ub = [max_dmui; max_dii; max_dnui; w_dop];

% Решение системы уравнений, 
[x,fval,exitflag,output,lambda] = linprog(f,A,b,[],[],lb, ub);

% В результате - x - вектор, содержащий, макс. отклонения от опорного
% режима, вычитаем их от номинального значения
%x,lambda.ineqlin,lambda.lower
mu_korr = mu_nom-x(1);
i_korr = i_nom-x(2);
nu_korr = nu_nom-x(3);

%x,lambda.ineqlin,lambda.lower
%fprintf('\n---------------\n');
%end
end