CC BY
47
УДК 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.
Ключевые слова: газовая печь, структурно-параметрическая идентификация, непрерывная дробь, структурная функция, идентифицирующая матрица, дискретная передаточная функция.
Keywords: gas furnace, structural and parametric identification, continued fraction, structural function, identifying matrix, discrete gear function.
При изучении любых объектов, процессов, явлений основной задачей является построение их моделей. Как результат познания модель представляет собой отображение в той или иной форме свойств, закономерностей, физических и других характеристик, присущих исследуемому объекту. Любая научная или инженерная деятельность в разной степени использует формальное или содержательное описание процессов или явлений в той или иной области науки и техники. В различных научных направлениях разрабатываются свои подходы, способы и методы построения и использования модели.
Идентификация технологического процесса - это важнейший этап создания любой автоматизированной системы управления этим технологическим процессом. Для контроля и управления технологическим процессом необходимо знать, как влияет то или иное входное воздействие на выходную переменную, характеризующую его протекание. Поэтому актуальность проблем моделирования и идентификации остается одной из центральных задач кибернетики, так как эффективное управление и контроль любым технологическим процессом возможны лишь тогда, когда известно математическое описание этого процесса.
В работе рассматривается газовая печь, в которую подается воздух и метан, на выходе получается смесь газов, содержащая двуокись углерода СО2. Скорость подачи воздуха постоянна, а скорость подачи метана произвольно изменяется. Результатом наблюдений являются N = 296 пар последовательных измерений случайных процессов {х(кД^), у(кД^); к = 1, ...296}, измеренных с шагом дискретизации At = 9 с [2], где х(кД^) - расход газа на входе (в кубических фунтах в минуту) (рис. 1), у(кД^) - процентное сдержание С02 в газе на выходе печи (рис. 2). Необходимо построить
динамическую модель газовой печи на основе измерений входного процесса х(кД^) и соответствующего выходного процесса у(кД^) данной системы.
Рис. 2. Содержание С02 в газе на выходе печи
Для решения задачи воспользуемся методикой структурно-параметрической идентификации стоха-
стических объектов на основе теории непрерывных дробей [5; 6]. Методика позволяет получить модель стохастического объекта в условиях отсутствия априорной информации о нем и его вход-выходных переменных и может быть использована при создании информационных систем управления и контроля технологическими процессами. Методика включает в себя проверку стационарности и однородности вход-выходных процессов объекта, получение дискретной модели идентифицируемого объекта, проверку адекватности полученной модели с помощью принципа вариации шага дискретизации. Приведем последовательность этапов, необходимых для осуществления данной методики.
1 этап. Проверка на стационарность вход-выходных процессов. Для проверки стационарности рассчитываются значения структурной функции процесса [6]
1 К-к
Сх (кД )=к-к £(х(Д*)- х(+к)Д1 )2. (1)
1=1
Таким же образом производится получение структурной функции для выходного процесса, т. е. для
{у ( к М )}} выполняется аналогичная последовательность шагов, приводимая далее.
Для построения модели каждой из структурных функций используется идентифицирующая матрица В. Висковатова [6]:
(-1) - строка Г 1 1 1 . . 1 ... ^
(0)- - строка Сх(0) Сх(А1) Сх(2А1) . . Сх(пАг) ...
1 - строка а1(0) а1 (Аі) < . аДтА) ...
а 2 (0) а 2 (а>) а 2 (2А>) . 2а п >
т - строка а т (0 ) а т (А ) < . (т а А> . (т а
... J
в которой элементы ат(пД1) последовательно определяются с помощью соотношения:
а т - 2 (п + ^) а т -1 (п + ^)
а т(пА1) = ■
(0)
(0)
(3)
где а-1 (иД) = 1(пД), а0 (пД1) = Сх (пД1), т=1,2,3,..., п=0,1,2,... Вычисление элементов ат(пД1) продолжается до появления нулевой строки. Элементы первого столбца матрицы (2) позволяют получить непрерывную дробь:
(4)
а2
1 + 2
(о )* 1
Сох(г) =
0.7567-1 + 0.6437-2
1 - 0.2907-1
ДПФ имеет нуль 7н =- 0.850 и полюс 7п= 0.290, на
основании значений которых можно сделать вывод об устойчивости модели и, следовательно, - о стационарности входного процесса.
1 +...
которая после сворачивания представляет собой модель структурной функции сигнала в форме дискретной передаточной функции (ДПФ) виртуального объекта [6]. Модель (4) позволяет вынести решение, подтверждающее стационарность (устойчивость модели) или нестационарность (неустойчивость модели) процесса [6]. Если данные измерений вход-выходных случайных процессов нестационарны, то производится процедура стационаризации с помощью взятия разностей требуемого порядка [2; 6].
Графики структурных функций входного и выходного процессов приведены на рис. 3 и 4 соответственно. На основании значений Сх (кД1) была составлена идентифицирующая матрица (2) и аппроксимирована непрерывной дробью модель структурной функции входного процесса:
Рис. 3. Структурная функция Сх (кА)
Су(к)
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25
Рис. 4. Структурная функция Су (кД)
Модель структурной функции Су (кЛ1) имеет вид:
^ 4.3627-1 + 4.1107-2
<7) = 1-0.4967-' .
Все нули и полюса находятся внутри единичного круга, следовательно, выходной процесс стационарен.
Таким образом, в результате выполнения 1 этапа методики идентификации установлено, что вход-выходные процессы являются стационарными.
2 этап. Проверка на однородностьпо корреляционной функциивход-выходных случайных процессов. Реализация стационарного входного процесса по принципу дихотомии разбивается на 2 интервала п1 и п2 с целью получения на каждом из них модели структурной функции с помощью модифицированного алгоритма В. Висковатова [6]. В случае совпадения мо-12
делей на интервалах п и п по структуре и параметрам с некоторой заданной вероятностью можно утверждать, что однородность входного случайного процесса не нарушается. В противном случае следует продолжить применение принципа дихотомии до совпадения моделей на интервалах. Так как в процессе решения задачи 8Р-идентификации стохастического объекта следует использовать причинно-следственный принцип системного анализа, это означает, что для проверки однородности достаточно исследовать только реализацию входного процесса. Если однородность входного процесса не нарушается, то выходной процесс также является однородным. В противном случае на интервалы неоднородности разбивается реализация и входного, и выходного процессов.
Для входного процесса на интервалах п1 = [1,148] и п2 = [149,296] были получены модели структурных функций (5) и (6) соответственно:
^ () 0.7417-' + 0.6487-2 (5)
°0х1(2) =----. _______,---, (5)
1 - 0.4027
0.6767-' + 0.6567-2
ной процесс работы газовой печи тоже является однородным.
В результате применения 2 этапа методики установлено, что вход-выходные процессы являются однородными.
3 этап. Вычисление статистических характеристик вход-выходных стационарных случайных процессов. Корреляционная функция Яхх (кА>) входного процесса представлена на рис. 5, а взаимная корреляционная функция вход-выходных сигналов Яху (кА>)
- на рис. 6. Интервал корреляции входного процесса х(кМ) согласно [1, 7] ткор >15М =15 9 = 135 с.
Рис. 5. Корреляционная функция Rxx(kAt)
1 (6)
1 - 0.3727
Модели (5) и (6) однозначно совпадают по структуре, однако несколько отличаются по параметрам. Для уточнения результата проводилась дополнительная проверка на независимость остатков моделей (5) и (6) с помощью критерия серий [1]. Для уровня значимости а = 0.05 гипотеза о независимости остатков моделей принимается. Следовательно, с доверительной вероятностью 95 % можно утверждать, что входной процесс работы газовой печи однороден. Так как входной процесс является однородным, то и выход-
0
-0,5 --1 -1,5 --2 -2,5 --3 -3,5 J
1 2 3 4 5 6 7
.13 14 15-16 17 1^-19 20
Кху(к)
Рис. 6. Взаимная корреляционная функция Rxy(kAt)
4 этап. Получение модели идентифицируемого объекта. При помощи модифицированного алгоритма В. Висковатова [6] до появления нулевой строки заполняется идентифицирующая матрица:
к
(-1) - строка (0) - сторока 1 - сторока
т - сторока
Кхх(0) < х х РЇ Яхх(2А1) .. . Кхх(пА1)
(0) £ Рб Яху(А1) Яху(2А1) .. !Л у( >
а1(0) а 1 (Аі) < а1 < а1
а 0 а 2 "А а 2 > а 2 >
а т (0) а 3 . "А : < : .т а < . (т а
(7)
элементы которой последовательно рассчитываются с помощью соотношения (3), где а-1 (пМ) = Яхх (пМ),
а0(пЛ1;) = Яху(пА1), т=1,2,3,..., а п=0,1,2,... Тогда
элементы первого столбца матрицы (7) позволяют получить непрерывную дробь
1 + -
1+-
а
(0);
1 + ...
(8)
нули и полюса непрерывных передаточных функций (НПФ), полученные до и после децимации, совпадают, то на основании условия SP-идентифицируемости [6] можно утверждать, что математическая модель объекта восстановлена достоверно. В противном случае необходимо уменьшить начальный шаг дискретизации технологических средств измерения и вернуться к этапу 1. Если же начальный шаг дискретизации не может быть уменьшен, то оценить истинные свойства непрерывного объекта не представляется возможным, и для оценки качества идентификации необходимо воспользоваться классическими критериями [1; 2].
После децимации модель стохастического объекта принимает вид:
^ ч -1.422 + 0.472z-1 - 0.873z -2
G(z) =-------------:---------2---.
1 - 1.022z -1 + 0.463z-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) адекватна идентифицируемому объекту и в дальнейшем может использоваться для прогнозирования, мониторинга, контроля и диагностики процессов газовой печи.
Литература
1. Бендат, Дж. Прикладной анализ случайных процессов / Дж. Бендат, А. Пирсол. - М.: Мир, 1989. - 464 с.
2. Бокс, Дж. Анализ временных рядов. Прогноз и управление / Дж. Бокс, Г. Дженкинс. - М.: Мир, 1974. -Вып. 1. - 406 с.
3. Дейч, А. М. Методы идентификации динамических объектов / А. М. Дейч. - М.: Наука, 1985. - 240 с.
4. Идентифицирующая функция и ее свойства / В. Я. Карташов [и др.] // Сборник научных трудов VII Всероссийской научно-практической конференции «Системы автоматизации в образовании, науке и производстве -AS2009». - Новокузнецк: СибГИУ, 2009. - С. 399 - 405.
5. Карташов, В. Я. Непрерывные дроби (определения и свойства) / В. Я. Карташов. - Кемерово: Изд-во Кемеровского госуниверситета, 1999. - 88 с.
6. Карташов, В. Я. Идентификация стохастических объектов: учебное пособие / В. Я. Карташов, М. А. Новосельцева. - Кемерово, 2010. - 108 с.
7. Романенко, А. Ф. Вопросы прикладного анализа случайных процессов / А. Ф. Романенко, Г. А. Сергеев. -М.: Советское радио, 1968. - 247 с.
Информация об авторе:
Новосельцева Марина Александровна - кандидат технических наук, доцент математического факультета кафедры автоматизации исследований и технической кибернетики КемГУ, 8-906-926-70-00, [email protected].
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]
8 = —ЫЫ + ^-1аГ27 . (9)
А1 А1
Поскольку Яху (0) = 1.148 ф 0, то идентифицирующая функция на основе (7) принимает вид:
^ ч -1.422 + 2.2747-1 -1.2887-2 пч
0(7) =-------------:----------2--.(10)
1 -1.8957 -1 +1.0557-2
Нули идентифицирующей функции равны 7^ = 0.800 ± 0.5161, а полюса - 7^ = 0.948 ± 0.3971.
На основании (9) полюса в 7-плоскости соответствуют полюсам в 8-плоскости = -0.005 ± 0.0441. Нули
идентифицирующей функции в 8-плоскости соответствуют значениям = -0.011 ± 0.0641.
5 этап. Реализация принципа вариации шага дискретизации. Принцип вариации шага дискретизации является достаточным условием взаимно однозначного соответствия между непрерывным объектом и дискретными моделями [6]. Принцип вариации определяет качество идентификации и осуществляется в форме децимации значений Яхх(кА1) и Я (кА1) [6],
после этого осуществляется возврат к этапу 4. Если
