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

clear;%очищаем рабочее пространство
Ts = 10e-6;%шаг интегрирования
Decim = 20;%интервал прореживания
Limit = 2000000;%кол-во запоминаемых точек
% Scope-ы запоминают последнии tt секунд расчёта
tt = Ts*Decim*Limit;


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

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

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

% номинальные величины ГЭД
% ... мощность на валу, Вт
Pnom = 5000e3*2;
% ... линейное напряжение, В (rms)
Unom = 3000;
% ... механическая скорость, об/мин
NmNom = 165;
% ... число пар полюсов
Pp = 6;
% ... коэффициент мощности
CosFi = 0.87;
% ... КПД
Eff = 0.962;
% ... приведенный к валу момент инерции, кг*м^2
J = 87e3*0.5;
% ... полная мощность, ВА
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*1.18;
% ... ток, А (ampl)
Inom = Snom/(Unom*SQRT3)*SQRT2*0.5;%0.5 - т.к. обмоток две

% параметры ГЭД
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.87 * 1.5;%2.87;%Ом
Fe = 12;%Гц
Lls = Xls/(Fe*PI2);%Гн
Llr = Xlr/(Fe*PI2);%Гн
Lm = Xm/(Fe*PI2);%Гн

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

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