% Cкрипт для задания параметров модели

clear;%очищаем рабочее пространство

%% ПАРАМЕТРЫ МОДЕЛИ

addpath('MCU_Wrapper');
addpath('motor');
Ts = 10e-6;%шаг интегрирования
Decim = 1;%интервал прореживания
DisableScope = {
    "Idc";
    "Udc";
    };

GED = "23550";
% GED = "22220";

% начальная скорость ГЭД, доля от NmNom
w0 = 0;%0.5;%-0.75;%
% пусковой момент, о.е.
Mst = 0.6;%0.6;

% разрешаем/запрещаем сбросы/набросы момента нагрузки
changingLoadEnable = 0;%1
% разрешаем/запрещаем шум в измеренном токе
noiseEnable = 0;%1;%
% ... мощность шума
NP = 0.08;


%% НОМИНАЛЬНЫЕ ВЕЛИЧИНЫ ГЭД
% ... мощность на валу, Вт
Pnom = 6300e3;
% ... линейное напряжение, В (rms)
Unom = 3300;
% ... механическая скорость, об/мин
NmNom = 180;
% ... число пар полюсов
Pp = 6;
% ... коэффициент мощности
CosFi = 0.87;
% ... КПД
Eff = 0.968;
% ... приведенный к валу момент инерции, кг*м^2
J = 87e3*0.1;


%% РАСЧЕТЫ
% разкомментирование всех блоков

modelName = [bdroot '/Measurements'];
blocks = find_system(modelName,  ...
    'IncludeCommented', 'on', ...
    'FollowLinks', 'on', ...
    'LookUnderMasks', 'all', ...
    'BlockType', 'Scope'); 
for i = 1:length(blocks)
    set_param(blocks{i}, 'Commented', 'off');
end
% отключение графиков для ускорения смуляции
for i = 1:length(DisableScope)
    set_param([modelName '/'] + DisableScope{i}, 'Commented', 'on');
end

% для упрощения записи
SQRT2 = sqrt(2);
SQRT3 = sqrt(3);
PI2 = pi*2;

% ... полная мощность, ВА
Snom = Pnom/CosFi/Eff;
% ... механическая скорость, рад/с
WmNom = NmNom/60*PI2;
% ... момент на валу, Н*м
Mnom = Pnom/WmNom;
% ... эл. скорость, рад/с
WeNom = WmNom*Pp;
% ... эл. скорость, Гц
FeNom = WeNom/PI2;
% ... потокосцепление статора, Вб
PsiNom = Unom*SQRT2/(WeNom*SQRT3);
% ... напряжение на входе инвертора, B
UdcNom = Unom*SQRT2;
% ... ток, А (ampl)
Inom = Snom/(Unom*SQRT3)*SQRT2*0.5;%0.5 - т.к. обмоток две

% схема замещения ГЭД
if GED == "22220"
    GED
    Rs = 11.8e-3;%Ом
    Xls = 72.7e-3;%72.7e-3;%Ом
    Rr = 11.1e-3*2.0;%*0.8;%Ом
    Xlr = 85.5e-3;%Ом
    Xm  = 2.9322;%2.87;%Ом
    Fe = 18;%Гц
    Lls = Xls/(Fe*PI2);%Гн
    Llr = Xlr/(Fe*PI2);%Гн
    Lm = Xm/(Fe*PI2);%Гн
elseif GED == "23550"
    GED
    Rs = 0.0282;%Ом
    Xls = 0.4016;%Ом
    Rr = 0.139;%Ом
    Xlr = 0.2006;%Ом
    Xm  = 5.2796;%Ом
    Fe = 18.2;%Гц
    Lls = Xls/(Fe*PI2);%Гн
    Llr = Xlr/(Fe*PI2);%Гн
    Lm = Xm/(Fe*PI2);%Гн
end

% ёмкость на входе INU, Ф
Cdc = 50e-3;
% снаберы в INU
Csn = Pnom/(1000*WeNom*Unom^2)/10;%Ф (0.5 - т.к. преобразователей два)
Rsn = 2*Ts/Csn*10;%Ом

% постоянная времени фильтра для тока ГЭД, c
Tiac = 30e-6;