Идентификация модифицированного объекта Гаммерштейна
Рассмотрение задачи компенсации влияния входной инерционности при активной идентификации структуры, представляющей собой последовательное соединение линейного динамического звена, нелинейного статического звена и второго линейного динамического звена.
Рубрика | Физика и энергетика |
Вид | статья |
Язык | русский |
Дата добавления | 31.08.2018 |
Размер файла | 409,9 K |
Отправить свою хорошую работу в базу знаний просто. Используйте форму, расположенную ниже
Студенты, аспиранты, молодые ученые, использующие базу знаний в своей учебе и работе, будут вам очень благодарны.
Размещено на http://www.allbest.ru/
ИДЕНТИФИКАЦИЯ МОДИФИЦИРОВАННОГО ОБЪЕКТА ГАММЕРШТЕЙНА
Р.Ш. Галиуллин Ренат Шагитович Галиуллин - аспирант.
Самарский государственный технический университет
443100, г. Самара, ул. Молодогвардейская, 244
E-mail: galren@mail.ru
В статье рассматривается задача компенсации влияния входной инерционности при активной идентификации структуры, представляющей собой последовательное соединение линейного динамического звена, нелинейного статического звена и второго линейного динамического звена. Задача решается путем адаптации известного метода идентификации объекта Гаммерштейна. Приводится несколько способов такой адаптации.
Ключевые слова: идентификация, компенсация, динамическое звено, нелинейность, Гаммерштейн.
компенсация инерционность идентификация структура
Многие технологические объекты управления могут быть описаны как объекты со структурой Гаммерштейна [1], однако управляющее воздействие подается на вход данного объекта через исполнительное устройство, которое может быть описано как линейное динамическое звено, и его передаточная функция известна.
Таким образом, идентифицируемая структура представляет собой последовательное соединение линейного динамического звена, нелинейного статического звена и второго линейного динамического звена.
Рассматривается задача активной идентификации данной структуры при адаптации известного метода идентификации объекта Гаммерштейна [2] путем компенсации влияния входной инерционности [3] и использования тестового сигнала в виде двоичного белого шума.
На рис. 1 приведена структурная схема описанной модели, где - тестовый сигнал в виде двоичного белого шума, - передаточная функция линейного динамического звена объекта Гаммерштейна, - передаточная функция входного динамического звена, например исполнительного устройства, - функция, описывающая входную нелинейную характеристику объекта Гаммерштейна, на выходе снимается сигнал .
Рис. 1. Структурная схема модели
Характеристики модели , неизвестны и должны быть оценены по результатам экспериментальных исследований, причем желательна оптимизация условий данного эксперимента, обеспечивающая адекватность модели.
Пусть входное динамическое звено описывается передаточной функцией вида
(1) |
Линейное динамическое звено объекта Гаммерштейна представлено инерционным звеном первого порядка с передаточной функцией:
(2) |
Нелинейность имеет вид
(3) |
При конкретном виде тестового воздействия на вход линейной подструктуры объекта поступает воздействие вида .
Для оценки возможности идентификации описанного объекта без компенсации входной инерционности проведена идентификация данного объекта с использованием известного алгоритма [2], когда тестовый сигнал имеет характер двоичного белого шума.
Оценка импульсной переходной функции можно определить из соотношения [2]
(4) |
где - нечетная функция; - оценка взаимной корреляционной функции выходного и входного сигналов в матричной форме; - оценка взаимной корреляционной функции выходного и входного сигналов.
Оценка взаимной корреляционной функции входного и выходного сигнала определяется следующим образом:
R, |
(5) |
где - количество отсчетов реализации наблюдаемых сигналов; - количество отсчетов памяти объекта; - интервал квантования по времени.
Оценка нелинейности может быть представлена как
, |
(6) |
где
(7) |
оценка математического ожидания на разных амплитудах;
(8) |
оценка взаимной корреляционной функции выходного и входного сигналов на разных амплитудах.
Результаты идентификации показаны на рис. 2; для оценки точности идентификации введем такие оценки, как максимальное значение ошибки (9) (обозначим как ) и среднеквадратическая ошибка (10).
, |
(9) |
где - импульсная переходная функция линейного динамического звена модели объекта Гаммерштейна.
. |
(10) |
Рис. 2. Идентификация модифицированного объекта типа Гаммерштейна: 1 - импульсная переходная функция линейного динамического звена модели объекта Гаммерштейна; 2 - оценка импульсной переходной функции модели объекта Гаммерштейна; 3 - оценка импульсной переходной модели объекта Гаммерштейна с одним компенсационным звеном; 4 - оценка импульсной переходной модели объекта Гаммерштейна с тремя компенсационными звеньями
Как видно из рис. 2, без использования компенсационного звена оценка импульсной переходной функции совершенно не отвечает виду импульсной переходной функции линейного динамического звена модели объекта. При описании точности идентификации величинами максимальной ошибки и среднеквадратической ошибки становится очевидным, что значения ошибок становятся велики по сравнению с экспериментальными значениями идентифицируемой характеристики: =0,3657, =0,0533.
Для получения адекватной модели объекта данного типа необходимо подать на вход объекта воздействие , вид которого привел бы к компенсации влияния входного динамического звена. При такой компенсации воздействие должно иметь вид двоичного белого шума.
Очевидно, что вводимое компенсирующее звено должно решать обратную задачу и в идеальном случае его передаточная функция должна иметь вид .
Исследуем различные виды компенсирующих цепей.
Первый вид компенсационного звена - использование замкнутой структуры, в которой модель входной инерционности включена в прямую и обратную связи и тем самым образует цепочечную структуру [3] (рис. 3). Варьируя параметр в пределах и контролируя устойчивость системы, можно добиться выполнения условия инвариантности в заданном частотном диапазоне [3]. Результат идентификации приведен на рис. 2.
Рис. 3. Структурная схема компенсирующей цепи с цепочечной структурой
Включение одного компенсационного звена позволило нам получить более точную оценку импульсной переходной функции, чем без компенсации (=0,2247, =0,0328).
На рис. 2 показаны результаты идентификации при цепочечной структуре из 3 звеньев. Как видно, с увеличением количества компенсационных звеньев результат идентификации улучшается: =0,0985 =0,0197.
Дальнейшее же увеличение количества звеньев приводит к незначительному его улучшению, например при 5 звеньях: =0,0924, =0,0180.
Для увеличения эффективности процедуры идентификации предложен следующий вид структуры компенсационного звена (рис. 4).
Рис. 4. Структурная схема компенсационного звена с объектом в обратной связи
В данной структуре - звено с пропорционально-интегральным законом преобразования:
(11) |
Структура модели объекта, на входе которого расположено компенсационное звено данного вида, приведена на рис. 5, а результат идентификации - на рис. 6.
Рис. 5. Структурная схема модели объекта с компенсационным звеном: - компенсационное звено
Рис. 6. Оценка импульсной характеристики:
1 - импульсная переходная функция линейного динамического звена модели объекта Гаммерштейна; 2 - оценка импульсной переходной модели объекта Гаммерштейна
Использование интегральной составляющей в компенсационном звене позволяет получить адекватную оценку импульсной переходной функции при =0,1356, =0,0231.
Сравнение результатов идентификации, показанных на рис. 2 и 6, показывает большую адекватность модели, полученной при идентификации с компенсационным звеном со структурой, приведенной на рис. 4. Хотя значения ошибок несколько выше, но, как видно из рис. 6, характер оценки импульсной переходной функции в большей мере соответствует характеру данной функции самого объекта. В том числе запаздывание по времени, наблюдаемое на рис. 2, в результате идентификации при компенсирующей цепи с цепочечной структурой полностью исчезло.
Таким образом, применение компенсационных звеньев с разной структурой позволяет проводить идентификацию объекта с дополнительной входной инерционностью с помощью наиболее эффективного метода, использующего тестовый сигнал в виде двоичного белого шума.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Буштрук А.Д. Структурная идентификация нелинейных динамических объектов // АиТ. - 1989. - №10. - С. 84-96.
2. Чостковский Б.К., Юдашкин А.А. Активная идентификация нелинейных динамических объектов типа Гаммерштейна // АиТ. - №1. - 1992. - С. 96-103.
3. Рапопорт Э.Я., Тян В.К. Достижение заданной инвариантности в стохастических системах комбинированного управления // Куйбыш. политех. ин-т. Деп. в ВИНИТИ 1989.
Размещено на Allbest.ru
...Подобные документы
Общее описание интегро-дифферецирующее звена. Дифференцирующие и интегрирующие свойства звена. Дифференциальное уравнение для цепи. Уравнение и график переходной и импульсной функции. Интегро-дифферецирующее звено с преобладанием дифференцирующих свойств.
реферат [438,3 K], добавлен 17.12.2011Электропривод звена промышленного робота. Типовой технологический процесс и выбор манипулятора. Выбор и проверка электродвигателя. Расчет динамических параметров привода, определение его основных характеристик. Расчет расхода энергии и КПД за цикл.
курсовая работа [1,4 M], добавлен 30.04.2012Определение инерционных свойств средств измерений. Построение временных (переходных) характеристик СИ. Конструкция и динамические свойства термометра сопротивлений. Экспериментальное определение динамических характеристик звена первого и второго порядка.
контрольная работа [106,4 K], добавлен 01.02.2013Определение свойств объекта, подлежащего исследованию. Изменение сопротивления медного проводника. Процессы распространения тепловой энергии. Идентификация типа дифференциального уравнения. Входной и выходной параметры. Размерность входного возмущения.
курсовая работа [190,5 K], добавлен 13.03.2014Исследование динамического поведения механической системы с использованием теорем и уравнений теоретической механики. Дифференциальное уравнение движения механической системы. Законы движения первого груза, скорость и ускорение в зависимости от времени.
реферат [107,8 K], добавлен 27.07.2010Исследование формы и расчётов характеристики динамического торможения. Расчет эквивалентного момента торможения, критического скольжения и момента, электромеханической характеристики ЭД. Схема динамического торможения АД с короткозамкнутым ротором.
лабораторная работа [15,6 K], добавлен 12.01.2010Задачи на применение первого закона Кирхгофа. Параллельное соединение элементов. Второй закон Кирхгофа, его применение. Последовательное соединение конденсаторов, их эквивалентная емкость. Обратная емкость конденсаторов, соединенных последовательно.
реферат [85,5 K], добавлен 15.01.2012Построение и расчет зубчатого зацепления и кулачкового механизма. Проектирование и кинематическое исследование зубчатой передачи и планетарного редуктора. Определение уравновешенной силы методом Жуковского. Построение диаграмм движения выходного звена.
курсовая работа [400,8 K], добавлен 23.10.2014Исследование процессов, происходящих в простейших электрических цепях переменного тока, содержащих последовательное соединение активных и индуктивных сопротивлений. Измерение общей силы тока, активной и реактивной мощности; векторная диаграмма напряжений.
лабораторная работа [79,2 K], добавлен 11.05.2013Построение планов положений и кинематических диаграмм. Определение скорости и ускорения ведомого звена в исследуемом положении двигателя при помощи диаграмм. Определение сил приложенных к звеньям механизма. Определение потребной мощности двигателя.
контрольная работа [240,2 K], добавлен 10.08.2012Определение положения мгновенного центра скоростей для каждого звена механизма и угловые скорости всех звеньев и колес. Плоскопараллельное движение стержня. Расчет скорости обозначенных буквами точек кривошипа, приводящего в движение последующие звенья.
контрольная работа [66,5 K], добавлен 21.05.2015Специфика измерения силы тока амперметром и напряжения вольтметром. Методика расчета падения напряжения на приемниках по закону Ома и по второму закону Кирхгофа на различных участках цепи. Сравнительный анализ расчетных и измерительных параметров цепи.
лабораторная работа [22,9 K], добавлен 12.01.2010Построение рациональных эксплуатационных режимов асинхронного двигателя, выбор системы управления. Исследование двухмассового динамического стенда на базе математической модели. Техническая разработка лабораторного стенда на базе асинхронного двигателя.
магистерская работа [2,0 M], добавлен 20.10.2015Определение ионосферы и линейного слоя, расчёт диалектической проницаемости ионосферы без учёта магнитного поля. Распределение магнитного поля в точке попадания на Землю отражённого луча. Закон изменения электронной концентрации для линейного слоя.
курсовая работа [321,8 K], добавлен 14.07.2012Составление баланса активной и реактивной мощностей генератора и нагрузки. Проверка его выполнимости для симметричного и несимметричного режимов. Расчет фазного и линейного напряжения и мощности генератора. Построение топографической диаграммы токов.
контрольная работа [374,5 K], добавлен 16.05.2015Динамический, структурный, кинематический и силовой анализ механизма, построение плана скоростей и ускорений. Выбор расчетной схемы и проектный расчет механизма на прочность. Построение эпюр и подбор сечений звена механизма для разных видов сечений.
курсовая работа [118,9 K], добавлен 18.09.2010Порядок построения кинематической схемы рычажного механизма по структурной схеме, коэффициенту изменения скорости выходного звена и величине его полного перемещения. Число подвижных звеньев механизма, построение диаграммы перемещения и плана скоростей.
курсовая работа [63,4 K], добавлен 11.11.2010Расчет средней скорости и среднего ускорения в интервале заданного времени. Поиск силы, действующей на тело, движущееся с ускорением. Потенциальная энергия груза, расчет его ускорения. Поиск линейного ускорения с использованием второго закона Ньютона.
контрольная работа [207,3 K], добавлен 23.09.2013Уравнение теплопроводности: физический смысл и выводы на примере линейного случая. Постановка краевой задачи остывания нагретых тел, коэффициент теплопроводности. Схема метода разделения переменных Фурье применительно к уравнению теплопроводности.
курсовая работа [245,8 K], добавлен 25.11.2011Решение задачи идентификации коэффициента температуропроводности непрерывнолитого стального цилиндрического слитка. Математическая модель теплового процесса. Методы поиска градиента функции с помощью сопряженной задачи и численного дифференцирования.
практическая работа [96,8 K], добавлен 02.07.2012