Корни характеристического уравнения апериодического звена

3.4 Апериодическое звено второго порядка

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

У нас есть модель механического демпфера. Это поршень на пружине, он движется внутри цилиндра, может перемещается вверх-вниз. Его положение – это интересующая нас функция Y(t), сверху на него воздействует возмущающая сила (U(t)), на стенках поршня действует сила вязкого трения. (См. рис. 3.4.1)

Рисунок 3.4.1. Расчетная схема амортизатора.

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

— масса поршня;

— положение поршня (выходная переменная);

— приложенная сила (входное воздействие);

— сила тяжести;

– сила сопротивления пружины;

– сила вязкого трения (пропорциональная скорости движения поршня).

Считаем, что в нулевой момент времени поршень находится в равновесии. Тогда начальное положение поршня — y₀ в равновесии, где скорость и ускорения равны 0, можно посчитать из уравнения 2.

Перепишем уравнение равновесия в отклонениях от нулевого состояния:

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

Приведем данное уравнение к классическому виду:

Уравнение динамики апериодического звена 2−го порядка имеет следующий вид:

Если D Рисунок 3.4.2 Апериодическое звено 2-го порядка (два варианта)

Амплитудно-фазовая частотная характеристика (АФЧХ):

Домножив числитель и знаменатель формулы (3.4.5) на комплексно-сопряженные скобки и , получаем:

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

Анализируя поведение и при и при , получаем:

Модуль АФЧХ (амплитуда), то есть mod(W(i·ω)) = |W(i·ω)| из формулы 3.4.5:

Подставляя в формулы (3.4.6) или в формулу (3.4.5) различные значения ω можно построить векторы, соответствующие различным значениям ω:

Рисунок 3.4.3 Годограф АФЧХ апериодического звена 2-го порядка

Из формул 3.4.6 очевидно, что на рисунке годографа 3.4.3 :

omega_5>omega_4>omega_3>omega_2>omega_1>0\ 2) 0 >varphi_1>varphi_2>varphi_3>varphi_4>varphi_5>varphi_6″ alt=»1) omega_6>omega_5>omega_4>omega_3>omega_2>omega_1>0\ 2) 0 >varphi_1>varphi_2>varphi_3>varphi_4>varphi_5>varphi_6″ src=»https://habrastorage.org/getpro/habr/upload_files/663/02f/0f1/66302f0f16e854680e603612571bc0af.svg» width=»733″ height=»44″/>

Используя формулу 3.4.6 можно показать что при

Из рисунка видно, что .

Формула фазового сдвига:

omega_3 Rightarrow j=1.» alt=» omegaleq omega_3 Rightarrow j = 0;\ omega>omega_3 Rightarrow j=1.» src=»https://habrastorage.org/getpro/habr/upload_files/ee5/e65/0f3/ee5e650f39001b205a1e75c0ca3f70a8.svg» width=»733″ height=»40″/>

Для фазового сдвига удобно представить апериодическое звено в виде последовательного соединения двух звеньев (см. рис. 3.4.2). Известно, что при последовательном соединении звеньев общий сдвиг фазы равен сумме фазовых сдвигов:

Логарифмическая амплитудная характеристика (ЛАХ)

Графики А(ω), φ (ω), Lm(ω) имеют вид:

Рисунок 3.4.4 АЧХ и ФЧХ апериодического звена 2-го порядка Рисунок 3.4.5 ЛАХ и ЛФЧХ апериодического звена 2-го порядка

В инженерных расчетах часто график Lm(ω) представляют виде отрезков ломаных, тогда:

при — звено близко к идеальному усилительному звену

при — звено близко к идеальному интегрирующему звену

при 1/T_3″ alt=»omega>1/T_3″ src=»https://habrastorage.org/getpro/habr/upload_files/230/304/ad3/230304ad37f36d023cf3233756aa804a.svg»/>- звено близко к дважды интегрирующему звену

В граничном случае или отмеченные на графике Lm(ω) (см. рис. 3.4.5 выше) точки «излома» совпадают:

Рисунок 3.4.6 ЛАХ и ЛФЧХ апериодического звена 2-го порядка в граничном случае

Если звено “переходит” в разряд колебательных звеньев. Поэтому постоянная Т₁ в уравнении динамики (3.4.1) играет роль демпфирующего фактора, увеличение Т₁ (в колебательном звене) приводит к уменьшению или к полному исчезновению колебаний.

Найдем переходную функцию звена — реакцию на воздействие единичное воздействие 1(t).

Для нахождения функции по формуле Хэвисайда (см. раздел 2.8 Некоторые способы нахождения оригинала по известному изображению), запишем корни полюса изображения, т.е. те значения «s» при которых обращается в ноль:

Тогда по формуле Хэвисайда:

Вычисляя пределы получим формулу для переходной функции звена:

Весовая функция получается дифференцированием :

Рисунок 3.4.7 Переходная функция апериодического звена 2-го порядка

Примерами апериодического звена 2-го порядка являются:

1) двигатель постоянного тока при учете инерционности самого якоря (механической) и цепи якоря (электрической);

2) электрический усилитель руля автомобиля с учетом инерционности (механической и электрической) ротора;

3) двойные R − C или R – L цепочки

Рисунок 3.4.9 Пример апериодического звена 2-го порядка

Если звено представлено в переменных состояния в матричной форме таким образом:

то звено будет апериодическим 2-го порядка, если:

Видео:Теория автоматического управления. Лекция 7. Типовые звенья САУСкачать

Пример

В качестве примера возьмём модель демпфера, которую мы уже использовали в лекциях. (см. Рисунок 3.4.10) Структурная схема модели описывает уравнения динамики, описанные в начале статьи. Свойства системы заданы в списке общих сигналов проекта (см. рис. 3.4.11). Для получения из демпфера апериодического звена 2-го порядка необходимо увеличить силу трения таким образом, чтобы (как показано выше) коэффициент T₁ был больше, чем 2 х T₂. В этом случае D>0 и из колебательного звена мы получим апериодическое 2-го порядка.

Рисунок 3.4.10 Структурная схема модели демпфера. Рисунок 3.4.11 Параметры модели

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

Рисунок 3.4.12. Параметры для модели демпфера в виде звена

Выполним моделирование переходного процесса при ступенчатом изменении приложенной силы и сравним переходные процессы в двух вариантах модели демпфера. График переходного процесса (см. рис. 3.4.13) показывает, что переходные процессы в двух моделях полностью идентичны:

Рисунок 3.4.13 Переходные процессы в двух моделях.

График частотных характеристик звена (ЛАХ и ФЧХ) представлен на рисунке 3.4.14 На графике видно две точки излома характеристики ЛАХ в которых наклон последовательно меняется с 0, до 20дБ/дек и с 20дБ/дек до 40 дБ/дек.

Рисунок 3.4.14 Частотные характеристика ЛАХ и ФЧХ

Для демонстрации влияния изменения Т₁ на свойства звена выполним моделирование, в котором структурная схема является эталонной, а в модели звена будем уменьшать коэффициент силы трения (коэффициент T₁).

Источником воздействия будет меандр, с периодом 3 секунды.

Для изменения свойств звена создадим блок на языке программирования. Данный блок, в процессе моделирования, постепенно уменьшает коэффициент Т₁ для модели в виде звена. Этот же блок готовит данные для отображения на 3D графике переходного процесса.

Общая схема модели приведена на рисунке 3.4.15.

Рисунок 3.4.15 Схема демпфера с изменения свойств блока

Меандр задает изменение приложенной силы 0 – 30 Н (входного воздействия) с полупериодом 1.5 сек. График изменения положения приведен на рисунке 3.4.16 Видно, что на первом изменении графики совпадают, но потом по мере накопления отличий в параметрах динамика изменения положения начинает меняться.

Рисунок 3.4.16 Графики положения демпферов.

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

Рисунок 3.4.17 Начальная часть графика

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

Рисунок 3.4.18 Конечная часть моделировани

ЗD поверхность отображает переходный процесс при ступенчатом увеличении воздействия в блоке меандр. По оси Z отражается положение демпфера, по оси Y – время после увеличения входного воздействия в блоки меандр, по оси X – изменений T₁ (уменьшение силы трения).

Рисунок 3.4.19 Поверхность переходного процесса при снижении трения

В заключение, сравним переходные процессы для разных параметров T₁ (разных коэффициентов трения). Поскольку все основные блоки в SimInTech являются векторными, создадим модели 7-ми демпферов из одного звена. Для этого в главном окне программы подготовим 7 векторов значений с разными коэффициентами трения. Скрипт приведен на рисунке 3.4.20.

Рисунок 3.4.20 Скрипт модели для задания параметров 7 демпферов

Четвертый вектор содержит переходное значение T₁. Как было показано выше, переходное значение T₁, при котором апериодическое звено второго порядка превращается в колебательное рассчитывается по формуле T₁ = 2хT_2.

В модели, в свойствах блока указываем эти векторы в столбце «формулы», и теперь блок может рассчитывать одновременно 7 демпферов одним блоком. (см. рис. 3.4.21)

Рисунок 3.4.21 Настройка параметров блока для векторного расчета

Общая схема модели в этом случае будет выглядеть как показано на рисунке 3.4.22 Ступенчатое изменение силы передается в блок «Размножитель», где преобразуется в вектор из 7 воздействий. Данный вектор передается в блок, где и происходит расчёт семи вариантов демпфера.

Рисунок 3.4.22 Схема модели 7-и демпферов

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

Рисунок 3.4.23 Перемещение 7 демпферов при ступенчатом воздействии

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

Рисунок 3.4.25 Частотные характеристики 7-и демпферов

Модели с примерами для самостоятельного изучения можно взять по ссылке.

Видео:Характеристическое уравнение в ДУСкачать

Апериодическое звено первого порядка

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

Инерционное звено описывается дифференциальным уравнением первого порядка

,

где T – постоянная времени звена (T>0);

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

К апериодическим звеньям можно отнести: R-L и R-C цепи, генераторы постоянного тока, фильтры, термисторы, механические устройства, имеющие массу и силу трения (без пружин) и другие подобные устройства, в которых возможно накопление какого-либо вида энергии и её рассеивания.

Операторное уравнение апериодического звена

Передаточная функция звена

На структурных схемах графически инерционное звено изображается следующим образом:

Временная характеристика, представляющая реакцию звена на ступенчатое воздействие x_вх(t)=1(t), определяется зависимостью .

Выходная величина в переходном режиме определяется

где вынужденная составляющая выходной величины

Cвободная составляющая выходной величины x_вых(t) определяется из следующего выражения

где P_k – корни характеристического уравнения звена

т.е.

Отсюда

Начальное значение для переходной функции найдется

при t=0,

т.е.

или

Окончательно получаем следующее выражение для переходной функции:

(3.11)

или (3.12)

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

при , при .

представлена на рис.3.5,б.

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

Амплитудно-фазовая характеристика апериодического звена

где модуль вектора W (jω);

аргумент вектора W (jω).

АФХ представляет собой окружность радиусом с центром в точке 0, лежащей на оси абсцисс на расстоянии от начала координат.

Уравнения вещественной и мнимой характеристик

Логарифмическая амплитудная частотная характеристика (ЛАЧХ) может быть получена путём логарифмирования выражения для A(ω).

В этом выражении слагаемое представляет постоянную величину, не зависящую от частоты. Рассмотрим вторую составляющую ЛАЧХ

.

Полагая, что

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

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

Дата добавления: 2016-03-20 ; просмотров: 10714 ; ЗАКАЗАТЬ НАПИСАНИЕ РАБОТЫ

Видео:10) ТАУ для чайников Части 4.1. и 4.2. Типовые динамические звенья. Усилитель. Апериодическое звено.Скачать

Звено второго порядка (колебательное звено)

Дата добавления: 2013-12-23 ; просмотров: 28109 ; Нарушение авторских прав

Такие звенья описываются дифференциальным уравнением вида:

Если на вход звена подать единичную функцию Хэвисайда от времени 1[t], при нулевых начальных условиях системы, то реакция на выходе будет называться переходной функцией (или переходной характеристикой), которую часто обозначают как h(t). Сигнал 1[t] — это, в некотором смысле, эталонный испытательный сигнал. Существуют и другие эталонные испытательные сигналы. Например, бесконечный импульс нулевой длины (дельта-функция Дирака), гармонический сигнал, периодические прямоугольные импульсы.

Преобразуем по Лапласу это уравнение:

Определим передаточную функцию звена:

Если записать уравнение без входного воздействия (нулевые входные воздействия U = 0) и сократить Y, то есть: T 2 p 2 + 2ξTp + 1 = 0, то такое уравнение будет называться характеристическим, поскольку характеризует исключительно внутренние свойства звена. Обратите внимание, что в записи звена содержатся три параметра:

T — постоянная времени (в секундах); ξ — коэффициент затухания (безразмерная величина); k — передаточный коэффициент. В зависимости от величины ξ звенья второго порядка классифицируются по видам: ξ = 0 — консервативное звено второго порядка; 0 2 p 2 + 2ξTp + 1 = 0. И оно имеет действительные отрицательные корни: Данное звено можно представить в виде последовательно соединенных звеньев с различными постоянными времени: Tогда при T1 > T2 переходная характеристика звена имеет вид: То есть в решении присутствуют затухающие экспоненты. Типичное поведение звена с такими параметрами показано на рис. 4.6. Рис. 4.6. Реакция апериодического звена на единичный входной сигнал В частном случае, когда ξ = 1, оба корня будут одинаковыми, отрицательными: 15 Колебательное звено 2-го порядка (0 2 p 2 + 2ξTp + 1 = 0. Корни разные, комплексно-сопряженные, с отрицательной вещественной частью: Так как корни мнимые, то в поведении звена присутствует колебательная составляющая. Именно за эту особенность поведения звено получило название колебательного (см. рис. 4.7 и рис. 4.8). Рис. 4.7. Реакция колебательного звена на входной единичный сигнал (ξ = 0.5) Рис. 4.8. Реакция колебательного звена на входной единичный сигнал (ξ = 0.2) Из графиков видно, что с ростом ξ колебательность звена уменьшается, исчезая при ξ ≥ 1 Переходная функция звена имеет вид: При малых ξ значение A приближается к 1, а значение φ — к 90°. По физическому смыслу ω0 представляет собой собственную частоту колебаний. 16 Консервативное звено 2-го порядка (ξ = 0) Характеристическое уравнение звена следующее: Корни одинаковые, комплексно-сопряженные, с нулевой вещественной частью: Так как корни чисто мнимые, то поведением звена являются незатухающие колебания (ξ = 0), см. рис. 4.9. Рис. 4.9. Реакция колебательного звена на входной единичный сигнал (ξ = 0) Переходная функция звена имеет вид: h(t) = k · (1 – cos(t/T)). Из графика экспериментальным путем можно определить единственный параметр T = T0/(2 · π). 5. Лекция 05. Динамические регрессионные модели, заданные в виде передаточной функции Построим регрессионную модель динамической системы на примере. Зададим модель в виде передаточной функции (см. рис. 5.1). Рис. 5.1. Модель черного ящика в виде передаточной функции Примерный вид динамических сигналов на входе и на выходе показан на рис. 5.2. Ограничимся временем рассмотрения сигналов, равным T. Рис. 5.2. Возможный вид зависимости входного сигнала X от времени t и зависимости выходного сигнала Y от t для случая с непрерывным временем После дискретизации, связанной с обработкой информации на цифровых машинах, эти сигналы будут выглядеть так, как показано на рис. 5.3. Обратите внимание на то, что отдельные отсчеты отстоят друг от друга на расстоянии Δt. Важно, что отсчеты стоят достаточно часто. Всего этих отсчетов — n, то есть T = n · Δt, Рис. 5.3. Возможный вид зависимости входного сигнала X от времени t и зависимости выходного сигнала Y от t для случая с дискретным временем Допустим, что зависимости, представленные на рис. 5.2 и рис. 5.3, описываются передаточной функцией следующего вида (заметим, что, как и в лекции 02, вид зависимости выдвигается гипотетически и гипотеза должна быть, в конце концов, подтверждена или опровергнута): Заменяя значок «p» на «d/dt» (обращаем ваше внимание, что такую замену можно производить только для случая нулевых начальных условий: X(0) = 0 и Y(0) = 0) и учитывая, что передаточная функция — это, по определению, отношение выхода к входу, то есть W = Y/X, получаем дифференциальное уравнение 2-го порядка: Дважды проинтегрируем это выражение и получим для некоторого произвольного момента времени t: Коэффициенты A1, A2, A3, A4 требуется определить. Для этого выразим уравнение в разностном виде через суммы: где n — число экспериментальных точек. Заметим, что мы заменили интегралы и непрерывное течение времени — на суммы и дискретное представление времени. Чтобы вычислить суммы и двойные суммы экспериментально заданных зависимостей x и y, воспользуемся для удобства табл. 5.1 (чтобы излишне не загромождать таблицу, мы не стали дописывать в суммах Δti и Δτj). Таблица 5.1. Таблица исходных данных и вспомогательных расчетов i Xi Yi X1 Y1 X1 Y1 X1 Y1 X2 Y2 X1 + X2 Y1 + Y2 2X1 + X2 2Y1 + Y2 X3 Y3 X1 + X2 + X3 Y1 + Y2 + Y3 3X1 + 2X2 + X3 3Y1 + 2Y2 + Y3 … … … … … … … m Xm Ym X1 + X2 + … + Xm Y1 + Y2 + … + Ym mX1 + … + Xm mY1 + … + Ym … … … … … … … n Xn Yn X1 + X2 + … + Xn Y1 + Y2 + … + Yn nX1 + … + Xn nY1 + … + Yn Ошибку в некоторой m-ой точке можно записать так: Как и ранее, ошибка показывает, насколько отходит теоретическое значение Ym от экспериментального значения. Суммарная ошибка (вносимая всеми точками), которую надо минимизировать, будет: Величина ошибки зависит от значений параметров A1, A2, A3, A4. Поэтому F является функцией от четырех переменных: F(A1, A2, A3, A4). Чтобы найти минимум функции F, доставляемый за счет параметров A1, A2, A3, A4, надо взять частные производные F по каждому из параметров и приравнять каждую производную к нулю. В результате получаем систему линейных уравнений: Получено четыре уравнения с четырьмя неизвестными A1, A2, A3, A4. Из решения системы уравнений вычисляем неизвестные коэффициенты и дополняем ими модель, где коэффициенты уже определены как числа: Задача определения коэффициентов модели решена. Разумеется, как и ранее, необходимо сравнить получаемое из этой модели решение Y теоретическое с Y, заданным экспериментально, и вычислить ошибку F. Далее проверить ее значение по критерию — допустимо ли значение вычисленной ошибки, или гипотезу о виде модели требуется сменить на более точную. Считая что коэффициенты модели теперь нам известны, построим для заданного примера реализацию, имитирующую поведение системы, описанной передаточной функцией. Для этого воспользуемся уже однажды полученной формулой: Реализация модели представлена на рис. 5.4. Рис. 5.4. Техническая реализация передаточного звена после определения коэффициентов регрессионной модели При переходе от интеграла к численному суммированию мы воспользовались методом прямоугольников. Разбив площадь под кривой y на ряд прямоугольников одинаковой ширины Δt (см. рис. 5.5), получаем, что площадь i-го прямоугольника равна yi · Δt, а S — сумма площадей всех n прямоугольников — будет приблизительно равна площади под кривой (интегралу от функции y). Очевидно, что приближение тем точнее, чем меньше значение Δt. Рис. 5.5. Применение метода прямоугольников для численного вычисления интегралов 6. Лекция 06. Модель в виде фильтра Каллмана Каллманом была доказана теорема о том, что любой динамический сигнал может быть представлен в виде: Рис. 6.1. Графическое представление фильтра Каллмана на схемах Идея фильтра Каллмана заключается в том, что выход системы в i-ый момент времени определяется входным сигналом, его предысторией и предысторией самого состояния системы. Чем больше имеется членов ряда, то есть чем больше переменных Y учитывается в записи модели, тем глубже память системы. Заметим, что наличие члена Yi – 1 в модели динамической системы соответствует наличию первой производной, Yi – 2 — второй производной и т. д. Допустим, известны следующие экспериментальные данные: состояния сигналов Xi и Yi в n временных точках (табл. 6.1). Таблица 6.1. Таблица экспериментальных данных i Xi Yi X1 Y1 X2 Y2 … … … n – 1 Xn – 1 Yn – 1 n Xn Yn Поскольку для каждой экспериментальной точки Xi надо указать ее соседей, задаваемых рядом, то удобно отсчеты представить в расширенной таблице, используемой для расчета (см. табл. 6.2). Таблица 6.2. Таблица экспериментальных данных и промежуточных расчетов i Xi Xi – 1 … Yi Yi – 1 Yi – 2 m Xm Xm – 1 … Ym Ym – 1 Ym – 2 m + 1 Xm + 1 Xm … Ym + 1 Ym Ym – 1 m + 2 Xm + 2 Xm + 1 … Ym + 2 Ym + 1 Ym … … … … … … … Находим ошибку между значением экспериментально снятой точки и теоретическим ее значением (гипотезой): Суммарная ошибка F (сумма берется по всем экспериментальным точкам) должна быть минимизирована относительно определяемых переменных A1, A2, …, B1, B2, …, C: После взятия частных производных от F по A1, A2, …, B1, B2, …, C, приравнивания их к нулю и составления системы уравнений получается линейная множественная регрессионная модель, из которой определяются неизвестные коэффициенты A1, A2, …, B1, B2, …, C модели. Поскольку коэффициенты модели определены, построим реализацию (см. рис. 6.2), имитирующую поведение системы, описанной фильтром Каллмана. Рис. 6.2. Вариант технической реализации фильтра Каллмана «Блок задержки» в представленной реализации необходим для того, чтобы сдвинуть сигнал на такт и получить соседний отсчет для следующей переменной ряда модели. В зависимости от среды реализации блок задержки можно организовать разными способами. Например, в случае реализации блока задержки в среде моделирования Stratum-2000, первый способ может быть основан на перезаписи информации из одной переменной (ячейки) в другую, на что требуется один такт. Таким образом, можно организовать задержку сигнала на любое число тактов. Например, задержка сигнала X относительно Y будет составлять 3 такта, если выполнить следующую последовательность операций: A1 := X; A2 := A1; Y := A2. Во втором способе задержка организуется при помощи массива: на каждом такте нужно, чтобы цифры были перемещены в соседние ячейки. На рис. 6.3 приведена схема настройки (автоматического нахождения коэффициентов). Рис. 6.3. Схема автоматической настройки коэффициентов модели «на ходу» На рис. 6.4 приведена схема проверки фильтра Каллмана. Рис. 6.4. Схема проверки работы модели фильтра Каллмана 7. Лекция 07. Модель динамической системы в виде Фурье представления (модель сигнала) Этот способ моделирования динамических систем основывается на том, что в любом сигнале присутствуют гармонические составляющие. В зависимости от частоты, составляющие называются гармониками (первая, вторая и так далее). Сумма гармоник с соответствующими весами составляет модель сигнала. Пусть, например, в некотором сигнале присутствует сумма трех гармоник: 3 · cos(t) + 2 · cos(3t) + 0.5 · cos(5t). Это значит, что в сигнале присутствует первая гармоника с амплитудой 3, третья гармоника с амплитудой 2, пятая гармоника с амплитудой 0.5. Сам суммарный сигнал выглядит так, как показано на рис. 7.1. Рис. 7.1. Пример гармонического сигнала Спектр этого сигнала показан на рис. 7.2. Ясно, что в нашем примере больший вес (амплитуду) в сигнале имеет (более других представлена) первая гармоника, наименьший вес имеет пятая гармоника. Рис. 7.2. Пример спектра гармонического сигнала Любой сигнал, сколь сложен бы он ни был, может быть представлен суммой гармоник. Более простой сигнал представляется меньшим числом гармоник, более сложный — большим. Быстро меняющийся сигнал, содержащий резкие пики, имеет в своем составе гармоники высоких порядков. Чем больше гармоник представлено в модели сигнала, тем точнее, в общем случае, модель отражает реальный сигнал. Пусть задан некий сигнал X(t) (рис. 7.3). Рис. 7.3. Временной сигнал на входе преобразования Фурье (возможный вид) Определимся со временем рассмотрения сигнала: если сигнал периодический, то время рассмотрения равно периоду p сигнала; если сигнал непериодический, то периодом сигнала считается все время его рассмотрения. … … … … Ai и Bi — это веса соответствующих гармоник, присутствующих в сигнале; i — номер гармоники. Формулы их расчета называются прямым преобразованием Фурье. Значение 2π · i/p = ωi — это частота i-ой гармоники. Отметим также, что частота i-ой гармоники связана с частотой первой гармоники простым соотношением: ωi = i · ω1. Отметим важную особенность данного способа представления: вместо всего сигнала во всех его подробностях достаточно хранить вектор чисел, представляющих весовые коэффициенты составляющих его гармоник: (A0, A1, A2, …, B1, B2, …). То есть эти числа полностью характеризуют исходный сигнал, так как по ним сигнал можно полностью восстановить формулой обратного преобразования Фурье: Именно эти числа используются также при обработке сигнала в модели динамической системы. Изображение этих чисел на графике в зависимости от номера гармоники (частоты) называется спектром сигнала (рис. 7.4). Спектр показывает, насколько присутствует в сигнале соответствующая составляющая. Спектр — это частотная характеристика сигнала. Рис. 7.4. Сигнал, представленный в частотной области на выходе преобразования Фурье, спектр сигнала (возможный вид) Здесь сигнал представлен в частотной области. Всегда по формулам прямого преобразования Фурье можно перейти из временной области в частотную, а по формулам обратного преобразования Фурье перейти из частотной области во временную. В какой области (частотной или временной) работать с сигналом в отдельный момент, решают из соображений удобства, наглядности и экономии вычислений. Заметим, что емкие с точки зрения вычислений операции интегрирования и дифференцирования сигнала во временной области заменяются на операции алгебраического сложения и умножения в частотной области, что с вычислительной точки зрения реализуется намного точнее и быстрее. Система чисел Ai и Bi является полной характеристикой сигнала. Такой же полной характеристикой сигнала является система чисел S и φ, которые также образуют спектр (рис. 7.5). S — это амплитудно-частотная характеристика (АЧХ), φ — фазо-частотная характеристика (ФЧХ). Рис. 7.5. Сигнал, представленный в частотной области, амплитудно-частотная и фазо-частотная характеристика сигнала (возможный вид) Системы «A и B» и «S и φ» являются полностью равнозначными. Переход из системы «A и B» в систему «S и φ» производится по следующим формулам: Si = sqrt(Ai 2 + Bi 2 ) — абсолютная амплитуда сигнала; φi = arctg(Bi/Ai) — фаза сигнала, при сложении гармоник нужно учитывать сдвиг фаз (сдвиг фаз проиллюстрирован на рис. 7.8). В случае с системой «S и φ» обратное преобразование Фурье имеет вид: Рис. 7.6 и рис. 7.7 разъясняют смысл коэффициентов A и B разных гармоник. Эти коэффициенты — амплитуды синусов и косинусов соответствующих частот (гармоник). Во временной области графически они соответствуют размаху гармонических колебаний (рис. 7.6 и рис. 7.7); в частотной — высоте спектральной полоски на соответствующей частоте (рис. 7.4). Рис. 7.6. Геометрическая иллюстрация параметров А и ω для косинусной составляющей гармонического сигнала Рис. 7.7. Геометрическая иллюстрация параметров В и ω для синусной составляющей гармонического сигнала Смысл чисел Si и φi разъяснен на рис. 7.8. Рис. 7.8. Геометрическая иллюстрация параметров S и φ для составляющей гармонического сигнала 8. Лекция 08. Модель динамической системы в виде Фурье представления (модель объекта) Пусть имеется входной динамический сигнал X(t) и объект F, преобразующий этот сигнал в выходной Y(t) (см. рис. 8.1). Если объект описывается дифференциальными уравнениями, то таким преобразованием является интегрирование входного сигнала и вычисление Y(t). Интегрирование, как было ранее показано, — операция, требующая значительных вычислительных ресурсов и имеющая значительную погрешность при реализации на цифровых машинах. Рис. 8.1. Схема моделирования динамического объекта при переходе из временной области представления в частотную Если перейти от описания входного сигнала во временной области к описанию в частотной области (см. рис. 8.1), а от дифференциальных уравнений перейти к частотной характеристике объекта, — то есть, фактически, заменить сигнал на частотную модель сигнала, а объект на частотную модель объекта, — то с вычислительной точки зрения процесс преобразования сигнала упростится. Конечно, полученный результат тоже будет частотной моделью выходного сигнала, которую для получения окончательного ответа придется сконвертировать во временную область Y(t). Процесс такой конвертации из частотной области во временную и обратно называется преобразованием Фурье (есть и другие преобразования). Для тех объектов, для которых известна их модель в частотной области, такой подход достаточно просто реализуется на компьютере и позволяет достичь любой наперед заданной точности. Модель объекта в частотном виде называется передаточной функцией или АЧХ (амплитудно-частотной характеристикой). Объекты, для которых известны АЧХ, обычно называют типовыми звеньями (усилительное звено, апериодическое, колебательное и т. д.). Пусть, для примера, характеристика объекта в частотной области следующая (см. рис. 8.2). Рис. 8.2. АЧХ (возможный вид) Амплитудно-частотная характеристика (АЧХ) показывает, насколько пропускается объектом на выход соответствующая гармоника. Значение ki характеризует коэффициент усилeния гармонического сигнала на определенной частоте ωi. Моделирование прохождения сигнала через объект в этом виде заключается в умножении коэффициента Ai гармоники с частотой ωi входного сигнала X(t) на коэффициент усиления ki при той же гармонике с частотой ωi в АЧХ: Ai * = Ai(ωi) · ki(ωi). (Для коэффициента B преобразование аналогично.) В результате получается коэффициент Ai * выходной гармоники данной частоты ωi. Процедура выполняется для всех частот, представленных во входном сигнале и АЧХ. После получения спектра выходного сигнала можно восстановить сигнал как временную зависимость с помощью формулы обратного преобразования Фурье. Заметим главное: моделирование прохождения сигнала через динамический объект свелось к операции умножения двух переменных, точнее, к операции поэлементного умножения вектора одних переменных на вектор других переменных. Схема преобразования показана на рис. 8.3. Рис. 8.3. Схема процедуры преобразования сигнала при использовании метода Фурье Если бы временной сигнал проходил через звено, которое во временной области представлено дифференциальным уравнением, то пришлось бы его интегрировать, что, конечно, приводит к погрешностям результата. В частотной области достаточно перемножить значения коэффициентов ряда Фурье сигнала и звена при одинаковых частотах. Очевидно, что достоинством метода является замена дифференциальных уравнений модели на алгебраические. Разумеется, данный подход может быть использован только для объектов, у которых известен вид передаточной функции. (Кстати, для неизвестных случаев АЧХ может быть получена численным разложением в ряд.) В процессе моделирования набора объектов для преобразования сигнала (например, протяженных трактов радиоэлектронных устройств) иногда приходится применять прямое и обратное преобразование Фурье неоднократно. На практике последовательные блоки часто называют каскадами. Пусть мы имеем радио-электронное устройство (РЭУ), состоящее из 5 блоков (см. рис. 8.4). Блоки 1, 2, 4, 5 — линейные и представлены соответствующими известными АЧХ; блок 3 — нелинейный, поэтому АЧХ для него неизвестна. Примером линейного блока может служить апериодическое звено, колебательное звено и т. д. (см. Лекцию 05. Динамические регрессионные модели, заданные в виде передаточной функции). Примером нелинейного блока может служить устройство ограничения сигнала (срез) по амплитуде. Как видно из рис. 8.4, сначала входной сигнал X(t) прямым преобразованием Фурье переводится в частотную область и проходит в виде спектра через АЧХ 1 и 2 линейного блока, затем обратным преобразованием Фурье сигнал после 2 блока переводится во временную область. Проходим нелинейный блок 3 во временном представлении. Результат работы блока 3 снова преобразуем прямым преобразованием Фурье в частотную область и проходим через АЧХ блоков 4 и 5. В конце полученный спектр преобразуется с помощью обратного преобразования Фурье во временную область, — вид сигнала, Z(t), является результатом моделирования. Рис. 8.4. Пример моделирования тракта, содержащего нелинейные блоки, с использованием метода Фурье Метод, который мы рассмотрели, является одним из самых быстродействующих. Это связано с заменой операций интегрирования и дифференцирования, встречающихся в моделях динамических звеньев, на операции сложения и умножения при переходе в частотную область. Такая процедура обеспечивает точность и быстродействие модели. Для метода важно, с какой частотой вы дискретизируете сигнал при разложении в ряд Фурье. Если частота дискретизации мала, то есть отсчеты в сигнале следуют редко, с большими интервалами, то часть сигнала остается потерянной, так как между отсчетами может оказаться резко возросший и опавший пик, информация о котором пропадет. То есть говорят, что малая частота дискретизации срезает высокие частоты в сигнале. (Пик — это и есть высокочастотная составляющая, которая может быть потеряна). По теореме Котельникова, чтобы не потерять соответствующую гармонику, требуется дискретизировать сигнал с частотой не менее чем в 2 раза большей, чем самая высокая частота из представленных в аналоговом сигнале: где Wдискр. = 1/Δtдискр. — частота дискретизации, Wmax — максимальная частота, присутствующая в сигнале 9. Лекция 09. Оценка качества модели Оценка качества показывает, насколько теоретические вычисления по построенной модели отклоняются от экспериментальных данных. Наличие связи двух переменных называется корреляцией. Если оценка качества применяется до исследования, то она решает задачу: есть ли связь между входом X и выходом Y и оценивает силу этой связи. 📹 Видео ЛОДУ 2 порядка c постоянными коэффициентамиСкачать Линейное неоднородное дифференциальное уравнение второго порядка с постоянными коэффициентамиСкачать Лекция №3. Апериодическое звено. Прокопенко В. А.Скачать 23) Построение Л.А.Ч.Х. и Л.Ф.Ч.Х. системы по её передаточной функцииСкачать 16. Линейные неоднородные дифференциальные уравнения 2-го порядка с постоянными коэффициентамиСкачать Теория автоматического регулирования. Лекция 3. Временные характеристики САУСкачать Частотные характеристики | Утро с теорией управления, лекция 5Скачать 15. Линейные однородные дифференциальные уравнения второго порядка с постоянными коэффициентамиСкачать Урок 228 (осн). ЛР Определение коэффициента ослабления гамма-излучения различными материаламиСкачать #78. Роль иудейского эгрегора и ответственность иудаизма в развитии земной цивилизации.Скачать Типовые динамические звенья | Вечер с теорией управления, вебинар 3Скачать Реакция на результаты ЕГЭ 2022 по русскому языкуСкачать Линейное неоднородное дифференциальное уравнение с постоянными коэффициентами 4y''-y=x^3-24x #1Скачать c04 3, Динамические звенья 2: типовые звенья и их АФЧХСкачать 12) ТАУ для чайников. Часть 4.4. Интегрирующее звено.Скачать Задача Коши ➜ Частное решение линейного однородного дифференциального уравненияСкачать Поделиться или сохранить к себе: Смотрите также Сера — химические свойства, получение, соединения. Нитрат кальция: способы получения и химические свойства Кальций: способы получения и химические свойства Гидроксид натрия: способы получения и химические свойства Гидроксид кальция: способы получения и химические свойства Получение бензола Разница между глюкозой и сахарозой Гидроксид калия: способы получения и химические свойства