Построение оптимальной модели измерений для линейных динамических систем тема диссертации и автореферата по ВАК РФ 05.13.17, кандидат технических наук Самочернов, Игорь Валентинович

• Самочернов, Игорь Валентинович
• кандидат технических наук
• 2004

Оглавление диссертации кандидат технических наук Самочернов, Игорь Валентинович

Введение

Основные определения и постановка задачи

1.1. Математическая модель динамической системы.

1.2. Исследования рассматриваемых динамических систем

1.3. Постановка задачи

• * г

1.4. Выводы

Алгоритмы решения задач оптимизации наблюдателя

2.1. Оптимизация модели измерений для улучшения оценок состо

2.1.1. Информационная матрица в качестве основы критерия оптимальности

2.1.2. Использование ковариационйой матрицы

2.1.3. Связь между ковариационной и информационной матрицами

2.1.4. Выбор между моделями измерений.

2.1.5. Независимая от времени оптимальная модель.

2.1.6. Планирование моментов измерений

2.2. Оптимизация измерений для оценки неизвестных параметров системы.

2.3. Выводы.

Моделирование векторов состояний и наблюдений

3.1. Постановка задачи.

3.2. Общее решение уравнения состояний.

3.2.1. Вычисление переходной матрицы состояний.

3.2.2. Вычисление ковариационной матрицы.

3.2.3. Моделирование нормально распределенных векторов

3.2.4. Алгоритм моделирования.

3.3. Методы Эйлера решения стохастических дифференциальных уравнений. . . '

3.3.1. Описание методов

3.3.2. Алгоритм моделирования, основанный на методах Эй

3.4. Методы Мильштейна

3.5. Исследования на тестовом примере.

3.6. Выводы .51'

Алгоритмы поиска оптимального наблюдателя

4.1. Постановка задачи.

4.2. Полностью стохастический и смешанный подходы.

4.2.1. Выбор начального приближения для метода локального поиска.

4.2.2. Локальный поиск

4.3. Детерминированный подход .:.

4.3.1. Алгоритм отображения отрезка на многомерный гипер

4.3.2. Алгоритм одномерной глобальной оптимизации

4.4. Продолжение примера

4.5. Выводы

Исследование динамических систем

5.1. Модель системы чандлеровских колебаний

5.1.1. Моделирование вектора состояний.

5.1.2. Оптимальные модели измерений.

5.1.3. Проверка результатов

5.2. Следящая система управления электроприводом постоянного тока.

5.2.1. Получение уравнения состояний.

5.2.2. Переходная матрица состояний.

5.2.3. Моделирование вектора состояний.

5.2.4. Выбор оптимальной модели наблюдений.

5.2.5. Выбор частоты проведения измерений

5.3. Система стабилизации самолета по тангажу.

5.3.1. Построение оптимальной модели наблюдений

5.4. Выводы

Описание программной системы

6.1. Ядро комплекса

6.1.1. covcalc.exe.

6.1.2. dssolver.exe

6.1.3. kalman.exe.

6.1.4. optimizer.exe

6.2. Вспомогательные программы.

6.3. Поддержка проведения вычислительных экспериментов

Введение диссертации (часть автореферата) на тему «Построение оптимальной модели измерений для линейных динамических систем»

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

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

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

Задача в такой постановке встречалась в работах Mehra R.K и ранее у Meier L., Jonson C.D.? причем Mehra R.K. впервые сделал попытку применить методы, используемые в планирования регрессионных экспериментов, для исследования динамических систем. О подобных исследованиях в нашей стране нам ничего не известно, поэтому в данной работе сделана попытка восполнить этот пробел.

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

Задачи исследования. Для достижения цели диссертационной работы решены следующие задачи: предложены и исследованы алгоритмы и критерии оптимальности моделей наблюдений с точки зрения качества оценок состояний; исследованы методы компьютерного моделирования реализаций наблюдений и состояний динамических систем; в применении к основной задаче реализованы и исследованы некоторые методы глобальной оптимизации.

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

Положения, выносимые на защиту.

1. Алгоритмы построения и критерии оптимальности моделей наблюдений.

2. Рекомендации по использованию пяти методов моделирования динамических систем.

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

4. Результаты построения оптимальных моделей наблюдений для исследованных систем: следящей системы управления электроприводом постоянного тока, системы управления самолетом по тангажу, системы чандлеровских колебаний.

5. Разработанное программное обеспечение по моделированию и оптимизации моделей наблюдений динамических систем.

Научная новизна. Разработаны алгоритмы построения оптимальных моделей наблюдений для оценки состояний в дискретных и непрерывнодискретных линейных динамических системах. Предложены критерии оптимальности модели наблюдений в виде функционалов от ковариационной матрицы ошибок оценивания и информационной матрицы Фишера о состояниях системы. Проведено сравнение ряда алгоритмов моделирования реализаций векторов состояний и измерений: стандартного и основанных на разложении Тейлора-Ито уравнения состояний. Сделано сравнение нескольких алгоритмов поиска экстремумов в применении к задаче оптимизации проведения измерений. Разработано соответствующее программное обеспечение. Получены оптимальные модели измерений для системы стабилизации самолета по тангажу, системы управления электроприводом постоянного тока, системы чандлеровских колебаний.

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

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

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

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

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

Апробация работы. Основные результаты работы докладывались и обсуждались на Новосибирской межвузовской научной студенческой конференции "Интеллектуальный потенциал Сибири" (2000), конференции, посвященной дням науки НГТУ-2000 (исследования были поддержаны грантом университета), 6-й международной российско-корейской конференции KORUS-2002.

Публикации. По теме диссертационной работы опубликовано 8 работ. Из них 1 — в трудах международной конференции, 1 — в сборнике тезисов межвузовской конференции, 2 — в Научном вестнике НГТУ, 4 — в сборнике научных трудов НГТУ.

Структура и объем диссертации. Работа состоит из введения, шести глав основного текста и заключения. Объем работы — 125 страниц. Список литературы содержит 67 источников. Рисунков 23, таблиц 11.

Заключение диссертации по теме «Теоретические основы информатики», Самочернов, Игорь Валентинович

5.4. Выводы

По результатам исследований динамических систем можно сделать следующие выводы:

• Разработанные и реализованные алгоритмы позволяют построить модель измерений, использование которой значительно улучшает величину критерия оптимальности измерений динамической системы (от 1 до тысяч раз), величина абсолютной точности оценивания вектора состояний увеличивается до 3-х раз.'

• Линейная комбинация измерений оказывается эффективнее, чем несмешанные измерения с диагональной матрицей наблюдений if .

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

• Оптимальные модели наблюдений зависят не только от структуры динамической системы, но и от количества моментов измерений и интервала наблюдений.

• Оптимальные модели, полученные при помощи различных критериев оптимальности в общем случае различны.

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

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

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

Полностью стохастический метод поиска работает достаточно плохо на рассмотренных системах. Главным образом за сложности выбора между невысокой точностью результата и большим количеством вычислений. Двухэтапный и смешанный методы в большинстве случаев дают одинаковые результаты по критерию и моделям. При этом смешанный метод проигрывает двухэтапному по времени при малых размерностях пространства параметров (меньше 10) до 2-х раз, но начинает выигрывать до 3 и более раз при более высоких размерностях.

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

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

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

0. 0. 0. 0.

0.

13.00 26.00 39.00 52.00 Время, с а)

Ы ,2,4

6.

4.

2.

13.00 26.00 39.00 52.00 Время, с б)

Рис. 5.20. Параметры нестационарной модели наблюдений в зависимости от времени. Оптимальные модели по следу информационной матрицы (а) и ковариационной матрицы (б). Часть компонент

6. Описание программной системы

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

Комплекс состоит из трех частей:

1. Ядро, содержащее программные реализации алгоритмов моделирования, оптимизации моделей; •

2. Утилиты для отрисовки графиков, преобразования исходных данных и результатов, помощи в составлении отчетов;

3. Скрипты для поддержки автоматического выполнения групп вычислительных экспериментов над динамическими системами.

Все программы работают из командной строки операционной системы. Написаны на языках С++ [59] с поддержкой объектно-ориентированного подхода, на языке Perl [64] и языке командного интерпретатора. Для проверки работы разработанных алгоритмов существуют дублирующие программы для среды математического программирования Maple [45, 53].

Достоинством выбранного подхода является отличная масштабируемость. Новые программы в систему могут быть добавлены: самым простейшим образом: просто добавлением в нужный каталог. При этом новые программы могут использовать все предыдущие наработки.

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

Объем 40-ка исходных файлов на С++ составляет 280 Кб. На скрипто-вых языках — более 30-ти файлов объемом свыше 15 Кб. Программы собраны для операционной системы Microsoft Windows1, но с незначительными изменениями (кроме модуля построения графиков) могут быть скомпилированы для работы в операционной системе Unix1.

6.1. Ядро комплекса

Ядро состоит из следующих программ:

1. covcalc.exe — программа вычисления значений различных критериев оптимальности для заданной модели системы;

2. dssolver.exe — моделирование реализаций динамической системы различными методами;

3. kalman.exe — оценка состояний, а также вычисление абсолютного отклонения оценки состояний по набору реализаций с усреднением результата;

4. optimizer.exe — оптимизация модели наблюдений различными методами и с использованием различных критериев.

Все эти программы пользуются одним и тем же классом динамической системы и используют его методы для достижения своих целей. Текст заголовочного файла с описанием приведен на листинге 6.2. Диаграмма основных классов ядра изображена на рис. 6.1.

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

Описание модели хранится в файле формата XML [23]. Пример описания модели простой системы приведен на листинге 6.1. Файл параметров также хранится в формате XML и его содержимое зависит от программы. Имена и значения параметров для разных методов не пересекаются, если только они не обозначают одно и то же, поэтому информация о них хранится в одном файле из которого каждый метод выбирает необходимые ему данные.

Все зарегистрированные торговые марки являются собственностью их владельцев

Пользовательский уровень

Смешанный метод оптимизации А utility» covcalc «utility» dssolver «utility» kalman

Средний уровень ч 1 ' Ч X / Ч Х ' Ч * 1 Ч * ' / \ X / / X \ ' ' Ч * ' ' Ч \ / / Ч(Л/ ч/и/ utility» optimizer л \ t I / I

Нижний

Powel 11

I 1 \ \ jjL

MIxedMlnimlzerl

Стохастический метод оптимизации

-о CompiexMIn

Двухэтапный метод оптимизации

Rav Xmllnterface TVector TMatrix 1

•

1 V 1

ExprList

Expr

1 •

Генератор случайных чисел

Тч

Вычисление ^ отдельных выражений

Поддержка работы ^ \ с файлами параметров / \ /

1 I

Управление вычислением выражений

Ъ,

Операции с матрицами и векторами К

Рис. 6.1. Диаграмма основных классов программной системы

Листинг 6.1: Пример файла описания динамической системы xml version»" 1.0" encoding="windows-1251 "?> <dynasystem tag="Осциллятор с шумом" id=" oscillator "> <matrix name="F" cols="2" rows="2"> 0-1 10 </matrix> eigenvalues matrix="F"> I

-I eigenvalues> matrix name="Phi" cols="2" rows="2" numeric=" false "> cos (t—tau) sin(tau—t) sin(t—tau ) cos(t—tau) </matrix> matrix name="C" rows="2" cols="2"> 0 0 0 0 </matrix> matrix name="G" cols="l" rows="2"> 0.17320508075688772935274463415059 0 matrix> matrix name="P0" cols="2" rows="2"> .10 0 .01-</matrix> vector name="X0" dim="2"> 1.5 0 vector> matrix name="R" cols="l" rows="l">

03 </matrix> matrix name="H" cols="2" rows="r'> 11 </matrix> </ dynasystem>

6.1.1. covcalc.exe

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

Аргументы командной строки: <имя файла с описанием системы> <имя файла параметров> <критерий>.

Пример файла с описанием системы приведен в листинге 6.1. Из файла параметров могут выбираться следующие параметры:

• План измерений; •

• Точность интегрирования/шаг по времени;

• Параметры исследования частоты проведения измерений.

Параметры критериев могут быть следующими: SpP, SpM, DetP, DetM, SpP(dt), SpM(dt), DetP(dt), DetM(dt). Первые 4 параметра служат для вычисления значений критериев, а последние позволяют определить зависимость значений критериев от частоты проведения измерений.

Также данная программа умеет определять требуемый шаг численного интегрирования для достижения заданной точности вычисления интегралов вида (3.3).

6.1.2. dssolver.exe

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

Аргументы командной строки: <имя файла с описанием системы> <имя файла параметров> <метод моделирования> <файл для векторов состоящий [файл для векторов наблюдений].

Из файла параметров берется план измерений и точность интегрирования/шаг по времени.

Метод моделирования может быть одним из следующих:

• equationSolve — общее решение уравнения системы;

• explicitEuler.— явный стохастический метод Эйлера;

• implicitEuler — неявный стохастический метод Эйлера;

• explicitMilstein — явный стохастический метод Милыптейна;

• implicitMilstein — неявный стохастический метод Милыптейна.

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

6.1.3. kalman.exe

Вычисляет оценки состояний и значение в соответствии с формулой (5.1).

Аргументы командной строки: <имя файла с описанием системы> <имя файла параметров> [выходной файл оценок] [входной файл состояний] [входной файл наблюдений].

Если входные файлы состояний и наблюдений не указаны, то вычисляется величина по количеству выборок 7V; указанных в файле параметров. Иначе в выходной файл записываются оценки состояний.

6.1.4. optimizer.exe

Заключение

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

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

2. Разработаны и исследованы алгоритмы выбора оптимальных моделей: последовательный алгоритм построения нестационарной модели, алгоритм получения стационарной оптимальной модели, алгоритм выбора модели из конечного множества моделей, алгоритм выбора оптимальной частоты проведения измерений. Полученные при помощи этих алгоритмов модели могут выигрывать у неоптимальных по значениям критерия до нескольких тысяч раз.

3. Исследованы алгоритмы моделирования реализаций динамических систем, основанные на классическом подходе, а также на методах Эйлера и Мильштейна для численной аппроксимации решения стохастических уравнений. Даны рекомендации для использования и улучшения качества алгоритмов моделирования.

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

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

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

7. Показано, что основными величинами одномерной оптимальной модели наблюдений следящей системы являются переменные напряжения тиристорного преобразователя и тока якорной цепи. Для моделей измерений более высоких размерностей требуется измерять различные комбинации по 3-4 компоненты состояний. Чем больше размерность модели, тем ближе отдельные параметры в модели наблюдений к границам допустимой области.

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

Список литературы диссертационного исследования кандидат технических наук Самочернов, Игорь Валентинович, 2004 год

1. Tom М. Apostol. Calculus. Volume II., Multivariable Calculus and Linear Algebra, with Applications to Differential Equations and Probability. John Wiley & Sons. 1969.

2. Atkinson A.C., Fedorov V.V. Optimal. Design: Experiments for Discriminating between Several Models. //Biometrika,' 1975, 62, №2, pp. 289-303.

3. Becareva N.D., Denisov V.I., Samochernov I N. About an algorithm of the observation matrix optimization. Abstracts of the 6-th Russian-Korean International Symposium on Science And Technology. Vol. 1, pp. 236-239, June 24-30, 2002.

4. T.J. Broida, S. Chandeashekhar, R. Chellappa. Recursive 3D Motion Estimation from a Monocular Image Sequence. 1996.

5. Richard S. Bucy, Peter D. Joseph. Filtering for Stochastic Processes with Applications to Guidance. New York: John Wiley & Sons, Inc., 1968.

6. Madalena Chaves, Eduardo D. Sontag. State-Estimators for Chemical Reaction Networks of Feinberg-Hom-Jackson Zero Deficiency Type. arXiv:math.OC/0012130 vl, 2000.

7. Deniz Erdogmus, Jose C. Principe, Geetha Thampi. Adaptive linear observer for nonlinear systems. Controlo, 2002.

8. William Fong, Simon Godsill, Arnaud Doucet, Mike West. Monte Carlo Smoothing with Application to Audio Signal Enhancement. 2001.

9. Handbook of Global Optimization II ed. by Panos M. Pardalos, H. Edwin Romeijn. Kluwer Academic Publishers. 2001.

10. Robert Shaw Harvey. Development of A Precision Pointing System Using an Integrated Multi-Sensor Approach. UCGE Reports Number 20117. 1998.

11. Joachim Heel. Dynamical Systems and Motion Vision. Massachusetts Institute of Technology, 1988.

12. Leslie Lamport. ИГ^Х: A Document Preparation System. Addison-Wesley, 1986.

13. Tine Lefebvre, Herman Bruyninckx, Joris De Shutter. Kalman filters for nonlinear systems: a comparison of performance. 2001.

14. L.W.B. Leite, M.P.C. da Rocha. Deconvolution of non-stationary seismic process. Brazilian Journal of Geophysics, Vol. 18(1), 2000.

15. X. Rong Li, Chongzhao Han, Jie Wang. Discrete-Time Linear Filtering in Arbitrary Noise. CDC2000, 2000.

16. Mehra R.K. Optimal input signals for parameter estimation in dynamic systems — Survey and new results. IEEE Trans. Automat. Contr. vol. AC-19, 1974.

17. Mehra R.K. Optimization of Measurement Schedules and Sensor Designs for Linear Dynamic Systems. IEEE, 1975.

18. William H. Press. Numerical recipes in C: the art of scientific computing. Cambridge University Press 1988, 1992.

19. Ali H. Sayed. A Framework for State-Space Estimation with Uncertain Models. IEEE Transactions on Automatic Control, vol. 46, no. 7, pp. 998.1013, July 2001.

20. R.G. Strongin, Y.D. Sergeyev, Global Optimization with Non-Convex Constraints: Sequential and Parallel Algorithms, Nonconvex Optimization and Its Applications. Kluwer Academic Publishers, Dordrecht, 2000, 728 pp.

21. Hisashi Tanizaki. Nonlinear and Non-Gaussian State-Space Modeling with Monte Carlo Techniques: A Survey and Comparative Study. 2000.

22. F. Wang, V. Balakrushnan. Robust Adaptive Kalman Filters for Linear Time-Varying Systems with Stochastic Parametric Uncertainties.23. http://www.w3.org/TR/1998/REC-xml-19980210.

23. Zarrop M.B. Optimal Experiment Design for Dynamic System Identification. New York:Springer-Verlag, 1979.

24. Абденов А.Ж., Денисов В.И., Чубич B.M. Введение в планирование экспериментов для стохастических динамических систем: Учеб. пособие / Новосиб. гос. техн. ун-т. — Новосибирск, 1993, 45 с.

25. Андреев Ю. Н. Управление конечномерными линейными объектами. М., "Наука", 1976.

26. Башарин А.В., Постников Ю.В. Примеры расчета автоматизированного электропривода на ЭВМ: Учебное пособие для вузов. — 3-е изд. — Л.: Энергоатомиздат. Ленингр. отд-ние, 1990. — 512 с.

27. Бекарева Н.Д., Парлюк А.В. Моделирование случайного процесса, заданного стохастическим дифференциальным уравнением. Сборник научных труцов НГТУ. — 1999. — №1.

28. Бекарева Н.Д., Самочернов И.В. Оптимизация модели измерений для линейных динамических систем. Сборник научных трудов НГТУ. — 2000. №5.

29. Бекарева Н.Д., Самочернов И.В. Исследование процедуры оптимизации модели измерений. Сборник научных трудов НГТУ. — 2001. — № 1.

30. Бекарева Н.Д., Самочернов И.В. О выборе между моделями наблюдений для динамических систем. Сборник научных трудов НГТУ. — 2003.-№1.

31. Бекарева Н.Д., Самочернов И.В. Выбор модели наблюдений для системы стабилизации самолета по тангажу. Сборник научных трудов НГТУ. 2003. - №2.

32. Березин И.С., Жидков Н.П. Методы вычислений. В 2-х томах. Физ-матгиз, 1959.

33. Браммер К., Зиффлинг Г. Фильтр Калмана-Бьюси. Пер. с нем. — М.: Наука. Главная редакция физико-математической литературы. 1982.

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

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

36. Денисов В .И. Математическое обеспечение системы ЭВМ — экспериментатор. — М.: Наука, 1977. — 251 с.

37. Д'Анжело Г. Линейные системы с переменными параметрами. Анализ и синтез. Под ред. Н. Т. Кузовкова. (Пер. с англ.). М., "Машиностроение", 1974.

38. Денисов В.И., Бекарева Н.Д., Самочернов И.В. Исследование свойств уравнения Риккати. Научный вестник НГТУ. — 2001. — № 1(10).

39. Денисов В.И., Бекарева Н.Д., Самочернов И.В. Сравнение нескольких методов моделирования состояний линейных динамических систем. Научный вестник НГТУ. — 2004. — № 1.

40. Денисов В.И., Чубич В.М. Алгоритмы синтеза планов экспериментов для стохастических динамических систем: Учеб. пособие. — Новосибирск: НГТУ, 1996. 36 с.

41. Дорф Р., Бишоп Р. Современные системы управления. — М.: Лаборатория базовых знаний, 2002, 832 с.т

42. Дьяконов В. Maple 6: Учебный курс. — СПб.: Питер 2001. — 608 с.

43. Ермаков С.М., Жиглявский А.А. Математическая теория оптимального эксперимента. — М.: Наука, 1987. — 320 с.

44. Ермаков С.М., Михайлов Г.А. Статистическое моделирование. 2-е изд., дополн. М.: Наука. Главная редакция физико-математической литературы, 1982. — 296 с.

45. Ивченко Г.И.,'Медведев Ю.И. Математическая статистика: Учеб. пособие для втузов. — М.: Высш. шк., 1984. — 248 с.

46. Климов А.С. Форматы графических файлов. К: НИПФ "ДиаСофт Лтд.", 1995' 480 с.

47. Котельников И., Чеботаев П. Издательская система LaTeX2e. Новосибирск: "Сибирский хронограф", 1998.

48. Крылов В.И., Бобков В.И., Монастырный П. И. Вычислительные методы. В 2-х т. Главная редакция физико-математической литературы изд-ва "Наука", М , 1977.

49. Кузнецов Д.Ф. Численное интегрирование стохастических дифференциальных уравнений. — Спб. Издательство С.-Петербургского Государственного университета, 2001.

50. Манзон Б.М. Maple V. Power Edition. — М.: Информационно-издательский дом "Филинъ" 1998. — 240 с.

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

52. Мороз А.И. Курс теории систем: Учеб. пособие для вузов по спец. у "Прикладная математика". — М.: Высш. шк., 1987. — 304 с.

53. Мудров А.Е. Численные методы для ПЭВМ на языках Бейсик, Фортран и Паскаль. — Томск: МП "Раско", 1991. — 272 с.

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

55. Парлюк А.В. Планирование дискриминирующего эксперимента для стохастических динамических систем: дисс. на соиск. учен. степ, канд. техн. наук (05.13.17). Новосиб. гос. техн. ун-т. 2001.

56. Ирэ Пол. Объектно-ориентированное программирование с использованием С++. К.: НИПФ "ДиаСофт Лтд:", 1995, — 480 с.

57. Пугачёв B.C., Синицын И.Н. Теория стохастических систем: Учеб. пособие для вузов по спец. "Прикладная математика". — М.: Логос, 2000,999с. '

58. Самочернов И.В. Оптимизация проведения измерений в линейных динамических системах. Сборник тезисов докладов Новосибирской межвузовской научной студенческой конференции "Интеллектуальный потенциал Сибири". — 2000.

59. Федоров В. В. Теория оптимального эксперимента (планирование регрессионных экспериментов). М., 1971, 312 с.t

60. Химмельблау Д. Прикладное нелинейное программирование. — М.: Мир, 1975.

61. Стивен Холзнер. Perl: специальный справочник. — СПб: Питер, 2000. -496с.

62. Чубич В.М. Планирование D-оптимальных входных сигналов для стохастических линейных дискретных систем: дисс. на соиск. учен. степ, канд. техн. наук (05.13.16). Новосиб. гос. техн. ун-т. 1995.

63. Штейнберг Ш.Е. Идентификация в системах управления. — М.: Энергоатомиздат, 1987.

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