Модификация метода расщепления для моделирования многосердцевинных волноводов с учётом насыщающегося усиления

Представлен эффективный безытерационный численный алгоритм для решения системы линейно связанных нелинейных уравнений Шрёдингера (НУШ), описывающих распространение световых импульсов в многосердцевинных оптических волокнах с учётом оптических потерь и насыщающегося усиления. Метод также может быть использован и в случае одного скалярного НУШ. Предлагаемый алгоритм основан на методе Фурье расщепления по физическим процессам, обладает вторым порядком точности по эволюционной переменной и экспоненциальным порядком точности по времени. Для возможности решения системы связанных НУШ предлагается использовать аппроксимацию матричной экспоненты с помощью разложения Паде.