Выполнено численное моделирование односкачкового и двухскачкового распространения частотно-модулированных сигналов в анизотропной ионосферной плазме на основе численного решения пространственно-временной бихарактеристической системы Гамильтона, неизвестными в которой являются компоненты волнового вектора, координаты, а также частота и время. Предполагалось, что источник излучения точечный, находится вне магнитоактивной плазмы, частотная модуляция декаметровых радиоволн линейная, а амплитуда магнитного поля постоянная. Рассмотрена модель двухслойной ионосферной плазмы с волнообразным возмущением. Проанализированы особенности выхода обыкновенной и необыкновенной волн из плоскости распространения, образование каустик пространственно временных геометрооптических лучей и доплеровское изменение частоты.