Коэффициенты характеристического полинома замкнутой скорректированной системы имеют следующие значения: d0 = 0,002625, d1 = 0,089175, d2 = 0,91025, d3 = 3,01497, d4 = 1,5499.
Матрица Гурвица имеет вид:
Диагональные миноры матрицы Гурвица:
Д1 = 0,089175 > 0;
;
.
Вывод: замкнутая скорректированная система устойчива, так как не имеет отрицательных миноров.
4.3 Исследование скорректированной АСР по критерию Михайлова
В пункте 4.2 был получен характеристический полином такой замкнутой скорректированной системы DЗАМ(p):
Тогда получим частотный годограф Михайлова путем перевода характеристического полинома (4.15) в частотную область:
Получим действительную и мнимую часть частотного годографа Михайлова (4.16), возведя частотный оператор j в соответствующую степень:
Вычисляем значения действительной и мнимой части частотного годографа Михайлова при изменении часты от 0 до требуемого значения, при котором можно сделать вывод об устойчивости системы. Шаг изменения частоты принимаем произвольный с учетом удобства восприятия графика. Вычисления проводим в MS Excel.
Таблица 4.1 Значения действительной и мнимой части частотного годографа Михайлова
|
Частота сигнала |
Действительная часть UD() |
Мнимая часть VD() |
|
|
0 |
1,55 |
0 |
|
|
2 |
-2,049 |
5,31654 |
|
|
4 |
-12,342 |
6,35268 |
|
|
6 |
-27,817 |
-1,17198 |
|
|
8 |
-45,954 |
-21,5378 |
|
|
10 |
-63,225 |
-59,0253 |
|
|
12 |
-75,094 |
-117,915 |
|
|
14 |
-76,017 |
-202,487 |
|
|
16 |
-59,442 |
-317,021 |
|
|
18 |
-17,809 |
-465,799 |
|
|
20 |
57,45 |
-653,101 |
По таблице 4.1 строим годограф Михайлова (рис. 4.1).
Рис. 4.1 Годограф Михайлова скорректированной системы
Вывод: замкнутая скорректированная система устойчива, т.к. годограф Михайлова начинается на положительной действительной полуоси и огибает против часовой стрелки начало координат, проходя последовательно 4 квадранта, где 4 - порядок характеристического полинома замкнутой системы.
4.4 Исследование скорректированной АСР по критерию Найквиста
Для построения АФЧХ разомкнутой скорректированной системы необходимо перейти к частотной форме записи передаточной функции.
Передаточная функция разомкнутой системы:
Проведем преобразования, выделив общий коэффициент усиления разомкнутой системы kРАЗ:
Передаточная функция разомкнутой системы в частотной форме (p = jщ):
Запишем выражения амплитуд А() и фаз ц() частотных передаточных функции для каждого звена (табл. 4.2). Вычисляем значения амплитуд А() и фаз ц() при изменении часты от 0 до значения, при котором UРАЗ(щ) и VРАЗ(щ) станут равны нулю. Также вычислим амплитуду АРАЗ() и фазу цРАЗ() передаточной функции разомкнутой системы, а также вещественную UРАЗ(щ) и мнимую VРАЗ(щ) составляющую частотной передаточной функции разомкнутой системы WРАЗ(jщ) (3.18 - 3.23). Шаг изменения частоты принимаем произвольный с учетом удобства восприятия графика. Вычисления проводим в MS Excel (табл. 4.3).
Таблица 4.2 Выражения амплитуд А() и фаз ц() передаточных функции звеньев системы
|
Частотная передаточная функция |
Амплитуда А() |
Фаза ц() |
|
Таблица 4.3 Значения амплитуд А() и фаз ц() передаточных функции звеньев системы
|
щ |
Амплитуды звеньев |
Фазы звеньев |
Амплитуда WРАЗ(jщ) |
Фаза WРАЗ(jщ) |
Действит. часть WРАЗ(jщ) |
Мнимая часть WРАЗ(jщ) |
|||||
|
A1(щ) |
A2(щ) |
A3(щ) |
ц1(щ) |
ц2(щ) |
ц3(щ) |
AРАЗ(щ) |
цРАЗ(щ) |
UРАЗ(щ) |
VРАЗ(щ) |
||
|
0,0 |
1,000 |
0,550 |
1,000 |
0,000 |
0,000 |
0,000 |
0,550 |
0,000 |
0,55000 |
0,00000 |
|
|
0,1 |
0,971 |
0,550 |
1,000 |
-0,215 |
-0,020 |
-0,015 |
0,534 |
-0,250 |
0,51709 |
-0,13203 |
|
|
0,2 |
0,896 |
0,550 |
1,000 |
-0,404 |
-0,040 |
-0,030 |
0,492 |
-0,474 |
0,43810 |
-0,22458 |
|
|
0,3 |
0,803 |
0,549 |
0,999 |
-0,554 |
-0,060 |
-0,045 |
0,441 |
-0,659 |
0,34855 |
-0,26976 |
|
|
0,4 |
0,712 |
0,548 |
0,999 |
-0,666 |
-0,080 |
-0,060 |
0,390 |
-0,806 |
0,27015 |
-0,28139 |
|
|
0,5 |
0,632 |
0,547 |
0,999 |
-0,747 |
-0,100 |
-0,075 |
0,345 |
-0,922 |
0,20863 |
-0,27500 |
|
|
0,6 |
0,564 |
0,546 |
0,998 |
-0,805 |
-0,119 |
-0,090 |
0,307 |
-1,014 |
0,16229 |
-0,26074 |
|
|
0,7 |
0,507 |
0,545 |
0,997 |
-0,845 |
-0,139 |
-0,105 |
0,275 |
-1,089 |
0,12767 |
-0,24394 |
|
|
0,8 |
0,460 |
0,543 |
0,996 |
-0,872 |
-0,159 |
-0,120 |
0,249 |
-1,150 |
0,10162 |
-0,22712 |
|
|
0,9 |
0,421 |
0,541 |
0,995 |
-0,889 |
-0,178 |
-0,135 |
0,227 |
-1,202 |
0,08175 |
-0,21135 |
|
|
1,0 |
0,388 |
0,539 |
0,994 |
-0,899 |
-0,197 |
-0,150 |
0,208 |
-1,246 |
0,06636 |
-0,19700 |
|
|
2,0 |
0,229 |
0,511 |
0,977 |
-0,833 |
-0,381 |
-0,297 |
0,114 |
-1,511 |
0,00683 |
-0,11386 |
|
|
3,0 |
0,178 |
0,472 |
0,949 |
-0,705 |
-0,540 |
-0,441 |
0,080 |
-1,687 |
-0,00924 |
-0,07905 |
|
|
4,0 |
0,155 |
0,429 |
0,913 |
-0,595 |
-0,675 |
-0,580 |
0,061 |
-1,850 |
-0,01678 |
-0,05861 |
|
|
5,0 |
0,144 |
0,389 |
0,871 |
-0,508 |
-0,785 |
-0,712 |
0,049 |
-2,006 |
-0,02053 |
-0,04418 |
|
|
6,0 |
0,137 |
0,352 |
0,825 |
-0,441 |
-0,876 |
-0,837 |
0,040 |
-2,154 |
-0,02192 |
-0,03323 |
|
|
7,0 |
0,133 |
0,320 |
0,778 |
-0,387 |
-0,951 |
-0,955 |
0,033 |
-2,293 |
-0,02180 |
-0,02475 |
|
|
8,0 |
0,130 |
0,291 |
0,729 |
-0,345 |
-1,012 |
-1,065 |
0,028 |
-2,422 |
-0,02076 |
-0,01818 |
|
|
9,0 |
0,128 |
0,267 |
0,682 |
-0,310 |
-1,064 |
-1,168 |
0,023 |
-2,542 |
-0,01922 |
-0,01313 |
|
|
10,0 |
0,126 |
0,246 |
0,636 |
-0,282 |
-1,107 |
-1,264 |
0,020 |
-2,653 |
-0,01745 |
-0,00927 |
|
|
12,0 |
0,124 |
0,212 |
0,551 |
-0,238 |
-1,176 |
-1,436 |
0,014 |
-2,850 |
-0,01388 |
-0,00417 |
|
|
14,0 |
0,123 |
0,185 |
0,476 |
-0,205 |
-1,228 |
-1,585 |
0,011 |
-3,018 |
-0,01078 |
-0,00134 |
|
|
16,0 |
0,123 |
0,164 |
0,412 |
-0,180 |
-1,268 |
-1,713 |
0,008 |
-3,161 |
-0,00829 |
0,00016 |
|
|
18,0 |
0,122 |
0,147 |
0,358 |
-0,161 |
-1,300 |
-1,825 |
0,006 |
-3,286 |
-0,00637 |
0,00092 |
|
|
20,0 |
0,122 |
0,133 |
0,313 |
-0,145 |
-1,326 |
-1,922 |
0,005 |
-3,393 |
-0,00492 |
0,00126 |
|
|
22,0 |
0,121 |
0,122 |
0,275 |
-0,132 |
-1,347 |
-2,008 |
0,004 |
-3,487 |
-0,00382 |
0,00138 |
|
|
24,0 |
0,121 |
0,112 |
0,242 |
-0,121 |
-1,365 |
-2,083 |
0,003 |
-3,570 |
-0,00299 |
0,00137 |
|
|
26,0 |
0,121 |
0,104 |
0,215 |
-0,112 |
-1,381 |
-2,150 |
0,003 |
-3,643 |
-0,00237 |
0,00130 |
|
|
28,0 |
0,121 |
0,097 |
0,191 |
-0,104 |
-1,394 |
-2,209 |
0,002 |
-3,707 |
-0,00189 |
0,00120 |
|
|
30,0 |
0,121 |
0,090 |
0,171 |
-0,097 |
-1,406 |
-2,262 |
0,002 |
-3,765 |
-0,00152 |
0,00109 |
|
|
32,0 |
0,121 |
0,085 |
0,154 |
-0,091 |
-1,416 |
-2,310 |
0,002 |
-3,817 |
-0,00123 |
0,00099 |
|
|
34,0 |
0,121 |
0,080 |
0,139 |
-0,086 |
-1,425 |
-2,353 |
0,001 |
-3,864 |
-0,00101 |
0,00089 |
|
|
36,0 |
0,121 |
0,076 |
0,126 |
-0,081 |
-1,433 |
-2,392 |
0,001 |
-3,906 |
-0,00083 |
0,00080 |
|
|
38,0 |
0,120 |
0,072 |
0,115 |
-0,077 |
-1,440 |
-2,428 |
0,001 |
-3,945 |
-0,00069 |
0,00071 |
|
|
40,0 |
0,120 |
0,068 |
0,105 |
-0,073 |
-1,446 |
-2,460 |
0,001 |
-3,980 |
-0,00058 |
0,00064 |
По значениям вещественной UРАЗ(щ) и мнимой VРАЗ(щ) части частотной передаточной функции разомкнутой скорректированной системы из таблицы 4.3 строим годограф Найквиста (рис. 4.2).
Рис. 4.2. Годограф Найквиста (АФЧХ) скорректированной системы
Вывод: замкнутая скорректированная система устойчива, так как АФЧХ разомкнутой системы не охватывает точку с координатами (-1; j0).
4.5 Исследование скорректированной АСР по логарифмическому критерию
Для построения ЛАХ и ЛФХ системы необходимо разложить передаточную функцию разомкнутой системы на элементарные звенья, амплитуды А() и фазы ц() которых приведены в пункте 4.3.
Затем на плоскости строятся ЛАХ и ЛФХ каждого отдельного звена и методом графического суммирования находятся результирующие ЛАХ и ЛФХ разомкнутой системы [**]:
Зная выражения амплитуд А() и фаз ц() частотных передаточных функции для каждого звена (табл. 4.2), представим их в логарифмическом масштабе (табл. 4.4).
Таблица 4.4 Выражения логарифмических амплитуд L() и фаз ц() передаточных функции звеньев
|
Частотная передаточная функция |
Амплитуда L() |
Фаза ц() |
|
Вычисляем значения логарифмических амплитуд L() и фаз ц() для каждого звена при изменении часты от 0 до значения, при котором результирующая ЛФХ пересекает значение - р. Также вычислим амплитуду LРАЗ() и фазу цРАЗ() передаточной функции разомкнутой системы WРАЗ(jщ). Шаг изменения частоты принимаем произвольный с учетом удобства восприятия графика. Вычисления проводим в MS Excel (табл. 4.5).
Таблица 4.5 Значения логарифмических амплитуд L() и фаз ц() передаточных функции звеньев
|
Частота щ |
Логарифмические амплитуды |
Общая ЛАХ LРАЗ(щ) |
Фазы звеньев |
Общая ЛФХ цРАЗ(щ) |
|||||
|
L1(щ) |
L2(щ) |
L3(щ) |
ц1(щ) |
ц2(щ) |
ц3(щ) |
||||
|
0,100 |
-0,259 |
-5,194 |
-0,001 |
-5,454 |
-0,215 |
-0,020 |
-0,015 |
-0,250 |
|
|
0,500 |
-3,990 |
-5,236 |
-0,013 |
-9,239 |
-0,747 |
-0,100 |
-0,075 |
-0,922 |
|
|
0,900 |
-7,521 |
-5,331 |
-0,042 |
-12,894 |
-0,889 |
-0,178 |
-0,135 |
-1,202 |
|
|
1,300 |
-10,016 |
-5,477 |
-0,088 |
-15,580 |
-0,900 |
-0,254 |
-0,194 |
-1,349 |
|
|
1,700 |
-11,798 |
-5,668 |
-0,149 |
-17,615 |
-0,868 |
-0,328 |
-0,253 |
-1,449 |
|
|
2,100 |
-13,106 |
-5,898 |
-0,226 |
-19,231 |
-0,820 |
-0,398 |
-0,312 |
-1,530 |
|
|
2,900 |
-14,841 |
-6,452 |
-0,425 |
-21,718 |
-0,718 |
-0,526 |
-0,427 |
-1,670 |
|
|
3,300 |
-15,426 |
-6,763 |
-0,546 |
-22,735 |
-0,670 |
-0,583 |
-0,483 |
-1,737 |
|
|
4,100 |
-16,254 |
-7,426 |
-0,826 |
-24,506 |
-0,585 |
-0,687 |
-0,593 |
-1,866 |
|
|
4,500 |
-16,551 |
-7,770 |
-0,984 |
-25,304 |
-0,549 |
-0,733 |
-0,647 |
-1,929 |
|
|
5,300 |
-16,994 |
-8,463 |
-1,332 |
-26,789 |
-0,486 |
-0,815 |
-0,750 |
-2,051 |
|
|
5,700 |
-17,160 |
-8,809 |
-1,520 |
-27,489 |
-0,459 |
-0,851 |
-0,801 |
-2,110 |
|
|
6,100 |
-17,300 |
-9,152 |
-1,716 |
-28,169 |
-0,435 |
-0,884 |
-0,849 |
-2,168 |
|
|
6,900 |
-17,520 |
-9,823 |
-2,132 |
-29,475 |
-0,392 |
-0,944 |
-0,944 |
-2,280 |
|
|
7,300 |
-17,607 |
-10,150 |
-2,350 |
-30,107 |
-0,374 |
-0,970 |
-0,989 |
-2,333 |
|
|
8,100 |
-17,748 |
-10,785 |
-2,801 |
-31,334 |
-0,341 |
-1,018 |
-1,076 |
-2,435 |
|
|
8,500 |
-17,805 |
-11,092 |
-3,034 |
-31,931 |
-0,327 |
-1,039 |
-1,118 |
-2,484 |
|
|
9,300 |
-17,900 |
-11,686 |
-3,510 |
-33,096 |
-0,301 |
-1,077 |
-1,198 |
-2,576 |
|
|
9,700 |
-17,939 |
-11,972 |
-3,753 |
-33,664 |
-0,290 |
-1,095 |
-1,236 |
-2,621 |
|
|
10,100 |
-17,974 |
-12,252 |
-3,998 |
-34,224 |
-0,279 |
-1,111 |
-1,273 |
-2,664 |
|
|
12,100 |
-18,103 |
-13,554 |
-5,248 |
-36,905 |
-0,236 |
-1,179 |
-1,444 |
-2,859 |
|
|
14,100 |
-18,184 |
-14,712 |
-6,508 |
-39,404 |
-0,204 |
-1,230 |
-1,591 |
-3,025 |
|
|
16,100 |
-18,237 |
-15,750 |
-7,754 |
-41,741 |
-0,179 |
-1,270 |
-1,719 |
-3,168 |
|
|
18,100 |
-18,274 |
-16,686 |
-8,971 |
-43,930 |
-0,160 |
-1,301 |
-1,830 |
-3,291 |
По таблице 4.5 строим логарифмические частотные характеристики системы (рис. 4.3).
Рис. 4.3 Логарифмические частотные характеристики разомкнутой системы (ЛАХ и ЛФХ)
По частотным характеристикам (рис. 4.3) графически определим запасы устойчивости системы: ДL = **** дБ и Дц = **** рад.
Вывод: замкнутая скорректированная система устойчива, так как ЛАХ разомкнутой системы пересекает ось абсцисс раньше, чем ЛФХ переходит через значение - р.
Раздел 5. Исследование АСР в среде Simulink
Программный комплекс MATLAB поставляется вместе с пакетом расширения Simulink, предназначенным для моделирования динамических систем, модели которых составляются из отдельных блоков (компонентов) [**].
Для исследования скорректированной системы в среде Simulink составим структурную схему модели АСР по управляющему воздействию (рис. 5.1).
Для получения линейных характеристик системы необходимо задать точки входа и выхода для линеаризации: «Input point» и «Output point».
Для получения переходных и частотных характеристик в пакете Simulink с помощью LTI Viewer необходимо выполнить следующие действия:
1. Выполнить команду меню Tools > Control Design > Linear Analysis, после чего появится окно менеджера Control and Estimation Tools Manager, в котором показываются активные точки, кнопка запуска линеаризации «Linearize Model» и окно выбора типа характеристик.
2. Выбрать нужную характеристику и нажать кнопку «Linearize Model», после чего на экран выводится окно LTI Viewer с построенным графиком.
При получении частотных характеристик необходимо поставить галочку «Open Loop» на активной точке выхода «Output point», так как частотные характеристики строятся только для разомкнутых систем.
На рисунке 5.2 представлена переходная характеристика АСР стабилизации угловой скорости ДПТ с корректирующим звеном, полученная в Simulink.
На рисунках 5.3 и 5.4 представлены логарифмические частотные характеристика и АФЧХ АСР угловой скорости ДПТ с корректирующим звеном.
Таким образом, …………..(вывод).
Рис. 3.3 Логарифмические частотные характеристики разомкнутой системы (ЛАХ и ЛФХ)
По частотным характеристикам (рис. 3.3) графически определим запасы устойчивости системы: ДL = **** дБ и Дц = **** рад.
Вывод: замкнутая система устойчива, так как ЛАХ разомкнутой системы пересекает ось абсцисс раньше, чем ЛФХ переходит через значение - р.
3.2 Корректирующее звено:
Передаточная функция корректирующего звена:
Передаточная функция разомкнутой системы с корректирующим звеном
Где:
Находим такое значение k0, чтобы статическая ошибка соответствовала разрешенной:
3.3 Критерий устойчивости Гурвица:
Составим характеристическое уравнение:
Определитель Гурвица:
Вывод: Система устойчива, так как не имеет отрицательных миноров.
3.4 Критерий устойчивости Михайлова:
|
? |
Ud(?) |
j*Vd(?) |
? |
Ud(?) |
j*Vd(?) |
|
|
0 |
1 |
0 |
11 |
-12,13026 |
-40,1687 |
|
|
0,5 |
0,973845 |
4,7369125 |
12 |
-14,62832 |
-73,8216 |
|
|
1 |
0,89334 |
9,3923 |
13 |
-17,34378 |
-115,3009 |
|
|
2 |
0,56928 |
18,1324 |
14 |
-20,27664 |
-165,2588 |
|
|
3 |
0,02782 |
25,5681 |
15 |
-23,4269 |
-224,3475 |
|
|
4 |
-0,73104 |
31,0472 |
16 |
-26,79456 |
-293,2192 |
|
|
5 |
-1,7073 |
33,9175 |
17 |
-30,37962 |
-372,5261 |
|
|
6 |
-2,90096 |
33,5268 |
18 |
-34,18208 |
-462,9204 |
|
|
7 |
-4,31202 |
29,2229 |
19 |
-38,20194 |
-565,0543 |
|
|
8 |
-5,94048 |
20,3536 |
20 |
-42,4392 |
-679,58 |
|
|
9 |
-7,78634 |
6,2667 |
21 |
-46,89386 |
-807,1497 |
|
|
10 |
-9,8496 |
-13,69 |
22 |
-51,56592 |
-948,4156 |
Список использованной литературы
1. Бесекерский, В. А. Теория автоматического управления [Текст]: учебное пособие / В.А. Бесекерский, Е.П. Попов. - СПб.: Изд-во «Профессия», 2003. - 751 c.
2. Мирошник, И. В. Теория автоматического управления. Линейные системы [Текст]: учебное пособие / И. В. Мирошник. - СПб.: Изд-во «Питер», 2005. - 336 с.
3. Щербаков, В.С. Теория автоматического управления. Линейные непрерывные системы [Текст]: учебное пособие / В.С. Щербаков, И.В. Лазута. - Омск: СибАДИ, 2013. - 142 с.
4. Щербаков, В.С. Основы моделирования систем автоматического регулирования и электротехнических систем в среде MatLab и Simulink [Текст]: учебное пособие / В.С. Щербаков, А.А. Руппель, В.А. Глушец. - Омск: Изд-во СибАДИ, 2003. - 160 с.
5. Щербаков В.С. Теория автоматического управления [Текст]: методические указания по выполнению курсовых работ (для студентов специальностей 220301, 140607) / В.С. Щербаков, Р.Ю. Сухарев. - Омск: СибАДИ, 2011. - 36 с.