%% 
%mu=0.566;
%i=1.087;
%nu=0.527;

% Опорный режим
R_opor = 2;

% Время регулирования
t_reg_min = 2;
t_reg_max = 5;

% Коэффициенты получать так:
% [ a1, a2, a3, a4, a5, a6, a7, b1, b2, b3, b4, b5, b6, b7 ]=make_koefs(12);

%%
% --------- что сделать сейчас: 1 - показать, 0 - не показывать
% Построить АФЧХ опорного режима для Tmin, Tmax
make_AF_R_opor = 0;

% Построить амплитудные функции чувствительности на опорном режиме
make_dA_dz_graphics = 1;



%% -----------------
% Построить АФЧХ опорного режима
%

% получаем коэф-ты для него
[ a1, a2, a3, a4, a5, a6, a7, b1, b2, b3, b4, b5, b6, b7 ]=make_koefs(R_opor);

% Числа, соотв. времени регулирования.. мин и макс
mu_min = (18-b1*t_reg_min)/(b3*t_reg_min);
mu_max = (18-b1*t_reg_max)/(b3*t_reg_max);

i_min=108/(b3*t_reg_min^2);
i_max=108/(b3*t_reg_max^2);

nu_min=216/(b3*t_reg_min^3);
nu_max=216/(b3*t_reg_max^3);

%mu=mu_min; nu=nu_min; i=i_min; % min
%mu=mu_max; nu=nu_max; i=i_max; % max

% получаем передаточные функции АП на опорном режиме при Тmin и Tmax
% зависят от переменных: мю, и, ню, b1, b3   (см. структурную схему АП)
G1 = get_tf_kren_zamkn(mu_min, i_min, nu_min, b1, b3);   % t min
G2 = get_tf_kren_zamkn(mu_max, i_max, nu_max, b1, b3);   % t max

P = bodeoptions; % Set phase visiblity to off and frequency units to Hz in options
P.MagUnits = 'abs';
P.MagScale = 'linear';
P.FreqScale = 'linear';
P.Xlim = [0, 1.8];
P.Ylim = {[0, 1],[-90, 0]};
P.MagLowerLimMode = 'manual';
P.MagLowerLim = -100;
P.Grid = 'on';
% АФЧХ при Тмин и Тмакс
%h = bodeplot(G1,'r',G2,'b--',P)
if (make_AF_R_opor==1)
    bodeplot(G1,'r',G2,'b--',P)
    %hold on
end

%% 
% Построить амплитудные функции чувствительности на опорном режиме
%make_dA_dz_graphics = 1;
%[ X1,Y1,X2,Y2,X3,Y3 ] = get_dAdz_arrays( R_opor, t_reg_min );
[ X1,Y1,X2,Y2,X3,Y3 ] = get_dAdz_arrays( 2, 2 );
%plot(X1,Y1,X2,Y2,X3,Y3)
%plot(X1,Y1)

grid on
% Построить фазовые функции чувствительности на опорном режиме
%[ X1,Y1,X2,Y2,X3,Y3 ] = get_dFdz_arrays( R_opor, t_reg_min );
plot(X1,Y1,X2,Y2,X3,Y3)
grid on

asd
%%
%[a,b]=make_koefs(1);
%[ a1, a2, a3, a4, a5, a6, a7, b1, b2, b3, b4, b5, b6, b7 ]=make_koefs(1);


%sdf=aa

%for k = 1:1:12
%[ a1, a2, a3, a4, a5, a6, a7, b1, b2, b3, b4, b5, b6, b7 ]=make_koefs(k);
%fprintf('R=%d\n',k);
%mu_min = (18-b1*t_reg_min)/(b3*t_reg_min)
%i=108/(b3*t_reg_min^2);
%fprintf('i=%0.3f\n\n',i);
%end
%mu_min = (18-b1*t_reg_min)/(b3*t_reg_min)
%mu_max = (18-b1*t_reg_max)/(b3*t_reg_max)



%%%%%%%%%%%%%%%%%%%%%%%%%


%% all R
A_dop = 0.05;
Fi_dop= 0.1;

[ a1, a2, a3, a4, a5, a6, a7, b1, b2, b3, b4, b5, b6, b7 ]=make_koefs(2);
mu = (18-b1*t_reg_min)/(b3*t_reg_min);
i=108/(b3*t_reg_min^2);
nu=216/(b3*t_reg_min^3);
%mu=0.566;
%i=1.087;
%nu=0.527;
%[ a1, a2, a3, a4, a5, a6, a7, b1, b2, b3, b4, b5, b6, b7 ]=make_koefs(12);

%[mu, i, nu] = get_korr(1, 2)


fprintf('Скорректированные коэф-ты:\n');
%for k = 12:1:12
%[ a1, a2, a3, a4, a5, a6, a7, b1, b2, b3, b4, b5, b6, b7 ]=make_koefs(k);
%fprintf('Режим=%d:\n',k);
[ MU, I, NU ] = getCorrectedVectors( )