«Моделирование систем»

1 Методы параметрической идентификации объектов управления.

2 Методы структурной идентификации объектов управления.

3 Методы математической обработки экспериментальной информации (регрессионный анализ).

4 Методы планирования эксперимента (полный факторный эксперимент).

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

1 Методы параметрической идентификации объектов управления.

структурную и параметрическую идентификацию.

На этапе параметрической идентификации выполняется экспериментальная проверка модели.

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

Используют критерии: модульный, квадратичный, показательный, минимаксный, взвешенный критерии. Задача сводится к оценке суммарной невязки, которая служит основным критерием, по нему проводится идентификация модели.

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

Методы параметрической идентификации

Методы различают в зависимости от модели.

Модели бывают:

1. Статические и динамические.

2. Детерминированные и стохастические.

3. Линейные и нелинейные.

4. Непрерывные и дискретные.

Идентификация делится:

1. Активные и пассивные методы.

2. Непрерывные и дискретные.

Параметрическая идентификация для статической детерминированной модели y = F (x )

Модель объекта линейная, имеет n входов, m выходов и структуру, описываемую системой уравнений, которая в векторной форме имеет вид:

Y = B 0 + BX .

Допустим, модель имеет несколько входов и один выход, содержит число k = n + 1 неизвестных параметров.

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

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

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

Рассмотрим адаптивный шаговый метод. Суть метода: значение параметров модели связываются на двух следующих друг за другом шагах:

где J – алгоритм адаптации.

В качестве такого алгоритма часто используют метод наискорейшего спуска.

Достоинства метода: возможность использования текущей информации.

Недостаток: возникают проблемы сходимости процесса адаптации.

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

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

Объект отражается в виде функции F (X , B ) с неизвестными параметрами B .

Неизвестная функция объекта F 0 (X ) представлена в виде известной функции с неизвестными параметрами Y = F (X , B ). Чтобы определить неизвестные параметры B , приравнивают состояние модели и объекта для каждого из наблюдений. Решение сводится к задаче минимизации суммарной невязки:

2 Методы структурной идентификации объектов управления.

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

Задача идентификации: количественная оценка степени идентичности модели реальному объекту.

В зависимости от априорной (исходной) информации об объекте различают структурную и параметрическую идентификацию.

Предметом структурной идентификации является определение вида функции Y теор связывающей входные переменные Х . Структурная идентификация включает в себя: постановку задачи; выбор структуры модели и её математическое описание; исследование модели.

Задачи вскрытия структуры объекта:1) выделение объекта из среды; 2) ранжирование входов и выходов объекта по степени их влияния на конечный целевой показатель; 3) определение рационального числа входов и выходов объекта, учитываемых в модели; 4) определение характера связи между входом и выходом модели объекта.

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

Идентификация систем ЛА

Формирование моделей на основе результатов наблюдений и исследование их свойств - вот, по существу, основное содержание науки. Модели ("гипотезы", "за­коны природы", "парадигмы" и т.п.) могут быть более или менее формализованны­ми, но все обладают той главной особенностью, что связывают наблюдения в некую общую картину. Решение задачи построения математических моделей динамических систем по данным наблюдений за их поведением составляет предмет теории иденти­фикации, которая тем самым становится элементом общей научной методологии. А поскольку мы "окружены" динамическими системами, методы идентификации систем имеют широкие приложения. Цель данного раздела заключается в том. чтобы дать минимальное представление об имеющихся методах идентификации, их обосновании, свой­ствах и применении.

Динамические системы

Говоря нестрого, система - это объект, в котором происходит взаимодействие между разнотипными переменными и формируются наблюдаемые сигналы.

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

Рис 3.2 Движение судна в горизонтальной Рис. 3.3 Система динамики рулевого

плоскости (δ-команда на руль, управления (δ- входной сигнал, ψ-выходной

ψ- курсовой угол) сигнал, υ-неизмеряемая помеха)

Рис. 3.4. Входо-выходные данные для системы рулевой динамики судна (интервал между замерами -10с.)

Пример Динамика управления судном.

Движение судна происходит под действием тяговой силы винта и зависит от положения рулей, силы и направления ветра и волн. См. рис. 3.2. В качестве подпроблемы можно рассмотреть частную задачу о зависимости курса судна (направление движения носовой части) от положения рулей при постоянном тяговом усилии. Эта система изображена на рис. 3.3. Записи данных наблюдений показаны на рис. 3.4. Длительность интервала наблюдений составила 25 мин, замеры осуществлялись каждые 10 с.

Процедура идентификации системы. Три основных компонента

Конструирование моделей по данным наблюдений включает три основных ком­понента.

1. Данные.

2. Множество моделей-кандидатов.

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

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

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

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

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

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

Контур идентификации системы. Процедура идентификации системы порождает следующую естественную логику действия: (1) собрать данные; (2) выбрать мно­жество

моделей; (3) выбрать наилучшую в этом множестве модель. Однако вполне

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

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

Численный метод не позволяет найти наилучшую по выбранному критерию модели;

Критерий выбран неудачно;

Множество моделей оказалось неполноценным в том смысле, что в этом мно-­
жестве вообще нет "достаточно хорошего" описания системы;

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

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

Параметрическая идентификация объектов.

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

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

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

Описание линейных систем.

Поскольку обработка сигналов в вычислительной машине производится дискретно, то целесообразным является описание линейных систем и сигналов на основе Z – преобразования. При этом непрерывные процессы и отклик системы дискретезируются с тактовым шагом T 0 . (См. рис 3.6).

k = t / T 0

Переход к дискретному времени k=t/T 0 позволяет описывать поведение линейной сиситемы с помощью разностного уравнения.

Используя понятие Z – оператора, где , достаточно просто представляется непрерывное звено.

Общий вид:

или (обратно) в области времени:

Обратно в области времени:

Дифференциальное уравнение системы:

Где τ – чистое запаздывание.

Передаточная функция, следовательно, имеет вид:

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

Рис. 9.1. Общая схема адаптивной системы управления (r, u - входы, - параметры).

1-й блок – это реальный объект, точное описание которого неизвестно. 2-й блок – более или менее точное его описание, 3-й – закон управления.

Модель (2-й блок) – приближённая и объект изменяется во времени. Улучшение модели называется идентификацией (повышение её точности). Идентификация имеет 2 стороны:

структурная;

параметрическая.

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

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

5-й блок (рисунок 9.1) осуществляет адаптацию параметров регулятора.

Адаптация означает подстройку параметров регулятора, его самообучение. Параметры регулятора зависят от параметров объекта и обычно выражаются через них. Изначально параметры объекта не всегда известны с достаточной степенью точности или “плывут” во времени, что и вынуждает прибегать к адаптации. Поэтому адаптации – изменению характеристик регулятора – обычно предшествует процедура уточнения характеристик объекта по результатам измерения входных и выходных величин, которая называется идентификацией. Вопросы идентификации и адаптации получили резкое развитие в последние 25-30 лет, благодаря быстрому росту ресурсов цифровых контроллеров с одновременным снижением их стоимости.

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

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

Рис. 9.2. (А - адаптер, Р – регулятор, О – объект).

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

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

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

В этой системе можно выделить этапы:

1. Идентификация объекта или системы в целом.

2. Расчет коррекции регулятора.

3. Коррекция (настройка) регулятора, изменение его структуры.

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

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

Параметрическая идентификация

Параметрическая идентификация моделей объектов позволяет сразу находить значения коэффициентов модели объекта по измеряемым значениям управляемого y и управляющего u сигналов объекта. При этом предполагается, что структура и порядок модели объекта уже известен. Измеряемые значения y и u представляются в виде временного ряда, поэтому в результате идентификации оцениваются параметры АРСС - модели объекта, или параметры его дискретной передаточной функции. Зная коэффициенты АРСС - модели и ее структуру можно перейти к непрерывным структурированным моделям и моделям в пространстве состояний.

В задачах параметрической идентификации используются модели объекта с шумом измерений, задаваемые передаточными функциями - и структурой. Считая порядки моделей заданными, задачей параметрической идентификации стохастической системы считается определение оценок коэффициентов полиномов модели A,B,C и D по результатам измерений входа u(t) и выхода y(t) . Свойства получаемых оценок (состоятельность, несмещенность и эффективность) зависят от характеристик внешних возмущений и метода идентификации, при этом существенную роль играет вид закона распределения внешних возмущений.

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