Синтез астатических систем
Составление структурной схемы замкнутой астатической системы автоматического управления. Определение минимальной установившейся ошибки, построение области устойчивости и моделирование в программе MatLab. Компенсация действия неконтролируемых возмущений.
Рубрика | Программирование, компьютеры и кибернетика |
Вид | курсовая работа |
Язык | русский |
Дата добавления | 15.02.2011 |
Размер файла | 523,2 K |
Отправить свою хорошую работу в базу знаний просто. Используйте форму, расположенную ниже
Студенты, аспиранты, молодые ученые, использующие базу знаний в своей учебе и работе, будут вам очень благодарны.
Размещено на http://www.allbest.ru/
Федеральное агентство по образованию РФ
ГОУ ВПО
"Московский Государственный Текстильный Университет им. А.Н. Косыгина"
ФИТАЭ
Кафедра автоматики и промышленной электроники
Курсовая работа
Тема: "Синтез астатических систем"
Дисциплина: "Теория Автоматического Управления"
Москва 2009
Исходные данные:
ПИ-закон
1. По заданной системе уравнений получить передаточную функцию объекта управления и составить структурную схему замкнутой САУ, считая регулятор звеном с входной величиной Е, выходной U и передаточной функцией
2. Выбрать регулятор таким образом, чтобы система обладала свойством астатизма по отношению к постоянному задающему g(t) и возмущающему f(t) воздействию.
Для того чтобы система обладала астатизмом необходимо выполнение условия - величина установившейся ошибки по заданию и возмущению должна быть равна 0. Проверим это условие на регуляторах:
· Выбираем П - регулятор
Т.к. величина установившейся ошибки не равна нулю, то система не обладает астатизмом,, а следовательно считать ошибку по каналу f-E нет необходимости.
· Выбираем И - регулятор
Обе ошибки равны нулю следовательно система с И - регулятором является астатической. Данный регулятор нам подходит
· Проверим наш ПИ - регулятор
Система так же является астатической, что является тем, что выбранный нами ПИ - регулятор подходит к условию задания.
3. Найти область значений постоянной времени регулятора для И - закона управления, обеспечивающих устойчивость системы.
Найдем передаточную функцию замкнутой системы по каналу g-y и затем выделим ее характеристическое уравнение.
Характеристическое уравнение Q(p):
Уравнение 3 порядка, следовательно, для устойчивости системы необходимо чтобы произведение внутренних коэффициентов характеристического уравнения было больше произведения внешних коэффициентов.
Следовательно область значений для И - закона можно определить интервалом .
4. Для И - закона регулирования найти минимальную возможную величину установившейся ошибки, если g(t)=2t и f(t)=-3t.
Для определения минимальной установившейся ошибки нам необходимо узнать при каком значении постоянной времени система находится на границе устойчивости. Для этого отыщем передаточную функцию и характеристическое уравнение системы с И - регулятором. И затем найдем значение
Находим значение постоянной времени на границе устойчивости:
Найдем величину установившейся ошибки при g(t)=2t и f(t)=-3t
Поскольку задающее воздействие у нас g(t)=2t, то используя преобразование Лапласа получаем:
тогда величина установившейся ошибки будет
Подставляем полученное значение и получаем
Найдем ошибку по каналу f-E
Подставляем и получаем
Тогда
5. Построить, с использованием ЭВМ, область устойчивости, на плоскости, параметров регулятора при использовании ПИ - закона, обосновать возможность и путь получения допустимой установившейся ошибки при и .
Для построение области устойчивости необходимо найти характеристическое уравнение передаточной функции для данного регулятора.
Отсюда:
Запишем условие, при котором система находится на границе устойчивости:
Выразим зависимость
от
Строим область устойчивости по по лученной зависимости:
Ти |
Кп |
|
0,01 |
4,00 |
|
0,02 |
1,50 |
|
0,03 |
0,67 |
|
0,04 |
0,25 |
|
0,05 |
0,00 |
|
0,06 |
-0,17 |
|
0,07 |
-0,29 |
|
0,08 |
-0,38 |
|
0,09 |
-0,44 |
Таким образом выбирая значения параметров регулятора над границей устойчивости - мы получаем устойчивую систему, и наоборот.
6. Вычислить значения параметров ПИ - регулятора, обеспечивающих устойчивость и установившуюся ошибку в системе = 0,06 при g(t)=2t и f(t)=0
Поскольку возмущение f(t)=0, то . Найдем :
для этого найдем передаточную функцию замкнутой системы по каналу g-E
По условию , тогда подставим это значение в получившееся выражение:
Таким образом для получения в системе установившейся ошибке равной 0,06 необходимо задать параметру постоянной времени значение 0,03.
7. Для интегрального закона регулирования и начальных условий , выбрать оптимальное значение постоянной времени регулятора по критерию: (Рассматривается движение в системе при g(t)=f(t)=0 и ненулевых начальных условиях).
Для решения мы будем использовать метод Мондельштама. Для этого нам необходимо найти передаточную функцию замкнутой системы и взять характеристическое уравнение:
Получаем уравнение:
Поочередно умножаем его на и на Е
Интегрируем полученное уравнение по частям
Получаем:
Интегрируем полученное уравнение по частям
Получаем:
Из полученных уравнений составим систему уравнений:
Выбираем =1,054
8. Для найденного в пункте 7 значения постоянной времени регулятора построить с помощью ЭВМ вещественную частотную характеристику P(щ), приняв входным воздействием g(t) и входной координатой E(t)/
Для нахождения вещественной характеристики нам понадобится передаточная функция замкнутой системы по каналу g-y.
Перейдем в частотную область p=jщ:
Домножаем на сопряженное знаменателю число и получаем:
Отделяем действительную часть U(щ):
При =1,054
w |
p(w) |
|
0 |
1 |
|
0,1 |
1,010058 |
|
0,2 |
1,039191 |
|
0,3 |
1,079471 |
|
0,4 |
1,099022 |
|
0,5 |
0,997092 |
|
0,6 |
0,58593 |
|
0,7 |
-0,06976 |
|
0,8 |
-0,48243 |
|
0,9 |
-0,56794 |
|
1 |
-0,5208 |
|
1,1 |
-0,44696 |
|
1,2 |
-0,3782 |
|
1,3 |
-0,32081 |
|
1,4 |
-0,27428 |
|
1,5 |
-0,23666 |
|
1,6 |
-0,20606 |
|
1,7 |
-0,18095 |
|
1,8 |
-0,16013 |
|
1,9 |
-0,14269 |
|
2 |
-0,12796 |
|
2,5 |
-0,08003 |
|
3 |
-0,05481 |
|
3,5 |
-0,03991 |
|
4 |
-0,03037 |
|
4,5 |
-0,02389 |
|
5 |
-0,01929 |
|
5,5 |
-0,0159 |
|
6 |
-0,01334 |
9. По вещественной характеристике P(щ) пункта 8 построить переходной процесс E(t) при единичном ступенчатом изменении g(t) и нулевых начальных условиях методом трапециидальных частотных характеристик.
Для построения переходного процесса нам необходимо разбить получившуюся вещественную характеристику на трапеции и построить переходный процесс для каждой из полученных трапеций.
R(0) |
Wo |
Wd |
? |
|||
I |
- |
0,099021688 |
0,38 |
0,1 |
0,263158 |
|
II |
+ |
1,666965285 |
0,88 |
0,43 |
0,488636 |
|
III |
- |
0,567943597 |
6 |
0,95 |
0,158333 |
1-я трапеция
t табл |
h(?) |
t=t табл/Wo |
h=R(0)*h(?) |
|
0 |
0 |
0 |
0 |
|
0,5 |
0,199 |
1,315789474 |
-0,019705316 |
|
1 |
0,386 |
2,631578947 |
-0,038222372 |
|
1,5 |
0,56 |
3,947368421 |
-0,055452145 |
|
2 |
0,709 |
5,263157895 |
-0,070206377 |
|
2,5 |
0,833 |
6,578947368 |
-0,082485066 |
|
3 |
0,928 |
7,894736842 |
-0,091892127 |
|
3,5 |
0,994 |
9,210526316 |
-0,098427558 |
|
4 |
1,039 |
10,52631579 |
-0,102883534 |
|
4,5 |
1,057 |
11,84210526 |
-0,104665924 |
|
5 |
1,067 |
13,15789474 |
-0,105656141 |
|
5,5 |
1,067 |
14,47368421 |
-0,105656141 |
|
6 |
1,054 |
15,78947368 |
-0,104368859 |
|
6,5 |
1,043 |
17,10526316 |
-0,103279621 |
|
7 |
1,035 |
18,42105263 |
-0,102487447 |
|
7,5 |
1,025 |
19,73684211 |
-0,10149723 |
|
8 |
1,024 |
21,05263158 |
-0,101398209 |
|
8,5 |
1,022 |
22,36842105 |
-0,101200165 |
|
9 |
1,025 |
23,68421053 |
-0,10149723 |
|
9,5 |
1,027 |
25 |
-0,101695274 |
|
10 |
1,027 |
26,31578947 |
-0,101695274 |
|
10,5 |
1,028 |
27,63157895 |
-0,101794295 |
|
11 |
1,025 |
28,94736842 |
-0,10149723 |
|
11,5 |
1,021 |
30,26315789 |
-0,101101144 |
|
12 |
1,015 |
31,57894737 |
-0,100507013 |
|
12,5 |
1,01 |
32,89473684 |
-0,100011905 |
|
13 |
1,005 |
34,21052632 |
-0,099516797 |
|
13,5 |
1 |
35,52631579 |
-0,099021688 |
|
14 |
0,997 |
36,84210526 |
-0,098724623 |
|
14,5 |
0,996 |
38,15789474 |
-0,098625601 |
|
15 |
0,995 |
39,47368421 |
-0,09852658 |
|
15,5 |
0,995 |
40,78947368 |
-0,09852658 |
|
16 |
0,995 |
42,10526316 |
-0,09852658 |
|
16,5 |
0,995 |
43,42105263 |
-0,09852658 |
|
17 |
0,995 |
44,73684211 |
-0,09852658 |
|
17,5 |
0,995 |
46,05263158 |
-0,09852658 |
|
18 |
0,995 |
47,36842105 |
-0,09852658 |
|
18,5 |
0,994 |
48,68421053 |
-0,098427558 |
|
19 |
0,992 |
50 |
-0,098229515 |
|
19,5 |
0,991 |
51,31578947 |
-0,098130493 |
|
20 |
0,991 |
52,63157895 |
-0,098130493 |
2-я трапеция
t табл |
h(?) |
t=t табл/Wo |
h=R(0)*h(?) |
|
0 |
0 |
0 |
0 |
|
0,5 |
0,24 |
0,568181818 |
0,400071669 |
|
1 |
0,461 |
1,136363636 |
0,768470997 |
|
1,5 |
0,665 |
1,704545455 |
1,108531915 |
|
2 |
0,833 |
2,272727273 |
1,388582083 |
|
2,5 |
0,967 |
2,840909091 |
1,611955431 |
|
3 |
1,061 |
3,409090909 |
1,768650168 |
|
3,5 |
1,115 |
3,977272727 |
1,858666293 |
|
4 |
1,142 |
4,545454545 |
1,903674356 |
|
4,5 |
1,138 |
5,113636364 |
1,897006495 |
|
5 |
1,118 |
5,681818182 |
1,863667189 |
|
5,5 |
1,092 |
6,25 |
1,820326092 |
|
6 |
1,051 |
6,818181818 |
1,751980515 |
|
6,5 |
1,018 |
7,386363636 |
1,696970661 |
|
7 |
0,993 |
7,954545455 |
1,655296528 |
|
7,5 |
0,974 |
8,522727273 |
1,623624188 |
|
8 |
0,966 |
9,090909091 |
1,610288466 |
|
8,5 |
0,966 |
9,659090909 |
1,610288466 |
|
9 |
0,97 |
10,22727273 |
1,616956327 |
|
9,5 |
0,975 |
10,79545455 |
1,625291153 |
|
10 |
0,982 |
11,36363636 |
1,63695991 |
|
10,5 |
0,987 |
11,93181818 |
1,645294737 |
|
11 |
0,993 |
12,5 |
1,655296528 |
|
11,5 |
0,997 |
13,06818182 |
1,66196439 |
|
12 |
0,997 |
13,63636364 |
1,66196439 |
|
12,5 |
0,997 |
14,20454545 |
1,66196439 |
|
13 |
0,997 |
14,77272727 |
1,66196439 |
|
13,5 |
0,998 |
15,34090909 |
1,663631355 |
|
14 |
1 |
15,90909091 |
1,666965285 |
|
14,5 |
1,002 |
16,47727273 |
1,670299216 |
|
15 |
1,005 |
17,04545455 |
1,675300112 |
|
15,5 |
1,008 |
17,61363636 |
1,680301008 |
|
16 |
1,011 |
18,18181818 |
1,685301904 |
|
16,5 |
1,011 |
18,75 |
1,685301904 |
|
17 |
1,012 |
19,31818182 |
1,686968869 |
|
17,5 |
1,009 |
19,88636364 |
1,681967973 |
|
18 |
1,008 |
20,45454545 |
1,680301008 |
|
18,5 |
1,006 |
21,02272727 |
1,676967077 |
|
19 |
1,001 |
21,59090909 |
1,668632251 |
|
19,5 |
0,998 |
22,15909091 |
1,663631355 |
|
20 |
0,996 |
22,72727273 |
1,660297424 |
3-я трапеция
t табл |
h(?) |
t=t табл/Wo |
h=R(0)*h(?) |
|
0 |
0 |
0,0000 |
0,0000 |
|
0,5 |
0,184 |
0,0833 |
-0,1045 |
|
1 |
0,256 |
0,1667 |
-0,1454 |
|
1,5 |
0,516 |
0,2500 |
-0,2931 |
|
2 |
0,655 |
0,3333 |
-0,3720 |
|
2,5 |
0,833 |
0,4167 |
-0,4731 |
|
3 |
0,863 |
0,5000 |
-0,4901 |
|
3,5 |
0,928 |
0,5833 |
-0,5271 |
|
4 |
0,974 |
0,6667 |
-0,5532 |
|
4,5 |
0,977 |
0,7500 |
-0,5549 |
|
5 |
1,012 |
0,8333 |
-0,5748 |
|
5,5 |
1,019 |
0,9167 |
-0,5787 |
|
6 |
1,013 |
1,0000 |
-0,5753 |
|
6,5 |
1,009 |
1,0833 |
-0,5731 |
|
7 |
1,006 |
1,1667 |
-0,5714 |
|
7,5 |
1,006 |
1,2500 |
-0,5714 |
|
8 |
1,008 |
1,3333 |
-0,5725 |
|
8,5 |
1,01 |
1,4167 |
-0,5736 |
|
9 |
1,016 |
1,5000 |
-0,5770 |
|
9,5 |
1,022 |
1,5833 |
-0,5804 |
|
10 |
1,025 |
1,6667 |
-0,5821 |
|
10,5 |
1,028 |
1,7500 |
-0,5838 |
|
11 |
1,029 |
1,8333 |
-0,5844 |
|
11,5 |
1,027 |
1,9167 |
-0,5833 |
|
12 |
1,025 |
2,0000 |
-0,5821 |
|
12,5 |
1,022 |
2,0833 |
-0,5804 |
|
13 |
1,019 |
2,1667 |
-0,5787 |
|
13,5 |
1,017 |
2,2500 |
-0,5776 |
|
14 |
1,016 |
2,3333 |
-0,5770 |
|
14,5 |
1,015 |
2,4167 |
-0,5765 |
|
15 |
1,014 |
2,5000 |
-0,5759 |
|
15,5 |
1,014 |
2,5833 |
-0,5759 |
|
16 |
1,014 |
2,6667 |
-0,5759 |
|
16,5 |
1,014 |
2,7500 |
-0,5759 |
|
17 |
1,013 |
2,8333 |
-0,5753 |
|
17,5 |
1,012 |
2,9167 |
-0,5748 |
|
18 |
1,011 |
3,0000 |
-0,5742 |
|
18,5 |
1,009 |
3,0833 |
-0,5731 |
|
19 |
1,008 |
3,1667 |
-0,5725 |
|
19,5 |
1,006 |
3,2500 |
-0,5714 |
|
20 |
1,005 |
3,3333 |
-0,5708 |
|
20,5 |
1,004 |
3,4167 |
-0,5702 |
|
21 |
1,003 |
3,5000 |
-0,5696 |
|
21,5 |
1,003 |
3,5833 |
-0,5696 |
|
22 |
1,002 |
3,6667 |
-0,5691 |
|
22,5 |
1,002 |
3,7500 |
-0,5691 |
|
23 |
1,002 |
3,8333 |
-0,5691 |
|
23,5 |
1,002 |
3,9167 |
-0,5691 |
|
24 |
1,001 |
4,0000 |
-0,5685 |
|
24,5 |
1 |
4,0833 |
-0,5679 |
|
25 |
1 |
4,1667 |
-0,5679 |
|
25,5 |
0,999 |
4,2500 |
-0,5674 |
|
26 |
0,999 |
4,3333 |
-0,5674 |
Суммируем графически полученные процесы и получаем
10. Определить время регулирования и максимальное перерегулирование в системе.
11. Разработать и начертить структурную схему комбинированной САУ, инвариантной по отношению к контролируемому возмущению .
Привести передаточную функцию устройства управления.
Проверить выполнение условия инвариантности.
Условия инвариантности:
, если
, если
При выборе корректирующих звеньев в виде обратных передаточных функций возникает проблема. Она обычно связана с тем, что порядок числителя корректирующего звена становится больше порядка знаменателя. Это означает, что частотные характеристики таких звеньев являются расходящимися, что говорит о том, что физически такие звенья нереализуемы. В тех случаях, когда корректирующие звенья применять необходимо, порядок числителя этих звеньев искусственно приравнивают к порядку знаменателя, отбрасывая в числители высшие порядки.
12. Предложить и обосновать методы компенсации действия неконтролируемых возмущений, если известен класс функций, которыми они описываются.
астатический автоматический управление моделирование программа
Решить проблему инвариантности можно, если известна предварительная информация о классе возмущающих воздействий. В частности, если известен математический аппарат, описывающий функцию f(t), заданный в виде решения некоторого дифференциального уравнения.
Процедура синтеза предусматривает:
1. восстановление вида дифференциального уравнения по заданному решению;
2. получение характеристического уравнения;
3. выбор передаточной функции регулятора, в которой знаменатель совпадает с видом полученного характеристического уравнения. Числитель передаточной функции регулятора выбирается того же порядка, что и знаменатель;
4. неизвестные коэффициенты числителя передаточной функции регулятора определяются из условий устойчивости замкнутой системы.
13. Провести моделирование в программе MatLab. Определить настройки регулятора, обеспечивающего минимизацию времени регулирования.
При нулевом задающем воздействии со значением регулятора, полученными в 7 пункте:
При единичном задающем воздействии:
Для снижения времени регулирования можно немного увеличить значения регулятора примерно до 1,085.
Так же значительно уменьшает время регулирование и введение пропорциональной составляющей, т.е. использование ПИ - регулятора. С его помощью легко можно уменьшить время регулирования примерно в 1,7 раза.
Список используемой литературы
· Лекции по курсу ТАУ, Румянцев Ю.Д.
· "Теория автоматического управления", Воронов А.А.
· "Теория систем автоматического управления", Бесекерский В.А.
Приложение
Размещено на Allbest.ru
Подобные документы
Составление и анализ математической модели объекта управления и структурной схемы системы. Построение областей устойчивости, требуемой точности и быстродействия статического регулятора. Анализ замкнутой системы управления с непрерывным регулятором.
курсовая работа [1,2 M], добавлен 12.04.2012Синтез системы автоматического управления корневым методом, разработанным Т. Соколовым. Определение передаточных функций по задающему и возмущающему воздействиям. Оценка устойчивости замкнутой нескорректированной системы регулирования по критерию Гурвица.
курсовая работа [1,3 M], добавлен 26.01.2015Особенности структурной и функциональной схем систем автоматического управления, характеристика и определение запаса ее устойчивости. Принцип управления по замкнутому циклу и ошибки переходного процесса. Использование регулятора для коррекции системы.
контрольная работа [827,6 K], добавлен 09.12.2011Характеристика системы управления двигателем постоянного тока. Моделирование системы управления в среде Matlab 6.1. Подбор параметров регуляторов структурной схемы в соответствии с предъявляемыми требованиями. Исследование электрической схемы системы.
курсовая работа [1,4 M], добавлен 29.11.2010Обзор методов составления математических моделей систем автоматического управления. Математические модели системы в векторно-матричной форме записи. Моделирование в пакете программы Simulink. Оценка устойчивости системы, рекомендации по ее применению.
курсовая работа [514,5 K], добавлен 10.11.2011Порядок оценки точности системы автоматического управления по величине установившейся ошибки при типовых воздействиях, механизм ее повышения. Разновидности ошибок и методика их вычисления. Определение ошибок по виду частотных характеристик системы.
реферат [103,3 K], добавлен 11.08.2009Исследование системы автоматического управления при помощи программного обеспечения MATLAB и пакета Simulink. Изучение замкнутой системы согласно критериям устойчивости Гурвица, Михайлова и Найквиста. Реализация модели "жесткого" спутника Земли.
методичка [911,6 K], добавлен 10.10.2010Исследование линейных динамических моделей в программном пакете Matlab и ознакомление с временными и частотными характеристиками систем автоматического управления. Поиск полюса и нуля передаточной функции с использованием команд pole, zero в Matlab.
лабораторная работа [53,1 K], добавлен 11.03.2012Схемотехнический синтез системы автоматического управления. Анализ заданной системы автоматического управления, оценка ее эффективности и функциональности, описание устройства и работы каждого элемента. Расчет характеристик системы путем моделирования.
курсовая работа [3,4 M], добавлен 21.11.2012Математические процессы, происходящие в системах автоматического управления. Определение передаточных функций разомкнутой и замкнутой систем, критерии устойчивости. Физический смысл логарифмических асимптотических амплитудных частотных характеристик.
курсовая работа [2,3 M], добавлен 12.05.2014