
- •2. Некоторые методы анализа систем
- •2.1. Построение и упрощение моделей объекта (Лекция 3)
- •2.1.1. Упрощение моделей объекта
- •2.1.1.1. Замена нестационарной модели набором стационарных
- •2.1.1.2. Линеаризация модели
- •2.1.2. Пример построения модели объекта
- •2.1.2.1. Характеристика объекта
- •2.1.2.2. Разработка математической модели
- •Линеаризация
- •2.1.2.4. Линейная модель объекта
- •2.2. Временные характеристики систем управления (Лекция 4)
- •2.2.1. Линейные нестационарные системы
- •2.2.1.1. Общий вид описания системы
- •2.2.1.2. Решение однородного уравнения
- •2.2.1.3. Решение неоднородного уравнения
- •2.2.1.4. Фундаментальная матрица системы и её свойства
- •2.2.1.5. Весовая и переходная матрицы системы
- •2.2.2. Линейные стационарные системы
- •2.2.2.1. Фундаментальная матрица стационарной линейной системы
- •2.2.2.2. Весовая и переходная матрицы стационарной системы
- •2.3. Передаточная функция (Лекция 5)
- •2.3.1. Некоторые операторы
- •2.3.1.1. Дифференциальный оператор
- •2.3.1.2. Оператор Лапласа
- •2.3.2. Передаточная функция и резольвента
- •2.3.2.1. Фундаментальная матрица (резольвента)
- •2.3.2.2. Матричная передаточная функция
- •2.3.2.3.Простейший пример определения резольвенты (Пример п2.2)
- •2.3.2.4. Метод Фаддеевой (Сурье )определения резольвенты
- •2.3.2.5. Пример применения метода Фаддеевой (Пример 2.3)
- •2.4. Частотные характеристики (функции) систем (Лекция 6)
- •2.4.1. Частотная переходная функция
- •2.4.1.1. Напоминание о представлении комплексных чисел
- •2.4.1.2. Частотная передаточная функция
- •2.4.1.3. Смысл компонент матричной частотной функции (Пример2.4)
- •2.4.1.4. Виды используемых частотных характеристик
- •2.4.1.5. Пример вычисления характеристик (Пример2.5)
- •2.4.2. Логарифмические частотные характеристики
- •2.4.2.1. Смысл логарифмических частотных характеристик
- •2.4.2.2. Определение логарифмических частотных характеристик
- •2.4.2.3. Асимптотические логарифмические частотные характеристики
- •2.4.2.4. Иллюстрация построения асимптотических характеристик (Пример 2.6)
- •2.5. Структурные схемы систем (Лекция 7)
- •2.5.1. Схемы соединения звеньев
- •2.5.1.1. Представление звеньев и связей в виде структурных схем
- •2.5.1.2. Последовательное соединение звеньев
- •2.5.1.3. Параллельное соединение
- •2.5.1.4. Соединение с обратной связью
- •2.5.1.5. Пример определения матрицы возвратной разности (Пример 2.7)
- •2.5.2. Структурные преобразования линейных систем
- •2.5.2.1. Назначение и содержание структурных преобразований
- •2.5.2.2. Правила структурных преобразований линейных систем
- •2.5.2.3. Дополнительные правила для стационарных линейных систем
- •2.5.2.4. Иллюстративный пример (Пример2.8)
- •2.6. Анализ устойчивости систем управления (Лекция 8)
- •2.6.1. Определение устойчивости систем
- •2.6.1.1. Номинальное состояние и понятие устойчивости
- •2.6.1.2. Определения устойчивости решений
- •2.6.1.3. Устойчивость линейных дифференциальных систем
- •2.6.1.4. Пример смесительного бака (Пример2. 9)
- •2.6.2. Устойчивость линейных стационарных систем
- •2.6.2.1. Представление реакции системы с различными собственными числами
- •2.6.2.2. Представление реакции системы с кратными собственными числами
- •2.7.1.2. Необходимые условия
- •2.7.1.3. Достаточные условия
- •2.7.1.4. Пример применения алгоритма Раусса (Пример 2.11)
- •2.7.2. Частотные критерии
- •2.7.2.1. Нестрогое обоснование частотных критериев
- •2.7.2.2. Критерий Михайлова
- •2.7.2.3. Критерий Найквиста
2.3.2.5. Пример применения метода Фаддеевой (Пример 2.3)
Пример
П2.3: дана
система
Найти Ф(p).
Решение:
Способ №1 (непосредственное вычисление).
,
т.к.
,
тогда
,
,
;
Способ №2 (метод Сурье).
n=2,
k=1,
α2
,
α1,
.
α0,
.
.
2.4. Частотные характеристики (функции) систем (Лекция 6)
2.4.1. Частотная переходная функция
2.4.1.1. Напоминание о представлении комплексных чисел
Частотные характеристики характеризуют зависимость реакции системы от частоты периодического входного воздействия. Зависимость от частоты может характеризоваться изменением амплитуды и фазы входного сигнала. Для представления этих двух важнейших характеристик удобно использовать комплексные числа, которые, как известно, содержат в себе обе эти характеристики.
Поэтому вспомним формы представления комплексных чисел и используемые при этом термины. Из курса математики известно, что любое комплексное число имеет вид:
z = a + jb,
где a=Re z – его действительная часть (от латинского Realis -действительный),
b=Im z – его мнимая часть (от латинского Imaginarias - мнимый).
На комплексной плоскости комплексное число представляется в виде, показанном на рис. 2.2, где a и b – действительная и мнимая части числа, A – его модуль, а φ – аргумент.
Рис. 2.2. Представление числа на комплексной плоскости
Комплексное число можно представить также в показательной и тригонометрической форме:
z = A(z)ejarg z = A(z) (cos j + j sin j),
где
A(z)
=
- модуль комплексного числа, который в
теории управления называется
амплитудой;
arg z = arctg b/a = φ – аргумент, который в теории управления называется фазой.
Выражения тригонометрических функций через показательные можно получить с помощью формул Эйлера:
cos φ = (ejφ + e- jφ) /2, sin φ = (ejφ - e- jφ) /2j.
2.4.1.2. Частотная передаточная функция
Периодическое входное воздействие можно представить в показательной форме:
w(t) = wm ejwt , t ³0. (2.25)
Пусть оно действует на систему
=
Ax(t)
+ Bw(t),
(2.26)
y(t) = Cx(t). (2.27)
Так как уравнение (2.26) неоднородное, то его решение состоит из суммы общего решения однородного уравнения и частного решения неоднородного. Общее решение было ранее определено. Частное решение неоднородного уравнения может быть определено в виде:
xч(t) = xm ejwt, xm = const. (2.28)
Подставив (2.25) и (2.28) в (2.26), получим:
xm jw ejwt = A xm ejwt + B wm ejwt,
(jwI -A)xm ejwt = B wm ejwt ,
xч (t) =(jwI -A)-1 B wm ejwt . (2.29)
Общее решение однородного уравнения в показательной форме имеет следующий вид:
xo (t) = eAtx(0). (2.30)
Сумма (2.29) и (2.30) представляет собой соответствующее воздействию (2.25) решение уравнения (2.26):
x(t) = eAtx(0) + (jwI -A)-1B wm ejwt . (2.31)
Первое слагаемое этого решения характеризует переходную погрешность, а второе установившуюся.
При x(0) = 0 получим:
x(t) = (jwI -A)-1B wm ejwt.
Подставим это в (2.27) и получим для вектора выходных величин:
y(t) = Cx(t) = C(jwI -A)-1 B wm ejwt . (2.32)
Функция, позволяющая определить значения компонент вектора выходных величин по значениям компонент вектора входных величин при разных значениях частоты, называется частотной переходной характеристикой системы. Из (2.32) следует её вид:
H(jw) = C(jwI -A)-1 B (2.33)
Её можно получить заменой в передаточной функции аргумента p на jw. В рассматриваемом случае она представляет собой матрицу, количество строк в которой совпадает с размерностью вектора y, а количество столбцов - с размерностью вектора w.
Если H(jw) известно и известно входное воздействие, то можно найти выходные величины, используя известные правила перемножения комплексных чисел.