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

• Черникова, Оксана Сергеевна
• кандидат технических наук
• 2007

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

ВВЕДЕНИЕ.

1. ПОСТАНОВКА ЗАДАЧИ ИССЛЕДОВАНИЯ.

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

2.1. ОЦЕНИВАНИЕ НЕИЗВЕСТНЫХ ПАРАМЕТРОВ.

2.1.1. Критерии идентификации для стохастических линейных дискретных систем.

2.1.2. Алгоритмы вычисления значения и градиентов критериев идентификации.

2.1.3. Выводы.

2.2. ПЛАНИРОВАНИЕ ОПТИМАЛЬНЫХ D- И А- ВХОДНЫХ СИГНАЛОВ.

2.2.1. Основные свойства информационной матрицы и D- и А- оптимальных планов.

2.2.2. Прямая градиентная процедура построения непрерывных оптимальных планов.

2.2.3. Двойственная градиентная процедура построения непрерывных планов.

2.2.4. Выводы.

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

3.1. ОЦЕНИВАНИЕ НЕИЗВЕСТНЫХ ПАРАМЕТРОВ.

3.1.1. Критерии идентификации для стохастических линейных дискретных систем.

3.1.1. Выводы.

3.2. ПОСТРОЕНИЕ ИНФОРМАЦИОННОЙ МАТРИЦЫ ОДНОТОЧЕЧНОГО ПЛАНА.

3.2.1. Вывод информационной матрицы одноточечного плана.

3.2.2. Нахождение установившихся значений ковариационной матрицы ошибки одношагового предсказания, ее первой и второй производной.

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

3.2.3. Алгоритм вычисления информационной матрицы одноточечного плана.

3.2.4. Выводы.

3.3. ПЛАНИРОВАНИЕ ОПТИМАЛЬНЫХ D- И А- ВХОДНЫХ

СИГНАЛОВ.

3.3.1. Основные свойства информационной матрицы и D- и А- оптимальных планов.

3.3.2. Прямая градиентная процедура построения непрерывных планов.

3.3.3. Двойственная градиентная процедура построения непрерывных планов.

3.3.4. Выводы.

4. АКТИВНАЯ ИДЕНТИФИКАЦИЯ НА ПРИМЕРЕ РАЗЛИЧНЫХ ДИНАМИЧЕСКИХ МОДЕЛЕЙ.

4.1. МОДЕЛИРОВАНИЕ ВЫБОРКИ ИЗМЕРЕНИЙ.

4.2. НЕКОТОРЫЕ МОДЕЛИ ДИНАМИЧЕСКИХ СИСТЕМ.

4.2.1. Модель бурильной машины.

4.2.2. Модель электрокардиостимулятора.И

4.3. Выводы.

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

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

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

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

На практике различают задачи идентификации в узком и широком смыслах [1]. При решении задач идентификации в широком смысле априорная информация о системе либо незначительна, либо вообще отсутствует. Для идентификации подобной системы необходимо решение ряда дополнительных задач, к которым можно отнести, например, выбор класса модели или определение её типа.

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

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

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

Активная идентификация предполагает подачу на вход изучаемой системы специально синтезированного пробного сигнала. Несмотря на трудности, связанные с необходимостью нарушения технологического режима, методы активной идентификации систем предоставляют экспериментатору г значительно большие возможности в получении качественной модели по сравнению с методами пассивной идентификации. Эти возможности обусловлены самой идеологией активной идентификации, базирующиеся на сочетании приемов параметрического оценивания с концепцией планирования эксперимента (см., например, [6-11]).

Собственно процедура активной идентификации состоит из следующих этапов:

• вычисление оценок параметров по измерительным данным, соответствующим некоторому пробному сигналу (пассивная идентификация);

• синтез на основе полученных оценок оптимального по некоторому критерию сигнала (планирование);

• регистрация выходных данных системы, отвечающих поданному синтезированному сигналу;

• пересчет оценок неизвестных параметров.

В случае необходимости этапы 2-4 повторяются.

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

Рис В. 1. Контур идентификации системы

Таким образом, процедура активной идентификации предполагает выполнение следующих основных этапов:

1. Подготовка данных наблюдений

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

2. Определение структуры математической модели

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

Существуют различные подходы к форме задания идентификационных моделей, а также ситуации, требующие соответствующих методов исследования. Математические модели можно классифицировать по разным признакам, например, в зависимости от учета фактора времени, уровня неопределенности, по типу используемого математического аппарата и т.д. Более определённо, выделяют линейные и нелинейные, стационарные и нестационарные, дискретные, непрерывные и непрерывно-дискретные системы. Другой вариант классификации предусматривает возможность идентификации в терминах детерминированных или стохастических процессов. И, наконец, модели можно задавать в пространстве состояний, в виде ARMAX - моделей (моделей авторегрессии с детерминированными управляющими сигналами, трендами и случайными возмущениями), моделей в виде передаточных функций.

Описание динамической системы в терминах пространства состояний [12-14] позволяет учесть имеющиеся физические представления о механизмах работы системы. В отличие от передаточных функций, которые используются, в основном, при описании моделей линейных стационарных систем, методы пространства состояний позволяют создать компактную форму представления любых систем: линейных и нелинейных, стационарных и нестационарных, непрерывных и дискретных. Используемый при этом математический аппарат позволяет создавать мощные программные средства для анализа и синтеза динамических систем.

3. Оценивание параметров, входящих в модель

На этом этапе требуется найти значения параметров по имеющимся экспериментальным данным, что является задачей параметрической идентификации [11, 15-18]. При решении задачи оценивания методами параметрической идентификации выбирается некоторый критерий качества, явно зависящий от вектора неизвестных параметров.

Для нахождения оценок неизвестных параметров используют различные методы оценивания, среди которых наиболее распространёнными являются метод максимального правдоподобия [18-21], метод наименьших квадратов и его обобщения [22-26], метод стохастической аппроксимации [1,2] и другие.

4. Планирование эксперимента

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

5. Проверка адекватности модели

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

С начала 70-х годов начинает интенсивно развиваться теория идентификации динамических моделей. Современный этап развития теории и практики параметрической идентификации характеризуется все более широким распространением идей и методов активной идентификации. Этому способствовали, в частности, труды А.Ж. Абденова, Ю.П. Адлера, В.Г. Горского, В.И. Денисова, И.Л. Еланцевой, Г.К. Круга, А.А. Попова, Ю.А. Сосулина, A.M. Талалая, В.М. Чубича, ВА.Фатуева в нашей стране и Дж.К. Гудвина (G.C. Goodwin), М.Б. Зейропа (М.В. Zarrop), Р.К. Мехры (R.K. Mehra ), P.JI. Пейна (R.L. Payne) за рубежом.

Р.К. Мехра [28-30], по-видимому, первым стал использовать идеи и методы теории планирования оптимального эксперимента для идентификации динамических систем, описываемых моделями в пространстве состояний, во временной и частотной областях. Для многомерных стохастических линейных дискретных систем он получил выражение для информационной матрицы Фишера, соответствующее случаю, когда подлежащие оцениванию параметры входили только в матрицы состояния и управления. Мехра доказал, что при переходе от статических детерминированных систем к стохастическим линейным дискретным системам свойства нормированных информационных матриц планов не изменяются и остаются справедливыми теоремы эквивалентности для D- и А- оптимальных планов. В частотной области Р.К. Мехра рассмотрел в асимптотике основные свойства информационной матрицы и передоказал для этого случая теоремы эквивалентности.

В дальнейшем А.Ж. Абденов в [31] продолжил исследование вопросов активной идентификации многомерных стохастических динамических систем, описываемых моделями в пространстве состояний во временной и частотной областях. В соавторстве А.А. Поповым он уточнил в [32] выражение для информационной матрицы, полученное Р.К.Мехрой в [28], и затем в статье с

В.И. Денисовым в [33] обобщил полученный в [32] результат на случай, когда неизвестные параметры входили не только в матрицы состояния и управления, но и в матрицу измерения. Отметим, что в [28,32,33] выражения информационных матриц Фишера были получены с использованием фильтра Калмана в форме нормированной обновляющей последовательности. Другой подход к вычислению информационной матрицы Фишера был предложен А.А. Поповым в статье [34]. Как и в [32], он рассмотрел случай вхождения неизвестных параметров в матрицы состояния и управления, но запись уравнений фильтра Калмана в форме ненормированной обновляющей последовательности позволила А.А. Попову оптимизировать процедуру вычисления информационной матрицы за счет исключения необходимости извлечения квадратного корня из матрицы для каждого дискретного момента времени.

В работах [35,36] А.Ж. Абденов исследовал теоретические и прикладные аспекты применения D- оптимального планирования для многомерных стохастических линейных дискретных и непрерывно-дискретных систем, при этом как и в [28-30] он рассмотрел входной сигнал случайной природы и вхождение параметров в матрицы состояния и управления.

В.М. Чубич в [37,38] исследовал структуру информационной матрицы из [32] и на этой основе разработал прямые и двойственные градиентные процедуры планирования D- и А- оптимальных планов. В последствии В.М. Чубич в работе [39] обобщил выражение для информационной матрицы во временной области [34] на случай, когда подлежащие оцениванию параметры в различных комбинациях содержались в матрицах состояния, управления, возмущения, измерения, в начальных условиях и в ковариационных матрицах помех динамики и ошибок измерений и предложил соответствующий алгоритм ее вычисления. Изложение вопросов активной идентификации для стохастических линейных дискретных систем, описываемых моделями в пространстве состояний и ARMAX-моделями с единой методологической точки зрения, было приведено В.И. Денисовым, совместно с В.И. Чубичем и И.Л. Еланцевой в [40].

Вопросы активной идентификации математических моделей на основе передаточных функций обсуждались в работах [8,9]. Отметим также, что из публикаций по идентификации моделей стохастических динамических систем в форме передаточных функций к работам [28-30,35,36], пожалуй, наиболее близки в методологическом отношении работы [41, 42]. Здесь дан практически исчерпывающий анализ задачи синтеза оптимальных сигналов и приведены алгоритмы ее численного решения.

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

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

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

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

1. Получено выражение для информационной матрицы одноточечного плана в частотной области.

2. Разработан алгоритм вычисления информационной матрицы одноточечного плана в частотной области.

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

4. Сформулированы и доказаны основные свойства информационных матриц и теоремы эквивалентности D- и А- оптимальных планов во временной и частотной областях.

5. На основе прямого и двойственного подхода разработаны оригинальные градиентные алгоритмы синтеза D- и А- оптимальных входных сигналов во временной и частотной областях.

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

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

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

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

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

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

Апробация работы. Основные результаты диссертационной работы опубликованы в ряде статей, докладывались и обсуждались на следующих международных и всероссийских конференциях: 7th Korea-Russian International

Symposium on Science and Technology at the University of Ulsan (Корея, Ульсан, th • 2003); 8 Korea-Russian International Symposium on Science and Technology

Россия, Томск, 2004); 9th Russian-Korea International Symposium on Science and

Technology (Россия, Новосибирск, 2005); всероссийская научная конференция молодых ученых «Наука. Технологии. Инновации» (Россия, Новосибирск,

2006), а также на научных семинарах кафедры прикладной математики

Новосибирского государственного технического университета (НГТУ).

Публикации. По материалам диссертации опубликовано 10 печатных работ, в том числе: 4 - в ведущих научных журналах и изданиях, входящих в перечень, рекомендованный ВАК РФ, 2 - в сборниках научных трудов, 4 - в трудах и материалах Международных и Российской конференций.

Диссертационная работа выполнялась в соответствии с планами НИР, проводимых в 2002-2006 гг. (Отчеты по НИР: «Математическое моделирование многофакторных объектов на основе наблюдений» и «Моделирование статических и динамических многофакторных объектов стохастической природы с исследованием вероятностных закономерностей» под руководством д.т.н., проф. В.И. Денисова), а также при финансовой поддержке Министерства образования и науки РФ в рамках программы «Развитие научного потенциала высшей школы» (код проекта РНП.2.1.2.43).

Структура диссертации. Диссертационная работа состоит из введения, четырех глав основного содержания, заключения, списка использованных источников из 79 наименований и приложений. Общий объем диссертации составляет 150 страниц, включая 17 таблиц и 16 рисунков.

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

4.3 Выводы

1. В зависимости от выбранных модельных структур, критериев идентификации и оптимальности, качество оценивания удалось улучшить от 2.6 (см., например, табл. 4.8) до 10.2 (см., например, табл. 4.14) раз в пространстве параметров и от 1.7 (см., например, табл. 4.8) до 6.3 (см., например, табл. 4.3) раз в пространстве откликов, что позволяет говорить о целесообразности и эффективности применения процедур активной идентификации как во временной, так и в частотной областях.

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

3. В силу того, что для одних модельных структур предпочтительнее оказался критерий А- оптимальности (см., например, табл.4.13, табл.4.14 ), а для других - критерий D-оптимальности (см., например, табл. 4.3, табл. 4.4), экспериментатору можно порекомендовать использовать процедуры планирования по различным критериям оптимальности.

ЗАКЛЮЧЕНИЕ

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

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

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

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

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

Список литературы диссертационного исследования кандидат технических наук Черникова, Оксана Сергеевна, 2007 год

1. Справочник по теории автоматического управления / Красовский А.А. -М.: Наука, 1987.-712 с.

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

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

4. Гроп Д. Методы идентификации систем. М.: Мир, 1979. - 302 с.

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

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

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

8. Круг Г.К., Сосулин Ю.А., Фатуев В.А. Планирование эксперимента в задачах идентификации и экстраполяции. М.: Наука, 1977. - 207 с.

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

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

11. Льюнг Л. Идентификация систем: теория для пользователей М.: Наука, 1991.- 431 с.

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

13. Казаков И.Е. Статистическая теория систем управления в пространстве состояний. М.: Наука, 1975. - 432 с.

14. Заде Л., Дезоер Ч. Теория линейных систем. Метод пространства состояний. М: - Наука, 1970. - 204 с.

15. Рао С. Линейные статистические методы и их применения. М.: Наука,1968.-547 с.

16. Дейч A.M. Методы идентификации динамических объектов. М.: Энергия, 1974.-240 с.

17. Крамер Г. Математические методы статистики. М.: Мир, 1975. - 648 с.

18. Спиди К., Браун Р., Гудвин Дж. Теория управления: идентификация и оптимальное управление. М.:Мир, 1973. - 248 с.

19. Боровков А.А. Математическая статистика: оценка параметров. Проверка гипотез: учеб. пособие для математических и физических спец. вузов. -М.: Наука, 1984.-472 с.

20. Клепиков Н.П., Соколов С.Н. Анализ и планирование экспериментов методом максимума правдоподобия. М.: Физматгиз, 1964. - 184 с.

21. Уилкс С. Математическая статистика. М.: Наука, 1967. - 632 с.

22. Линник Ю.В. Метод наименьших квадратов и основы математико -статистической теории обработки наблюдений. М.: Физматгиз, 1958. -334 с.

23. Айвазян С.А., Енюков И.С., Мешалкин Л.Д. Прикладная статистика. Исследование зависимостей. М.: Финансы и статистика, 1985. - 487 с.

24. Дрейпер Н., Смит Г. Прикладной регрессионный анализ. М.: Финансы и статистика, 1986. - 366 с.

25. Справочник по прикладной статистике: в 2 томах. T.l/Э.Ллойд, У. Ледерман, Ю. Тюрин. М.: Финансы и статистика, 1989. - 510 с.

26. Лоусон Ч., Хенсон Р. Численное решение задач методом наименьших квадратов. М.: Наука, 1986. - 232 с.

27. Денисов В.И., Попов А.А. Пакет программ оптимального планирования эксперимента. М.: Финансы и статистика, 1986. - 159 с.

28. Mehra R.K. Optimal Input Signals for Parameter Estimation in Dynamic Systems Survey and New Results // IEEE Trans. Automat. Control. - 1974. -v.AC - 19. - No. 6. - P. 753 - 768.

29. Абденов А.Ж. Активная идентификация для стохастических динамических систем, описываемых моделями в пространстве состояний: дис.д-ра техн.наук: 05.13.01 / Абденов Амирза Жакенович. -Новосибирск, 1999,- 377 с.

30. Абденов А.Ж., Попов А.А. Планирование D оптимальных входных воздействий при идентификации линейных систем / А.Ж. Абденов, А.А. Попов; Новосиб. электротехн. ин - т. - Новосибирск, 1981. - 12 с. - Деп. в ВИНИТИ 19.02.82, №771 - 82.

31. Попов А.А. Вычисление информационной матрицы Фишера в задаче планирования входного сигнала для динамических систем //Сб. науч. тр. НГТУ. 1998. - №2(11). - С. 8 - 16.

32. Абденов А.Ж. Активная идентификация линейных динамических систем для решения задач калмановской фильтрации. I. Теоретические и алгоритмические аспекты // Науч.вест. НГТУ. 1998. - №4. - С. 3-15.

33. Абденов А.Ж. Активная идентификация линейных динамических систем для решения задач калмановской фильтрации. I. Практические аспекты // Науч.вест. НГТУ. 1998. - №4. - С. 16-21.

34. Чубич В.М. Планирование D- оптимальных входных сигналов для стохастических линейных дискретных систем: дис.канд. техн.наук: Спец.05.13.16 / Чубич Владимир Михайлович; Новосиб. гос .техн. ун-т.-Новосибирск, 1995. 97 с.

35. Денисов В.И., Чубич В.М. А-оптимальное планирование на примере одной стохастической линейной дискретной системы: материалы междунар. научно-технической конференции «Информатика и проблемы телекоммуникаций». Новосибирск, 1997. - С. 112-117.

36. Рябых О.С. Новое обобщенное выражение для информационной матрицы Фишера в задаче планирования входного сигнала для динамических систем/ Денисов В.И., Чубич В.М., Рябых О. С. // Науч. вестник НГТУ. -2001.-№2(11).-С. 29-42.

37. Денисов В.И., Еланцева И.Л., Чубич В.М. Активная идентификация стохастических линейных в пространстве состояний и ARMAX -моделями // Сиб.журн.индустр. матем. 2000. - №1(5). - С. 87-100.

38. Goodwin G.C., Payne R.L. Dynamic system identification: experimental design and data analysis. New York: Academic Press, 1977. - 273 p.

39. Zarrop M.B. Optimal experimental design for dynamic system identification. -New York : Springer-Verlag., 1979 197 p.

40. Gupta N.K., Mehra R.K. Computational Aspects of Maximum Likelihood Estimation and Reduction in Sensitivity Function Calculation // IEEE Trans. Automat.Contr. 1974.-vol. AC-19 - P. 774-783.

41. Astrom К.J. Maximum Likelihood and Prediction Error Methods // Automatica. 1980. - v. 16. - P. 551 -574.

42. Огарков M.A. Методы статистического оценивания параметров случайных процессов. М.:Энергоатомиздат, 1980.- 206 с.

43. Абденов А.Ж., Денисов В.И., Чубич В.М. Введение в оценивание и планирование для стохастических динамических систем. Новосибирск: НГТУ, 1993.-45 с.

44. Налимов В.В., Чернова Н.А. Статистические методы планирования экстремальных экспериментов. М.:Наука, 1965. -340 с.

45. Налимов В.В. Теория эксперимента. М.: Наука, 1971. - 208 с.

46. Адлер Ю.П., Маркова Е.В., Грановский Ю.В. Планирование эксперимента при поиске оптимальных условий. М.: Наука, 1976. -279 с.

47. Налимов В.В., Голикова Т.Н. Логические основания планирования эксперимента. М: Металлургия, 1976. - 128 с.

48. Базара М., Шетти К. Нелинейное программирование. М.: Мир, 1982. -583 с.

49. Васильев В.П. Численные методы решения экстремальных задач. -М.:Мир, 1982.-372 с.

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

51. Моисеев Н.Н., Иванилов Ю.П., Столярова Е.М. Методы оптимизации. -М.: Наука, 1980.-351 с.

52. Поляк Б.Т. Введение в оптимизацию . М.: Наука, 1983. - 384 с.

53. Сухарев А.Г., Тимохов В.В., Федоров В.В. Курс методов оптимизации. -М.: Наука, 1986.-328 с.

54. Денисов В.И., Чубич В.М. Планирование D-оптимальных управляющих сигналов для стохастических линейных дискретных систем // Научный вестник НГТУ. 1995. - №1. - С.17-31.

55. Черникова О.С. Активная идентификация стохастических линейных дискретных систем в частотной области/ Денисов В.И., Чубич В.М., Черникова О.С.// Сиб.журн.индустр. матем.- Том 10Х.- №1(29). 2007. -С.71-89.

56. Кули Дж. У., Льюис П.А.У., Уэлч П. Д. Быстрое преобразование Фурье и его применения к анализу временных рядов // Стат. методы для ЭВМ/ Энслейн К., Рэлстон Э., Уилф Г.С. М.:Наука, 1986. - С. 373-434.

57. Деммель Дж. Вычислительная линейная алгебра. Теория и приложения. -М.:Мир, 2001. 430 с.

58. Васюков В.Н. Цифровая обработка сигналов и сигнальные процессоры в системах подвижной радиосвязи. Новосибирск : НГТУ, 2006. - 292 с.

59. Круглов В.В. Математические пакеты расширения MATLAB: Специальный справочник. СПб: Питер, 2001. - 480 с.

60. Дьяконов В.П., Абраменкова И.В. MATLAB 5.0/5.3. Система символьной математики. М.:Нолидж, 1999. - 640 с.

61. Черникова О.С. Эквивалентность выражений для информационных матриц Фишера в задаче активной идентификации стохастических линейных дискретных систем в частотной области //Сборник научных трудов НГТУ. 2006. - №4(46). - с. 11-18.

62. Гилл Ф., Мюррей Ч., Райт М. Практическая оптимизация. М.: Мир, 1982

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

64. Черникова О.С. Моделирование выборки измерений по описанию стохастической линейной стационарной системы с дискретным временем в пространстве состояний/ Денисов В.И., Чубич В.М., Черникова О.С. //Науч. вестник НГТУ. 2002. - №2(13). - С. 64-70.

65. Мюллер П., Нойман П., Шторм Р. Таблицы по математической статистике. М.: Наука, 1982. - 272 с.

66. Ермаков С.М., Михайлов Г.А. Статистическое моделирование. М.: Наука, 1982.-296 с.

67. Деммель Дж. Вычислительная линейная алгебра. Теория и приложения. -М.: Мир, 2001.-430 с.

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

69. Бегун П.И., Шукейло Ю.Л. Биомеханика. Спб.: Политехника, 2000. -463 с.

70. Баранович В.Ю., Таричко Ю.В. Лечение аритмий: учебное пособие. -РУДН, 2006- 121 с.