Доклады о будущих и современных технологиях
ПАКЕТЫ КОМПЬЮТЕРНОЙ МАТЕМАТИКИ В ЗАДАЧАХ УПРАВЛЕНИЯ
К. А. Ракутин, М. П. Цыганков
Научный руководитель - М. П. Цыганков, д-р техн. наук, профессор Ярославский государственный технический университет
Для моделирования систем управления технологическими процессами все более широкое применение находят пакеты компьютерной математики (ПКМ). Они включают мощные инструментальные средства, позволяющие оперировать сложными математическими объектами без привлечения программирования.
Тем не менее, в необходимых случаях средства программирования могут быть привлечены для повышения эффективности процедуры решения задачи.
Очень важна возможность работы в ПКМ с группами данных. Это достигается за счет применения инструментария матричных и векторизованных (поэлементных) вычислений. Такой возможности лишены пакеты классического программирования. Групповая обработка позволяет рассматривать группы скалярных переменных как объекты с определенными для них методами такой обработки.
Векторизованные операции делают возможным поэлементное выполнение вычислений, без использования циклов, что повышает наглядность программного кода и его компактность. Это значительно упрощает процесс подготовки модели к выполнению имитационных экспериментов. Оказывается удобной ее коррекция по результатам их выполнения.
В работе решается задача имитации рециркуляции тепловой энергии в высокотоннажном реакторном процессе производства технического углерода. Целью имитации является анализ эффективности управления рециркулирующим тепловым потоком с точки зрения экономии топлива (природного газа) и повышения выхода продукта на пропущенное сырье.
Использован ПКМ Ма1ЬаЬ, допускающий возможности построения п - мерных матриц, применение инструментария матричных и векторизованных вычислений, численного и аналитического интегрирования для систем дифференциальных уравнений без привлечения программирования. Например, вычисление текущего значения вектора С мольных теплоемкостей путем аппроксимации их зависимостей от температуры полиномом третьего порядка выполняется по формуле С=Кс*Т (матричная операция), а их перевод в вектор массовых теплоемкостей Ст - согласно выражению: Ст=С./Мт - векторизованная операция (нотация Ма1ЬаЬ). Здесь Т - вектор степеней значений температуры (1 Т Т2 Т3)т ; Кс - матрица коэффициентов аппроксимирующих полиномов; Мт - вектор молекулярных масс.