По таблице 3.1 строим годограф Михайлова
Рис. 3.1 Годограф Михайлова
Вывод: замкнутая система устойчива, т.к. годограф Михайлова начинается на положительной действительной полуоси и огибает против часовой стрелки начало координат, проходя последовательно 3 квадранта, где 3 - порядок характеристического полинома замкнутой системы.
3.3 Исследование АСР по критерию Найквиста
Критерий Г. Найквиста позволяет по амплитудно-фазовой частотной характеристике (АФЧХ) разомкнутой системы оценить устойчивость замкнутой системы с отрицательной обратной связью.
АФЧХ разомкнутой системы можно построить на комплексной плоскости [+1; j], откладывая вещественную UРАЗ(щ) и мнимую VРАЗ(щ) составляющую частотной передаточной функции разомкнутой системы WРАЗ(jщ), или в полярной системе координат, откладывая угол фазы ц() разомкнутой системы относительно оси +1 и в этом направлении откладывать вектор амплитуды разомкнутой системы А().
Вещественную UРАЗ(щ) и мнимую VРАЗ(щ) функции WРАЗ(jщ) можно определить по АЧХ АРАЗ() и ФЧХ цРАЗ() разомкнутой системы:
;
.
АЧХ передаточной функции разомкнутой системы АРАЗ() равна произведению АЧХ отдельных звеньев, а ФЧХ цРАЗ() - сумме ФЧХ звеньев:
;
.
Найти АЧХ А() и ФЧХ ц() можно по вещественной U(щ) и мнимой V(щ) составляющим частотной передаточной функции W(jщ) звена.
Амплитуда А() и фаза ц() частотной передаточной функции W(jщ):
;
.
Формулировка критерия Найквиста
Если разомкнутая система устойчива, то для устойчивости системы в замкнутом состоянии необходимо и достаточно, чтобы АФЧХ разомкнутой системы на комплексной плоскости [+1; j] при изменении частоты от 0 до не охватывала точку с координатами (-1; j0). Если АФЧХ разомкнутой системы проходит через точку с координатами (-1; j0), то система находится на границе устойчивости.
Для построения АФЧХ разомкнутой системы необходимо перейти к частотной форме записи передаточной функции разомкнутой системы.
Передаточная функция разомкнутой системы:
;
.
Проведем преобразования, выделив общий коэффициент усиления разомкнутой системы kРАЗ:
;
.
;
Передаточная функция разомкнутой системы в частотной форме (p = jщ):
.
Запишем выражения амплитуд А() и фаз ц() частотных передаточных функции для каждого звена (табл. 3.2). Вычисляем значения амплитуд А() и фаз ц() при изменении часты от 0 до значения, при котором UРАЗ(щ) и VРАЗ(щ) станут равны нулю. Также вычислим амплитуду АРАЗ() и фазу цРАЗ() передаточной функции разомкнутой системы, а также вещественную UРАЗ(щ) и мнимую VРАЗ(щ) составляющую частотной передаточной функции разомкнутой системы WРАЗ(jщ). Шаг изменения частоты принимаем произвольный с учетом удобства восприятия графика. Вычисления проводим в MS Excel (табл. 3.3).
Таблица 3.2 Выражения амплитуд А() и фаз ц() передаточных функции звеньев системы
|
Частотная передаточная функция |
АЧХ А() |
ФЧХ ц() |
|
Таблица 3.3 Значения амплитуд А() и фаз ц() передаточных функции звеньев системы
|
Частота щ |
Амплитуды звеньев |
Фазы звеньев |
Амплитуда WРАЗ(jщ) |
Фаза WРАЗ(jщ) |
Действит. часть WРАЗ(jщ) |
Мнимая часть WРАЗ(jщ) |
|||
|
A1(щ) |
A2(щ) |
ц1(щ) |
ц2(щ) |
AРАЗ(щ) |
цРАЗ(щ) |
UРАЗ(щ) |
VРАЗ(щ) |
||
|
0,0 |
8,190 |
1,000 |
0,000 |
0,000 |
8,190 |
0,000 |
8,19000 |
0,00000 |
|
|
1,0 |
7,496 |
0,893 |
-0,415 |
-0,655 |
6,692 |
-1,069 |
3,21843 |
-5,86721 |
|
|
2,0 |
6,148 |
0,679 |
-0,722 |
-1,185 |
4,175 |
-1,906 |
-1,37434 |
-3,94283 |
|
|
4,0 |
4,046 |
0,352 |
-1,054 |
1,292 |
1,426 |
0,238 |
1,38568 |
0,33590 |
|
|
6,0 |
2,901 |
0,197 |
-1,209 |
0,937 |
0,571 |
-0,272 |
0,55031 |
-0,15333 |
|
|
8,0 |
2,238 |
0,122 |
-1,294 |
0,728 |
0,273 |
-0,566 |
0,23028 |
-0,14639 |
|
|
10,0 |
1,815 |
0,082 |
-1,347 |
0,592 |
0,149 |
-0,755 |
0,10824 |
-0,10183 |
|
|
12,0 |
1,524 |
0,058 |
-1,384 |
0,499 |
0,089 |
-0,885 |
0,05640 |
-0,06892 |
|
|
14,0 |
1,312 |
0,044 |
-1,410 |
-2,712 |
0,057 |
-4,121 |
-0,03192 |
0,04758 |
|
|
16,0 |
1,152 |
0,034 |
-1,430 |
-2,764 |
0,039 |
-4,194 |
-0,01931 |
0,03381 |
|
|
18,0 |
1,026 |
0,027 |
-1,445 |
-2,805 |
0,028 |
-4,250 |
-0,01232 |
0,02471 |
|
|
20,0 |
0,925 |
0,022 |
-1,458 |
-2,838 |
0,020 |
-4,296 |
-0,00821 |
0,01854 |
|
|
22,0 |
0,842 |
0,018 |
-1,468 |
-2,865 |
0,015 |
-4,333 |
-0,00567 |
0,01422 |
|
|
24,0 |
0,772 |
0,015 |
-1,476 |
-2,888 |
0,012 |
-4,364 |
-0,00404 |
0,01113 |
|
|
26,0 |
0,713 |
0,013 |
-1,484 |
-2,907 |
0,009 |
-4,391 |
-0,00295 |
0,00886 |
|
|
28,0 |
0,663 |
0,011 |
-1,490 |
-2,924 |
0,007 |
-4,414 |
-0,00221 |
0,00717 |
|
|
30,0 |
0,619 |
0,010 |
-1,495 |
-2,938 |
0,006 |
-4,433 |
-0,00168 |
0,00587 |
|
|
32,0 |
0,580 |
0,009 |
-1,500 |
-2,951 |
0,005 |
-4,451 |
-0,00130 |
0,00487 |
|
|
34,0 |
0,546 |
0,008 |
-1,504 |
-2,962 |
0,004 |
-4,466 |
-0,00103 |
0,00408 |
|
|
36,0 |
0,516 |
0,007 |
-1,508 |
-2,972 |
0,004 |
-4,480 |
-0,00082 |
0,00345 |
|
|
38,0 |
0,489 |
0,006 |
-1,511 |
-2,981 |
0,003 |
-4,492 |
-0,00066 |
0,00295 |
|
|
40,0 |
0,465 |
0,006 |
-1,514 |
-2,989 |
0,003 |
-4,503 |
-0,00054 |
0,00254 |
|
|
42,0 |
0,443 |
0,005 |
-1,517 |
-2,996 |
0,002 |
-4,513 |
-0,00044 |
0,00220 |
|
|
44,0 |
0,422 |
0,005 |
-1,519 |
-3,003 |
0,002 |
-4,522 |
-0,00037 |
0,00192 |
|
|
46,0 |
0,404 |
0,004 |
-1,521 |
-3,009 |
0,002 |
-4,530 |
-0,00031 |
0,00168 |
|
|
48,0 |
0,387 |
0,004 |
-1,523 |
-3,014 |
0,002 |
-4,538 |
-0,00026 |
0,00148 |
|
|
50,0 |
0,372 |
0,004 |
-1,525 |
-3,019 |
0,001 |
-4,545 |
-0,00022 |
0,00131 |
|
|
52,0 |
0,358 |
0,003 |
-1,527 |
-3,024 |
0,001 |
-4,551 |
-0,00019 |
0,00117 |
|
|
54,0 |
0,344 |
0,003 |
-1,529 |
-3,028 |
0,001 |
-4,557 |
-0,00016 |
0,00104 |
|
|
56,0 |
0,332 |
0,003 |
-1,530 |
-3,032 |
0,001 |
-4,563 |
-0,00014 |
0,00094 |
|
|
58,0 |
0,321 |
0,003 |
-1,532 |
-3,036 |
0,001 |
-4,568 |
-0,00012 |
0,00084 |
|
|
60,0 |
0,310 |
0,002 |
-1,533 |
-3,040 |
0,001 |
-4,573 |
-0,00011 |
0,00076 |
По значениям вещественной UРАЗ(щ) и мнимой VРАЗ(щ) части частотной передаточной функции разомкнутой системы из таблицы 3.3 строим в MS Excel годограф Найквиста (рис. 3.2).
Рис. 3.2.Годограф Найквиста (АФЧХ)
Вывод: замкнутая система устойчива, так как АФЧХ разомкнутой системы не охватывает точку с координатами (-1; j0).
3.4 Исследование АСР по логарифмическому критерию
Амплитудная частотная характеристика (АЧХ) - это зависимость отношения амплитуды сигнала на выходе звена к амплитуде на входе A = ymax/xmax от частоты входного сигнала .
Фазовая частотная характеристика (ФЧХ) - зависимость угла сдвига по фазе сигнала на выходе звена от частоты входного сигнала .
Логарифмические частотные характеристики - АЧХ и ФЧХ, представленные в логарифмическом масштабе.
Параметры L и lgщ определяются следующим образом:
где ;
.
Ордината ЛАХ L измеряется в децибелах [дБ].
Децибел - логарифмическая единица уровней, затуханий и усилений.
1 Б = 10 дБ - это увеличение мощности сигнала в 10 раз.
2 Б = 20 дБ - это увеличение мощности сигнала в 100 раз;
3 Б = 30 дБ - это увеличение мощности сигнала в 1000 раз и т.д.
Абсцисса ЛАХ и ЛФХ lgщ измеряется в декадах [дек].
Декада - логарифмическая единица частот, соответствующая изменению частоты щ в 10 раз.
Логарифмические критерии, так же как и критерий Найквиста, позволяют судить об устойчивости замкнутой системы по виду ЛАХ и ЛФХ разомкнутой системы.
Формулировка логарифмического критерия
Если разомкнутая система устойчива, то для устойчивости системы в замкнутом состоянии, ЛАХ разомкнутой системы должна пересечь ось абсцисс раньше, чем ЛФХ, спадая окончательно, перейдёт через значение - р. То есть, на частоте среза щср величина фазы ц должна быть меньше значения | - р |.
Запас устойчивости по амплитуде ДL - это величина допустимого увеличения общего коэффициента усиления разомкнутой системы, при котором замкнутая система окажется на границе устойчивости .
Запас устойчивости по фазе Дц - это величина допустимого увеличения запаздывания по фазе разомкнутой системы на частоте среза щср, при котором замкнутая система окажется на границе устойчивости .
Для построения ЛАХ и ЛФХ системы необходимо разложить передаточную функцию разомкнутой системы на элементарные звенья, амплитуды А() и фазы ц() которых приведены в пункте 3.3. Допускается использовать асимптотические ЛАХ, которые графически представляют собой ломаные прямые линии (табл. 3.4).
Затем на плоскости строятся ЛАХ и ЛФХ каждого отдельного звена и методом графического суммирования находятся результирующие ЛАХ и ЛФХ разомкнутой системы:
;
.
Таблица 3.4 Асимптотические логарифмические частотные характеристики типовых динамических звеньев
Зная выражения амплитуд А() и фаз ц() частотных передаточных функции для каждого звена (табл. 3.2), представим их в логарифмическом масштабе (табл. 3.5).
Таблица 3.5 Выражения логарифмических амплитуд L() и фаз ц() передаточных функции звеньев
|
Частотная передаточная функция |
Амплитуда L() |
Фаза ц() |
|
Вычисляем значения логарифмических амплитуд L() и фаз ц() для каждого звена при изменении часты от 0 до значения, при котором результирующая ЛФХ пересекает значение - р. Также вычислим амплитуду LРАЗ() и фазу цРАЗ() передаточной функции разомкнутой системы WРАЗ(jщ). Шаг изменения частоты принимаем произвольный с учетом удобства восприятия графика. Вычисления проводим в MS Excel (табл. 3.6).
Таблица 3.6 Значения логарифмических амплитуд L() и фаз ц() передаточных функции звеньев
|
Частота щ |
Логарифмические амплитуды |
Общая ЛАХ LРАЗ(щ) |
Фазы звеньев |
Общая ЛФХ цРАЗ(щ) |
|||
|
L1(щ) |
L2(щ) |
ц1(щ) |
ц2(щ) |
||||
|
0,100 |
18,257 |
-5,194 |
-0,001 |
13,062 |
-0,215 |
-0,020 |
|
|
0,500 |
18,060 |
-5,236 |
-0,013 |
12,811 |
-0,747 |
-0,100 |
|
|
0,900 |
17,633 |
-5,331 |
-0,042 |
12,260 |
-0,889 |
-0,178 |
|
|
1,300 |
17,036 |
-5,477 |
-0,088 |
11,472 |
-0,900 |
-0,254 |
|
|
1,700 |
16,336 |
-5,668 |
-0,149 |
10,519 |
-0,868 |
-0,328 |
|
|
2,100 |
15,585 |
-5,898 |
-0,226 |
9,461 |
-0,820 |
-0,398 |
|
|
2,500 |
14,822 |
-6,162 |
-0,318 |
8,341 |
-0,769 |
-0,464 |
|
|
2,900 |
14,069 |
-6,452 |
-0,425 |
7,192 |
-0,718 |
-0,526 |
|
|
3,300 |
13,340 |
-6,763 |
-0,546 |
6,031 |
-0,670 |
-0,583 |
|
|
3,700 |
12,642 |
-7,089 |
-0,680 |
4,873 |
-0,626 |
-0,637 |
|
|
4,100 |
11,977 |
-7,426 |
-0,826 |
3,725 |
-0,585 |
-0,687 |
|
|
4,500 |
11,346 |
-7,770 |
-0,984 |
2,592 |
-0,549 |
-0,733 |
|
|
4,900 |
10,746 |
-8,116 |
-1,153 |
1,477 |
-0,516 |
-0,775 |
|
|
5,300 |
10,178 |
-8,463 |
-1,332 |
0,383 |
-0,486 |
-0,815 |
|
|
5,700 |
9,638 |
-8,809 |
-1,520 |
-0,691 |
-0,459 |
-0,851 |
|
|
6,100 |
9,125 |
-9,152 |
-1,716 |
-1,743 |
-0,435 |
-0,884 |
|
|
6,500 |
8,637 |
-9,490 |
-1,921 |
-2,773 |
-0,412 |
-0,915 |
|
|
6,900 |
8,172 |
-9,823 |
-2,132 |
-3,783 |
-0,392 |
-0,944 |
|
|
7,300 |
7,728 |
-10,150 |
-2,350 |
-4,772 |
-0,374 |
-0,970 |
|
|
7,700 |
7,304 |
-10,471 |
-2,573 |
-5,740 |
-0,357 |
-0,995 |
|
|
8,100 |
6,898 |
-10,785 |
-2,801 |
-6,688 |
-0,341 |
-1,018 |
|
|
8,500 |
6,508 |
-11,092 |
-3,034 |
-7,618 |
-0,327 |
-1,039 |
|
|
8,900 |
6,134 |
-11,392 |
-3,270 |
-8,528 |
-0,313 |
-1,059 |
|
|
9,300 |
5,775 |
-11,686 |
-3,510 |
-9,421 |
-0,301 |
-1,077 |
|
|
9,700 |
5,429 |
-11,972 |
-3,753 |
-10,296 |
-0,290 |
-1,095 |
|
|
10,100 |
5,096 |
-12,252 |
-3,998 |
-11,154 |
-0,279 |
-1,111 |
|
|
12,100 |
3,590 |
-13,554 |
-5,248 |
-15,211 |
-0,236 |
-1,179 |
|
|
14,100 |
2,301 |
-14,712 |
-6,508 |
-18,919 |
-0,204 |
-1,230 |
|
|
16,100 |
1,174 |
-15,750 |
-7,754 |
-22,330 |
-0,179 |
-1,270 |
|
|
18,100 |
0,175 |
-16,686 |
-8,971 |
-25,482 |
-0,160 |
-1,301 |
|
|
20,100 |
-0,722 |
-17,538 |
-10,148 |
-28,409 |
-0,144 |
-1,327 |
|
|
22,100 |
-1,537 |
-18,318 |
-11,283 |
-31,138 |
-0,132 |
-1,348 |
|
|
24,100 |
-2,282 |
-19,037 |
-12,372 |
-33,691 |
-0,121 |
-1,366 |
|
|
26,100 |
-2,969 |
-19,703 |
-13,417 |
-36,089 |
-0,112 |
-1,382 |
|
|
30,000 |
-4,171 |
-20,875 |
-15,331 |
-40,376 |
-0,097 |
-1,406 |
По таблице 3.6 строим логарифмические частотные характеристики системы (рис. 3.3)
Рис. 3.3 Логарифмические частотные характеристики разомкнутой системы (ЛАХ и ЛФХ)
3.5 Проверка статистической ошибки
Статическая ошибка:
Раздел 4. Анализ автоматической системы регулирования, c корректирующим звеном.
4.1 Расчет передаточной функции корректирующего звена
Коррекция динамических свойств линейных систем осуществляется для выполнения рассмотренных ранее требований по точности, устойчивости и качеству переходных процессов. Осуществляется коррекция путем введения в систему специальных корректирующих звеньев с особо подобранной передаточной функцией WКЗ(p).
В данном случае корректирующее звено включается в систему относительно прямой ветви последовательно, образуя новую структуру системы (рис. 4.1).
Рис. 4.1 Блок-схема стабилизации угловой скорости ДПТ с корректирующим звеном
Корректирующие звенья аналоговых линейных электрических систем обычно выполняются в виде пассивных или активных четырехполюсников (АЧ и ПЧ) постоянного тока.
Рассмотрим принципиальную электрическую схему заданного корректирующего звена в виде пассивного четырехполюсника (рис. 4.2).
Рис.4.2.Пассивное корректирующее звено.
Параметры корректирующей цепи:
Передаточная функция корректирующего звена:
Передаточная функция разомкнутой системы с корректирующим звеном:
Где:
Находим такое значение k0, чтобы статическая ошибка соответствовала разрешенной:
4.2 Исследование скорректированной АСР по критерию Гурвица
Для исследования устойчивости замкнутой АСР по критерию устойчивости Гурвица понадобится характеристический полином замкнутой скорректированной системы.
Передаточная функция прямой ветви скорректированной системы:
Передаточная функция обратной связи системы:
Получим эквивалентную передаточную функцию прямой ветви системы, охваченную обратной связью или передаточную функцию замкнутой скорректированной системы:
; (4.9)
Передаточная функция замкнутой скорректированной системы после упрощения:
Выражение, находящееся в знаменателе передаточной функции замкнутой скорректированной системы, является характеристическим полиномом замкнутой скорректированной системы:
Подставим численные значения коэффициентов и преобразуем выражение:
%Вводим в командную строку Matlab:%
>> syms p;
expand((2.5*p+1)*(0.2*p+1)*(0.035*0.15*p^2+0.15*p+1)+0.235*(0.3*p+1)*10*0.5*1.3*0.36)
Поскольку степень полинома замкнутой системы n равна 4, то матрица Гурвица будет иметь размер 4х4.