Вывод: Данная система в замкнутом состоянии неустойчива, т.к. не все коэффициенты первого столбца таблицы Рауса положительны.
2.3 Критерий Михайлова.
Для устойчивости АСР n-го порядка необходимо и достаточно, чтобы годограф Михайлова, начав движение от вещественной положительной оси комплексной плоскости, обошел в положительном направлении (против часовой стрелки) последовательно n квадрантов, нигде не обращаясь в нуль.
Выпишем характеристическое уравнение замкнутой АСР:
0,0042P³+0,072P²+1,012P+447,68=0
| лист | 
| 8 | 
Произведем подстановку P=jw:
0,0042(jw)³+0,072(jw)²+1,012jw+447,68=0.
Выделим вещественную и мнимую часть многочлена:
P(w)+jQ(w)=(447,68-0,072w²)+j(1,012w-0,0042w³).
Зададим значения 0<w<¥ и вычислим P(w) и Q(w):
Таблица 3 - Значения для построения АФЧХ
| W | 0 | 4 | 8 | 12 | 16 | 20 | 24 | 28 | 32 | 
| P(w) | 447,68 | 446,53 | 443,07 | 437,31 | 429,24 | 418,25 | 406,21 | 391,23 | 373,95 | 
| Q(w) | 0 | 3,779 | 5,946 | 4,886 | -1,011 | -13,36 | -33,773 | -63,862 | -105,24 | 
| 36 | 40 | 44 | 48 | 52 | 56 | 60 | 64 | 70 | 84 | ... | ∞ | 
| 354,4 | 332,4 | 308 | 281,8 | 252,9 | 221,8 | 188,5 | 153 | 94,8 | -60,3 | - ∞ | |
| -159,5 | -228,3 | -313 | -415,9 | -537,9 | -680,9 | -846,5 | -1000 | -1370 | -2400 | - ∞ | 
По данным таблицы 3 построим годограф Михайлова. (рисунок 1)
Вывод: Данная система в замкнутом состоянии неустойчива, т.к. годограф Михайлова, начав движение от вещественной положительной оси комплексной плоскости, не обошел в положительном направлении последовательно 3 квадранта, нигде не обращаясь в нуль.
2.4 Критерий Найквиста.
Критерий Найквиста позволяет судить об устойчивости замкнутой системы по АФЧХ разомкнутого контура.
Запишем характеристическое уравнение разомкнутой АСР:
446,68
 Wр(p)=    0,0042P³+0,072P²+1,012P+1
Wр(p)=    0,0042P³+0,072P²+1,012P+1
Сделаем замену P=jw и подставим в данное уравнение:
446,68
 Wр(jw)=    -0,0042jw³-0,072w²+1,012jw+1
Wр(jw)=    -0,0042jw³-0,072w²+1,012jw+1
Выделим в знаменатели действительную и мнимую часть и домножим на сопряженное:
| лист | 
| 9 | 
 446,68
                                       
446,68
Wр(p)= (-0,072w²+1)+j(1,012w-0,0042w³)
 446,68·((-0,072w²+1)-j(1,012w-0,0042w³))
                              
446,68·((-0,072w²+1)-j(1,012w-0,0042w³))
Wр(p)= ((-0,072w²+1)+j(1,012w-0,0042w³))·((-0,072w²+1)-

 446,68·((-0,072w²+1)-j(1,012w-0,0042w³))
                                           
446,68·((-0,072w²+1)-j(1,012w-0,0042w³)) 
-j(1,012w-0,0042w³)) = (-0,072w²+1)²+ (1,012w-0,0042w³)²
 32,161w²+446,68
                         
32,161w²+446,68
P(w) = (-0,072w²+1)²+(1,012w-0,0042w³)²
1,876w³-452,04
 Q(w) = (-0,072w²+1)²+(1,012w-0,0042w³)²
Q(w) = (-0,072w²+1)²+(1,012w-0,0042w³)²
Таблица 4 - Значения для построения АФЧХ
| W | 0 | 4 | 8 | 12 | 16 | 20 | 24 | 28 | 32 | 
| P(w) | 446,68 | 67,195 | 51,79 | 45,486 | 28,468 | 13,992 | 6,828 | 3,587 | 2,04 | 
| Q(w) | -452 | -23,2 | 10,523 | 25,005 | 23,733 | 15,309 | 9,176 | 5,697 | 3,731 | 
| 38 | 44 | 48 | 52 | 56 | 60 | 64 | 68 | 72 | 76 ∞ | 
| 0,988 | 0,53 | 0,372 | 0,267 | 0,197 | 0,148 | 0,114 | 0,089 | 0,07 | 0,057 ... 0 | 
| 2,16 | 1,36 | 1,055 | 0,806 | 0,64 | 0,517 | 0,424 | 0,352 | 0,295 | 0,25 ... 0 | 
По данным таблицы 4 построим годограф Найквиста (рисунок 2)
Вывод: Система, будучи нейтральной в разомкнутом состоянии, неустойчива в замкнутом состоянии, т.к. расширенный годограф Найквиста при 0<w<∞ охватывает точку с координатой (-1,j0).
2.5 Логарифмический критерий:
Запишем передаточную функцию разомкнутой системы:
446,68
 Wр(p)= 0,0042P3
+ 0,072P2 + 1,012P + 1
  Wр(p)= 0,0042P3
+ 0,072P2 + 1,012P + 1Определим значения сопрягающих частот:
| лист | 
| 10 | 


 1                 1                1
           1                 1                1
w3= T3 w4= T4 w5= T5
1 1 1


 w3=0,94  = 1,06 1/c   w4=0,072 = 13,88 1/c      
w5=0,062=16,129
1/c
w3=0,94  = 1,06 1/c   w4=0,072 = 13,88 1/c      
w5=0,062=16,129
1/c
Запишем выражение суммарных ЛАЧХ и ЛФЧХ:
Lp(w)=L1(w)+L2(w)+L3(w)+L4(w)+L5(w)+L6(w);
Φp(w)=Ф1(w)+Ф2(w)+Ф3(w)+Φ4(w)+Φ5(w)+Φ6(w);
Вычислим значения Фазочастотных характеристик для каждого звена:
 Q(w)
        
          Q(w)   
Φ(w)=arctg P(w).
Ki Ki·(1-Tijw) Ki -Tijw·Ki



 Wi(jw)=      Tijw+1   =        (Tiw)²+1         =    (Tiw)²+1   +    (Tiw)²+1
Wi(jw)=      Tijw+1   =        (Tiw)²+1         =    (Tiw)²+1   +    (Tiw)²+1
-Tiw·Ki
 Φi=arctg       Ki  
   =arctg(-Tiw).
Φi=arctg       Ki  
   =arctg(-Tiw).
Φ7=0; Φ5=arctg(-0,051w); Φ6=-Π/2;
Вычислим таблицу для Φ3:
Таблица 5 – Значения ВЧХ для первого периодического звена
| W | 0 | 15 | 30 | 60 | 105 | 135 | 150 | 175 | 230 | ... | ∞ | 
| Φ3 | 0 | -0,65 | -0,99 | -1,25 | -1,38 | -1,42 | -1,44 | -1,46 | -1,5 | ... | -∞ | 
Вычислим таблицу для Φ4:
Таблица 6 – Значения ВЧХ для второго периодического звена
| W | 0 | 15 | 30 | 60 | 105 | 135 | 150 | 175 | 230 | ... | ∞ | 
| Φ4 | 0 | -0,31 | -0,56 | -1,89 | -1,14 | -1,23 | -1,26 | -1,31 | -1,37 | ... | -∞ | 
Вычислим таблицу для Φ5:
| лист | 
| 11 | 
Таблица 7 – Значения ВЧХ для третьего периодического звена
| W | 0 | 15 | 30 | 60 | 105 | 135 | 150 | 175 | 230 | ... | ∞ | 
| Φ5 | 0 | -0,17 | -0,32 | -0,58 | -0,86 | -0,98 | -1,03 | -1,09 | -1,19 | ... | -∞ | 
Совместим графики суммарных ЛФЧХ и ЛАЧХ. (рисунок 3)
Вывод: Система находится на грани устойчивости.
| лист | 
| 11 | 
Уважаемый посетитель!
Чтобы распечатать файл, скачайте его (в формате Word).
Ссылка на скачивание - внизу страницы.