УДК 6S1.5.G15
СТРУКТУРНО-ПАРАМЕТРИЧЕСКАЯ ИДЕНТИФИКАЦИЯ МОДЕЛИ ГАЗОВОЙ ПЕЧИ НА ОСНОВЕ ИЗМЕРЕНИЙ ЕЕ ВХОД-ВЫХОДНЫХ ПРОЦЕССОВ
М. А. Новосельцева
STRUCTURAL AND PARAMETRIC IDENTIFICATION OF THE GAS FURNACE MODEL BASED ON MEASUREMENTS OF ITS INPUT AND OUTPUT PROCESSES
M. A. Novoseltseva
В работе описывается получение модели газовой печи на основе метода структурно-параметрической идентификации по дискретным измерениям случайных входных и выходных процессов. Метод основан на теории непрерывных дробей. Метод включает в себя проверку стационарности и однородности процессов печи, получение дискретной модели печи, которая в дальнейшем может использоваться для прогнозирования, мониторинга, контроля и диагностики процессов газовой печи.
The paper describes buildinga gas furnace model basedon the method of structural and parametric identification by discrete measurements of accidental input and output processes. The method is based on the theory of continued fractions. The method includes checkingthe stationarity and homogeneity of the furnaceprocesses, receiving the discrete model of the furnace, which can be further used for predicting, monitoring, control and diagnostics of the gas furnace processes.
При изучении любых объектов, процессов, явлений основной задачей является построение их моделей. Как результат познания модель представляет собой отображение в той или иной форме свойств, закономерностей, физических и других характеристик, присущих исследуемому объекту. Любая научная или инженерная деятельность в разной степени использует формальное или содержательное описание процессов или явлений в той или иной области науки и техники. В различных научных направлениях разрабатываются свои подходы, способы и методы построения и использования модели.
Идентификация технологического процесса - это важнейший этап создания любой автоматизированной системы управления этим технологическим процессом. Для контроля и управления технологическим процессом необходимо знать, как влияет то или иное входное воздействие на выходную переменную, характеризующую его протекание. Поэтому актуальность проблем моделирования и идентификации остается одной из центральных задач кибернетики, так как эффективное управление и контроль любым технологическим процессом возможны лишь тогда, когда известно математическое описание этого процесса.
В работе рассматривается газовая печь, в которую подается воздух и метан, на выходе получается смесь газов, содержащая двуокись углерода СО2. Скорость подачи воздуха постоянна, а скорость подачи метана произвольно изменяется. Результатом наблюдений являются N = 296 пар последовательных измерений случайных процессов {х(кД^), у(кД^); к = 1, ...296}, измеренных с шагом дискретизации At = 9 с [2], где х(кД^) - расход газа на входе (в кубических фунтах в минуту) (рис. 1), у(кД^) - процентное сдержание С02 в газе на выходе печи (рис. 2). Необходимо построить
динамическую модель газовой печи на основе измерений входного процесса х(кД^) и соответствующего выходного процесса у(кД^) данной системы.
Рис. 2. Содержание С02 в газе на выходе печи
Для решения задачи воспользуемся методикой структурно-параметрической идентификации стохастических объектов на основе теории непрерывных дробей [5; 6]. Методика позволяет получить модель стохастического объекта в условиях отсутствия априорной информации о нем и его вход-выходных переменных и может быть использована при создании информационных систем управления и контроля технологическими процессами. Методика включает в себя проверку стационарности и однородности вход-выходных процессов объекта, получение дискретной модели идентифицируемого объекта, проверку адекватности полученной модели с помощью принципа вариации шага дискретизации. Приведем последовательность этапов, необходимых для осуществления данной методики.
Сх (кД )=к-к £(х(Д*)- х(+к)Д1 )2. (1)
Таким же образом производится получение структурной функции для выходного процесса, т. е. для
{у ( к М )}} выполняется аналогичная последовательность шагов, приводимая далее.
Для построения модели каждой из структурных функций используется идентифицирующая матрица В. Висковатова [6]:
(-1) - строка Г 1 1 1 . . 1 ... ^
(0)- - строка Сх(0) Сх(А1) Сх(2А1) . . Сх(пАг) ...
а 2 (0) а 2 (а>) а 2 (2А>) . 2а п >
т - строка а т (0 ) а т (А ) < . (т а А> . (т а
в которой элементы ат(пД1) последовательно определяются с помощью соотношения:
а т - 2 (п + ^) а т -1 (п + ^)
а т(пА1) = ■
где а-1 (иД) = 1(пД), а0 (пД1) = Сх (пД1), т=1,2,3,..., п=0,1,2,... Вычисление элементов ат(пД1) продолжается до появления нулевой строки. Элементы первого столбца матрицы (2) позволяют получить непрерывную дробь:
(о )* 1
Сох(г) =
ДПФ имеет нуль 7н =- 0.850 и полюс 7п= 0.290, на
основании значений которых можно сделать вывод об устойчивости модели и, следовательно, - о стационарности входного процесса.
которая после сворачивания представляет собой модель структурной функции сигнала в форме дискретной передаточной функции (ДПФ) виртуального объекта [6]. Модель (4) позволяет вынести решение, подтверждающее стационарность (устойчивость модели) или нестационарность (неустойчивость модели) процесса [6]. Если данные измерений вход-выходных случайных процессов нестационарны, то производится процедура стационаризации с помощью взятия разностей требуемого порядка [2; 6].
Графики структурных функций входного и выходного процессов приведены на рис. 3 и 4 соответственно. На основании значений Сх (кД1) была составлена идентифицирующая матрица (2) и аппроксимирована непрерывной дробью модель структурной функции входного процесса:
Рис. 3. Структурная функция Сх (кА)
Су(к)
Рис. 4. Структурная функция Су (кД)
Модель структурной функции Су (кЛ1) имеет вид:
^ 4.3627-1 + 4.1107-2
<7) = 1-0.4967-& .
Все нули и полюса находятся внутри единичного круга, следовательно, выходной процесс стационарен.
Таким образом, в результате выполнения 1 этапа методики идентификации установлено, что вход-выходные процессы являются стационарными.
делей на интервалах п и п по структуре и параметрам с некоторой заданной вероятностью можно утверждать, что однородность входного случайного процесса не нарушается. В противном случае следует продолжить применение принципа дихотомии до совпадения моделей на интервалах. Так как в процессе решения задачи 8Р-идентификации стохастического объекта следует использовать причинно-следственный принцип системного анализа, это означает, что для проверки однородности достаточно исследовать только реализацию входного процесса. Если однородность входного процесса не нарушается, то выходной процесс также является однородным. В противном случае на интервалы неоднородности разбивается реализация и входного, и выходного процессов.
Для входного процесса на интервалах п1 = [1,148] и п2 = [149,296] были получены модели структурных функций (5) и (6) соответственно:
^ () 0.7417-& + 0.6487-2 (5)
°0х1(2) =----. _______,---, (5)
ной процесс работы газовой печи тоже является однородным.
В результате применения 2 этапа методики установлено, что вход-выходные процессы являются однородными.
- на рис. 6. Интервал корреляции входного процесса х(кМ) согласно [1, 7] ткор >15М =15 9 = 135 с.
Рис. 5. Корреляционная функция Rxx(kAt)
Модели (5) и (6) однозначно совпадают по структуре, однако несколько отличаются по параметрам. Для уточнения результата проводилась дополнительная проверка на независимость остатков моделей (5) и (6) с помощью критерия серий [1]. Для уровня значимости а = 0.05 гипотеза о независимости остатков моделей принимается. Следовательно, с доверительной вероятностью 95 % можно утверждать, что входной процесс работы газовой печи однороден. Так как входной процесс является однородным, то и выход0
-0,5 --1 -1,5 --2 -2,5 --3 -3,5 J
.13 14 15-16 17 1^-19 20
Кху(к)
Рис. 6. Взаимная корреляционная функция Rxy(kAt)
(-1) - строка (0) - сторока 1 - сторока
т - сторока
Кхх(0) < х х РЇ Яхх(2А1) .. . Кхх(пА1)
(0) £ Рб Яху(А1) Яху(2А1) .. !Л у( >
а1(0) а 1 (Аі) < а1 < а1
а 0 а 2 "А а 2 > а 2 >
а т (0) а 3 . "А : < : .т а < . (т а
элементы которой последовательно рассчитываются с помощью соотношения (3), где а-1 (пМ) = Яхх (пМ),
а0(пЛ1;) = Яху(пА1), т=1,2,3,..., а п=0,1,2,... Тогда
элементы первого столбца матрицы (7) позволяют получить непрерывную дробь
нули и полюса непрерывных передаточных функций (НПФ), полученные до и после децимации, совпадают, то на основании условия SP-идентифицируемости [6] можно утверждать, что математическая модель объекта восстановлена достоверно. В противном случае необходимо уменьшить начальный шаг дискретизации технологических средств измерения и вернуться к этапу 1. Если же начальный шаг дискретизации не может быть уменьшен, то оценить истинные свойства непрерывного объекта не представляется возможным, и для оценки качества идентификации необходимо воспользоваться классическими критериями [1; 2].
После децимации модель стохастического объекта принимает вид:
^ ч -1.422 + 0.472z-1 - 0.873z -2
G(z) =-------------:---------2---.
Идентифицирующая функция имеет нули, равные zf2 = 0.166 ± 0.766i, и полюса -z|I2 = 0.511 ± 0.449i,
которые соответствуют нулям и полюсам в s-плоскости: sf2 = -0.027 ± 0.075i, s1I2 = -0.043 ± 0.040i.
На основании этих результатов можно сделать вывод, что полюса и нули в s-плоскости для At = 9 с и At = 18 с совпадают с некоторой погрешностью. Следовательно, полученная модель (10) адекватна идентифицируемому объекту и в дальнейшем может использоваться для прогнозирования, мониторинга, контроля и диагностики процессов газовой печи.
Информация об авторе:
Новосельцева Марина Александровна - кандидат технических наук, доцент математического факультета кафедры автоматизации исследований и технической кибернетики КемГУ, 8-906-926-70-00, aanov@pochta.ru.
Marina A. Novoseltseva - Candidate of Technical Science, Assistant Professor at the Department of Research Automation and Technical Cybernetics, Kemerovo State University
Статья поступила в редколлегию 30.01.2014 г.
свернув которую можно получить ДПФ объекта. В случае ненулевых начальных условий выражение (8) будет представлять собой идентифицирующую функцию [4].
Непрерывную модель объекта получают на основе взаимно однозначного отображения 8- и 7-плоскостей [6]
Поскольку Яху (0) = 1.148 ф 0, то идентифицирующая функция на основе (7) принимает вид:
^ ч -1.422 + 2.2747-1 -1.2887-2 пч
Нули идентифицирующей функции равны 7^ = 0.800 ± 0.5161, а полюса - 7^ = 0.948 ± 0.3971.
На основании (9) полюса в 7-плоскости соответствуют полюсам в 8-плоскости = -0.005 ± 0.0441. Нули
идентифицирующей функции в 8-плоскости соответствуют значениям = -0.011 ± 0.0641.
после этого осуществляется возврат к этапу 4. Если