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