3. ЧАСТОТНЫЕ ХАРАКТЕРИСТИКИ СИСТЕМ АВТОМАТИЧЕСКОГО УПРАВЛЕНИЯ (АФЧХ, ЛАХ, ФЧХ) ч. 3.1

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

Лекции по курсу «Управление Техническими Системами», читает Козлов Олег Степанович на кафедре «Ядерные реакторы и энергетические установки», факультета «Энергомашиностроения» МГТУ им. Н.Э. Баумана. За что ему огромная благодарность.


Данные лекции только готовятся к публикации в виде книги, а поскольку здесь есть специалисты по ТАУ, студенты и просто интересующиеся предметом, то любая критика приветствуется.


В предыдущих сериях:
1. Введение в теорию автоматического управления.
2. Математическое описание систем автоматического управления 2.1 — 2.3, 2.3 — 2.8, 2.9 — 2.13


В этом разделе мы будем изучать частотные характеристики, тема сегодняшней статьи:
3.1. Амплитудно-фазовая частотная характеристика: годограф АФЧХ, ЛАХ, ФЧХ


Будет интересно познавательно и жестко.




3.1. Амплитудно-фазовая частотная характеристика: годограф АФЧХ, ЛАХ, ФЧХ


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



Рис. 3.1.1 – Схематическое представление синусаидального воздействия

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

$sin(\omega \cdot t+\phi) = sin \left[ \omega(t + \frac{\phi}{\omega}) \right] = sin \left[ \omega(t + \Delta t) \right];$



Рисунок 3.1.2 – График представления синусаидального воздействия

$\phi $ — сдвиг фазы (не редконазывают — фаза);
$A$ — амплитуда;
$A \equiv A(\omega); \phi \equiv \phi( \omega) $ т.е. амплитуда на выходе звена(системы) и сдвиг фазы зависят от частоты входного воздействия x(t).


Используем показательную форму записи функции единичного гармонического воздействия и отклика на это воздействие (рис. 3.1.1):

$sin(\omega \cdot t) = e^{i \cdot \omega \cdot t}$


$$display$$\left \{ \begin{align} x(t)&= e^{i \cdot \omega \cdot t} = cos (\omega \cdot t) + i \cdot sin(\omega \cdot t)\\ y(t)&= A\cdot sin(\omega \cdot t+ \phi) \Rightarrow A \cdot e^{i(\omega \cdot t+\phi)}=A \cdot e^{i\cdot \omega \cdot t} \cdot e^{i \cdot\phi} = A(\omega) \cdot e^{i \cdot \omega \cdot t} \cdot e^{i \cdot \phi(\omega)} \end{align} \right.\ \ \ \ \ \mathbf{(3.1.1)} $$display$$


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

$T^2_2y''(t)+T_1y'(t) = K[\tau \cdot x'(t)+x(t)].\ \ \ \ \ \mathbf{(3.1.2)} $


В показательной форме:

$(T^2_2\cdot s^2+T_1 \cdot s+1)\cdot Y(s) = K[\tau \cdot s+1] \cdot X(s)$


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

$W(s)= \frac{K(\tau \cdot s +1)}{T^2_2 \cdot s^2+T_1 \cdot s+1}\ \ \ \ \ \mathbf{(3.1.3)} $


Запишем в показательной $x, x', y, y', y''$ форме используя соотношения 3.1.1:

$x(t)=e^{i \cdot \omega \cdot t};\\ x'(t)=i \cdot \omega \cdot e^{i \cdot \omega \cdot t};\\ y(t) = A \cdot e^{i \cdot \omega \cdot t} \cdot e^{i \cdot \phi};\\ y'(t)=A \cdot i \cdot \omega \cdot e^{i \cdot \omega \cdot t} \cdot e^{i \cdot \phi};\\ y''(t)=A \cdot (i \cdot \omega)^2 \cdot e^{i \cdot \omega \cdot t} \cdot e^{i \cdot \phi};$


Подставим эти соотношения в (3.1.1) получим:

$T^2_2\cdot A \cdot (i \cdot \omega)^2 \cdot e^{i \cdot \omega \cdot t} \cdot e^{i \cdot \phi} +T_1\cdot A \cdot i \cdot \omega \cdot e^{i \cdot \omega \cdot t} \cdot e^{i \cdot \phi} + A \cdot e^{i \cdot \omega \cdot t} \cdot e^{i \cdot \phi}= \\ =K[\tau \cdot i \cdot \omega \cdot e^{i \cdot \omega \cdot t}+e^{i \cdot \omega \cdot t}] \Rightarrow \\ A \cdot e^{i \cdot \phi} \cdot [T^2_2\cdot (i \cdot \omega)^2+ T_1 \cdot (i \cdot \omega)+1] = K[\tau \cdot (i \cdot \omega)+1] \\ $


Поскольку $A \equiv A(\omega); \phi \equiv \phi( \omega) $ (амплитуда на выходе звена(системы) и сдвиг фазы зависят от частоты входного воздействия), то можно записать:

$A(\omega) \cdot e^{i \cdot \phi(\omega)}=\frac{K[\tau \cdot i \cdot \omega+1]} {T^2_2(i \cdot \omega)^2+ T_1(i \cdot \omega)+1}\ \ \ \ \ \mathbf{(3.1.4)}$


если вспомнить, что в преобразования Лапаласа $s = i \cdot \omega$, то:

$\frac{K[\tau \cdot i \cdot \omega+1]} {T^2_2(i \cdot \omega)^2+ T_1(i \cdot \omega)+1} = \frac{K[\tau \cdot s+1]} {T^2_2 \cdot s^2+ T_1\cdot s+1} =W(s);$


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

$A(\omega) \cdot e^{i \cdot \phi(\omega)} = W(i \cdot \omega) =W(s)\ \ \ \ \ \mathbf{(3.1.5)}$


$W(i \cdot \omega) $ — Амплитудно-фазовая частотная характеистика (АФЧХ)
Иногда $W(i \cdot \omega) $ называют частотной передаточной функцией.
Модуль АФЧХ=$mod W(i \cdot \omega) $ тождественно равен амплитуде выходного сигнала:

$A(\omega) = |W(i \cdot \omega)|\ \ \ \ \ \mathbf{(3.1.6)}$


Сдвиг фазы выходного сигнала:

$\phi(\omega) = arg( W(i \cdot \omega))\ \ \ \ \ \mathbf{(3.1.7)} $


Обычно АФЧХ $W(i \cdot w)$ изображается на комплексной плоскости. Формулы (3.1.6) и (3.1.7) позволяют изобразить $W(i \cdot w)$ в полярных координатах $(r, \phi).()$
Так же можно изображать в традиционных декартовых координатах:

$W(i \cdot \omega) = \underbrace { U(\omega)}_{Re} + i \cdot \underbrace { V(\omega)}_{Im} \ \ \ \ \ \mathbf{(3.1.8)} $


Если использовать для представления W(s) форму W(s)=K·N(s)/L(s), где L(s)- полиномы по степеням s, (причем свободные члены равны 1), а К – общий коэффициент усиления звена (системы), то

$W(i \cdot\omega) = \frac{K \cdot N(i \cdot\omega)}{L(i \cdot \omega)} \Rightarrow | W (i \cdot \omega)| = \frac{K|N(i \cdot \omega)|}{|L(i \cdot \omega)|} = A(\omega)\ \ \ \ \ \mathbf{(3.1.9)} $


Сдвиг фазы $\phi(\omega)$ можно определить по виду многочленов $N(i \cdot \omega)$ и $L(i \cdot \omega)$ (см. формулу (3.1.9)) т.е. как разность фаз (аргументов) числителя и знаменателя:

$\phi(\omega) = arg(N(i \cdot\omega))-arg(L(i \cdot \omega))\ \ \ \ \ \mathbf{(3.1.10)}$


Постоим АФЧХ для «абстрактного» звена (системы) с передаточной функцией:

$W(s) = W(i \cdot \omega) = W(s) |_{s =i \cdot \omega} = \frac{K \cdot N(i \cdot \omega)}{L(i \cdot \omega))}$


Подставляя в формулу различные значения $\omega$ получаем набор векторов, на комплексной плоскости


Рисунок 3.1.3 – Годограф абстрактного звена.

Рассмотрим действительную и мнимую части полученных векторов Из рисунка 3.1.3 видно, что:

$$display$$\left \{ \begin{align} u( \omega)&= A( \omega) \cdot cos( \phi(\omega)) \\ v(\omega)&= A( \omega) \cdot sin( \phi(\omega)) \end{align} \right.$$display$$


Амплитуда и сдвиг фазы рассчитываются, для векторов соответсвующих положительным частотам и лежащих в 4 квадранте $\omega_1,\omega_2, \omega_3$ по формулам:

$$display$$\left \{ \begin{align} A( \omega)&= \sqrt {u^2( \omega)+v^2( \omega)} ,\\ \phi(\omega)&= arctg \frac{v(\omega)}{u(\omega)} .\end{align} \right.$$display$$


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

$\phi(\omega) = -\pi \cdot j + arctg \frac{v(\omega)}{u(\omega)}\ \ \ \ \ \mathbf{(3.1.11)}$


где:
j = 0, 2, 3, 4..., если вектор в I и IV квадрант;
j = 1, 3, 4, 4..., если вектор в II и III квадранте.

Во всех технических системах отклик системы как правило отстает от входного воздействия, то есть сдвиг фазы всегда отрицательный. Исходя из формулы 3.1.10 степень полинома L(s) выше, чем полинома N(s). Поскольку обычно степень полинома L(s) выше, чем полинома N(s), то с увеличением частоты на входе в звено (в систему) сдвиг фазы обычно отрицателен, т.е. сигнал на выходе звена еще больше отстает по фазе от входного сигнала при увеличении частоты.
В предельном случае, если частота растет до бесконечности, мы можем вообще не получить выходного воздействий. Обычно при ω→ ∞ величина амплитуды на выходе звена стремится к 0, то есть lim A(ω→∞) = 0.


$W(i \cdot \omega)$ при замене $\omega$ на $-\omega$ имеет зеркальное изображение.


Анализируя годографы АФЧХ при $\omega$ > 0 (сплошная линия на рисунке 3.1.3) и при $\omega$ < 0 (пунктирная линия) видим, что:
$u(\omega) = u(-\omega)$ – четная функция, следовательно график симметричен относительно оси ординат, а
$v(\omega) = -v(-\omega)$ – нечетная функция и ее график центрально-симметричен относительно начала координат.



Рисунок 3.1.4 – «Зеркальная» симметрия относительно оси ординат.


Рисунок 3.1.5 – «Центральная» симметрия относительно начала координат.

Кроме анализа свойств звена (системы) по годографу АФЧХ широкое распространение получили анализ логарифмической амплитудной характеристики (ЛАХ) и фазочастотной характеристики (ФЧХ).


ЛАХ определяется как Lm(ω)=20lgA(ω).


Поскольку зачастую удобнее использовать десятичные логарифмы (lg), чем натуральные(ln), в теории управления (также и в акустике) значительно чаще используется специальная единица – децибел (1/10 часть Бела):
+1Бел – единица, характеризующая увеличение в 10 раз.
+1дБ (децибел) – соответствует увеличению в $\sqrt[10]{10}$ раз.


В формуле Lm(ω)=20lgA(ω) величина Lm(ω) измеряется также в децибелах. Происхождение множителя 20 таково: A(ω) – амплитуда, линейная величина, а мощность — квадратичная величина (например, напряжение в сети измеряется в Вольтах, а мощность ($N = \frac{U^2}{R}$) пропорциональна квадрату напряжения, поэтому в формуле для Lm(ω) стоит множитель 20 (чтобы привести ЛАХ (Lm(ω)) к традиционной мощностной характеристике).


Если $Lm(ω_1)$ больше $Lm(ω_2)$ на 20 дБ, то это означает, амплитуда $А(ω_1)$ больше амплитуды $А(ω_2)$ в 10 раз, $\frac{А(ω_1)}{ А(ω_2)} =10$

Окончательно: Lm(ω)=20lg│W(iω)│= 20lgA(ω)


Из этого следует, что +1 децибел (+1 дБ) соответствует увеличению амплитуды в $\sqrt[20]{10}$ раз (очень малая величина); -1 дБ – уменьшение амплитуды в $\sqrt[20]{10}$ раз.


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



Рисунок 3.1.6 – пример графика АЧХ


Рисунок 3.1.7 – пример графика ФЧХ

Учитывая, что “ω” обычно изменяется на порядки и значение A(ω) – также на порядки, график Lm(ω) строится, фактически, в логарифмических координатах, т.е. Lm(ω) =Lm(lg(ω)), например:



Рисунок 3.1.8 – пример графика ЛАХ

Наклон (– 40 дБ/дек) соответствует уменьшению амплитуды в 100 раз при увеличении частоты в 10 раз.



Рисунок 3.1.9 – пример графика ЛФЧХ

Рассмотренные характеристики Lm(ω), то есть ЛАХ и ФЧХ имеют широкое распространение при анализе динамических свойств звена (системы), например, при анализе устойчивости САР (см. раздел “Устойчивость систем автоматического управления”).



Рисунок 3.1.10 – пример ЛАХ и ФЧХ для сложной системы

Пример 1


В качестве примера построим АФЧХ для демпфера, модель которого, разобрана в этой статье.... Добавим на схему блок «Построение частотных характеристик» качестве входа возьмем возмущающее воздействие, в качестве выхода положение положение груза. Для наглядности иллюстрации примем в качестве выхода положение в миллиметрах (х1000), поскольку модель у нас размерная и результат получается в метрах уже достаточно маленьким примерно 0.004 метра. см. рис. 3.11



Рисунок 3.1.11 Схема для построения частотных характеристик.

Параметры блока «Построение частотных характеристик» приведены на рисунке 3.1.12, для илюстации зависимости АЧХ и ЛАХ. Результат работы блока — график с выбранными параметрам изображен на рисунке 3.1.13:



Рисунок 3.1.12 – Парамеметры блока «Частотные характеристики».


Рисунок 3.1.13 – Частотные харатктеристики в АЧХ, ЛАХ, ФЧХ в линейном масштабе по ω.

Анализ графика в линейном масштабе по ω чаще всего не очень удобен, поскольку весь график собиается в узкой области, а дальше график абсолютной амплитуды практически сливается с 0. Если мы хотим исследовать частоты хотя бы до 1000 Гц, мы увидем практически вертикальные и горизонтальные прямые. Изменения масштаба шакалы АЧХ и ω на логарифмический позволяет более удобно исследовать частотные характеристики (см. рис. 3.1.14).


На рисунке 3.1.14 представлены частотные характеристики демпфера в логарифмическом масштабе и иллюстарция соотношения между абсолютной величиной амплитуды АФЧХ и ЛАХ в децибелах.



Рисунок 3.1.14 – Частотные харатктеристики в АЧХ, ЛАХ, ФЧХ в логарифмитическом масштабе по ω.

Пример 2


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

Для начала посмотрим модель в виде блоков.


Модель подготовленная для анализа представлена на рисунке 3.1.15. В отличие от исходной модели, описанной ранее, входное воздействие задается блоком ступенька с скачком с 0 до 1 на 10 секунде расчтеа. В блоке линейная функции происходит пересчет сигнала ступенька:
0 — соответствует 200 бар в камере, (конечное состояние в предыдущем примере),
1 — соответствует 400 бар в камере.
Это сделано для того, что бы можно было подавать синусоидальный сигнал и не получать отрицательное давление в камере плунжера. Так же для наглядности графика, мы усиливаем выходное перемещение переводя его из метров в миллиметры.



Рисунок 3.1.15 – Модель гидравлического демпфера.

Частотные характеристики рассчитанные в конце рассчета приведены на рисунке 3.1.16. Видно что характеристики отличаются от простого пружинного демпфера сравните с 3.1.14



Рисунок 3.1.16 – Частотные хараткеристики гидравлического депфера

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


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


Модель гидравлического демпфера собранного из библиотечных блоков SimInTech представлена на рисунке 3.1.7


Рисунок 3.1.17 – Модель гидравлическогго демпфера для рассчета частотных характеристик.

Расчеты с моделью показывают, что при сохранении общего вида графиков, значения полученные при для «подробной модели» отличаются от линеаризованной модели, (см. рис. 3.18 — 3.19)



Рисунок 3.1.18 – АЧХ подробной модели привода, полученная прямым моделированием.


Рисунок 3.1.19 – ЛАХ подробной модели привода, полученная прямым моделированием.

Использование прямого моделирования, для получения характеристик, является более надежным способом и работает не только с линейными моделями, но может быть использован для построения характеристик некоторых реальных объектов, если их можно подключить к среде моделирования и воздействовать в реальном режиме времени. Однако затраты на вычисления значительно больше. Например для получения характеристик демпфера пришлось выполнить процесса 40 000 секунд модельного времени, на обычном компьютере это заняло порядка 35 минут. График процесса перемещения плунжера в процессе вычисления характеристик приведен на рисунке 3.1.20.


Рисунок 3.1.20 – Перемещения плунжера в процессем моделирования.

Блок «Гармонический анализатор»имеет выходы:
Re(w*t) – текущее значение действительной части амплитудно-фазовой частотной характеристики исследуемой системы;
Im(w*t) – текущее значение мнимой части амплитудно-фазовой частотной характеристики.
Это позволяет построить годограф исследуемой системы с помощью фазового портрета. (см. рис. 3.1.21)



Рисунок 3.1.21 – Годограф системы гидравличесого демпфера.

Модели, использованные для илюстарции в лекции можно взять здесь...

Продолжение следует.

Источник: https://habr.com/ru/post/526190/


Интересные статьи

Интересные статьи

Вновь обращаемся к экспертизе, собранной в нашем «Мире Hi-Fi». Ранее мы посмотрели, где производят аудиотехнику, рассказали о теории и практике DIY-акустики и настройке д...
С декабря 2020 года автобусный билет будет снабжаться метаданными, которые можно будет проверить через МКПП — это единая система, действующая по всей стране. То есть пассажир будет уверен...
В данной статье я бы хотел познакомить читателей с одним из проектов Apache Software Foundation сообщества — NlpCraft. NlpCraft — библиотека с открытым исходным кодом, предназначенная для...
Давайте поговорим о high-load, high-extensible, ООП, DDD & EDA. Речь пойдет о разнице между high load & high extensibility — это разные цели и задачи. Для которых нужны принц...
Я глубоко погружен в WSL (Windows Subsystem для Linux), и теперь, когда WSL2 доступен в Windows Insiders, это отличное время, чтобы по-настоящему изучить доступные опции. Очень интересная фича, к...