Логарифм десятичный с защитой 0

• a, ω, f – амплитуда, круговая частота и фаза
• x(t), y(t) – входной и выходной сигналы. b - константа, определяющая ограничение "снизу". Таким образом, блок работоспособен при любых изменениях входного сигнала, а расчет не останавливается при пересечении входным сигналом нуля.

//Описание блока, возвращающего выходное значение, вычисленное
//по алгоритму:
//y(t) = a·lg(ω·x(t) + f), если ω·x(t) + f > 0
//y(t) = b, если ω·x(t) + f ≤ 0
//y(t) = b, если ω·x(t) + f > 0 и a·lg(ω·x(t) + f) < b
// a, ω, f - амплитуда, круговая частота и фаза
// x(t), y(t) - входной и выходной сигналы
// b - константа, определяющая ограничение "снизу"
//Таким образом, блок работоспособен при любых изменениях входного сигнала,
//а расчёт не останавливается при пересечении входным сигналом нуля
//Блок векторизован
// a - амплитуда, в функции a·lg[w·x(t)+f]. Коэффициент a. По умолчанию [1]
// w - частота. Коэффициент ω. По умолчанию [1]
// f - фаза. Коэффициент f. По умолчанию [0]
// b - значение выхода при x <= 0. На выходе будет сформировано это значение
//в те моменты, когда вход не положительный, либо когда выход стал меньше значения b.
//Фактически, это ограничение "снизу" для выходного сигнала. По умолчанию -2
//Свойства могут быть векторами (кроме b), размерности которых должны быть равны
//размерности входа, в этом случае на выходе будет также вектор
//Входы: x - порт ввода аргумента функции
//Выходы: outport - порт вывода значения функции

block: (
  type = "Логарифм десятичный с защитой 0",
  points=[(-200 , 0),(-184 , 0),(-200 , -16),(408 , 308)],
  a = [1],
  w = [1],
  f = [0],
  b = -2
  ),