function [ mu_korr, i_korr, nu_korr ] = get_korr( rez, t ) % Получить скорректированные кэффициенты по номеру режима и %GET_KORR Summary of this function goes here % Detailed explanation goes here A_dop = 0.05; Fi_dop= 0.1; %dA_tek = 0; dW_tek = 0.08; 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_zamkn(mu_o, i_o, nu_o, b1, b3); % t min awo = get_Awo(G1, w_min, dW_tek); fwo = get_Fwo(G1, w_min, dW_tek); [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); [ a1, a2, a3, a4, a5, a6, a7, b1, b2, b3, b4, b5, b6, b7 ]=make_koefs(rez); mu_nom = (18-b1*t)/(b3*t); i_nom=108/(b3*t^2); nu_nom=216/(b3*t^3); dmui = mu_nom - mu_o; dii = i_nom - i_o; dnui = nu_nom - nu_o; % dmu di dnu dw f = [1; 1; 1; 0]; f = -f; 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]; b = [ A_dop; -A_dop+0.001; Fi_dop; -Fi_dop+0.0001]; lb = [0; 0; 0; 0.1]; ub = [dmui; dii; dnui; w_dop]; [x,fval,exitflag,output,lambda] = linprog(f,A,b,[],[],lb, ub); %x 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