Подякувати Студентському архіву довільною сумою
Admin
26.02.2023 12:38
Дякуємо, що користуєтесь нашим архівом!
Інформація про навчальний заклад
ВУЗ:
Національний університет Львівська політехніка
Інститут:
Не вказано
Факультет:
Не вказано
Кафедра:
Кафедра автоматизованих систем управління
Інформація про роботу
Рік:
2012
Тип роботи:
Розрахункова робота
Предмет:
Теорія управління
Частина тексту файла (без зображень, графіків і формул):
Міністерство освіти і науки України
Національний університет “Львівська політехніка”
Кафедра автоматизованих систем управління
Розрахунково-графічна робота
З дисципліни
Теорія Управління
Теоретичні основи управління
k0
1
5
4
14
1
Завдання
1. Записати передатну функцію при даних коефіцієнтах:
2. Побудувати амплітудно - фазочастотну характеристику (АФЧХ).
3. Побудувати амплітудно - частотну характеристику (АЧХ).
4. Побудувати фазочастотну характеристику (ФЧХ).
5. Побудувати логарифмічну амплітудно - фазову частотну
характеристику (лафчх).
6. Оцінити систему на стійкість за двома критеріями:
а) алгебраїчним;
б) частотним.
7. Знайти вираз перехідної функції h(t). Побудувати її графік.
8. Побудувати розміщення нулів і полюсів функції.
9. Здійснити корекцію системи:
а) якщо система не стійка, то знайти коректуючі елементи, такі, щоб вона стала
на межу стійкості;
б) якщо система на межі стійкості, то внести такі елементи, щоб вона мала запас
стійкості за фазою φ= π/12 = 15°;
в) якщо система стійка, то внести коректуючі елементи, щоб запас стійкості
збільшився вдвічі. .
10. Знайти вираз перехідної функції h(t) відкоректованої системи.
11. Здійснити декомпозицію системи третього порядку на типові ланки.
12. Побудувати фазовий портрет відкоректованої системи.
1. Записати передатну функцію при заданих коефіцієнтах.
Передатною функцією називається відношення перетвореної за Лапласом вихідної дії до перетвореної за Лапласом вхідної дії при нульових умовах і відсутності збурень.
;
За даними коефіцієнтами передатна функція має такий вигляд:
;
2. Побудувати амплітудно - фазочастотну характеристику (АФЧХ).
Сам по собі вираз b(w) є дійсним і не повинен містити ніяких j.
Передатна функція набуде вигляду: ;
a()
b()
a()
b()
a()
b()
0
0,071429
0
0,65
0,081022
-0,01585
1,3
0,139375
-0,07541
0,05
0,071478
-0,00102
0,7
0,082832
-0,01762
1,35
0,15025
-0,09037
0,1
0,071626
-0,00205
0,75
0,084878
-0,01956
1,4
0,162811
-0,11071
0,15
0,071875
-0,00309
0,8
0,087191
-0,0217
1,45
0,176737
-0,13943
0,2
0,072226
-0,00415
0,85
0,08981
-0,02409
1,5
0,19027
-0,18162
0,25
0,072683
-0,00523
0,9
0,092778
-0,02677
1,55
0,197146
-0,24561
0,3
0,073252
-0,00634
0,95
0,096152
-0,02982
1,6
0,17782
-0,34141
0,35
0,073937
-0,0075
1
0,1
-0,03333
1,65
0,084362
-0,4589
0,4
0,074745
-0,0087
1,05
0,104405
-0,03742
1,7
-0,11958
-0,50143
0,45
0,075687
-0,00996
1,1
0,109472
-0,04226
1,75
-0,29733
-0,37166
0,5
0,076771
-0,01129
1,15
0,115328
-0,04807
3
-0,02614
0,012648
0,55
0,078011
-0,0127
1,2
0,122133
-0,05518
4
-0,00991
0,007207
0,6
0,079422
-0,01422
1,25
0,130075
-0,06405
7
-0,00151
0,002064
Рис.1. АФЧХ
3. Побудувати амплітудно - частотну характеристику (АЧХ).
Для побудови АЧХ у середовищі MATLAB виконуємо послідовність команд:
num=[1];
den=[1 5 4 14];
[mag,phase,w]=bode(num,den);
plot(w,mag(:)),grid
Після виконання таких команд отримаємо графік АЧХ, зображений на рис.2.
Рис.2. АЧХ
4. Побудувати фазочастотну характеристику (ФЧХ).
Для побудови ФЧХ у середовищі MATLAB виконуємо послідовність команд:
num=[1];
den=[1 5 4 14];
[mag,phase,w]=bode(num,den);
plot(w,phase(:)),grid on
Після виконання таких команд отримаємо графік ФЧХ, зображений на рис.3.
Рис.3. ФЧХ
5. Побудувати логарифмічну амплітудно–фазочастотну характеристику (лафчх).
Логарифмічна амплітудно - фазочастотна характеристика є сукупність двох характеристик - ЛАЧХ та ФЧХ, побудованих на одному графіку.
Зручніше користуватись десятковим логарифмом і будувати окремо логарифмічно-амплітудну і фазову характеристики.
– одиниці вимірювання дБ (децибел)
Бел представляє собою логарифмічну одиницю, що відповідає десятикратному збільшенню потужності.
Особливість логарифмічної характеристики, зображеної на рис. 4, полягає в тому, що вона відображається сукупністю відрізків прямої.
Для побудови ЛАФЧХ у середовищі MATLAB виконуємо послідовність команд:
num=[1];
den=[1 5 4 14];
[mag,phase,w]=bode(num,den);
semilogx(w,20*log10(mag(:))),grid on
Рис.4. ЛАФЧХ
6. Оцінити систему на стійкість за двома критеріями
Стійкість – здатність системи повертатись в стан рівноваги після припинення дії вимушуючих сил (вхідна дія, збурення).
Для того, щоб оцінити систему на стійкість необхідно знаменник передатної функції прирівняти до 0.
Критерії стійкості (кс) поділяються на дві групи:
а) алгебраїчні (кс Вишнєградського, кс Рауса, кс Гурвіца);
б)частотні (кс Михайлова, кс Найквіста, логарифмічний кс);
Стійкість нас цікавить тому, що нестійкі системи – непрацездатні. Невеличкі збурення виводять їх з планової траєкторії, до якої вони ніколи не зможуть повернутися.
Стійкість – внутрішня властивість системи, яка не залежить від величини вхідної дії чи від величини збурення.
Для великих систем ця характеристика вироджується в надійність чи живучість.
6.1.Оцінити систему на стійкість за алгебраїчним критерієм.
З алгебраїчних критеріїв для оцінки стійкості системи третього порядку доцільно використати критерій Вишнєградського.
Цей критерій використовується для систем не вище третього порядку.
Якщо різниця добутків середніх і крайніх коефіцієнтів поліному (характеристичного рівняння) більше від нуля і кожен з коефіцієнтів є додатним, то система стійка.
Тобто якщо виконується умова:
Т2*Т1-Т3*Т0>0;
Згідно завдання :
5*4-1*14>0 Отже, розімкнута система стійка.
6.2.Оцінити систему на стійкість за частотним критерієм.
Приймемо, що система незамкнута, тоді оцінити її на стійкість можна за критерієм Михайлова. Відокремлюючи дійсну і уявну частину, поліном D(p) приводимо до виду:
D(p)=a(()+jb(()
– парні степені
– непарні степені
Геометричне місце точок кінця вектора D(j() при зміні частоти 0<(<( називається годографом Михайлова.
Динамічна система, що описується лінійним диференційним рівнянням n-го порядку стійка, якщо при зміні частоти від 0 до ( годограф Михайлова послідовно проходить в напрямку проти годинникової стрілки n квадрантів комплексної площини і не перетворюється в 0.
Критерій Михайлова, використовується для розімкнених систем управління.
Два наступних критерія, а саме – критерій стійкості Найквіста та логарифмічний критерій стійкості використовуються для оцінки на стійкість замкнутих систем управління.
Оцінку стійкості почнемо з визначення коренів рівнянь:
Складемо таблицю значень:
a()
b()
0
14
0
1
9
3
2
-6
0
3
-31
-15
Рис.5.Годограф Михайлова
Як видно з рис. 5 Годограф Михайлова для даної системи проходить в напрямку проти годинникової стрілки і не перетворюється в нуль, тому це свідчить про те що система є стійкою.
7. Знайти вираз перехідної функції h(t). Побудувати її графік.
Перехідна функція – показує перехідний процес на виході ланки, якщо на вході діє одинична ступінчата функція. Одинична ступінчата функція має розмірність таку, як і вхідна величина.
Для побудови графіку перехідної функції h(t) у середовищі MATLAB виконуємо послідовність команд:
num=[1];
den=[1 5 4 14];
t=0:0.1:20;
y=step(num,den,t);
plot(t,y),grid
Рис.6. Графік перехідної функції
Визначити запас стійкості за фазою та частотою зрізу можна таким способом: сформувати m-файл і набрати в ньому послідовність команд
num=[1];
den=[1 5 4 14];
h=tf(num,den);
margin(h); grid on
Рис. 7.
8. Знайти розміщення нулів і полюсів системи
Нулі системи – це ті р, при яких чисельник передатної функції перетворюється в нуль, а полюси – це р, при яких знаменник перетворюється в нуль.
Розміщення нулів і полюсів системи можна отримати, користуючись засобами MATLAB:
num=[1];
den=[1 5 4 14];
[mag, phase, w]=bode(num, den);
t=0:0.1:16;
y=step(num, den, t);
[p, z]=pzmap(num, den);
hold; plot (p, ‘bx’); grid on;
Рис.8. Нулі і полюси системи
Для знаходження числових полюсів системи потрібно набрати в середовищі MATLAB послідовність команд:
w=tf([1],[1 5 4 14]);
pole(w)
Отримаємо:
ans =
-4.7762
-0.1119 + 1.7084i
-0.1119 - 1.7084i
9. Провести корекцію системи
Визначимо стійкість системи:
5*4-1*(14+1)>0 - система стійка.
Оскільки дана система є стійкою, необхідно внести коректуючий множник, щоб збільшити запас стійкості.
Приймемо k=1.7, тоді розімкнена система виглядає так :
Замкнена система має наступний вигляд:
Знаходимо kd :
5*(4+ kd)-1*15,7=0; 5*kd=6,7; kd=1,34;
Wk(p) – відкоректована система.
Для того, щоб побудувати Bode Diagram відкоректованої системи необхідно взяти значення з наступної перехідної функції :
Визначити запас стійкості за фазою та частотою зрізу можна таким способом: сформувати m-файл і набрати в ньому послідовність команд.
num=[1.34 1.7];
den=[1 5 4 14];
h=tf(num,den);
margin(h)
Рис.9.
10. Знайти вираз перехідної функції відкоректованої системи
Даний вираз обраховано за допомогою програми MathCad:
11. Здійснити декомпозицію системи третього порядку на типові ланки
При знаходженні нулів і полюсів функції виявилося, що система має один дійсний і два комплексних корені. Тому систему можна декомпонувати на дві типові ланки, що проілюсровано на Рис.10.
Рис.10. Декомпонована система в середовищі MATLAB Simulink
Рис.11. Результат моделювання роботи окремих компонентів системи
12. Побудувати фазовий портрет системи
Для наочного зображення складних процесів управління застосовують поняття фазового простору, яке полягає в наступному: диференційне рівняння замкнутої системи n-го порядку приводиться до системи рівнянь І порядку.
Система n-го порядку замкнута і описується системою n рівнянь І порядку, та вводиться поняття зображуваної точки М0, а траєкторія руху точки називається фазовим портом.
Побудуємо фазовий портрет за допомогою середовища MATLAB Simulink.
Рис.12. Схема для побудови фазового портрету
Рис.13. Фазовий портрет системи
Висновок: виконуючи цю розрахунково-графічну роботу, я засвоїв матеріал курсу теоретичні основи управління, виконав всі поставлені в роботі завдання та закріпив навички роботи в середовищі MATHLAB. В результаті виконання роботи я встановив, що задана система є на межі стійкості, тому я ввів коректуючу ланку для того щоб зробити систему стійкою, а також здійснив декомпозицію системи на типові ланки і побудував фазовий портер системи.
03.02.2013 12:02-
Ділись своїми роботами та отримуй миттєві бонуси!
0%
Оголошення від адміністратора