
- •«Теория автоматического управления»
- •Оценка точности в установившемся режиме
- •Проверка устойчивости исходной системы
- •Расчет корректирующего устройства
- •Построение области устойчивости скорректированной системы
- •5. Построение графика переходного процесса и оценка качества
- •5.1 Моделирование системы на авм
- •Моделирование системы на цвм
- •6. Вычисление и минимизация квадратичной интегральной оценки
- •Вычисление и минимизация дисперсии сигнала ошибки при случайных воздействиях
- •8. Анализ нелинейной системы
5. Построение графика переходного процесса и оценка качества
скорректированной системы
5.1 Моделирование системы на авм
Исходным выражением для моделирования является передаточная функция разомкнутой скорректированной системы Wск(р).
Составим аналоговую модель системы на основании моделей типовых динамических звеньев. Схема модели приведена на рисунке 7.
Рассчитаем машинные коэффициенты модели.
Машинные коэффициенты преобразующего элемента:
Машинные коэффициенты корректирующего устройства:
Машинные коэффициенты управляющего устройства:
Машинные коэффициенты исполнительного устройства:
Машинные коэффициенты объекта управления:
Рисунок
7 - Схема аналоговой модели системы
Моделирование системы на цвм
Для цифрового моделирования системы используем среду Matlab приложение Simulink. Схема цифровой модели системы, составленная из типовых звеньев библиотеки Simulink, показана на рисунке 8.
Рисунок 8 - Схема цифровой модели системы
Результатом моделирования является переходная характеристика по каналу ошибки xз – ε , регистрируемая прибором Scope1. Переходная характеристика системы по каналу ошибки xз – ε приведена на рисунке 9.
Рисунок 9 - Переходная характеристика системы по каналу ошибки xз – ε
По графику переходного процесса определим перерегулирование, время переходного процесса и сравним полученные значения с заданными.
ε(∞)=0,5 - что удовлетворяет заданию
6. Вычисление и минимизация квадратичной интегральной оценки
при типовом воздействии
Исходным выражением для вычисления квадратичной интегральной оценки является передаточная функция замкнутой скорректированной системы по каналу хз-ε при единичном ступенчатом воздействии (то есть принимаем задающее воздействие Хз(t)=1, а следовательно Хз(р)=1/р). Запишем эту передаточную функцию:
Раскроем скобки и приведем передаточную функцию к стандартному виду:
Подставим в формулу численные значения:
Запишем выражение для изображения переходной составляющей сигнала ошибки:
Так как передаточный коэффициент разомкнутого контура kрк≥10 допускается упростить выражение для изображения переходной составляющей с учетом следующих условий:
kрк+1=kрк и kрк-1=kрк
Преобразуем выражение для изображения переходной составляющей сигнала ошибки, используя вышеприведенные условия:
=
Для вычисления квадратичной оценки по изображению используют равенство Парсеваля, которое имеет вид:
Где
Таким образом коэффициенты di и vi равны:
do =0,0138 vo =0,00019
d1 =0,306 v1 =-0,143
d2 =1,786 v2 =3,19
d3 = kрк
Составим определитель Δ по правилу составления определителя Гурвица
Определитель v составляем из определителя путем замены коэффициентов верхней строки на коэффициенты 0, 1 и 2.
Вычисляем квадратичную интегральную оценку:
Задаваясь численными значениями kрк, составляем таблицу зависимости квадратичной интегральной оценки от коэффициента kрк, которая приведена ниже.
Таблица 3. – Расчетные данные для построения кривой зависимости интегральной оценки от передаточного коэффициента разомкнутого контура
kрк |
Qкв |
0 |
∞ |
1 |
1,039 |
2 |
0,63 |
5 |
0,37 |
8 |
0,325 |
10 |
0,317 |
10,5 |
0,366 |
10,56 |
0,3165 |
11 |
0,31,68 |
12 |
0,3174 |
20 |
0,39 |
24 |
0,47 |
Рисунок 9 - Кривая зависимости интегральной оценки от передаточного коэффициента разомкнутого контура
Из рисунка 9 видно, что оптимальным значением kрк является kрк=10,56. При данном значении коэффициента разомкнутого контура система будет работать в оптимальном режиме, обеспечивая минимальную площадь под графиком переходного процесса.
Теперь перейдем от коэффициента kрк к передаточному коэффициенту kу. Для этого воспользуемся следующей формулой:
При коэффициенте kрк=12,96 и квадратичной интегральной оценке равной Q=0.3218 передаточный коэффициент управляющего устройства kу.= 56,28.
Вывод: В этом разделе с помощью квадратичной интегральной оценки получили оптимальное значение передаточного коэффициента управляющего устройства (kу=45,8). Этот коэффициент получился меньше, чем тот, что был выбран в разделе 2 (kу=56,28). Следовательно, при выборе этого коэффициента точность системы в установившемся режиме увеличится, но могут получиться более колебательные переходные процессы.