Идентификация объектов, описываемых линейными разностными и дифференциальными уравнениями в форме Коши с вещественным аргументом

Абденова, Гаухар Амирзаевна

Системный анализ, управление и обработка информации (по отраслям)

автореферат диссертации по информатике, вычислительной технике и управлению, 05.13.01, диссертация на тему:Идентификация объектов, описываемых линейными разностными и дифференциальными уравнениями в форме Коши с вещественным аргументом

кандидата технических наук: Абденова, Гаухар Амирзаевна
город: Новосибирск
год: 2012
специальность ВАК РФ: 05.13.01
цена: 450 рублей

Диссертация по информатике, вычислительной технике и управлению на тему «Идентификация объектов, описываемых линейными разностными и дифференциальными уравнениями в форме Коши с вещественным аргументом»

Автореферат диссертации по теме "Идентификация объектов, описываемых линейными разностными и дифференциальными уравнениями в форме Коши с вещественным аргументом"

На правах рукописи

АБДЕНОВА Гаухар Амирзаевна

ИДЕНТИФИКАЦИЯ ОБЪЕКТОВ, ОПИСЫВАЕМЫХ ЛИНЕЙНЫМИ РАЗНОСТНЫМИ И ДИФФЕРЕНЦИАЛЬНЫМИ УРАВНЕНИЯМИ В ФОРМЕ КОШИ С ВЕЩЕСТВЕННЫМ АРГУМЕНТОМ

Специальность 05.13.01 - Системный анализ, управление и обработка

информации (промышленность)

005049498

АВТОРЕФЕРАТ

диссертации на соискание ученой степени кандидата технических наук

7 ФЕВ 2013

НОВОСИБИРСК-2013

005049498

Работа выполнена в Федеральном государственном бюджетном образовательном учреждении высшего профессионального образования «Новосибирский государственный технический университет»

Научный руководитель: Воевода Александр Александрович

доктор технических наук, профессор

Официальные оппоненты: Куцый Николай Николаевич,

доктор технических наук, профессор, Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования «Национальный исследовательский Иркутский государственный технический университет», кафедра автоматизированных систем, профессор

Толстиков Александр Сергеевич, доктор технических наук, доцент, Федеральное государственное унитарное предприятие «Сибирский государственный ордена Трудового Красного Знамени научно-исследовательский институт метрологии», отдел государственной службы времени, частоты и определения параметров вращения Земли, начальник

Ведущая организация: Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования «Томский государственный университет систем управления и радиоэлектроники»

Защита состоится «19» февраля 2013 г. в 10.00 часов на заседании диссертационного совета Д 212.173.05 при ФГБОУ ВПО «Новосибирский государственный технический университет» по адресу: 630092, г. Новосибирск, пр. К. Маркса, д. 20, корпус 1, конференц-зал.

С диссертацией можно ознакомиться в библиотеке Новосибирского государственного технического университета.

Автореферат разослан « » января 2013 г.

Ученый секретарь диссертационного совета

Шпилевая Ольга Яковлевна

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность темы исследования. На практике многие объекты в большинстве своем являются нестабильными, часто имеют малое количество наблюдений1, характеризуются наличием ряда циклов или этапов развития, на каждом из которых достаточно сильно изменяются тенденции поведения. Математически поведение этих объектов можно описывать разными видами моделей.

На сегодняшний день методика построения моделей достаточно хорошо проработана для стационарных объектов, менее - для нестационарных объектов, зависящих от временного параметра, и не рассмотрена вовсе в случае объектов с явно выраженными изменяющимися тенденциями в поведении, входные и выходные данные которых зависят от невременного аргумента2 и содержат в себе малое количество наблюдений. Проблемы возникают при построении модели на участке, содержащий изменения тенденций поведения. Наличие изменений тенденции поведения на всем интервале наблюдения приводит к построению ошибочной модели. Поэтому в процессе построения модели очень важно выделять участки данных наблюдений, на которых модель строится на основе данных наблюдений с неизменяющимися статистическими свойствами. Объединение локальных моделей и построение при малом объеме выборки единой адекватной модели на всем участке наблюдения, где переменные состояния зависят не обязательно от временного (вещественного) аргумента, является актуальной задачей.

Теоретико-методологической основой исследования по вопросам параметрической идентификации, по оцениванию состояния и управления в технических и других системах послужил обзор концепций, сформулированных в работах зарубежных теоретиков и практиков, авторами которых являются: П. Эйкхофф, К. Острем, Р. Мехра, Р. Бен Мрад, Е. Фараг, А. Кисикайя, X. Кай-ран Ахмет, JI.C. Льюнг, P.JL Кашьяп, А.Р. Pao, М.С. Гревел и др. Эти задачи решались и в трудах отечественных ученых: A.C. Анисимова, Ю.Е. Воскобойни-кова, В.Г. Горского, Н.Д. Егупова, A.A. Ломова, В.Н. Овчаренко, М.А. Огар-кова, В.Н. Подладчикова, И.Н. Синицына, В.А. Фурсова и др.

Целью диссертационной работы является повышение эффективности оценивания состояния объектов с малым объемом наблюдения и с явно выраженными изменяющимися тенденциями в поведении путем идентификации объекта математическими моделями в форме Коши с вещественным аргументом.

Задачи исследования. Для достижения заданной цели поставлены следующие задачи исследования.

1. При малом объеме наблюдений разработать методику построения линейной кусочно-разностной и кусочно-дифференциальной стохастической модели

1 Под малым количеством наблюдений в данной работе понимается, такое количество наблюдений, которое превышает число неизвестных параметров модели динамики в 3-5 раз.

2 Пример: результаты выхода некоторого сорта стали в зависимости от угара кремния, измеряемого в процентах.

в форме Коши, где переменные входа, состояния и выхода, а также переменные параметры в модели зависят не обязательно от временного аргумента.

2. Рассмотреть объект с нестабильной динамикой развития, состояние которого зависит от невременного аргумента и ранее описываемый линейными и нелинейными регрессионными зависимостями с целью перехода к описанию с помощью разностной модели в форме Коши.

3. Рассмотреть случай, преобразования модели с распределенными параметрами в форме Коши к стохастической линейной модели в форме Коши сосредоточенного типа.

4. Для малого количества фиксированных данных наблюдений разработать алгоритм расчета шумовых характеристик одномерной модели в форме Коши без учета оценок параметров модели динамики.

5. Для повышения точности оценок параметров модели динамики в непрерывно-дискретной модели в форме Коши необходимо исследовать вопрос уменьшения числа обусловленности информационной матрицы Грама с помощью процедуры масштабирования входных и выходных данных.

Объектом исследования являются объекты, с изменяющимися тенденциями в поведении, при малом объеме наблюдений и функционирующие в условиях неопределенности.

Предметом исследования выступают аналитический и численный анализы некоторых типов моделей, их модификации с целью перехода к описанию поведения исследуемого объекта линейной стохастической моделью в форме Коши.

Методы исследования. Поставленные в работе задачи решены с использованием методов и положений математического анализа, дифференциального исчисления, теории фильтрации, теории математической статистики, методов вычислительной математики.

Научная новизна работы и результаты, выносимые на защиту.

1. Алгоритм построения стохастической линейной модели в форме Коши для линейных и нелинейных регрессионных уравнений в случае зашумленных входных и выходных переменных.

2. Алгоритм преобразования модели объекта, описываемого стохастическими моделями распределенного типа к адекватной стохастической модели в форме Коши сосредоточенного типа.

3. Алгоритм уменьшения числа обусловленности информационной матрицы Грама на основе масштабирования исходных входных и выходных данных наблюдений в линейных непрерывно-дискретных моделях в форме Коши.

4. Алгоритм рекуррентного оценивания дисперсий шумов в моделях динамики объекта и измерительной системы на основе дискретных наблюдений.

5. Методика построения единой линейной кусочно-разностной и кусочно-дифференциальной модели в форме Коши для решений задач тестового и практического характеров.

Теоретическая значимость. Разработанные в диссертации методика и алгоритмы позволяют при малых выборках построить в виде одной линейной стохастической кусочно-разностной и/или кусочно-дифференциальной модели в форме Коши, в которых переменные параметры, состояния объекта и выход измерительной системы зависят от невременного аргумента.

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

Методика построения стохастической дифференциальной модели в форме Коши внедрена для расчета индуктивных и емкостных параметров элементов аттенюатора в ФГУП «РТРС», филиал «Липецкий областной радиотелевизионный передающий центр» (г. Липецк, РФ).

Методика построения стохастической кусочно-дифференциальной модели внедрена с целью использования полученной модели для расчета количества прогнозных поставок зернового сырья в мукомольно-комби-кормовые комбинаты при решении задач транспортной логистики. Методика апробирована на реальных данных объема поставок сельскохозяйственной продукции в течение периода с 2001 по 2004 год в ОАО «Восточно-Казахстанский мукомольно-комбикормовый комбинат» (г. Семей, РК).

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

Следует отметить, что диссертационная работа поддержана грантом по заданию Министерства образования и науки РФ, по проекту № 7.559.2011, «Исследование предельных точностей оптических методов измерения параметров движения и мехатронных методов управления движением и разработка новых робототехнических и электромеханических систем», Темплан, гос. per. номер НИР 01201255056.

(Томск, 2004); III Международной научно-практической конференции аспирантов и молодых ученых «Страны СНГ в условиях глобализации» (Москва, РУДН, 2004); VI Среднеевропейской научно-технической конференции «Компьютерные методы и системы в автоматике и электротехнике» (Ченстохова, Польша, 2005); Международной научно-технической конференции «ИКИ-2001», «Измерение, контроль, информатизация» (Барнаул, 2001, 2003, 2005, 2007, 2009, 2010); II Российско-индийском объединенном семинаре «Computational Intelligence and Modern Heuristics in Automation and Robotics», 10-13 сентября 2011 г.; XI Международной конференции «Актуальные проблемы электронного приборостроения», Новосибирск, 2-4 октября 2012 г.

Публикации. По теме диссертационной работы опубликовано 15 печатных работ, в том числе 6 работ в журналах из Перечня ведущих рецензируемых научных журналов и изданий, в которых опубликованы основные научные результаты диссертаций на соискание ученой степени кандидата наук, 6 статей в сборниках научных трудов и 3 работы в материалах международных конференций.

Личный вклад автора. Все результаты, составляющие основное содержание диссертации, получены автором лично. Из 15 опубликованных статей 8 работ написаны в соавторстве, в них автором сформулированы задачи и получены теоретические результаты; в работе, связанной с процедурой масштабирования, задача поставлена научным руководителем, а алгоритм решения и вычислительные расчеты получены автором самостоятельно.

Обоснованность и достоверность научных положений и выводов подтверждаются результатами теоретических, экспериментальных исследований и расчетов, публикациями в рецензируемых изданиях, апробацией на международных конференциях, а также использованием результатов диссертационной работы при решении задач производственного характера.

Структура и объем работы. Диссертация состоит из введения, четырех глав, заключения, списка литературы из 137 наименований и восьми приложений. Общий объем работы составляет 190 страниц, включая 33 таблицы и 36 рисунков. Основное содержание работы изложено на 139 страницах.

КРАТКОЕ СОДЕРЖАНИЕ РАБОТЫ

Во введении обоснована актуальность темы диссертации, сформулированы цель и задачи исследования, определены научная новизна и практическая ценность работы, дана характеристика полученных результатов.

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

Коши с целью уменьшения числа обусловленности матрицы Грама.

В главе дан краткий теоретический материал по инструментариям, которые применяются в исследовании вышеперечисленных задач. Так, для построенной стохастической модели в форме Коши используется фильтр Калмана. Известно, что этот аппарат содержит в себе много трудностей, связанных с наиболее точным оцениванием параметров уравнений динамики, характеристик шумов измерительной системы, начального состояния и динамики исследуемого объекта. Чтобы избежать проблем с расходимостью при калмановской оценке состояния, требуется, чтобы построенная модель прежде всего удовлетворяла четырем свойствам: устойчивости, наблюдаемости, управляемости, идентифицируемости. Поэтому в работе приведены известные критерии проверок этих свойств. Для применения методов параметрической идентификации, когда в поведении объекта происходят изменения тенденции, весь интервал наблюдения за объектом при малом объеме выборки визуально разделяется на подынтервалы, где объект не испытывает изменения тенденции в поведении. В связи с этим в главе изложен критерий Грегори Чоу, с помощью которого проверяются границы подынтервалов, где происходят изменения тенденции в поведении объекта.

К описанию динамических систем моделями в форме Коши исследователи прибегают все чаще. В данной работе под динамической системой будем понимать любой объект, для которого однозначно определено понятие состояния как совокупности некоторых величин, зависящих не только от времени, но и от любого другого независимого аргумента и задан закон, который описывает изменение начального состояния в зависимости от значений другого независимого аргумента. Поэтому такие ключевые понятия параметрической идентификации, как «модели в пространстве состояний», «стационарные», «нестационарные модели» и т. д., в данной работе рассматриваются как понятия, связанные с переменными, которые зависят необязательно от временного (вещественного) аргумента и заменяются терминами «модели в форме Коши», «модели с постоянными и переменными коэффициентами».

В модели в форме Коши динамика системы описывается стохастическими конечно-разностными или дифференциальными уравнениями, а измерительная система - в форме стохастических дискретных алгебраических уравнений. Исходя из этого, линейная система на подынтервалах может быть описана линейной стохастической дискретной моделью с постоянными параметрами в зависимости от вещественного аргумента к в стандартной форме вида (1),(2): х(£ + 1) = Г-х(£) + С-и(/с) + Ь-'№(/с); х(0) = х0, к = О.ЛГ-1, (1)

у(к +1) = И • х(к +1) + \(к +1), * = 0ГЛГЧ (2)

где у (к), \(к), и (к), лу (к), \(к) - соответственно векторы выхода, переменных состояния, управляющих входных сигналов, шума динамики объекта, ошибок измерений (наблюдений).

Линейная непрерывно-дискретная система на подынтервалах может быть описана моделью с постоянными коэффициентами в уравнениях в зависимости от вещественного аргумента I в виде:

-^р- = А • х(/) + В • и(/) + С • лу(/); х(/0) = х0, (3)

у(4) = Н-х(4) + у(4), * = 1,1ке[10,1Л (4)

... с!х(1)

где х(/) - и-вектор состояния;--«-вектор производных от переменных

сИ

состояния; и(/) — г- вектор входного сигнала; \у(/) - вектор белого гауссовского процесса с нулевым математическим ожиданием и ковариационной матрицей (3 размера х0 - вектор начального состояния; А, В, С - постоянные матрицы размеров (пхп), (п х г) и (п х ц) соответственно; /- независимый параметр, [/0, ] - весь интервал наблюдения, на протяжении которого ведется наблюдение за объектом; у(1к) - известный вектор измерений (наблюдений) размера т; Н - матрица наблюдения размера (шхп); {\(1к), к = \, Щ — белая гауссовская последовательность с нулевым математическим ожиданием и ковариационной матрицей Я размера (ихт); N - количество дискретных моментов наблюдения.

Дифференциальные уравнения в частных производных на практике нередко используются для описания систем, изменяющихся во времени и в пространстве. Например, при описании процесса теплопереноса используется уравнение в частных производных параболического типа. Пространственно-временная модель в форме Коши описывается в виде:

д2д(х, ?) и д^ф^) дд(х, г) , дд(х, I) 1 „

а,-;--Но,--Не,--на,--не,-а(х, Л +

1 дх2 дхд( 1 дх 1 Ы 1 (5)

+/,-и(х, 0 + р1 • и(х, /) = 0,

г(хк,^) = 11-д(хк,^) + £(хк,^), к = 1,п; я = 0 ,т, (6)

где а,, 6,, с,, с!\, е,, р1 - постоянные коэффициенты; х, ? - пространственная координата и время (0<х<Ь, />0) соответственно; ц{х,1) - функция состояния системы с распределенными параметрами с некоторыми граничными и начальными условиями; и(х,1) - входное управляющее воздействие распределенного типа, которое удовлетворяет определенным амплитудным ограничениям мт1п < < итах; - белый гауссовский процесс распределенного типа с нулевым математическим ожиданием и известной дисперсией (). Для (5) граничные условия имеют вид д(0,О) = д0; ц(Ь, /у) = £/£. В (6) г(хк, ) — дискретный распределенный выход измерительной системы, где индексы к и 5 означают, что пространственно-временная функция состояния может измеряться только в дискретных пространственных точках хк и в дискретные моменты времени ; И - заданный весовой коэффициент; £(хк,15) - белый гауссовский шум измерительной системы распределенного типа с нулевым математическим ожиданием и известной дисперсией Я. При выше

перечисленных условиях ставится задача оценивания коэффициентов уравнения (5) на основе входа и выхода системы.

Таким образом, в первой главе поставлена задача разработки методики построения математических моделей в форме Коши для объекта со следующими ограничениями и характеристиками: предположение гладкости функции входа и выхода системы, малое количество наблюдений на выходе измерительной системы, наличие явно выраженных изменений тенденций в поведении объекта и зависимость входных данных, переменных состояния, выходных данных и переменных параметров от вещественного аргумента. В ходе разработки методики были рассмотрены следующие вопросы: постановка задачи перехода от линейных и нелинейных регрессионных уравнений и модели систем распределенного типа к линейной модели в форме Коши сосредоточенного типа; описание критериев проверки свойств модели в форме Коши, таких как устойчивость, наблюдаемость, управляемость и идентифицируемость; визуальное разделение данных наблюдений на подынтервалы с подтверждением этого разбиения с помощью критерия Г. Чоу; постановка задачи исследования числа обусловленности информационной матрицы Грама для повышения точности оценок параметров модели.

Во второй главе излагаются алгоритмы перехода от линейных и нелинейных регрессионных уравнений и моделей распределенного типа к уравнениям в виде стохастических моделей в форме Коши сосредоточенного типа.

Пусть линейная регрессионная модель представлена в виде, когда зашумлены вход и выход:

></у) = о,+6,-•*(/,) + *((,), У(1])=Х(11) + Б(1]), у = й, [/,_„/,], 1 = где 3(1]) и £(/.) - случайные величины, характеризующие ошибки значений входа и выхода соответственно; х(/у) — истинные, но наблюдаемые с шумами значения входа системы. Предполагается, что шумы е(1}) и д(1]) распределены нормально с нулевым средним и дисперсиями ег2 (>',.) и <т2(у1), соответственно, и коэффициентом корреляции р, = 0. Для простоты опустим все индексы, соответствующие подынтервалам.

Оценки параметров а и Ь определяются на основе МНК, а истинные значения входа х, из условия минимума функционала:

Р = 1 /2■ ти(V/ - (г •+ Ь ■ х, I)2 /<т2(У,) + (у,- -Х1у/а2(V,.)].

Такой подход не всегда приемлем на практике из-за смещенности МНК-оценок параметров модели. В связи с этим в работе предлагается более эффективная модель в форме Коши.

Введем следующие обозначения: через *,(£) и х2(к)обозначим выходные и входные переменные объекта в момент значения невременного аргумента к соответственно; через у, (к) — наблюдение за состоянием х, (к), а через у2 (к) -наблюдение за состоянием х2(к).

С помощью процедуры регуляризирующего сплайна значения данных наблюдений {у,(£ + 1), к = 0,п-\} и {у2(к +1), к = 0,п-\} преобразовываются в последовательности {х,(А: + 1), к = 0, п-\) и {х2{к +1), к = 0,и-1} соответственно. Подставляя эти значения в переменные модели динамики объекта с учетом аддитивных помех, получим следующую систему двух конечно-разностных уравнений:

Гх, (к +1) = /„ • (к) + /12 • х2 (к) + с, ■ (к),

[х2(А: + 1) = /21 ■х1{к) + /22 ■х2(к) + с2-\02(к), к = 0,п-1.

С помощью МНК можно получить следующие оценки параметров переходной матрицы Р,где /п ~}\2, /21 ~/21, /22 ~)п-

С учетом полученных оценок можно записать соотношения вида (1), позволяющие осуществлять прогнозирование и фильтрацию значений переменных состояния с помощью аппарата дискретного фильтра Калмана:

1х,(А: + 1) = /11 -хх(к) + /п-х2{к) + сх -м>х(к),

[х2(к +1) = /21 -х^к) + /22 • х2(к) + с2 -у^2(к), к = 0,п — \.

С учетом измерительной системы и некоторых обозначений получим: х(к +1) = Б" ■ х(к) + С • \\(к), х(0) = х0, у(А: + 1) = Н-х(А: + 1) + у(А: + 1), к = 0,п-1,

где Ё- переходная матрица, С—диагональная матрица, Н-матрица наблюдений.

Методика построения модели в форме Коши может быть применена и на случай нелинейной регрессионной модели. В работе рассмотрен частный случай нелинейной регрессионной модели вида

Г = А-Ха-Г'-а, (7)

где А, а — некоторые положительные постоянные, 0 <а < 1.

Введем обозначения к = Х/У, /(к) = ПУ, которые позволят записать (7)

в виде /(к) = А-к", где невременной параметр к - это вход, функция /(к) — выход исследуемого объекта. Наблюдения за объектом дают возможность записать две последовательности данных, соответствующие входу и выходу: {*„,*„..., ¿V-, Дл'Ь {/(к0), /(к,),...,/(*„_,), /(*„)}.

Поставим задачу исследования скорости изменения функции /(к) в зависимости от входа к. Представим скорость изменения значений выхода в виде следующей линейной модели в форме Коши непрерывно-дискретного типа: х{к) = а-х{к) + Ъ + ™(к), х(к0) = х0 (а); /(к:)=х(к1) + у(к!), 1 = (б), (8) где к - непрерывное значение входа; х(к) — непрерывное значение состояния исследуемого объекта; х(к) — скорость изменения функции состояния; х(к0) = х0 — начальное значение входа при к0 и являющейся гауссовской величиной с математическим ожиданием х0 и дисперсией Р(к0 ) = Р0; \\>(к) — белый гауссовский шум с нулевым математическим ожиданием и дисперсией 2; к1 — дискретные значения входа; /(£,) — дискретные значения выхода измерительной системы; х{к1) — истинные значения состояния соответствующего

входа у(£,) - белая гауссовская последовательность с математическим ожиданием д и дисперсией Я; а, Ъ - неизвестные параметры; N - объем выборки.

Алгоритм построения модели в форме Коши состоит из следующих шагов: Шаг 1. Предположим {х(к1) я /(£,), «= 1,.ЛГ}.

Шаг 2. Обозначим среднюю скорость изменения значений выхода иссле-

Ах А/

дуемого объекта через Цср= — ~— на промежутке [£;; к) + Ак], где

Ак = км-к!, А/ == / (к1 + Ак) - / (/с;), / = 1, N-1. Тогда соотношение (8 б) с учетом обозначений шага 2 и предположения из шага 1 запишется в виде:

отношение в виде х(км) = (1 + а ■ (км - /с,)) • х(к1 ) + Ь- (км - к ¡) + Мк1) • (км - к ¡).

Шаг 3. Вводятся обозначения: а, = 1 + а ■ (км — к¡), и{к1) = км - к1, £(к1) = У1>(к1)-(км-к1), ¡ = 1, N-1. Тогда последнее соотношение в шаге 2 запишется в виде х(км) = а, • х{к) ) + Ь- и{к1) + е(к1), г = 1, N -1, х(ка) = х0. Используя МНК, можно получить соотношения для оценивания неизвестных параметров ах,Ь.

Шаг 4. Рассчитываются дисперсии ошибок динамики объекта Q, ошибок измерительной системы Я и ошибки начального состояния Р(к0) = Р(].

Шаг 5. Оценивается состояние объекта на основе уравнений дискретного фильтра Калмана, который позволяет получить фильтрационные оценки х(к11 к1) относительно переменной состояния {х(к1), / = И, N +1,...,}.

Таким образом, функция (7) преобразована в линейную модель, для которой предложен алгоритм построения линейной модели в форме Коши.

Далее предлагается алгоритм преобразования модели с распределенными параметрами в адекватную модель в форме Коши во временной области.

Алгоритм состоит из следующих шагов.

Шаг 1. Примем интервалы квантования для ¿их равными Д/ и Ах. Тогда из (5) получим соотношение в частных разностях:

Шаг 2. Проведем группировку значений функций д(х,1) по величинам соответствующих координат к и 5 узлов сетки и переобозначим коэффициенты при дискретных значениях функции д(х,1). Тогда выражение (9) в конечных разностях запишется в виде

а'Чк+\,з +ь-Чк+и+1 +с'Як,5 1 +е'Чк-1,5 +Р-Мк,5 =0- 0°)

х(км)-х(к1)

, что позволит записать последнее со-

(9)

Шаг 3. Запишем систему уравнений, получаемых из (10) для к = 1,2,затем умножим все соотношения этой системы на (-1), а затем слагаемые, содержащие индексы +1), оставим в левой части, а остальные слагаемые с индексами 5 перенесем в правую часть. Тогда получим

..............................(П)

1 ~ЪЧп+1,5+1 =еЧп-1,5 + +«^+1,5 + У»Л,5 = 0'

Граничные условия для (11) заданы в виде 0■Ах = 0, (я +1) • Ах = £.

Шаг 4. Вводя обозначения, (11) можно записать в виде

К.Х,+1=М-Х,+С-и,+В.ц,+(Ь2-01+1+а2-ъ+е2-0Л

Умножим слева соотношение (12) на й-1, получим

Х„, = Ы"1 М X, + И"1 с; и5 + И"1 Б + Я"1 X

Шаг 5. Соотношение (13) требует введения дополнительных обозначений и преобразований. Пусть для начала рл =(е2 ■ #л[0ихи|а2 ■ ?/;|Ь2 • ) - матрица размером (я х (и + 3)). Теперь и-вектор управления и, увеличим на три

компоненты, введя обозначения: Ф = И-1 • М; Г = I*"1 - Б; С = И"1 • С.

Все последние матрицы имеют размерность их и. Первую компоненту обновленного вектора и? заменим на 1, а матрица С увеличится одним столбцом матрицы рл. и дополнительными (и + 2)-м и (п + 3)-м столбцами р5: С = (е2|к-'-С|а2|Ь2); и,=(1 и, 7, в5)т.

Окончательно выражение (13) запишется в виде:

(12)

(13)

^X,+GuJ+D|iI.

(14)

Временной параметр 5 принимает значения 5 = 0,1,2,..., а начальные условия

хо=(?1,о 92,о ■•• Ч„,о Г- (15)

Контролируемые входные воздействия будут иметь размерность вектора (и + 3)х1 и 5 = 0 будет иметь форму = {в1 м10 «2,о •••"«,<) т]5 вг+1У, а неконтролируемое входное воздействие ц? будет иметь вид: ц., ~{м\,з Мг^ ■■■М„,!)Т■ После введенных обозначений матрица С преобразуется к структуре вида

G =

а с

0 а

а с

пх(н+3)

а матрица Г = Е, где Е - единичная матрица.

Шаг 6. Введем обозначения ys. = z, Х( = q и для уравнения (6) при условии h = const. Тогда соотношение (6) можно записать в векторной форме:

У,=й-Х,+Е,, 5 = 0,1,2,... . (16)

Уравнение (14) с начальным условием (15), а также уравнение (16) -модель стохастической системы в форме Коши во временной области, соответствующая пространственно-временному описанию в форме Коши (5), (6).

Для повышения точности параметров модели в форме Коши в работе предлагается осуществлять анализ числа обусловленности в информационной матрице Грама Б = Хг • X. Поэтому ставится задача исследования влияния масштабирования вход-выходных переменных на уменьшение числа обусловленности этой матрицы на частном примере. Пусть для компактности объект описывается в виде (3), (4) при Ж=11, с точностью до неизвестных параметров, входящих во второе и третье уравнения:

х1(1) = х2(1) + м>](1), *,(/„) = х1Н,

■ х2{1) = а1ххх{1) + агъхъ{1) + Ъ21щ{1) + м!2{1), х2(10) = х20, (17)

х3 (/) = а31х, (/) + а33х} (/) + ЬЪ2и2 (/) + и'3 (/), х3 (/„ ) = х30.

Ь2(/,) = х3(/,)^/Д / = ТДо.

Предполагается, что известны: начальные условия динамики поведения объекта х,(/0), х2(10), х3(/0); входные данные {и,(/,.), и2(/;) / = 0,Щ и выходные данные {у! (/,) у2(/, ) / = 0, Ы}. В этих условиях требуется рассчитать числа обусловленности для матриц Грама относительно входных и выходных наблюдений без масштабирования и после процедуры масштабирования. Решение задачи проводиться по следующему алгоритму:

Шаг 1. Проведем сглаживание последовательности данных выхода измерительной системы >>,(/,)={0.0572 ... 2.0242} и _у2(/,■) = {2.5072 ... 11.925} с помощью регуляризирующего сплайна с равномерной сеткой /г, =/,+1 -/, =1, / = 0,N -1 и коридором фильтрации ? = {0.011,..., 0.011}. Рекомендуется сделать не менее 10-и итераций процедуры сглаживания. В результате получаем следующие значения: у1(/,.)={0.046 ... 2.03}; у2(/,)={0.77 ... 11.93}.

Шаг 2. Приравняем последовательности {х1(/,) = х1 =у,(/,.), х3(/,) = х3 =

= у2 (/,), г =0Л0}, получим х, ={0.046 0.101... 2.029}; х3 = {0.77 0.453... 11.932};

Шаг 3. Проведем процедуру численного дифференцирования значений последовательности {х, =х, (/,.), г = 0, Щ кубическим сплайном вида

ад = Л 0 - О + Л+1 • I - ^ I ■ (1 - / )[(2 ■- 1)Мк + (1 + 1)МЫ ],

где / е [/,., 1М ], / = 0, ТУ — 1, М,. = ), г = 0, N. Вычислим производную этого соотношения. При этом для левых границ отрезка [/,., /,+1], 1 = 0 и производные будут вычисляться по формуле

а для правой границы отрезка [/ЛГ_1, 1М~\, по формуле

5'(/,) = /,+1 -/, -^(2-6./ + 3-/2)-М; -1(2-3-/2)-Л/,.+1. (20)

Используя формулы (19) и (20) получим последовательность значений производных х, = х2=[0.0157 0.1019 ... -0.4123] , которая позволила восстановить все значения переменных х,, х2, х3;

Шаг 4. На основе данных входа и,,и2 и выхода значений переменных х,, х2> х3 необходимо построить матрицу Грама Б = ХГХ, где матрица X размера 30 х 15, имеет следующую структуру:

1)0 0 ^ х2(/1) х3(/,) и (/.) и2{1{) 0 ООО 00000 0

о х2(/,) 0 0 0 0 0 0 X, х2 х3 щ и2 0 0 0 0 0

0 0 х3(/,) 0 0 0 0 0 0 0 0 0 0 Х[ х2 х3 щ и2

х,(/10) 0 0 х2(1ю)х3{1ю)и Сю) и2(/10) 0 0 0 ООО 0000

0 х2(/10)0 0 0 0 0 0 Х[ х2 х3 щ и2 0 0 ООО

^ 0 0*3(/10), , 0 0 0 0 0 0 0 0 0 0 X! х2 х3 щ и2

01 0000 00 000 00 00 00 0 0 0 а210 а23 й2100 00 0 0 00 0000 00 0 0а310а330й3^

V.) 0 0 Л

Х1(/1)= 0 ¿2(/,) 0 ; ......*,<,(/, о):

о о ¿з(/,)

Введем следующие обозначения:

*,(/,„) 0 0 0 ¿2(Ао) 0

о о ¿з(/10)

(*10(/>о))г )г;

Х = ((х, (/,))' (*2(/2))

'*,</,) х2(/,) х3(/,) «/,(/,)«,(/,) 0 0 0 0 0 0 0 0 0 0 ^

ХЦ= 0 0 0 0 0 х1(/|)^(/,)х3(/,)",(/,) «2(А) 0 0 0 0 0 0 0000 0000 0 *,(/,) х2(/,) х3(/,) «,(/,) м2(/,Х

Vй

ГЗДо) ^С.о) *з(Ло> "1('ш)«2(Ао) 000000000 О"

хи,0 =

х2 х3 м, и2 0 0 0 0 0

0 0 0 0 о

О О О О 0 0 0 0 0 0 х

2 У

хи = (хиг хи[ хи[ ... хи[0)г; в, = (0 1 0 0 0 0 0 0 0 0 0 0 0 0 0),

е2 = (00000 а210 а23Ь210 0 0 0 0 0), 03 = (0000000000а310а33 0Ь32),

в=(в[ в2 03 ) • 14

Применив все обозначения, запишем соотношение Х = Хи-0, для которого с помощью МНК получим

т т

Л Л Л

е = (хи -хиг-хи -х.

(21)

Шаг 5. Из соотношения (21) для матрицы Грама Б = Хи • XII, используя формулу АТф) = /!.тах//1тт, вычислим число обусловленности АГ(О) я 2.95 • 105;

Шаг 6. Для уменьшения числа обусловленности матрицы Грама, проведем процедуру масштабирования входных и выходных данных с помощью масштабов вида:

л /

Му =

ту1 О

у2)

10 о

м„ =

тих 0

2 У

10 о

Запишем (17) и (18) с использованием матриц масштабирования Мц, Му Х(/) = А • Х(/) + В • й(/) + !»(/), Х(/0 ) = Х0, у(/,) = Н.Х(/,) + у(/,), / = й7,

(22) (23)

где й(/) = Мц •!!(/), В =В-М~', у(/,.) = Му -у(/,), Н = Му-Н, у(/,.) = М -Щ),

М"1 =

тух 0

\-1

У 2

^1/10 о' 0 1

м:1 =

И», 0

"2 У

1/10 о 0 1

Для (22), (23) построим матрицу Б = Хи • Х1Г, для которой ЛГф) = 2.59 • 104.

Таким образом, расчеты показали уменьшение числа обусловленности матрицы Грама после процедур масштабирования входных и выходных данных, приблизительно в 10 раз.

В третьей главе описывается методика построения кусочно-разностных и кусочно-дифференциальных моделей в форме Коши на основе малого объема данных наблюдений за объектом с изменяющимися тенденциями в поведении. Особенность построения подобных моделей заключается в выявлении границ подынтервалов, внутри которых не изменяются тенденции поведения объекта. Далее, так как модели стохастические, для построения одной модели в форме Коши предлагается алгоритм расчета дисперсии шумов динамики и измерителя. В главе показан способ оценивания параметров непрерывной модели по дискретным наблюдениям. Алгоритмы и методика апробированы на тестовых примерах. Алгоритм оценивания дисперсии в модели в форме Коши состоит из следующих шагов.

Шаг 1. Представим модель наблюдений на подынтервале в виде х('ш) = /■ х(1к) + М1к), х(10) = х0, у(1к+х) = х(1к+1) + у(4+,), к = 0.ЛГ-1 . Для оценивания среднего значения и дисперсий шумов сформированы ряды по двум и трем смежным измерениям следующим образом:

)(2) = >'(4) -У0к-\). к = 2,3,... N; Ч4)(3) =у{1к)-1/2X4-,)-1/2><4_2), к=3,4,.. Ж

Шаг 2. Оценивание среднего значения прироста ряда д по формуле

<?(4 14) = <?(4-1 IД + ц"^'(4}(2)-<7(4-. 14-,)), = , ад 1/0 = 0.

Окончательно <7 ~ </(/у | /Л,).

/¿/яг 5. Построение последовательностей произведений центрированных значений невязок по формуле

=2№)(3)-|?(414)), к=зл-.,м.

Шаг 4. Оценивание дисперсии шума модели динамики сг2 в соответствии с выражением

<71(1к 11к) = а;(4_, 1Я/*)(и,) -^(4-1)14-,)), * = , ¿1(121 /2) = 0.

Окончательно

Шаг 5. Построение последовательностей значений невязок по формуле

Я4)(у) =\тк)т-Мк \1к)?-*1ок 14)]. *=2,3,...,*.

Шаг б. Оценивание дисперсии шума измерителя сг2 на основе выражения

о-2(4 14) = о"2(4-1 14-1)+|«М>'(4)м-^2(4-, 14-1)). * = |/,) = о.

Окончательно

Шаг 7. Предсказанные значения оценок состояния на один шаг \ 4)

и оценки фильтрации х(1м \1к+1) рассчитываются в соответствии с алгоритмом калмановской фильтрации.

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

Шаг 1. Изменение тенденции ряда наблюдений для моделей первого порядка определяется визуально и подтверждается на основе алгоритма Г. Чоу, описанного в главе 1. Это позволяет весь период наблюдений разбить на подынтервалы.

Шаг 2. Для каждого из полученных подынтервалов рассчитываются дисперсии шумов измерителей для дискретной модели по формуле

Д.0=—!— Х0#)-й№)2> = где - оценка математичес-

И—\к=\ N к=\

кого ожидания на г'-м подынтервале, М - число подынтервалов, либо по рекуррентному алгоритму, рассмотренному выше.

Шаг 3. При больших шумах необходимо данные сглаживать с помощью алгоритма регуляризирующего кубического сплайна. Суть данного алгоритма

состоит в определении дважды дифференцируемого сплайна вида $(х)=ц+/}(х-х,)+с;(х-х,)2 +с11(х-х1)3 на каждом интервале [х,.,х,+1)е[а,6], / = 1 ,п,

где [а, 6] - область определения сплайна, и в минимизации функционала

ь п

= |(5 (х))2с& + - 5(х, ))2, где р - весовые коэффициенты.

а ¡=1

Шаг 4. На основе данных наблюдений с помощью МНК оцениваются коэффициенты дискретной и/или непрерывно-дискретной модели на каждом подынтервале.

Шаг 5. Для каждого подынтервала оцениваются дисперсии помех состояния объекта. Для этого разностное уравнение записывается в виде

у=1 р=1

либо по предыдущему алгоритму, изложенному выше.

Для непрерывно-дискретной модели необходимо предварительно получить новые измерения в виде рассчитанных значений производных х(к). Так как и(к) - гауссовский белый процесс с нулевым средним, оценки дисперсии помех состояния объекта на подынтервалах вычисляются по формуле

1 "Н _

йо (0 = ——- X н'< (кр), ¡ = 1, М, где /V,- объем выборки /' -го подынтервала, л,- -1 р=о

Шаг 6. Вычисляются оценки дисперсии помех начального состояния Рю (0) на подынтервалах как рекурсивное решение для устойчивого состояния процесса с помощью соотношения

Р?о=0,г = 1М, ,= 0,1,2,....

Шаг 7. Переходим на следующий подынтервал / = / +1. Если / < М, переходим на шаг 2, иначе на шаг 8.

Шаг 8. Для получения единой модели оценки параметров на всех подынтервалах объединяются и представляются одной функцией от некоторого независимого параметра с помощью сглаживающего кубического сплайна. Окончательно для всего интервала наблюдений записывается кусочно-разностная или кусочно-непрерывная модель в форме Коши в виде

х,. (к +1) = ¿Д- • х0 (к) • и,р (к) + (к), I = Щ х,. (0) = хю,

у=1 Р=1

(к +1) = Н, ■ X,.(А: +1) + у,.(¿ + 1), £ = 0,^-1,

МО(0+^,(0, х,(о)=х,0,

М1ы) = Н, ■x¡(lk+t) + v,(lkit), к = Щ

¿¿/яг 9. Полученные модели подвергаются совокупной проверке свойств, связанных с устойчивостью, наблюдаемостью, управляемостью и идентифици-

руемостью при и{к) = 0 и у(к) = 0. В работе приводятся критерии этих свойств.

Шаг 10. Полученные кусочно-разностная или кусочно-непрерывная модели при необходимости используются для оценивания состояния объекта с помощью уравнений фильтра Калмана.

Для получения непрерывно-дискретной модели из дискретной в работе предлагаются два подхода оценивания параметров: 1) использование сглаживающего кубического сплайна с регуляризацией как инструмент, позволяющий осуществлять расчеты численного дифференцирования; 2) в решении двух задач: а) оценивание параметров дискретной модели по дискретным измерениям; б) оценивание коэффициентов непрерывной модели динамики объекта по оценкам параметров дискретной модели.

Таким образом, в главе изложены: алгоритм оценивания дисперсий шумов модели в форме Коши, методика построения кусочно-разностной и/или кусочно-дифференциальной модели в форме Коши на основе данных наблюдений, численные расчеты по апробированию предложенной методики и алгоритмов.

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

Согласно методике построения кусочно-разностной модели в форме Коши, для подынтервала № 1 получены следующие характеристики: шум измерителя Д10=9.45МПа2, шум динамики 01 о = 10.64 МПа2; шум начального состояния Рю = 10.64 МПа2. При этих значениях разностная модель имеет вид: [х, (т +1) = 0.45х, (т) + 6.22щ (т) + 0.59и2 (т +1) + (т), х1 (0) = х0,

|7,(/и + 1)=х,(/и + 1) + у1(/я + 1), т = 0,9, где х{т) - значение прочности цемента, г/, (от) - значения количества вводимого гипса, и2{т) - значение количества вводимой минеральной добавки, у{т) - значение прочности цемента на основе измерительной системы.

Аналогично для подынтервала № 2: Я20 = 22.5МПа2; @20 = 26.39 МПа2; Р20 = 26.39МПа2 и разностная модель в форме Коши имеет вид

х2 (т+1) =0.06х2 (т) - 2.94м, (т) +15.98и2 {т +1) + и>2 (/и), х2 (0) = х0, у2(т+Х) =х2(т+1) + У2(т+1), т = 9,\Ъ.

Для получения единой кусочно-разностной модели на всем интервале исследования параметры , /2, соответствующие моделям двух подынтервалов, можно «сшить» с помощью сглаживающего кубического сплайна и получить функцию /(да). Аналогично, можно построить и функции (т), §2 (,и) • Полученная кусочно-разностная модель в форме Коши позволит применить аппарат фильтра Калмана для решения задачи фильтрации состояния прочности цемента как достоверной оценки. Реальные и фильтрационные значения прочности цемента представлены на графике (рисунок 1).

—•—ТВО, реальные значения

■ -•- - ТВО, фильтрац. значения

—*— Норм, усл., 3 суток, реальные значения

- — - Норм, усл., 3 суток, фильтрац. значения

— Норм, усл., 7 суток, реальные значения

- -О- - Норм, усл., 7 суток, фильтрац. значения

—»—Норм, усл., 14 суток, реальные значения

- -Норм, усл., 14 суток,фильтрац. значения

—'•— Норм, усл., 28 суток, реальные значения

Рисунок 1 - Зависимость прочности цементного камня, изготовленного на основе свежеприготовленного клинкера, от введенного в него волластонита

Расчеты показали, что более достоверные оценки фильтрации прочности цемента достигают своего максимального значения добавки при 9 % от общей массы исходного количества клинкера и что разброс значений ошибок оценок фильтрации составляет ± 1,5 %.

Далее, в главе рассмотрена задача построения кусочно-дифференциальной модели в форме Коши для плотности и фунщии распределения, необходимой для определения показателя качества промышленной продукции. Для выявления качества выпускаемой предприятием продукции проведены измерения длин деталей {уъ,У\,-..,Уц', N = 50}. Для полученного ряда требуется построить математическую модель плотности и функции распределения вероятностей проведенных измерений. Допустим, что в ходе контроля качества продукции получены частоты попадания длин элементов из выборочной партии в контрольные интервалы. Для контроля качества изделий (оценки плотности распределения) строится непрерывно-дискретная модель в форме Коши вида:

х = + х(0) = х0, ук=хк+ек, (24)

где х - состояние объекта, характеризующее частоту попадания ряда измерений; ук - выход измерительной системы; ек— случайная погрешность к-то

Количество волластонита, %

измерения; х0 - начальное состояние; г = 1 ,М, М - число подынтервалов.

Интервал наблюдения визуально разбивается на два подынтервала [1, 4] и [4, 7], где граница наблюдения определяется существенным изменением поведения объекта. На каждом подынтервале вычисляются значения параметров а{ и Ь1 и строится модель вида (24). В результате вычислений найдены параметры на подынтервале [1, 4]- {а^-0.05, ^ =5.08}, а на подынтервале [4,7]-{¿2 =-0.21 , ¿2 =-2.77}. Для набора накопленных частот, то есть для подынтервала [1, 7], параметры модели равны {а = -0,09; ¿=10,58}. На рисунке 2 (а) представлены данные накопленных частот и график одной дифференциальной модели. Как видно, модель не обладает достаточной

точностью.

§ 30

л № 20

ши*

1 I

/ С!

^ !

а)

5 6 7 Длина детали, I

4 5 6 7 Длина детали, I

Рисунок 2 - Гистограмма накопленных частот и графики дифференциальных моделей: а - дифференциальная модель для всего интервала [1; 7]; б - кусочно-дифференциальная модель на подынтервалах [1; 4] и [4; 7].

С целью повышения точности модели интервал наблюдения разбивается на два подынтервала: [1, 4] и [4, 7]. На каждом подынтервале строится своя

дифференциальная модель с параметрами: {а, = -0.31; ¿,=6.0} и {а2=-0.42;

Ь2=23.8} соответственно. График полученной кусочно-дифференциальной модели представлен на рисунке 2 (б). Численные результаты показали, что уменьшение погрешности модели на двух подынтервалах определено вычислением значений среднеквадратичных отклонений исходных данных от результатов моделей, которые для случая на всем интервале равно 3.49, а для кусочно-дифференциальной модели равно 1.03 .

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

Результаты наблюдения и обработки статистических данных показывают,

что объемы поставок зернового сырья на предприятие имеют сезонный характер. Это позволило использовать разработанную методику и получить два дифференциальных уравнения в форме Коши, построенных на двух подынтервалах:

Ш) = 0,08-хО) + 56,12■ и(0 + s{t), t е[t0,tN], x(t0) = х0 = y(tt);

= + k = \,9~,

\x{t) = -0,42• x{t)-12,10• u(f) + e(t), te[tm,tm], x(t0) =x0 = y2(l);

\y{tk) = x{tk) + v{tk), k = \J .

К полученным моделям (25) и (26), применив аппарат фильтра Калмана, получить предсказанные и фильтрационные оценки (рисунок 3).

, 600 т

(25)

(26) можно

г 5оо ■■

то X

® 400

Штящ

1 зоо

о о

о 20D i-

03 J-'' о

| 100 - -

О I

- реальные предсказанные | фильтрационные I

Рисунок 3 - Графики поведения реальных, предсказанных и фильтрационных поставок зерновых культур за 12 месяцев 2004 г.

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

В работе решена задача оценки невязки аналитических и экспериментальных данных индуктивности микрополосковой линии на основе дифференциальной модели в форме Коши с использованием фильтра Калмана.

В заключении сформулированы основные результаты исследования.

В приложении приведены результаты обработки данных и акты о внедрениях результатов исследований.

ОСНОВНЫЕ РЕЗУЛЬТАТЫ ИССЛЕДОВАНИЯ

В диссертационной работе получены следующие основные результаты.

1. Алгоритмы построения линейных моделей в форме Коши вместо линейных и нелинейных регрессионных моделей, где входные данные, состояние

динамики и выходные данные измерительной системы зависят от невременного аргумента и содержат в себе малое число наблюдений.

2. Алгоритм преобразования стохастической модели с распределенными параметрами к адекватной модели в форме Коши во временной области.

3. Рекуррентные алгоритмы вычислений дисперсий шумов динамики и измерительной системы на основе дискретных данных наблюдений.

4. Алгоритм уменьшения числа обусловленности информационной матрицы Грама для непрерывно-дискретной модели в форме Коши на основе процедуры масштабирования входных и выходных данных.

5. Методика построения стохастической линейной кусочно-разностной и кусочно-дифференциальной модели в форме Коши для объектов с малым объемом наблюдений, где вход, состояние, выход измерительной системы и коэффициенты модели динамики зависят не обязательно от временного аргумента.

6. На основе разработанных методики и алгоритмов получены следующие модели в форме Коши: кусочно-разностная модель для оценивания прочности цементного камня; кусочно-дифференциальная модель для плотности и функции распределения длин деталей, необходимой для оценивания качества промышленных продукций; кусочно-дифференциальная модель для прогнозирования объема поставок зернового сырья в мукомольно-комбикормовые комбинаты; дифференциальная модель для оценки невязки аналитических и экспериментальных данных индуктивности микрополосковой линии.

ОСНОВНЫЕ ПУБЛИКАЦИИ ПО ТЕМЕ ДИССЕРТАЦИИ

Статьи в журналах, рекомендованных ВАК:

1. Абденова Г.А. Структурно-параметрическая идентификация систем с распределенными параметрами с использованием модели типа «вход-состояние-выход» // Научный вестник НГТУ. Новосибирск: НГТУ, 2010. № 1 (38). С. 9-16.

2. Абденова Г.А. Прогнозирование значений уровня временного ряда на основе уравнений фильтра Калмана // Ползуновский вестник. Барнаул: АлтГТУ, 2010. № 2. С. 4-6.

3. Абденова Г.А., Апсалямов H.A., Толысбаев Б.С. Прогнозирование уровней поставок зернового сырья на мукомольное производство // Известия международной академии наук высшей школы. Новосибирск: НГТУ, 2006. № 1 (35). С. 114-124.

4. Абденова Г.А., Ильина JI.B., Раков М.А. Применение модели в пространстве состояний для достоверного оценивания прочности цементного камня // Фундаментальные исследования, 2011. № 8. С. 589-592.

5. Абденова Г.А., Рубанович М.Г., Хрусталёв В.А. Оценка невязки экспериментальных данных и аналитических расчетов интегральной индуктивности микрополосковой линии // Научный вестник НГТУ. Новосибирск: НГТУ, 2011. №3 (44). С. 29-42.

6. Абденова Г.А. Воевода A.A. Обусловленность информационной матрицы Грама в задаче идентификации: масштабирование входных и выходных

сигналов // Доклады ТУ СУР. Томск: Из-во ТУ СУР, 2012. №1 (25). Часть 1. С.131-135.

Статьи, опубликованные в других изданиях:

7. Абденова Г.А. Стохастическая нестационарная модель для функции Кобба-Дугласа в пространстве состояний // Сб. научных трудов НГТУ. Новосибирск: Изд-во НГТУ, 2009. № 3 (57). С. 41-46.

8. Абденова Г.А. Построение кусочно-дифференциальной модели в пространстве состояний для контроля качества продукции // Сб. научных трудов НГТУ. Новосибирск: Изд-во НГТУ, 2009. № 3 (57). С. 47-56.

9. Абденова Г.А., Воевода A.A. Идентификация непрерывных многомерных систем по дискретным наблюдениям // Сб. научных трудов НГТУ. Новосибирск: Изд-во НГТУ, 2010. № 2 (60). С. 165-168.

10. Абденова Г.А., Воевода A.A. Оценивание параметров и характеристик шумов нестационарных процессов в стохастических системах, описываемых в пространстве состояний // Сб. научных трудов НГТУ. Новосибирск: Изд-во НГТУ, 2010. № 3 (61). С. 11-18.

11. Абденова Г.А. Методика построения кусочно-разностной модели в пространстве состояний на основе нестационарного ряда наблюдений // Сб. науч. трудов НГТУ. Новосибирск: Изд-во НГТУ, 2011. № 2(64). С. 3-12.

12. Абденова Г.А., Воевода A.A. Один алгоритм перехода от регрессионной модели к модели в форме пространства состояний // Сб. научных трудов НГТУ Новосибирск: Изд-во НГТУ, 2012. № 3 (69). С. 3-10.

13. Абденова Г.А., Толысбаев Б.С. Моделирование поставки зерна на мукомольное производство // Компьютерные методы и системы в автоматике и электротехнике: материалы VI Среднеевропейской научно-технической конференции, Ченстохова, Польша, 16-17 сентября 2005 г. Ченстохова: Изд-во Ченстоховского технологического университета, 2005. Т. 2. С. 57-59.

14. Abdenova G.A. The method of construction piece-differential model on the basis of nonstationary observations // The Second Russian-Indian Joint Workshop on Computational Intelligence and Modern Neuristics in Automation and Robotics: proceedings RFBR and DST Sponsored, Novosibirsk, 10-13 September 2011. Novosibirsk: Publ. of NSTU, 2011. Additional volume. P.7-10. [Методика построения кусочно-дифференциальной модели на основе нестационарных наблюдений]

15. Абденова Г.А. Масштабирование вход-выходных данных и число обусловленности матрицы Грама // Актуальные проблемы электронного приборостроения: материалы XI международной конференции, Новосибирск, 2-4 октября 2012 г. Новосибирск: Изд-во НГТУ, 2012. Т. 6. С.16-18.

Отпечатано в типографии Новосибирского государственного технического университета 630073, г. Новосибирск, пр. К. Маркса, 20, тел. 346-08-57 формат 60x84 1\16, объем 1,5 п. л., тираж 100 экз., заказ № 187 , подписано в печать 15.01.2013

Оглавление автор диссертации — кандидата технических наук Абденова, Гаухар Амирзаевна

ВВЕДЕНИЕ.

1. ПОСТАНОВКА ЗАДАЧ ПАРАМЕТРИЧЕСКОЙ ИДЕНТИФИКАЦИИ.

МОДЕЛИ В ФОРМЕ КОШИ С НЕВРЕМЕННЫМ АРГУМЕНТОМ

1.1. Проблемы и задачи идентификации.

1.2. Переход от регрессионных моделей к линейным моделям в форме Коши. в зависимости от невременного аргумента.

1.3. Задачи параметрической идентификации для линейных моделей. в форме Коши

1.4. Преобразование стохастической модели распределенного типа к. стохастической линейной модели в форме пространства состояний.

1.5. Проблема плохой обусловленности в задачах параметрической. идентификации моделей в форме Коши.

1.6. Совокупная проверка свойств модели в форме Коши.

1.7. Оптимальная фильтрация.

1.8. Аппроксимация регуляризирующим кубическим сплайном.

1.9. Анализ тенденции ряда наблюдений при структурных изменениях.

1.10. Постановка задач для исследования.

2. АЛГОРИТМЫ ПЕРЕХОДА ОТ РЕГРЕССИОННЫХ И ПРОСТРАНСТВЕННО-ВРЕМЕННЫХ МОДЕЛЕЙ К МОДЕЛЯМ В ФОРМЕ КОШИ СОСРЕДОТОЧЕННОГО ТИПА.

2.1. Построение модели в форме Коши для линейной регрессионной модели

2.1.1. Оценивание параметров в регрессионной модели на основе. данных наблюдений.

2.1.2. Алгоритм построения модели в форме Коши в случае, когда. вход и выход носят случайный характер.

2.2. Алгоритм построения линейной модели в форме Коши для нелинейной. регрессионной модели.

2.3. Параметрическая идентификация систем с распределенными. параметрами с использованием модели типа «Вход-Состояние-Выход».

2.4. Процедура масштабирования входных и выходных данных наблюдений

2.5. Выводы по главе.

3. ПОСТРОЕНИЕ КУСОЧНО-РАЗНОСТНЫХ И КУСОЧНО-ДИФФЕРЕН

ЦИАЛЬНЫХ МОДЕЛЕЙ В ФОРМЕ КОШИ.

3.1. Особенность построения кусочно-разностных и кусочно-. дифференциальных моделей.

3.2. Получение соотношений для оценок дисперсий шумов динамики. объекта и измерительной системы.

3.3. Методика построения кусочно-разностных и кусочно-. дифференциальных моделей в форме Коши.

3.3.1. Методика построения математических моделей в форме Коши.

3.3.2. Способы оценивания параметров непрерывно-дискретной модели

3.3.3. Исследование методики построения кусочно-дифференциальной. модели в форме Коши для тестовых примеров.

3.4. Выводы по главе.

4. ПОСТРОЕНИЕ МОДЕЛЕЙ В ФОРМЕ КОШИ С ПЕРЕМЕННЫМИ.

ПАРАМЕТРАМИ ДЛЯ РЕШЕНИЯ ПРОИЗВОДСТВЕННЫХ ЗАДАЧ.

4.1. Построение кусочно-разностной модели для оценивания прочности. цементного камня.

4.2. Построение кусочно-дифференциальной модели для плотности и. функции распределения частот длин промышленных деталей.

4.3. Оценивание индуктивности микрополосковой линии.

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

4.5. Выводы по главе.

Введение 2012 год, диссертация по информатике, вычислительной технике и управлению, Абденова, Гаухар Амирзаевна

Теоретико-методологической основой исследования по вопросам параметрической идентификации, по оцениванию состояния и управления в технических и других системах послужил обзор концепций, сформулированных в работах зарубежных теоретиков и практиков, авторами которых являются: П. Эйкхофф, К. Ост-рем, Р. Мехра, Р. Бен Мрад, Е. Фараг, А. Кисикайя, X. Кайран Ахмет, J1.C. Лыонг, Р.Л. Кашьяп, А.Р. Рао, Т. Кайлац, М.С. Гревел и др. [1, 3,10,17, 20, 26, 34, 82]. Эти задачи решались и в трудах отечественных ученых: A.C. Анисимова, Ю.Е. Воско-бойникова, В.Г. Горского, Н.Д. Егупова, A.A. Ломова, В.Н. Овчаренко, М.А. Огар

2 Пример: результаты выхода некоторого сорта стали в зависимости от угара кремния, измеряемого в процентах. кова, В.Н. Подладчикова, И.Н. Синицына, В.А. Фурсова и др.[2, 5, 12, 21, 38, 39, 49,50, 56, 74].

Целью диссертационном работы является повышение эффективности оценивания состояния объектов с малым объемом наблюдения и с явно выраженными изменяющимися тенденциями в поведении путем идентификации объекта математическими моделями в форме Коши с вещественным аргументом.

Задачи исследования. Для достижения заданной цели поставлены следующие задачи исследования.

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

2. Рассмотреть объект с нестабильной динамикой развития, состояние которого зависит от невременного аргумента и ранее описываемый линейными и нелинейными регрессионными моделями3 с целью перехода к описанию с помощью разностной модели в форме Коши.

5. Для повышения точности найденных оценок параметров в модели динамики, необходимо разработать алгоритм уменьшения числа обусловленности"* информационной матрицы Грама с помощью процедуры масштабирования входных и выходных данных.

3 Регрессионная модель — модель, основанная на уравнении регрессии или системе регрессионных уравнений, связывающих величины входных (объясняющих) и выходных (зависимых) переменных.

4 Понятие обусловленность употребляется для общей характеристики чувствительности параметров модели к ошибкам в исходных данных наблюдаемых величин.

Научная новизна работы и результаты, выносимые на защиту.

4. Алгоритм рекуррентного оценивания дисперсий шумов динамики объекта и измерительной системы на основе дискретных данных наблюдений.

5. Методика построения единой линейной кусочно-разностной или линейной кусочно-дифференциальной модели в форме Коши для решений задач тестового и практического характеров.

Теоретическая значимость. Разработанные в диссертации методика и алгоритмы позволяют при малых выборках построить в виде одной линейной стохастической кусочно-разностной и/или кусочно-дифференциальной модели в форме Коши, в которых переменные параметры состояния объекта и выход измерительной системы зависят от невременного аргумента.

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

Внедрение результатов исследования. Методика построения кусочно-разностной модели в форме Коши внедрена в ООО «Монтажстрой» (г. Новосибирск, РФ) для оценивания прочности цементного камня, хранившегося длительное время. Кусочно-разностная модель в форме Коши используется для оценивания активности цемента, прочность которого оценивается в зависимости от количества вводимой минеральной вяжущей добавки.

Методика построения кусочно-дифференциальной модели в форме Коши внедрена как инструмент, позволяющий осуществлять оперативный контроль качества продукции на АО «Востокмашзавод» (г. Усть-Каменогорск, РК) для изготовленных деталей металлургического оборудования, а именно редуктора поворота ковша, деталей для захвата автоматического и вакуумного ковша. Кусочно-дифференциальная модель описывает плотность распределения значений длин деталей, характеризующих точность как один из показателей контроля качества продукции.

Методика построения стохастической кусочно-дифференциальной модели в форме Коши внедрена для расчета индуктивных и емкостных параметров элементов аттенюатора в ФГУП «РТРС», филиал «Липецкий областной радиотелевизионный передающий центр» (г. Липецк, РФ).

Методика построения стохастической кусочно-дифференциальной модели внедрена для расчета количества прогнозных поставок зернового сырья в муко-мольно-комбикормовые комбинаты для решения задач транспортной логистики. Разработанная методика апробирована на тестовых примерах и на примерах с реальными данными относительно поставок сельскохозяйственной продукции в течение периода с 2001 по 2004 год в ОАО «Восточно-Казахстанский мукомоль-но-комбикормовый комбинат» (г. Семей, РК).

Следует отметить, что диссертационная работа поддержана грантом по заданию Министерства образования и науки РФ, по проекту № 7.559.2011, «Исследование предельных точностей оптических методов измерения параметров движения и мехатронных методов управления движением и разработка новых робото-технических и электромеханических систем», Темплан, гос. per. номер НИР 01201255056.

Апробация результатов исследования. Основные положения и выводы исследования докладывались и обсуждались на: 7-м Российско-Корейском международном симпозиуме науки и техники, 2003; Международном симпозиуме науки и техники (Новосибирск, 2003); Международной научно-практической конференции «Электронные средства и системы управления» (Томск, 2004); III Международной научно-практической конференции аспирантов и молодых ученых «Страны СНГ в условиях глобализации» (Москва, РУДН, 2004); VI Среднеевропейской научно-технической конференции «Компьютерные методы и системы в автоматике и электротехнике», (Ченстохова, Польша, 2005); Международной научно-технической конференции «ИКИ-2001», «Измерение, контроль, информатизация» (Барнаул, 2001, 2003, 2005, 2007, 2009, 2010); II Российско-индийском объединенном семинаре «Computational Intelligence and Modern Heuristics in Automation and Robotics», 10-13 сентября 2011 г.; XI Международной конференции «Актуальные проблемы электронного приборостроения», Новосибирск, 2-4 октября 2012 г.

Обоснованность н достоверность научных положений и выводов подтверждаются результатами теоретических, экспериментальных исследований и расчетов, публикациями в рецензируемых изданиях, апробацией на всероссийских и международных конференциях, а также использованием результатов диссертационной работы при решении задач производственного характера.

Заключение диссертация на тему "Идентификация объектов, описываемых линейными разностными и дифференциальными уравнениями в форме Коши с вещественным аргументом"

4.5. Выводы по главе

Разработанные во второй и третьей главах методика и алгоритмы построения математических моделей в форме Коши по данным наблюдений объекта с изменяющимися тенденциями в поведении, использованы для решения задач практического характера, в частности:

1. Получена кусочно-разностная модель в форме Коши для оценивания прочности цементного камня в зависимости от вводимой минеральной добавки. Полученная модель с помощью фильтра Калмана позволяет рассчитать оценки фильтрации как наиболее достоверные оценки прочности цементного камня, которые являются одной из существенных характеристик цементной продукции для ООО «Монтажстрой» (г. Новосибирск, РФ).

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

21 Предлагаемая методика была передана в мукомольно-комбикормовый комбинат г. Семей в 2005 году (Республика Казахстан). примере условных данных измерений длин деталей, изготовленных в промышленном предприятии АО «Востокмашзавод» (г. Усть-Каменогорск, РК).

3. Построена стохастическая модель в форме Коши для расчета индуктивности МПЛ в зависимости от ее ширины. Рассчитываются оценки невязки экспериментальных данных и аналитических расчетов индуктивности МПЛ, полученные методом токовых полос, с оценками фильтрации, рассчитанные по фильтру Кал-мана. Показано, что наилучшее согласование экспериментальных оценок и аналитических расчётов индуктивности достигается при определенном числе продольных разбиений МПЛ. На основе полученной модели были осуществлены расчеты индуктивных и емкостных параметров элементов аттенюатора для ФГУП «РТРС» филиал «Липецкий областной радиотелевизионный передающий центр» ( г. Липецк, РФ).

4. Построена кусочно-дифференциальная модель в форме Коши вместо регрессионной модели для прогнозирования поставок зерновых культур в муко-мольно-комбикормовый комбинат (г. Семей, РК) в течение 2004 года. Полученная кусочно-дифференциальная модель была использована для решения задачи ежемесячного прогнозирования поставок зернового сырья на 2005 год с последующей корректировкой в режиме реального времени.

22 Аттенюатор - это устройство для понижения напряжения или мощности электрических или электромагнитных колебаний.

125

ЗАКЛЮЧЕНИЕ

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

В данной диссертационной работе следующие результаты являются оригинальными:

1. Методика построения линейных стохастических кусочно-разностных и кусочно-дифференциальных моделей в форме Коши на основе данных наблюдений малого объема выборки. При этом переменные входа, состояния, выхода и переменные коэффициенты в модели динамики зависят от вещественного (не обязательно временного) аргумента;

2. Алгоритмы построения линейных моделей в форме Коши для объектов, описываемых линейными или нелинейными регрессионными моделями. Полученные модели позволят более достоверно оценивать состояние исследуемого объекта с помощью уравнений фильтра Калмана;

3. Алгоритм преобразования стохастической модели с распределенными параметрами в пространстве состояний к адекватной стохастической модели в форме пространства состояний во временной области. Это позволит решать задачи параметрической идентификации и прогнозирования;

4. Рекуррентный алгоритм вычислений дисперсий шумов динамики и измерительной системы для линейной дискретной модели в форме Коши на основе дискретных данных измерений малого объема выборки без учета оценок параметров модели динамики;

5. Алгоритм улучшения числа обусловленности для непрерывно-дискретной модели в форме Коши на основе процедуры масштабирования входных и выходных данных измерений. Данная процедура позволит повысить точность параметров модели динамики;

6. Построена кусочно-разностная модель в форме Коши для достоверного оценивания прочности цементного камня в зависимости от количества введенного минерального вяжущего вещества - волластонит;

7. Построена кусочно-дифференциальная модель в форме Коши для плотности и функции распределения частоты длин деталей, необходимой для решения задачи качества продукций;

8. Построена стохастическая модель в форме Коши для оценки точности аналитического расчета индуктивности микрополосковой линии, полученной методом токовых полос. Показано, что наилучшее согласование экспериментальных оценок и расчётов индуктивности, полученных по фильтру Калмана, достигается при определенном числе поперечных разбиений МПЛ, которое является более эффективным для аналитического метода;

9. Построена кусочно-дифференциальная модель в форме Коши для прогнозирования объема поставок зернового сырья на мукомольно-комбикормовые комбинаты, для решения задачи транспортной логистики.

127

Библиография Абденова, Гаухар Амирзаевна, диссертация по теме Системный анализ, управление и обработка информации (по отраслям)

1. Эйкхофф П. Основы идентификации систем управления. М.: Мир, 1975. 688 с.

2. Горский В.Г., Адлер Ю.П., Талалай А.М. Планирование промышленных экспериментов (модели динамики). М.: Металлургия, 1978.112с.

3. Astrom K.J. Máximum Likelihood and Prediction Error Methods // Automática. 1980. Vol. 16, No 5. P. 551-574.

4. Райбман H.C. Что такое идентификация? M.: Наука, 1970. 119с.

5. Анисимов А.С. Идентификация объектов управления: Учебное пособие. Новосибирск: Издание НЭТИ, 1985. 80с.

6. Береза В.Ю., Ясинский В.К. Фильтр Калмана-Бьюси для линейных стохастических динамических систем с разрывными траекториями // Кибернетика и системный анализ: международный научно-теоретический журнал. 2003. N 2. С.89-100.

7. Сейдж Э., Уайт Ч.С. Оптимальное управление системами. М.: Радио и связь, 1982. 292 с.

8. Драница Ю.П. Об одном методе моделирования нестационарных динамических систем и процессов // Вестник МГТУ, 2000. Том 3. №1. С.55-66.

9. Анисимов А.С., Толстиков А.С. О некоторых подходах к синтезу алгоритмов идентификации нестационарных объектов // В кн. «Проблема идентификации нестационарных объектов в измерительной технике». Труды СНИИМ. Новосибирск. 1971. Вып. 17. С. 36-56.

10. Mehra R.K. On the Identification of Variences and Adaptive Kalman filtering // IEEE Trans. Autom. Control, v. AC 15.1970, No.2. P. 175-184.

11. Каленова В. И., Морозов В. М., Соболевский П. М. К вопросу об исследовании линейных нестационарных систем // Вестник МГУ. Сер.1. 2009. С. 50-60.

12. Огарков М.А. Методы статистического оценивания параметров случайных процессов. М.: Энергоатомиздат, 1990. 208с.

13. Гинсберг К. С., Басанов Д. М. К вопросу об общей методологии идентификации // Труды 5 Международной конференции "Идентификация систем и задачи управления", Москва, 30 янв.-2 февр. 2006 г. М: Изд-во ИПУ РАН, 2006. С. 131-141.

14. Учебник в 5-ти томах: Т.1. Математические модели, динамические характеристики и анализ систем автоматического управления / Под ред. К.А. Пупкова, Н.Д. Егупова. М.: Изд-во МГТУ им. Н.Э. Баумана, 2004. 656 с.

15. Zhang Huanshui, Lu Xiao, Cheng Daizhan. Optimal estimation for continuous-time systems with delayed measurements // IEEE Trans. Autom. Contr. 2006. N 5. C. 823-827.

16. Острем К. Введение в стохастическую теорию управления. Пер. с англ. М.: Мир, 1973. 322 с.

17. Chou J.-H., Horng I.-R. Parameter identification of lumped time-varying systems via shifted Chebyshev series // Int. J. Syst. Sei. 1986.V.17. No. 3. P.459-464.

18. Шалаев Ю.Н. Моделирование нестационарных динамических систем методом изображающих векторов // Известия Томского политехнического университета. 2006. Т. 309. № 7. С.44-47.

19. Льюнг Л. Идентификация систем. Теория для пользователя. Под ред. Я.З. Цыпкина. М.: Наука, 1991.432 с.

20. Ломов A.A. Оценка трендов и идентификация динамики временных рядов на коротких интервалах наблюдения // Известия РАН. Теория и системы управления, 2009. №1. С. 5-17.

21. Ляпунов A.M. Общая задача об устойчивости движения. М.: Гостехиздат, 1950. 472с.

22. Крылов В.И., Бобков В.В., Мопастырный П.И. Вычислительные методы высшей математики. Минск: Вышэ. шк., 1972. Т1. 584 с.

23. Гребенюк Е.А. Методы анализа нестационарных временных рядов с неявными изменениями свойств (обзор) // Автоматика и телемеханика. 2005. №12.1. С.3-28.

24. Пучков А. 10., Шарпаев А. К. Прогнозирование рынка промышленной продукции на основе применения фильтра Калмана // Интеграл. 2009. № 2. С. 100-105.

25. Кашьяп PJL, Pao А.Р. Построение динамических стохастических моделей по экспериментальным данным. М.: Наука, 1983. 384 с.

26. Фурсов В.А., Чудилин А.Г. Проекционные методы оценивания параметров моделей прогнозирования по малому числу наблюдений // Искусственный интеллект: Симферополь, 2000. № 2. С. 212-218.

27. Фаддеев М.А., Марков К.А. Основные методы вычислительной математики: Учебное пособие. СПб: Издательство «Лань», 2008. 100с.

28. Бессонов A.A., Загашвили Ю.В., Маркелов A.C. Методы и средства идентификации динамических объектов. Л.:Энергоатомиздат, 1989. 280 с.

29. Coca Daniel, Billings Stephen A. Identification of finite dimensional model of infinite dimensional dynamical systems // Automatica: A Journal of IFAC the International Federation of Automatic Control. 2002. 38, N 11. C. 1851-1865.

30. Демиденко E.3. Линейная и нелинейная регрессии. M.: Финансы и статистика, 1981. 302 с.

31. Богословский C.B., Богословский B.C. Динамика нестационарных систем с равномерно изменяющимися во времени коэффициентами // Научные прибо-ростр. 2002. 12, N 3. С. 83-92.

32. Бутковский А.Г. Характеристики систем с распределенными параметрами (справочное пособие). М.: Наука, 1979. 224с.

33. Grewal M.S., Andrews А.Р. Kaiman filtering: theory and practice, using MATHLAB. N.Y.: Willey, 2001. 401p.

34. Кочетков С.А., Уткин В.А. Оценка характеристик шумов в линейных системах в режиме реального времени // Труды 8 Международной конференции «Идентификация систем и задачи управления», Москва, 26-30 января 2009 г. М.: ИПУ РАН, 2009. С. 1073-1093.

35. Сидоров Д. Н. Моделирование нелинейных нестационарных динамических систем рядами Вольтерра: идентификация и приложения // Сибирский журн.индустр. матем., 2000. Том 3. № 1. С. 182-194.

36. Nunez-Garcia Javier, Wolkenhauer laf. Random set system identification // IEEE Trans. Fuzzy Syst. 2002. Vol.10, N3. P. 287-296.

37. Згуровский M.3., Подладчиков B.H. Аналитические методы калмановской фильтрации для систем с априорной неопределенностью. К.: Наукова Думка, 1995. 584 с.

38. Гайшун И.В. Идентификация линейных нестационарных систем по реакции на обобщенные управления // Дифференциальные уравнения. 2008. 44, С. 301307.

39. Кедрин B.C. Методика прогнозирования нестационарных процессов на базе методов сингулярного разложения и искусственного интеллекта: автореф. дис. . канд. техн. наук по спец. 05.13.01. Иркутск. 2007 г. 22с.

40. Калинкина С.Ю. Методы представления интервальных динамических систем в пространстве состояний: автореф. дис. . канд. техн. наук по спец. 05.13.01. Барнаул. 2005 г. 19 с.

41. Тычинин А. В. Методы обратных задач динамики в задачах синтеза систем управления распределенным объектом // Вестник Самар. гос. техн. ун-та. Самара, 2005. N 39. С. 21-25.

42. Chang R.-Y., Yang S.-Y., Wang M.-L. Parameter identification of time-varying systems via generalized orthogonal polynomials // Int. J. Syst. Sei. 1988.V.19. No. 3. P.471-485.

43. Мышкис А.Д. Лекции по высшей математике: Учебное пособие. СПб.: Издательство «Лань», 2007. 688 с.

44. Кондратов Д. В. Прогнозирование временных рядов на основе использования полиномов Чебышева, наименее уклоняющихся от нуля // Вестник Самар. гос. техн. ун-та. Самара, 2005. № 32. С. 49-53.

45. Танеев P.M. Математические модели в задачах обработки сигналов. М.: Горячая линия Телеком, 2004. 80 с.

46. Методы анализа, синтеза и оптимизации нестационарных систем автоматического управления: Учебное пособие /К.А. Пупков, Н.Д. Егупов, В.Г. Коньков, JI.T. Милов, А И. Трофимов. Под ред. Н.Д. Егупова. М.: Изд-во МГТУ им. Н.Э. Баумана, 1999. 684с.

47. Воскобойников Ю.Е. Рекуррентное оценивание в динамических системах: Учеб. пособие. Новосибирск: Изд-во НГТУ, 2002. 92с.

48. Ахметжанов P.C. Никифоров А.Н. Применение вейвлет-анализа для исследования нестационарных процессов роторных систем // Проблемы машиностроения и автоматизации. 2005, № 2, с. 52-61.

49. Измайлов Р. А. Лебедев A.A. Вейвлет преобразования для анализа нестационарных процессов в центробежном компрессоре // Научно-технические ведомости СПбГПУ. 2008. т. 2. С. 179-182.

50. Li Kang, Peng Jian-Xun, Irwin George W. Быстрый метод идентификации нелинейных моделей. A fast nonlinear model identification method // IEEE Trans. Au-tom. Contr. 2005. Vol. 50, N 8. P. 1211-1216.

51. Фурсов В.А. Идентификация систем по малому числу наблюдений. Самара: Самар. гос. аэрокосм, ун-т, 2007. 80с.

52. Тырсин А. Н. Робастное построение линейных регрессионных моделей по экспериментальным данным // Заводская лаб.: Диагност, матер. 2005. вып. 71, № U.C. 53-57.

53. Згуровский М. 3., Панкратова Н. Д. Информационная платформа сценарного анализа в задачах технологического предвидения // Кибернетика и систем, анал.: Международный научно-теоретический журнал. 2003. № 4. С. 112-125.

54. Semper J. David Cabedo, demente Ismael Moya. Value at risk calculation through ARCH factor methodology: Proposal and comparative analysis // Eur. J. Oper. Res. 2003. Vol. 150, N 3. P. 516-528.

55. Стечкин C.B., Субботин Ю.Н. Сплайны в вычислительной математике. М.: Наука, 1976. 352 с.

56. Вержбицкий В.М. Основы численных методов. М.: Высш. шк., 2005. 840 с.

57. Ретинский B.C. Компьютерный мониторинг динамических процессов в сетях с распределенными параметрами // Информационные технологии. 2007. № 9. С.9-13.

58. Минаев Ю.Н. Филимонова О.Ю., Бенамеур JL Методы и алгоритмы решения задач идентификации и прогнозирования в условиях неопределенности в ней-росетевом логическом базисе. М.: Горячая линия Телеком, 2003. 205 с.

59. Шторм Р. Теория вероятностей. Математическая статистика. Статистический контроль качества: Перев. с нем. М.: Мир, 1970. 368с.

60. Костюк В.Н. Управление в переходных экономических процессах на основе анализа тенденций: автореферат . д-ра. экон. наук: 05.13.10 : Москва, 1999 г. 25 с.

61. Афанасьев В. Н., Крестникова Д. Г. Идентификация нестационарного объекта методом настраиваемой модели // Нейрокомпьютеры: разраб., применение. 2009. № 9. С. 39-46.

62. Чинков В. Н. Оценка методических погрешностей вычисления статистических характеристик эргодических случайных сигналов // Укр. метрол. ж. 2005. N3. С. 30-33.

63. Александров А. Г., Орлов Ю. Ф. Сравнение двух методов идентификации при неизвестных ограниченных возмущениях // Автоматика и телемеханика. 2005. № 10. С. 128-147.

64. Smidl Vaclav, Quinn Anthony. Mixture-based extension of the AR model and its recursive Bayesian identification // IEEE Trans. Signal Process. 2005. Vol. 53. № 9. P. 3530-3542.

65. Синицын И.Н. Фильтры Калмана и Пугачева. M.: Университетская книга, Логос, 2006. 640с.

66. Илюшин В. Б., Чадеев В. М. Алгоритм идентификации динамических объектов с учетом априорной информации о параметрах // Автоматика и телемеханика. 2006. № 7. С. 133-143.

67. Филипп В. Б., Щуров Н. И. Вейвлет-анализ как метод оценки нестационарных процессов в электротранспортном комплексе // Научные пробл. трансп. Сиб. и Дал. Вост., 2008. № 1. С. 267-270.

68. Абденов Х-Г.А., Абденова Г.А. О динамизации регрессионных моделей //

69. Проблемы развития малого и среднего бизнеса в РК: экономика, технология,финансы и менеджмент: материалы международной научно-технической конференции, Семипалатинск, 5-7 декабря 2001 г. Семипалатинск: Изд-во КазФЭИ, 2001. С. 403-405.

70. Конев В.В., Пергаменщиков С.М. Последовательные планы идентификации параметров динамических систем // Автоматика и телемеханика. 1981. №7. С.84-92.

71. Вальд А. Последовательный анализ. М.: Физматизд., 1960. 328с.

72. Абденова Г.А. Стохастическая нестационарная модель для функции Кобба-Дугласа в пространстве состояний // Сборник науч. трудов НГТУ. 2009. № 3(57). С. 41-46.

73. Kailath Т. An innovations approach to least-squares estimations, pt. 1: linear filtering in additive white noise // IEEE Trans. Autom. Control. V. AC-13. 1968. P.646-655.

74. Mehra R.K. Synthesis of optimal inputs for multiinput-multioutput (MIMO) Systems with process noise. Part I: frequency-domain synthesis. Part II: time-domain synthesis // Syst. Identif. Adv. and Case Stud. New York e.a., 1976. P.211-249.

75. Ли P. Оптимальные оценки, определение характеристик и управление. М.: Наука, 1966. 176 с.

76. Медич Дж. Статистически оптимальные линейные оценки и управление. М.: Энергия, 1973. 440 с.

77. Завьялов Ю.С., Квасов В.И., Мирошниченко В.Л. Методы сплайн-функций. М.: Наука, 1980. 352 с.

78. Абденова Г.А. Прогнозирование некоторых параметров технических систем на основе моделей регрессионного и динамического характеров // Сборник науч. трудов НГТУ. 2010. № 1(59). С. 3-12.

79. М. Де JIa Сен Идентификация одного класса непрерывных линейных нестационарных систем // Автоматика и телемеханика. 2002. № 2. С. 84-100.

80. Абденов А.Ж., Абденова Г.А., Снисаренко A.B. Построение и применение кубических сплайнов для сглаживания и дифференцирования данных наблюдений. (Методическое пособие для магистрантов). Новосибирск: Изд-во НГТУ, 2004. 31 с.

81. Mehra R. Identification and adaptive Kaiman filtering // Mechanics. 1971. № 3. P. 34-52.

82. Абденова Г.А. Адаптивное прогнозирование и фильтрация временного ряда с помощью моделей пространства состояний // Сборник науч. трудов НГТУ. 2010. № 1(59). С. 35-44.

83. Гельфанд A.M., Хмельник С.И. Дискретная фильтрация многомерных коррелированных нестационарных процессов // Доклады независимых авторов. 2008. Вып. № 7. С.6-18.

84. Передумов В.П. Об устойчивости моделирования на ЦВМ линейных динамических систем // Вестник Киев. Политехи. Ин-та.: Техн. кибернетика. 1983. Вып. 7. С. 39-41.

85. Абденова Г.А. Структурно-параметрическая идентификация систем с распределенными параметрами с использованием модели типа «вход-состояниевыход» // Научный вестник НГТУ. 2010. № 1. С. 9-16.

86. Лившиц К.И. Идентификация: Учебное пособие. Томск: изд-во Томск.ун-та, 1981. 132с.

87. Самарский A.A., Гулин A.B. Численные методы. М.: Наука, 1989. 430с.

88. Вапник В.Н. Восстановление зависимостей по эмпирическим данным. М.: Наука, 1979. 448с.

89. Абденова Г.А. Построение кусочно-дифференциальной модели в пространстве состояний для контроля качества продукции // Сборник науч. трудов НГТУ. 2009. № 3(57). С.47- 56.

90. Абденова Г.А., Рубанович М.Г., Хрусталев В.А. Оценка невязки экспериментальных данных и аналитических расчётов интегральной индуктивности мик-рополосковой линии // Научный вестник НГТУ, 2011. №3 (44). С.29-42

91. Абденова Г.А., Ильина Л.В., Раков М.А. Применение модели в пространстве состояний для достоверного оценивания прочности цементного камня // Фундаментальные науки. 2011. № 8. С. 589-592.

92. Абденова Г.А. Методика построения кусочно-разностной модели в пространстве состояний на основе нестационарного ряда наблюдений // Сб. науч. трудов НГТУ. 2011. № 2 (64). С. 3-12.

93. Mehra R.K. Optimal Input Signals for Parameter Estimation in Dynamic Systems -Survey and New Results // IEEE Trans. Automat. Control, 1974. Vol.19. №. 6. P.753-768.

94. Ченстоховского технологического университета, 2005. Т. 2. С. 57-59.

95. Воскобойников Ю.Е. Эконометрика в Excel: Учебное пособие. Новосибирск: НГАСУ. 2005.215с.

96. Соболь Б.В., Мишняков Н.Т., Поркшеян В.М. Практикум по высшей математике. Изд. 6-е. Ростов н/Д: Феникс, 2010. 630 с.

97. Тычинин А.В. Структурно-параметрический синтез системы граничного управления объектом с распределенными параметрами методами обратной динамики: автореф. дис. . канд. техн. наук по спец. 05.13.01. Самара. 2008. 22с.

98. Абденова Г.А., Апсалямов Н.А., Толысбаев Б.С. Прогнозирование уровней поставок зернового сырья на мукомольное производство // Известия международной академии наук высшей школы. Новосибирск, 2006. №1 (35). С.114-124.

99. Chow Gregory С. Tests of equality between sets of coefficients in two linear regressions //Econometrica. 1960.Vol. 28. № 3. P. 591-605.

100. Райбман H.C., Чадеев В.М. Построение моделей процессов производства. М.: Энергия, 1975. 375 с.

101. Деруссо П., Рой Р., Клоуз Ч. Пространство состояний в теории управления. М.: Наука, 1970. 620 с.

102. Сейдж Э.П., Мелса Дж.Л. Идентификация систем управления. М.: Наука, 1974.246 с.

103. Анжело Г.Д. Линейные системы с переменными параметрами. Пер. и ред. Н.Т. Кузовкова. М.: Машиностроение, 1974. 344 с.

104. Михайлов Ф.А., Теряев Е.Д., Булеков В.П. и др. Динамика нестационарных линейных систем. М: Наука, 1967. 344 с.

105. Фильтрация и стохастическое управление в динамических системах / Под ред. Леондеса К.Т. М.: Мир, 1980. 408 с.

106. Guo Lingzhong, Billings Stephen A. Identification of partial differential equation models for continuous spatio-temporal dynamical systems // IEEE Trans. Circuits and Syst. 2006. Vol.53. №8. P. 657-661.

107. Абденова Г.А. К проблеме анализа характеристик моделей систем в формепространства состояний // Сб. науч. трудов НГТУ. 2010. №3 (61). С. 3-10.

108. Абденова Г.А., Воевода A.A. Оценивание параметров и характеристик шумов нестационарных процессов в стохастических системах, описываемых в пространстве состояний // Сб. науч. трудов НГТУ. 2010. №3 (61). С.11-18.

109. Абденова Г.А. Прогнозирование значений уровня временного ряда на основе уравнений фильтра Калмана // Ползуновский вестник. Барнаул: Изд-во АлтГТУ, 2010. № 2. С. 4-6.

110. Абденова Г.А., Мальцев A.C. Кусочно-дифференциальная модель для полигона и кумулянты данных измерений относительно длин выпущенной продукции // Вестник Алт. гос. техн. ун-та им. И.И. Ползунова. Барнаул: Изд-во АлтГТУ, 2006. №2. С.157-161.

111. Абденова Г.А., Воевода A.A. Идентификация непрерывных многомерных систем по дискретным наблюдениям // Сборник науч. трудов НГТУ. Новосибирск: Изд-во НГТУ, 2010. № 2 (60). С. 165-168.

112. Воевода A.A., Трошина Г.В. О масштабировании данных «вход-выход» при идентификации объектов // Сборник науч. трудов НГТУ. Новосибирск: Изд-во НГТУ, 2010. №3. С.3-10.

113. Girshich М.А., Rubin H.A. Bayes approach to a quality control model // Ann. Math.Statist., 1952. Vol. 23. № 1. P. 183-192.

114. Ширяев A.H. Статистический последовательный анализ. M.: Наука,1972. 272 c.

115. Абденова Г.А., Воевода A.A. Обусловленность информационной матрицы Грама в задаче идентификации: масштабирование входных и выходных сигналов // Доклады ТУ СУР. Томск: Изд-во ТУ СУР, 2012. № 1(25). Часть 1. С.131-135.

116. Воробейчиков С.Э., Кабанова Т.В. Обнаружение момента разладки последо* f.вательности независимых случайных величин // Радиотехника и электроника, 2002. Т. 47. №10. С.1198-1203.

117. Page E.S. Continuous Inspection Schemes // Biometrika, 1954. №1. P. 141-154.

118. Форсайт Дж., Малькольм M., Маулер В. Машинные методы математических вычислений. М.: Мир, 1980. 279 с.

119. Бердов Г.И., Ильина Л.В., Машкин H.A. Влияние волластонита на прочность цементного камня из длительно хранившегося портландцемента // Строительные материалы. 2011. № 1. С. 48 49.

120. Бердов Г.И., Ильина Л.В. Активация цементов действием минеральных добавок // Международный журнал прикладных и фундаментальных исследований. 2010. №9. С.55-58.

121. Heinhold J., Gaede K.-W., Ingenieur-Statistic, München/Wien, 1964. p. 352.

122. Абденова Г.А., Воевода A.A. Один алгоритм перехода от регрессионной модели к модели в форме пространства состояний // Сборник науч. трудов НГТУ. Новосибирск: Изд-во НГТУ, 2012. № 3 (69). С.3-10.

123. АЛГОРИТМ РАСЧЕТА КОВАРИАЦИОННЫХ МАТРИЦ ДЛЯ СТОХАСТИЧЕСКОЙ МОДЕЛИ В ФОРМЕ КОШИ

124. АЛГОРИТМ РАСЧЕТА КОВАРИАЦИОННЫХ МАТРИЦ ДЛЯ СТОХАСТИЧЕСКОЙ МОДЕЛИ В ФОРМЕ КОШИ

125. Е{х(0) } = % Е{(х(0) х0 )(х(0) -х0)Т} = Р(0).

126. Рассматриваемая далее процедура идентификации параметров предполагает малый уровень требований к исходной информации.

127. Обозначив функционал для і -го уравнения через Jj, i = l,n, имеемj, = ріп Nf(tfi\(0 + tg{F«k(t)-x,(t+l) +1. Ji\8i\ t=Q J=і k=1finSir1. E44(04(0+SAÄ(0)27=1 k=l

128. После последовательных процедур приближений с помощью соотношения (П.1.5) и восстановления динамических параметров в, можно оценить дисперсии помех динамики. Для этого соотношение (П. 1.4) записывается в виде7=1 к=1

129. Далее, исходя из того, что гауссовский белый процесс, можно записатьследующую предварительную оценку для ковариационной матрицы помех динамики 1.:1. N-11. Qu=-rlLwKtk\ і=Хп.1. N к=о1. П.1.6)

130. Р0ш=рр0крг+д0, Р® =0, к = 1^20. (П. 1.7)

131. Надо полагать, что полученные оценки динамических параметров в и ковариационных матриц Я0,£>0,Р0 могут быть далеки от истинных значений. Поэтому необходимо провести некоторые повторные вычисления для уточнения всех начальных оценок.

132. На рисунке П. 1.1 показано множество значений для р к ., к > 0, і = 1,2 , и &=1,.,40 лежащих внутри и вне полосы оптимальности фильтра ±1.961{ЫХ12}. Для определения факта оптимальности или субоптимальности фильтра, объем выборки был взят равным N = 500.

133. Рисунок П. 1.1. Нормализованные оценки \р Др/ = 1,2 и «коридор»оптимальности фильтра

134. На рисунке П.1.2 и рисунке П.1.3 даны графики вектора наблюдения (штрих-пунктирная линия) и сглаживающего кубического сплайна (непрерывная линия) для объема выборки N = 100.

135. Рисунок П.1.2. Графики вектора наблюдения ух и сглаживающегокубического сплайна

136. Рисунок П. 1.3. Графики вектора наблюдения у2 и сглаживающегокубического сплайна уз 2

137. На рисунках П. 1.4 и П. 1.5 даны графики фильтрационной оценки вектора состояния системы (непрерывная линия) и соответствующий вектор измерения модели (штрих-пунктирная линия) для объема выборки N = 100.

138. Рисунок П. 1.4. Графики вектора наблюдения и фильтрационнойоценки вектора состояния х^

139. Рисунок П. 1.5. Графики вектора наблюдения у2 и фильтрационнойоценки вектора состояния х2

140. В таблицу П. 1.1 сведены численные результаты пассивной идентификации по вышеизложенному алгоритму.

Похожие работы

Информатика, вычислительная техника и управление
05.13.00