Примеры решения задач по теории автоматического управления

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

Если что-то непонятно вы всегда можете написать мне в воцап и я вам помогу!

Одномерные линейные непрерывные системы

К оглавлению…

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

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

Эквивалентная передаточная функция параллельно соединенных звеньев равна сумме передаточных функций этих звеньев (с учетом знака входа сумматора на пути сигнала).

Путь от входа к выходу системы называется прямой связью, от выхода ко входу — обратной связью. Если сигнал на пути меняет знак (обычно на инвертирующем входе сумматора), обратная связь называется отрицательной (ООС), если не меняет знак — положительной (ПОС). Замкнутый путь называется контуром, например, замкнутый контур обратной связи (ЗКОС). Эквивалентная передаточная функция соединения с обратной связью равна дроби, в числителе которой записана ПФ звена на прямом пути, а в знаменателе — единица минус произведение ПФ звеньев по замкнутому контуру обратной связи. Величина . называется определителем ЗКОС.

Особенности этого вида соединения звеньев:

  • если в системе есть хоть одна обратная связь, передаточная функция системы будет всегда представлять собой дробь;
  • знак перед произведением ПФ звеньев в знаменателе (в определителе ЗКОС) обычно противоположен знаку обратной связи.

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

По Мейсону структурная схема может быть описана целиком, без деления на звенья. Передаточная функция многоконтурной системы образует дробь, числитель которой равен сумме произведений передаточных функций прямых путей на совокупные определители ЗКОС, не касающихся этих путей, а знаменатель — единица минус сумма произведений определителей несоприкасающихся ЗКОС и передаточных функций общих ЗКОС. Следует внимательно относиться к ветвям, которые заходят извне в контур ОС, т.к. они могут образовывать неявные прямые пути по цепям обратных связей.

Предмет теория автоматического управления тау

Пример №1

Определить передаточную функцию схемы (рисунок 1.1,а).

Решение:

Видно, что без преобразований нельзя начинать сворачивать схему, в частности, нельзя объединить звенья и , как последовательно включенные, из-за связи в точке . Перенесем ветвь из узла в узел (рисунок 1.1,6).

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

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

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

Пример №2

Определить передаточную функцию схемы (рисунок 1.2).

Решение:

Применим правило Мейсона. В системе имеются обратные связи, поэтому ПФ представляет собой дробь. Прямой путь от входа к выходу только один, его касаются все пять ЗКОС, поэтому в числителе ПФ пишем просто произведение . Знаменатель начинаем описывать с несоприкасающихся контуров — контур I не имеет общих точек с контуром III и вложенным в него контуром И, поэтому записываем сначала произведение их определителей. Контур IV соприкасается с контурами I и III, поэтому просто добавляем произведение звеньев по нему , но умножаем его на определитель контура И, так как этот ЗКОС не имеет общих точек с IV. И в конце просто добавляем произведение звеньев контура V, поскольку он соприкасается со всеми остальными ЗКОС

Дифференциальное уравнение

К оглавлению…

Поведение линейных, непрерывных, стационарных систем с сосредоточенными параметрами описывается во времени обыкновенным дифференциальным уравнением (ОДУ) с постоянными коэффициентами

где слева — выходная функция и ее производные (результат), справа — входная функция и ее производные.

Для записи передаточной функции используется комплексная переменная Лапласа (иногда обозначаемая символом ). Чтобы получить ПФ, достаточно в ОДУ заменить производные на в соответствующей степени, отбросить символы функций и и разделить многочлен правой части дифференциального уравнения на многочлен левой части.

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

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

При переходе от разомкнутой системы к замкнутой, охваченной общей единичной отрицательной обратной связью (ООС), достаточно к знаменателю ПФ разомкнутой системы добавить ее числитель, чтобы получить ПФ замкнутой системы.

Пример №3

Определить передаточную функцию объекта регулирования, модель которого задана дифференциальным уравнением

Решение:

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

Решение задач по теории автоматического управления

Пример №4

При единичном скачке на входе реакция звена описывается функцией . Найти передаточную функцию звена.

Решение:

Преобразуем по Лапласу входной и выходной сигналы, пользуясь таблицей соответствия оригиналов и изображений (приложение ). Изображение входного воздействия равно , изображение реакции звена после приведения к общему знаменателю

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

Делим изображение реакции на изображение входного воздействия и получаем передаточную функцию звена

Пример №5

Система имеет нуль -3, комплексные сопряженные полюса и коэффициент усиления . Определить ПФ системы после её замыкания единичной ООС.

Решение:

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

Добавляя к знаменателю числитель, получаем ПФ замкнутой системы

Принципиальная схема

К оглавлению…

Если анализируется принципиальная электрическая схема, передаточная функция составляется с учетом известных закономерностей работы таких схем. Для индуктивных элементов (катушек, дросселей) операторное реактивное сопротивление равно для емкостных элементов , где — индуктивность (Генри), -емкость (Фарад), — комплексная переменная Лапласа.

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

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

Пример №6

Определить передаточную функцию схемы (рисунок 1.12).

Решение:

Схема представляет собой делитель напряжения с коэффициентом

поэтому передаточная функция равна

Пример №7

Определить передаточную функцию схемы (рисунок 1.13).

Решение:

Эквивалентное операторное сопротивление в цепи отрицательной обратной связи равно сумме

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

Пример №8

Составить структурную схему по дифференциальному уравнению объекта

Решение:

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

Затем составляют структурную схему, используя блоки интегрирования (т.е. деления на переменную Лапласа ), их число равно порядку системы (в данном случае трём). С выхода каждого интегратора организуют обратные связи к общему (входному) сумматору с инвертирующим входом, начиная с коэффициента при -1 производной. С выхода интеграторов организуют связи с коэффициентами из правой части ОДУ к выходному сумматору объекта (если производные здесь отсутствуют, то выходной сумматор не нужен, а блок с коэффициентом можно поместить и на выходе, и на входе системы, до главного сумматора). Полученная схема показана на рисунке 1.14.

Пример №9

Определить порядок объекта, записать его дифференциальное уравнение по передаточной функции

Решение:

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

Временные характеристики. Реакция на произвольное воздействие

К оглавлению…

Для решения дифференциального уравнения (нахождения реакции системы) с помощью преобразования Лапласа необходимо:

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

Рекомендуется:

а) перед вычислением корней обязательно нормировать ПФ по старшему коэффициенту при знаменателя;

б) не сокращать существующие нули и полюса с положительной действительной частью, ведущие к неустойчивости системы, если их части не являются целыми числами; остальные нули и полюса могут быть сокращены перед переходом во временную область;

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

г) комплексные сопряженные корни представлять одним общим квадратным трехчленом.

После разложения на простые дроби и вычисления вычетов полезно проверить правильность результата. Первое правило проверки -сумма дробей правой части должна быть равна изображению в левой части равенства. Второе правило проверки — сумма всех составляющих оригинала при (начальное значение оригинала) в соответствии со свойствами преобразования Лапласа должна быть равна .

Пример №10

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

Решение:

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

Пример №11

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

Решение:

Начальное значение оригинала (при ) вычисляется как предел

для производной по времени -го порядка от функции производится умножение изображения на т.е.

Поэтому

Конечное значение оригинала (при ) для устойчивых систем также вычисляется как предел

Для полной записи оригинала разлагаем изображение на простые дроби в соответствии с полюсами, находим вычеты и в полюсах методом подстановки полюсов (приложение Б)

По таблице соответствия оригиналов и изображений (приложение А) записываем оригинал в виде формулы . Проверка: при значение оригинала равно нулю, при соответственно 3.

Переходная и импульсная функции

К оглавлению…

К типовым функциям времени (реакциям системы) относятся переходная и импульсная переходная (весовая) функции.

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

Импульсной (весовой) функцией или называется реакция системы на единичный импульс при нулевых начальных условиях. Она является оригиналом передаточной функции.

Поскольку всегда то

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

Связь между импульсной и переходной функциями определяется соотношением

откуда

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

Решение задач по ТАУ

Пример №12

Для системы найти и .

Решение:

Поскольку порядок многочлена числителя ПФ равен порядку многочлена знаменателя , начальное значение переходной функции равно

Коэффициент усиления в установившемся режиме равен

Пример №13

Определить передаточную функцию объекта регулирования, если его весовая функция равна

Решение:

По таблице соответствия АЛ находим изображение весовой функции (а это уже и есть передаточная функция объекта)

Приведя все дроби к общему знаменателю, получим ПФ в стандартном виде

Пример №14

Найти весовую функцию системы, если переходная функция равна

Решение:

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

Другой путь решения — через преобразование Лапласа

убираем нулевой корень s в знаменателе, принадлежащий входному воздействию — скачку, получаем ПФ или изображение весовой функции

откуда весовая функция

Свободное движение системы

К оглавлению…

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

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

Используют два способа вычисления совокупного переходного процесса. В первом случае система обычно задается ОДУ, производят в соответствии со свойством дифференцирования преобразования Лапласа индивидуальное преобразование каждого члена дифференциального уравнения, вычисляются одновременно вынужденная и свободная составляющие.

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

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

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

Пример №15

Для системы, заданной ОДУ , найти реакцию на начальные условия .

Решение:

Преобразуем индивидуально каждый член ОДУ по Лапласу с учетом свойств дифференцирования оригинала при ненулевых начальных условиях

Группируем и переносим подобные члены, подставляем значения

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

записываем разложение на простые дроби, вычисляем вычеты в полюсах (смотри приложение Б), переходим к оригиналу по таблице А.1

При начальное значение , как и было задано.

Пример №16

Система задана ОДУ

Найти реакцию системы, если

Решение:

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

Определяем по характеристическому полиному числитель и свободную составляющую переходного процесса

Полное описание переходного процесса

Частотные характеристики. Основные частотные характеристики

К оглавлению…

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

получают амплитудную (АЧХ) , фазовую (ФЧХ) , действительную (ВЧХ) и мнимую (МЧХ) частотные характеристики, связанные соотношениями

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

то амплитудная характеристика будет равна

а фазовая характеристика

Обобщающей является амплитудно-фазовая частотная характеристика (АФЧХ или просто АФХ) — кривая (годограф), которую чертит на комплексной плоскости конец вектора при изменении частоты от 0 до .

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

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

Быстрая проверка правильности расчетов:

  • АФЧХ и АЧХ начинаются при значении ;
  • АФЧХ и АЧХ заканчиваются в нуле или при (для );
  • АФЧХ устойчивой системы, не имеющей нулей, проходит по часовой стрелке столько квадрантов, каков порядок характеристического полинома.

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

Примеры решения задач по ТАУ

Пример №17

Построить частотные характеристики системы с ПФ .

Решение:

Подставляем учитывая, что , снижаем порядок ( и т.п.), избавляемся от мнимости в знаменателе, умножая числитель и знаменатель дроби на комплексное выражение, сопряженное стоявшему в знаменателе, отделяем действительную и мнимую части, приводим в знаменателе подобные члены

В данном случае числители и знаменатели дробей (действительной и мнимой частей) на сократить нельзя. Составляем таблицу (таблица 1), используя обязательные значения частот (можно взять больше точек, но не меньше), и подставляем эти значения:

  • крайние частоты 0 и ;
  • частоты пересечения характеристик с осями (определяются путем приравнивания числителей дробей мнимой и действительной части к нулю и решения полученного уравнения);
  • частоты разрыва характеристики (находят, приравнивая знаменатель нулю и решая уравнение) и близкие к ним (чуть больше-чуть меньше) частоты;
  • прочие частоты для повышения точности расчета.

Приравнивая , получаем , откуда = 2,45.

Приравнивая , получаем 10 = 0, откуда = 0.

По виду биквадратного уравнения определяем, что частот разрыва (действительных корней) нет. Частоты 1 и 3 рад/с добавлены произвольно для более точного построения графика.

По одной таблице можно построить АФЧХ на комплексной плоскости (рисунок 1.25, а), индивидуально ВЧХ и МЧХ (рисунок 1.25, б), и после дополнительных расчетов АЧХ и ФЧХ (рисунок 1.25, В).

Пример №18

Записать аналитически реакцию системы с известными АЧХ и ФЧХ (рисунок 1.26) на воздействие

Решение:

Общий вид гармонического сигнала . Следовательно, входное воздействие характеризуется параметрами: амплитуда 3,5, фаза 0 рад, частота = 1 рад/с. Находим для этой частоты по графику

Отсюда амплитуда выходной величины равна 3,5 0,36 = 1,26; фаза выходной величины 0 — 0.785 рад и окончательный вид реакции

Пример №19

При воздействии найти сигнал на выходе системы с передаточной функцией .

Решение:

Получаем по ПФ аналитические выражения для АЧХ и ФЧХ

Для известной частоты 10 рад/с значения АЧХ и ФЧХ равны

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

Логарифмические частотные характеристики

К оглавлению…

Зависимость от называется логарифмической амплитудной частотной характеристикой (ЛАЧХ) или ЛАХ. Зависимость от называется логарифмической фазной частотной характеристикой (ЛФЧХ) или просто ЛФХ. Частоту откладывают либо в логарифмах (в декадах), либо в радианах, но с учетом логарифмического масштаба. Декада соответствует изменению частоты в 10 раз, откладывают в децибелах (дБ), в градусах.

Для упрощения при построении вручную действительную ЛАЧХ заменяют асимптотической, т.е. ломаной линией из прямых отрезков, имеющих стандартный наклон, кратный ±20дБ/дек.

Частоты пересечения отрезков со а называются частотами сопряжения, они соответствуют корням ПФ. Частоты пересечения ЛАЧХ с осью абсцисс называются частотой среза, они соответствуют значению или (усиление или ослабление сигнала на частоте среза отсутствует). Для удобства построения через значения сопрягающих частот проводят на графике вертикальные линии, а на свободном поле графика — вспомогательные линии со стандартными наклонами (-20) дБ/дек.

Частоты сопряжения находят по корням (постоянным времени ) простых дробей, на которые разбивают ПФ, или типовых звеньев, из которых состоит структурная схема системы регулирования.

Звено первого порядка (один действительный корень)

Звено второго порядка (комплексные сопряженные корни)

где — показатель затухания (коэффициент демпфирования), характеризует величину резонанса в звене. При = 1 резонанс отсутствует, при резонансный выброс стремится к бесконечности. При значениях < 0,6 асимптотическую ЛАЧХ корректируют на величину выброса , определяемого по формуле где — число одинаковых корней (кратность корня), либо по типовым характеристикам (таблица 2) и графикам.

Левую (начальную) часть ЛАЧХ (низкочастотную или НЧ-асимптоту) или ее продолжение проводят через точку с координатами и слева направо с наклоном (-20 дБ/дек) до первой (наименьшей) частоты сопряжения. Здесь это степень астатизма, — число нулевых корней знаменателя, — числителя; добротность — отношение свободных членов полиномов числителя и знаменателя ПФ после удаления нулевых корней.

Двигаясь вправо, на каждой частоте сопряжения продолжают ЛАЧХ с отклонением от предыдущего направления: для корня числителя вверх (+20 дб/дек); для корня знаменателя вниз (-20 дБ/дек). Если кратность корня , наклон асимптоты изменяется в раз. Общий наклон ЛАЧХ в конце равен -(-20 дБ/дек). Выбросы при комплексных корнях откладывают вверх для корней знаменателя, вниз для корней числителя, близкие выбросы суммируются графически.

ЛФЧХ устойчивых систем строят по шаблону, неустойчивых -по вычисляемым точкам. Приближенно считают, что участку ЛАЧХ с наклоном ±20 дБ/дек соответствует фазовый сдвиг около ±90°, а участку с наклоном ±40 дБ/дек сдвиг на ±180°; действительному корню знаменателя соответствует угол наклона ЛФЧХ на сопрягающей частоте , комплексной паре

У статических систем (степень астатизма = 0) НЧ-асимптота представляет собой прямую, параллельную оси частот, и значение в децибелах равно расстоянию этой прямой от оси частот . У астатических систем находят частоту пересечения НЧ-асимптоты или её продолжения с осью частот, откуда . Степень астатизма определяется по наклону НЧ-асимптоты относительно оси частот, частоты сопряжения находят по точкам пересечения асимптот — касательных, проведенных к линейным участкам реальной ЛАЧХ.

Пример №20

Построить ЛАЧХ системы, заданной структурной схемой (рисунок 1.31, а). Передаточная функция системы равна

Решение:

Определяем параметры НЧ-асимптоты:

  • порядок астатизма = 1 — 0 = 1 (имеется один нулевой корень в знаменателе);
  • добротность .

Нули в системе отсутствуют, полюс -5 имеется, отсюда частота сопряжения Строим график ЛАЧХ толстой сплошной линией, проводя слева вниз прямую линию с наклоном 1 х(-20 дБ/дек) через точку с координатами (20 дБ, 0) до первой частоты сопряжения (рисунок 1.31, б). Поскольку частота сопряжения соответствует полюсу, отклоняемся от текущего направления вниз на угол -20 дБ/дек, общий наклон ЛАЧХ в конце равен -40 дБ/дек. Корень действительный, поэтому резонанса нет, выбросы не учитываем.

Пример №21

Составить ПФ системы с заданной ЛАЧХ (рисунок 1.31, в), предполагая, что все корни имеют отрицательную действительную часть.

Решение:

На частотах сопряжения наблюдается отклонение характеристики от предыдущего направления вверх на +20 дБ/дек, на частотах сопряжения — вниз на -20 дБ/дек, поэтому передаточная функция будет иметь вид

Поскольку

и окончательно

Устойчивость непрерывных стационарных систем. Математический и физический признаки устойчивости

К оглавлению…

Устойчивость — это свойство системы возвращаться в исходное состояние равновесия после снятия воздействия, выведшего систему из этого состояния.

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

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

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

Пример №22

Оценить прямым методом устойчивость системы, описываемой дифференциальным уравнением

Решение:

Характеристическое уравнение системы

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

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

Пример №23

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

Решение:

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

Алгебраические критерии устойчивости. Критический коэффициент усиления

К оглавлению…

Критерий Гурвица: система устойчива, если все коэффициенты ее характеристического уравнения

и все диагональные миноры матрицы Гурвица положительны.

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

для системы четвертого порядка

Критерий Гурвица удобно использовать при устном счете для систем не выше четвертого порядка.

Критерий Рауса: система устойчива, если все коэффициенты ее характеристического уравнения и все элементы первого столбца таблицы Рауса положительны. Необходимое условие (положительность всех коэффициентов) совпадает с критерием Гурвица.

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

где — номер строки, — номер столбца, — вспомогательное число для /-той строки. Таблица содержит строку и ()/2 с округлением столбец.

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

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

Контрольная работа по теории автоматического управления ТАУ

Пример №24

Оценить по критерию Гурвица устойчивость системы

Решение:

Характеристическое уравнение

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

Проверяем достаточное условие по определителю Гурвица

Оба диагональных минора положительны. Так как необходимое и достаточное условия выполняются, система устойчива.

Пример №25

Оценить по Раусу устойчивость системы с характеристическим уравнением

Решение:

Необходимое условие выполняется.

Проверяем достаточное условие -составляем таблицу Рауса: число строк равно числу коэффициентов (шесть), число столбцов 6/2 = 3. Заполняем две первые строки попарно коэффициентами с четными и нечетными индексами. Последний коэффициент смещается вниз и влево ходом шахматного коня (три клетки вниз и одна влево), ниже него записываем нули. Вычисляем вспомогательное число и элементы третьей строки: ; откуда , затем элементы остальных строк.

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

Пример №26

Найти критическое значение коэффициента усиления системы с характеристическим уравнением

Решение:

Формулируем условия нахождения системы на границе устойчивости по критерию Гурвица (он наиболее удобен и нагляден для систем первого-третьего порядка):

  • на апериодической границе , откуда ;
  • на периодической границе , откуда следует . Учитывая опущенные знаки неравенств, делаем вывод, что система устойчива при значениях коэффициента усиления .

Оценить устойчивость по критерию Рауса системы с характеристическим уравнением

Частотные критерии устойчивости. Критерий Михайлова Согласно принципу аргумента, известному в теории комплексной переменной, если среди полюсов ПФ системы расположены справа от мнимой оси, а остальные — слева, то полное изменение аргумента комплексной функции равно

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

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

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

Система находится на апериодической границе устойчивости, если кривая при = 0 начинается в начале координат, и на периодической границе устойчивости, если кривая при проходит через начало координат. Частота незатухающих колебаний соответствует периодической границе устойчивости системы.

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

Действительная часть

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

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

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

Пример №27

Система пятого порядка с кривой Михайлова (рисунок 1.43) неустойчива, т.к. сначала вектор повернулся против часовой стрелки на три квадранта (три левых полюса), а затем по часовой стрелке на два квадранта (два правых полюса).

Решение:

Иначе: итоговый поворот равен одному квадранту, т.е. , тогда правых корней характеристического полинома (5-1 )/2 = 2.

Пример №28

Найти критическое значение коэффициента усиления системы с

по критерию Михайлова.

Решение:

Заменяя получим характеристическую функцию

Условия нахождения САУ на границе устойчивости

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

Тогда из второго уравнения определяем частоту

и подставляем ее значение в первое уравнение

Частота, соответствующая колебательной границе устойчивости

Курсовая работа по теории автоматического управления ТАУ

Пример №29

Используя вторую форму (следствие) критерия Михайлова, оценить устойчивость системы

Решение:

В характеристическом уравнении

заменяем , снижаем порядок и группируем

Здесь — это четная (действительная) функция , а — это нечетная (мнимая) функция .

Приравнивая поочередно четную и нечетную функции нулю, находим частоты 1,41 и 1,73, соответствующие пересечению кривой с осями координат, подставляем эти частоты в характеристическую функцию и заполняем таблицу. Строим графики четной и нечетной функций — они поочередно пересекают ось частот, т.е. их корни перемежаются, и общее число пересечений равно , следовательно, система устойчива (рисунок 1.44).

-разбиение по одному параметру

Областью устойчивости (0) называют область в пространстве изменяемых параметров, каждой точке которой соответствуют только левые корни характеристического уравнения. Остальные -области отличаются числом правых корней характеристического уравнения и обозначаются соответственно (1) — область с одним правым полюсом, (2) — с двумя и т.д.

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

Направление штриховки указывает на область с наибольшим числом левых корней. При каждом переходе через кривую навстречу штриховке один корень характеристического уравнения становится правым, в обратном направлении — левым. Выбранную по штриховке область-претендент (0) проверяют на устойчивость с помощью любого критерия, подставив значение параметра из этой области в характеристическое уравнение. Поскольку изменяемый параметр является действительной величиной, его допустимые значения лежат на отрезке действительной оси, заключенном внутри области устойчивости ( 0).

Пример №30

Найти методом -разбиения критические значения коэффициента усиления к системы, заданной передаточной функцией

Решение:

Разрешаем характеристическое уравнение системы

относительно исследуемого параметра

производим замену

снижаем порядок и группируем

Определяем частоты пересечения основной кривой с осями:

Строим основную и зеркальную кривые на комплексной плоскости, указывая направление возрастания частоты стрелкой на характеристике (рисунок 1.47). Наносим штриховку, обозначаем области с предполагаемым числом правых полюсов в скобках. Проверяем область-претендент (0) на устойчивость по критерию Гурвица при значении , выбранном на отрезке внутри области (0) между точками 0 и 60

Так как и необходимое, и достаточное условия устойчивости по Гурвицу при выполняются, то система будет устойчивой при любых значениях коэффициента усиления в интервале 0 < < 60. Критические значения коэффициента равны

Критерий Найквиста. Запасы устойчивости Упрощенная формулировка: система, устойчивая в разомкнутом состоянии или нейтральная, будет устойчивой в замкнутом состоянии, если АФЧХ разомкнутой системы при изменении частоты от нуля до плюс бесконечности не охватывает точку с координатами . Всегда подразумевается замыкание системы единичной ООС.

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

Оценка запасов устойчивости по АФЧХ. Запасы устойчивости по амплитуде в относительных единицах равны расстоянию от критической точки до ближайших точек пересечения АФЧХ с отрицательной действительной полуосью. В децибелах запас устойчивости по амплитуде находят как величину, обратную амплитуде вектора при угле -180° или , где — расстояние от точки пересечения АФЧХ с отрицательной действительной полуосью до начала координат. Норма или 6-12 дБ.

Запас устойчивости по фазе равен углу между отрицательной действительной полуосью и лучом, проведенным из начала координат в точку пересечения АФЧХ с дугой единичного радиуса. Запас по фазе находится в пределах от 0 до 180°, при проектировании обычно нормой является .

Система устойчива в замкнутом состоянии, если обратная АФЧХ разомкнутой системы охватывает точку .

Логарифмический критерий Найквиста (диаграмма Боде). Обычная формулировка: замкнутая система устойчива, если в момент пересечения ЛФЧХ разомкнутой системы линии -180° её ЛАЧХ отрицательна. Общая формулировка пригодна и для систем, неустойчивых в разомкнутом состоянии: замкнутая система устойчива, если на интервале положительности ЛАЧХ разомкнутой системы сумма переходов ее ЛФЧХ линии -180° равна , где — число правых корней характеристического уравнения разомкнутой системы.

Оценка запасов устойчивости по ЛЧХ. Запас устойчивости по амплитуде равен отклонению ЛАЧХ от нуля на ближайших к частоте среза частотах пересечения ЛФЧХ с линией минус 180°. Запас устойчивости по фазе равен отклонению ЛФЧХ на частоте среза от линии минус 180° к нулю.

Пример №31

Оценить устойчивость системы (рисунок 1.50) по Найквисту.

Решение:

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

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

Определяем частоты пересечения годографа с осями координат: мнимая часть отсутствует, из уравнения видно, что корни, т.е. частоты пересечения с мнимой осью, отсутствуют. Зато уравнение дает частоту разрыва характеристики . В подобном случае обычно берут еще две частоты (произвольно) — немного меньше частоты разрыва и немного больше, например, возьмём 0,1 и 10.

Замкнутая система также находится на колебательной границе устойчивости (рисунок 1.51), т.к. АФЧХ проходит через точку .

Пример №32

Оценить запасы устойчивости по АФЧХ после замыкания единичной ООС системы с

Решение:

Задача не требует построения АФЧХ. По критерию Гурвица следует, что в разомкнутом состоянии система устойчива, нулей нет, поэтому годограф Найквиста проходит два квадранта по часовой стрелке и не пересекает отрицательную действительную полуось. Таким образом, запас по амплитуде максимален . Полюса системы действительные -1 и -2, следовательно, резонанс в системе отсутствует и амплитуда вектора нигде не превышает величины = 1/2, запас устойчивости по фазе равен = 180°.

Качество непрерывных стационарных систем. Прямые оценки качества регулирования

К оглавлению…

Прямые оценки качества определяются по переходной характеристике, т.е. реакции системы на единичный скачок при нулевых начальных условиях (рисунок 1.55).

Время регулирования измеряется от начала переходного процесса до момента, после которого характеристика не отклоняется от установившегося значения более, чем на величину допустимой ошибки (обычно 5 %, реже 2 % от установившегося значения). Следует указывать, при какой зоне получено время регулирования.

Перерегулирование — величина максимального относительного заброса переходной характеристики от начальной величины за линию установившегося значения (в относительных единицах или %)

Если начальное и конечное значения характеристики равны нулю или одинаковы (и приняты условно за 0), возможны два способа оценки. При наличии разнополярных значений перерегулирование равно отношению величины второго экстремума к величине первого (рисунок 1.56, а), а если колебание одно (рисунок 1.56, б), то перерегулирование равно отношению величины максимального отклонения к величине входного воздействия (обычно это единица). Зону для оценки времени регулирования в первом случае определяют от значения первого максимума, во втором случае — от величины входного воздействия.

Время нарастания определяется: для процессов с перерегулированием как время от начала процесса до момента пересечения кривой линии установившегося значения; для любых процессов как время между моментами достижения заданных уровней установившегося значения (например, 10 и 90 %). Поэтому при оценке времени нарастания следует указывать, каким способом оно получено.

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

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

Степень затухания (демпфирования) — величина относительного уменьшения

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

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

Задачи теории автоматического управления ТАУ

Пример №33

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

Решение:

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

Пример №34

Определить величину перерегулирования и времени регулирования (рисунок 1.57)

Решение:

Перерегулирование = (1,5 — 1,0)/1,0 = 0,5 или 50 %. Для определения времени регулирования проводим параллельно линии установившегося значения две прямые на уровне . По точке последнего вхождения кривой в зону 2 получаем = 15 с.

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

При заданной зоне ошибки 2 % вместо коэффициента 3 берут приблизительно 4.

Найдя степень колебательности системы

определяют значение перерегулирования

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

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

Пример №35

Оценить показатели качества регулирования системы, имеющей нуль -0,125, полюса и коэффициент передачи 1,2.

Решение:

Коэффициент передачи на относительные показатели не влияет. Нуль -0,125, равный полюсу, взаимно с ним компенсируется. Следовательно, доминирующими являются комплексно-сопряженные полюса откуда , степень колебательности системы

и перерегулирование

или 45,6 %.

Пример №36

Оценить перерегулирование и время регулирования системы

с законом управления .

Решение:

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

Нули отсутствуют, из характеристического уравнения

находим полюса

Отсюда

а перерегулирование или 4,3 %.

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

Общие принципы оценки качества по вещественной частотной характеристике :

  • — конечное значение переходной характеристики численно равно начальному значению ВЧХ;
  • — начальное значение переходной характеристики численно равно конечному значению ВЧХ;
  • — кратность изменения масштаба ВЧХ и переходной характеристики одинакова;
  • — расширение полосы рабочих частот ведет к соразмерному повышению быстродействия системы;
  • время регулирования
  • перерегулирование а определяется по форме ВЧХ:

а) если ВЧХ монотонно убывает, то перерегулирование = 0;

б) если ВЧХ является положительной невозрастающей функцией, то перерегулирование < 18 %;

в) если ВЧХ имеет подъем от

г) если ВЧХ имеет отрицательный минимум со значением более 0,1 , то с его учетом

д) если ВЧХ терпит разрыв при , система совершает незатухающие колебания, и показатели качества не определяются.

При оценке качества регулирования по АЧХ обычно вычисляют значение частотного показателя колебательности, равное отношению максимума характеристики к ее начальному значению . При = 1 переходная характеристика системы не колебательна, при система находится на границе устойчивости, наблюдаются незатухающие колебания с частотой . Оптимальными считаются значения = 1,1..1,5, которым соответствует перерегулирование 10-30 % и запас по фазе 30-50°.

Пример №37

Оценить значение частотного показателя колебательности системы по её АЧХ (рисунок 1.68).

Решение:

Максимальное значение АЧХ равно 1,51, следовательно, показатель колебательности = 1,51/1,0 = 1,51, что ещё удовлетворяет минимальному запасу по фазе 30° и перерегулированию 30 %.

Пример №38

Найти значение перерегулирования и времени регулирования системы по заданной АФЧХ (рисунок 1.69)

Решение:

Частота = 1,45 рад/с, положительный максимум ВЧХ равен 1,09 при начальном значении = 1,0, отрицательный минимум 0,521. Отсюда получаем перерегулирование

и время регулирования не более

Интегральные оценки качества переходных процессов Интегральные показатели качества регулирования дают совокупную оценку быстродействия и колебательности без вычисления их значений. Они характеризуют отклонение реального переходного процесса от заданного идеального.

Интегральная линейная оценка (ИЛО) определяется площадью отклонения реального процесса от идеального ступенчатого. Для обеспечения требуемых динамических свойств САУ необходимо выразить величину через коэффициенты передаточной функции системы

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

Пример №39

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

линейная интегральная оценка

зависит от соотношения постоянных времени и . Минимум оценки достигается при их равенстве.

Точность в установившемся режиме

К оглавлению…

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

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

Передаточная функция ошибки воспроизведения задания определяется по ПФ разомкнутой системы как по передаточной функции замкнутой системы как .

Относительная величина установившейся ошибки называется коэффициентом статизма (статизмом) системы по соответствующему каналу:

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

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

Установившийся динамический режим имеет место при возмущенном движении системы с момента затухания свободной составляющей переходного процесса.

Если входное воздействие аппроксимируется полиномом от т.е. разлагается в степенной ряд

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

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

Обычно вычисляют не более трех первых коэффициентов ошибок. Коэффициенты передачи составляющих входного воздействия вычисляются по ПФ разомкнутой системы и называются:

Пример №40

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

Пример №41

Определить полную статическую ошибку для системы (рисунок 1.72), полагая, что

Выражение для суммарной ошибки в операторной форме

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

Пример №42

Определить три первых коэффициента ошибки, вынужденную составляющую ошибки от воздействия

и добротность по скорости для системы, имеющей в разомкнутом состоянии ПФ

Решение:

Находим передаточную функцию по каналу ошибки

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

В общем виде вынужденная составляющая ошибки воспроизведения задающего воздействия равна

Для задающего воздействия

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

Добротность по скорости вычисляем по ПФ разомкнутой системы

Многомерные системы регулирования. Переход к пространству состояний

К оглавлению…

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

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

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

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

По системе уравнений составляется матрица состояния (из коэффициентов при ) и матрица входа (из коэффициентов при входном воздействии ), по уравнению выхода составляется матрица выхода (из коэффициентов при )

Сопровождающая матрица (матрица Фробениуса) может быть записана прямо по ОДУ (по характеристическому полиному системы)

По уравнениям состояния или матрицам указанного вида легко восстановить ПФ или ОДУ, учитывая, что в последней строке сопровождающей матрицы записаны с конца, с обратным знаком коэффициенты нормированного характеристического многочлена, а в матрице — коэффициенты многочлена числителя передаточной функции в обратном порядке.

Пример №43

Дифференциальное уравнение объекта управления

Выбираем переменные состояния

Решение:

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

Пример №44

Пусть модель объекта управления имеет вид

тогда после нормирования (деления на 2), считая общий коэффициент перед правой частью уравнения равным единице, получим описание системы в пространстве состояний матрицами

Канонические представления

К оглавлению…

Стандартные формы описания систем в пространстве состояний с сопровождающей матрицей называются каноническими. Это каноническая управляемая форма (с упрощенной матрицей ) и каноническая наблюдаемая форма (с упрощенной матрицей ).

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

Пусть

тогда вычисления для перехода к канонической управляемой форме имеют вид

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

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

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

Пусть

тогда

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

Если матрица не сопровождающая, а произвольного вида, ее характеристический многочлен нужно вычислять как определитель , где — комплексная переменная Лапласа, 1 — единичная матрица.

Корни характеристического уравнения =0 являются собственными значениями матрицы . Матрицы подобны, если имеют одинаковые собственные значения (характеристические многочлены и их корни).

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

Пример №45

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

Каноническое управляемое представление (нормирование по не требуется, матрица имеет стандартный вид, всегда одинаковый)

Пример №46

По уравнению

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

матрица нулевая, поскольку , и окончательно (матрица с имеет стандартный вид, всегда одинаковый)

Пример №47

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

Решение:

Коэффициенты на главной диагонали матрицы равны её собственным значениям (полюсам системы) ; структурная схема соответствует рисунку 2.1. Матрицы и включены последовательно, поэтому, если вычеты 0,5 и 0,5 вписаны в матрицу (как показано), то в матрицу записываются единицы, и наоборот.

Пример №48

Оценить устойчивость системы, проверить подобие матрицы и матрицы

Система

матрица

Решение:

Характеристическая матрица

Характеристический многочлен (определитель характеристической матрицы)

По критерию Гурвица система устойчива, т.к. все коэффициенты характеристического многочлена положительны.

Характеристический многочлен матрицы

равен

Матрицы и подобны, поскольку равны их характеристические многочлены.

Описание по структурной схеме

К оглавлению…

На структурной схеме переменные состояния могут быть назначены разным образом, поэтому и описания системы в пространстве состояний будут отличаться. Все матрицы имеют нестандартный вид. Однако переменная всегда назначается на выходе блока с в знаменателе, а ОДУ первого порядка для каждого такого блока записывают в зависимости от вида знаменателя:

а) звено с нулевым корнем в знаменателе (рисунок 2.3, а)

б) звено с действительным корнем, две формы (рисунок 2.3, б)

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

Звено с комплексными сопряженными корнями (рисунок 2.3, в), не разлагается на два простых, поэтому вводят условно переменную состояния с промежуточным индексом и составляют два уравнения

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

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

Умножая матрицу на вектор и вектор на вход , получаем систему уравнений, которую затем совмещаем с уравнениями оставшейся части структурной схемы.

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

Пример №49

Описать систему (рисунок 2.4, а)

Решение:

Сначала рассматриваем сложный блок с переменной s в числителе, учитывая, что вектор с для него составлен единственной единицей и в вычислениях не нуждается, а переменная состояния на выходе блока имеет индекс 2:

Затем описываем всю систему, включая в нее этот блок:

и окончательно

Пример №50

Составляя уравнения состояния для случая, когда в цепи обратной связи есть звено дифференцирования с (рисунок 2.4, б) учитываем, что умножение на s в операторной области соответствует взятию производной во временной области.

Решение:

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

Синтез структурной схемы

К оглавлению…

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

Пример №51

Перейдем от матриц

к структурной схеме (рисунок 2.11), для чего выбираем число звеньев (равно порядку матрицы ), определяем корни знаменателей ПФ по диагональным элементам матрицы ( = -1 у блока с переменной на выходе и = -3 у блока с переменной ), находим коэффициенты прямых связей — числители ПФ блоков между и , между и (оба числителя равны 1). В схеме имеются две отрицательные обратные связи: единичная ООС от к и с коэффициентом 3 от к . На входе системы находится блок с коэффициентом 2, выход связан с системой через коэффициенты 1 матрицы .

Пример №52

Построить структурную схему объекта, заданного системой дифференциальных уравнений

Решение:

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

Пример №53

Построить структурную схему объекта по дифференциальному уравнению

Решение:

Поскольку порядок системы равен трем, используем три интегратора , включив их последовательно и установив сумматор на входе первого интегратора слева. К инвертирующему входу этого сумматора подключаем через согласующие сумматоры блоки с коэффициентами (по порядку): — с выхода первого интегратора, — с выхода второго интегратора, — с выхода третьего интегратора.

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

Основные матричные функции

К оглавлению…

— характеристическая матрица, аналог характеристического полинома одномерной системы .

— системная матрица (резольвента), называемая также передаточной матрицей или матрицей передаточных функций (МПФ) для переменных состояния, аналог системной функции .

— реальная МПФ для назначенных входов и выходов (передаточная матрица выходов), совпадает по виду с только в частном случае.

Пример №54

Система задана в пространстве состояний матрицами

Характеристическая матрица

Характеристический полином (определитель характеристической матрицы)

Присоединенная матрица

Решение:

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

матрица передаточных функций выходов

Решение уравнения движения

К оглавлению…

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

Реакция на выходе системы вычисляется с учетом матрицы

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

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

Пример №55

Найти при и начальных условиях уравнения движения системы

Решение:

Система задана в наблюдаемой форме с матрицей = [1 0], поэтому вектор начальных значений переменных формируем по выходу

Характеристическая матрица

Характеристический полином (определитель характеристической матрицы)

Резольвента

где присоединенная матрица

Заменяем по таблице соответствия изображения на оригиналы

Пример №56

Найти изображение реакции на системы

Решение:

Изображение входного воздействия

Вычисление фундаментальной матрицы

К оглавлению…

Поскольку

то фундаментальную матрицу определяют как матричную экспоненту от тремя способами:

а) разложением в бесконечный или конечный ряд

где — порядок системы.

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

б) по формуле Сильвестра

где — собственные значения матрицы (корни характеристического уравнения системы), или в развернутом виде

Здесь

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

в) Наконец, вычисляется и как обратное преобразование Лапласа от системной матрицы или .

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

Пример №57

Определим матричную экспоненту для системы с . Поскольку уже при получена нулевая матрица

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

Пример №58

Определить методом Сильвестра для системы

Решение:

Вычисляем характеристический полином, находим его корни

Вычисляем матрицы коэффициентов при собственных модах системы

Пример №59

Определить с помощью обратного преобразования Лапласа фундаментальную матрицу системы

Решение:

Находим корни характеристического полинома и адъюнкту -1

Общий вид разложения на простые дроби

Находим коэффициенты числителей простых дробей:

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

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

Управляемость и наблюдаемость систем

К оглавлению…

Для управляемости системы необходимо и достаточно, чтобы матрица управляемости вида

имела ранг, равный . При управляемости системы говорят также, что пара управляема.

Ранг матрицы равен порядку её наибольшего ненулевого минора. Матрица составляется присоединением справа к матрице произведения матриц , затем произведения и т.д. Размерность матрицы равна , где — число входов. Если ранг матрицы (обозначим его ) не равен единице, то вычисление матрицы можно закончить досрочно по формуле

Система полностью управляема при , полностью неуправляема при = 0, частично управляема при , порядок управляемости равен .

Для наблюдаемости системы необходимо и достаточно, чтобы матрица наблюдаемости

имела ранг, равный порядку системы . Символ означает транспонирование или перевод вектора-строки в вектор-столбец. Говорят иначе, что пара наблюдаема.

Система полностью наблюдаема при , полностью не-наблюдаема при = 0, частично наблюдаема при , порядок наблюдаемости равен .

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

Существует и иная форма составления матрицы наблюдаемости — по вертикали без транспонирования

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

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

Пример №60

Оценить управляемость системы (достаточно иметь пару и ).

Система

Решение:

Находим

Определитель матрицы управляемости

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

Задачу можно было не решать: числитель ПФ содержит только 1 (это видно из матрицы ), следовательно, сокращаемые пары нуль-полюс отсутствуют и система полностью управляема.

Система

Пример №61

Оценить управляемость системы.

Матрица диагональная (в каждой строке одна переменная с возрастающим индексом). Уже ясно, что система неуправляема по (по полюсу +1), поскольку в первом уравнении нет . Проверим вывод.

то

Система частично управляема, порядок управляемости равен двум.

Пример №62

Оценить наблюдаемость системы

С учетом того, что делаем вывод, что = 1 — система частично наблюдаема, порядок наблюдаемости равен 1.

Пример №63

Проверить управляемость системы

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

содержит сокращаемую пару (диполь) нуль -1/полюс -1, что ведет либо к неуправляемости, либо к ненаблюдаемости системы. От чего это будет зависеть? Составим описание системы в канонической управляемой форме и проверим управляемость

Система в таком представлении полностью управляема (но не вполне наблюдаема). Составим описание системы в канонической наблюдаемой форме и снова проверим управляемость

А теперь система управляема частично. Таким образом, если в ПФ системы обнаруживается сокращаемая пара, неуправляемость или ненаблюдаемость зависит от того, какое представление выбирается для перехода в пространство состояний. Если же в ПФ сокращаемые пары отсутствуют, система полностью управляема и наблюдаема.

Наблюдатели состояния

К оглавлению…

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

Наблюдатель в виде параллельного фильтра представляет собой модель объекта регулирования на интеграторах в каноническом управляемом представлении. Его вход подключается параллельно входу объекта регулирования, а с выходов интеграторов снимают идеальные значения переменных состояния объекта (оценки), которые обозначают значком «каре» А над символом переменной. Разница значений выходов объекта и наблюдателя называется невязкой (обозначается значком «тильда» ~ над символом сигнала), при совпадении модели с оригиналом невязка стремится к нулю.

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

Пример №64

Построить наблюдатель в виде параллельного фильтра к объекту с передаточной функцией

Решение:

Модель объекта (описание наблюдателя) соответствует канонической форме управляемости

Этому описанию отвечает структурная схема (рисунок 2.22)

Пример №65

Построим наблюдатель в виде фильтра Калмана для объекта, заданного системой дифференциальных уравнений

обеспечив показатели качества переходного процесса ошибки наблюдателя

Решение:

По матрицам коэффициентов объекта регулирования определяем его передаточную функцию (объект неустойчив)

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

компенсирующая добавка образуется обратными связями с коэффициентами (рисунок 2.23).

В соответствии с матрицей

характеристический полином наблюдателя имеет вид

или

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

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

Из равенства

находим неизвестные коэффициенты корректирующих обратных связей

Пример №66

Рассчитать параметры наблюдателя в виде фильтра Калмана (рисунок 2.24) с компенсирующим звеном, имеющим передаточную функцию

при тех же требованиях к качеству переходного процесса наблюдателя и параметрах ПФ модели объекта регулирования .

Решение:

Передаточная функция модели объекта регулирования равна

а характеристическое уравнение наблюдателя имеет вид

откуда, приравняв числитель нулю и нормируя, получаем

Желаемый характеристический полином третьего порядка формируем из корней с одинаковой действительной частью и -3, он равен

Приравнивая находим неизвестные коэффициенты

Проектирование модального регулятора

К оглавлению…

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

Закон управления для объекта второго порядка имеет вид

где — коэффициент усиления П-регулятора, — задание, — коэффициенты обратных связей регулятора по переменным состояния.

Подставив значение в уравнение состояния, получим систему уравнений, которая описывает замкнутую систему управления

и характеристический полином замкнутой системы

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

Приравняв коэффициенты полиномов при одинаковых степенях , получим

Расчет существенно упрощается, если объект представлен в канонической форме управляемости с .

Коэффициент усиления обычно находят из условия нулевой статической ошибки: либо по коэффициентам передаточной функции

откуда

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

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

Пример №67

Рассчитать параметры модального регулятора для объекта

при требованиях к качеству регулирования

Решение:

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

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

здесь — неизвестные коэффициенты, — коэффициент передачи объекта регулирования.

Тогда характеристическое уравнение замкнутой системы равно

Выберем распределение корней, обеспечивающее заданное качество процессов, например,

(все действительные полюса обеспечат нулевое перерегулирование и время регулирования не более 3/2 = 1,5 с). Сформируем желаемое характеристическое уравнение третьего порядка

Приравнивая коэффициенты при соответствующих степенях , получим расчетные соотношения

Отсюда находим параметры регулятора

Пример №68

ПФ объекта регулирования после нормирования имеет вид

заданные показатели качества: время регулирования 6 с, перерегулирование 0,02, выбрать параметры модального регулятора.

Решение:

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

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

При заданном времени регулирования с степень устойчивости для ошибки равна = 3/6 = 0,5, отсюда действительная часть комплексного корня будет равна -0,5. Действительный корень принимаем в 10 раз большим, т.е. -5, чтобы исключить его влияние на переходный процесс. По заданной величине перерегулирования = 0,02 вычисляем степень колебательности , после чего можно вычислить мнимую часть комплексного корня

По значениям корней -5 и находим вид желаемого характеристического полинома

Из условия нулевой ошибки регулирования значение коэффициента усиления регулятора

Значения коэффициентов обратной связи по переменным состояния равны

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

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

Расчет подтверждает, что установившаяся ошибка отсутствует, так как коэффициент передачи в установившемся режиме равен 2,05/2,05 = 1, а полученный характеристический полином системы регулирования равен желаемому. При единственной обратной связи

Кстати дополнительная теория из учебников тут.

Преобразования подобия

К оглавлению…

При анализе и синтезе многомерных систем необходимо уметь переходить от одной формы к другой — поскольку все эти системы подобные, такой переход называется преобразованием подобия или базиса.

Один из путей перехода, приемлемый для одномерной системы

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

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

откуда следует, что матрицы коэффициентов новой системы равны

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

  • исходная и преобразованная система должны иметь одинаковые собственные значения (характеристические многочлены и их корни);
  • преобразование базиса не меняет передаточную функцию системы.

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

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

Для перехода к управляемой форме должна быть полностью наблюдаема пара для перехода к наблюдаемой форме должна быть полностью наблюдаема пара .

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

где — матрица обратных связей замкнутой системы по переменным состояния.

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

Пример №69

Пусть преобразуемый к канонической управляемой форме объект третьего порядка описывается системой уравнений

Решение:

Характеристический полином объекта равен

матрица управляемости

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

Применяя формулы

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