Расчёт структурной схемы

Дана структурная схема:

Где:

W₁ = 10;

W₅ = K (1+10p)

W₆=10 / (1+2*10*0.2*p+10²p²)

1. Получить передаточную функцию разомкнутой системы W (p)

Вывод передаточной функции производится вручную любым из методов алгебраических и структурных преобразований блок — схемы.

Перенесём средний сумматор против хода сигнала, преобразуем при этом схема станет:

Заменим звено с единичной ООС на эквивалентное:

Перенесём правый сумматор против хода сигнала, переставим местами сумматоры и получим звено суммы и звено с отрицательной обратной связью их эквивалентная передаточная функция:

2. Исследовать устойчивость разомкнутой системы от буквенного параметра методами Гурвица и Михайлова Метод Михайлова:

Запишем характеристический полином системы:

D (p)=11000kp³+ (100+1540k)p²+ p (4+154k)+11k+1

D (p)=-11000jw³— (100+1540k)w²+ jw (4+154k)+11k+1

U (w)=-(100+1540)w²+11k+1

V (w)= -11000jw³ + jw (4+154k)

Для того, чтобы система находилась на границе устойчивости, необходимо чтобы:

Корень второго уравнения w=0 отбрасываем, т.к. для нахождения системы на границе устойчивости годограф Михайлова должен пройти через начало координат при w= 0.

Тогда из второго уравнения определяем Подставим в первое и получим

1452k² + 132k+5>0

тогда:

Метод Гурвица:

Запишем характеристический полином системы:

D (p)=11000kp³+ (100+1540k)p²+ p (4+154k)+11k+1

В общем виде

D (p) =a₃p³+a₂p²+a₁p+a₀

Так как система имеет третий порядок, то она будет находиться на границе устойчивости при равенстве нулю выражения:

a₁a₂-a₀a₃= (4+154k)*(100+1540k) -11000k*(11k+1)=0

или 1452k² + 132k+5>0

что одинаково с выше полученным уравнением,

3. Получить передаточную функцию W (p) системы, замкнутой единичной отрицательной обратной связью

4. Исследовать устойчивость замкнутой системы от буквенного параметра методам Гурвица. Получить области устойчивых и неустойчивых значений параметра в классе вещественных чисел Метод Михайлова:

Запишем характеристический полином системы:

Для того, чтобы система находилась на границе устойчивости, необходимо чтобы:

w² =

512k² + 1137k+5>0

Метод Гурвица:

Запишем характеристический полином системы:

В общем виде

D (p) =a₃p³+a₂p²+a₁p+a₀

Так как система имеет третий порядок, то она будет находиться на границе устойчивости при равенстве нулю выражения:

a₁a₂-a₀a₃= (4+1014k)*(100+140k) -1000k*(101k+11)=0

512k² + 1137k+5>0

5. Сформировать набор значений параметра, включающий все граничные и по одному из каждого интервала устойчивости и неустойчивости замкнутой системы

k₁ = -2.2163, k₂= - 0,0044, k3=1, k₄ = -10 k₅ = -1

6. Для каждого значения параметра из набора построить частотные характеристики, необходимые для исследования зависимости устойчивости замкнутой системы от параметра по критериям Найквиста и Михайлова В общем виде

D (p) =a₃p³+a₂p²+a₁p+a₀

Годограф Михайлова построим по формулам c помощью пакета MAPLE:

Из графика видно, что гадограф Михайлова, начавшись с положительной действительной оси обходит последовательно 3 квадранта против часовой стрелки, проходя через ноль, следовательно замкнутая система находится на границе устойчивости Проведем анализ при k₂= -0,0044 по критерию Найквиста с помощью пакета MatLab:

k1=tf ([44 6,16 -3,784 9,604],[-48,4 93,224 3,3224 0,9516])

subplot (121)

nyquist (k1,'b')

Из рисунка видно, что АФХ системы проходит через точку (-1;j0), следовательно, замкнутая система на границе устойчивости.

Теперь рассмотрим точку

Метод Михайлова:

Из графика видно, что годограф Михайлова, начавшись с положительной действительной оси обходит последовательно 3 квадранта против часовой стрелки, следовательно, замкнутая система устойчива.

Проведем анализ k₄ = -10 по критерию Найквиста с помощью пакета MatLab:

i1=tf ([100 000 14 000 -8600 -890],[-110 000 -15 300 -1536 -109])

subplot (211)

pzmap (i1,'b')

subplot (212)

nyquist (i1,'b')

Из расположения корней на комплексной плоскости видно, что система не имеет корней с положительной вещественной частью, а АФХ системы не охватывает точку (-1;j0), следовательно, замкнутая система устойчива.

Исследуем точку

Метод Михайлова:

Из графика видно, что годограф Михайлова, начавшись с положительной действительной оси не обходит последовательно 3 квадранта против часовой стрелки, следовательно, замкнутая система неустойчива.

7. Получить оценки качества временных характеристик разомкнутой системы

i1=tf ([-22 163 -3102,82 1906,018 189,467],[24 379,3 3313,102 337,3102 23,3793])

subplot (211)

step (i1,'b')

subplot (212)

pzmap (i1,'b')

i1=tf ([99.89 9.989 99.89 20],[99.89 11.989 100.9 2])

subplot (211)

step (i1,'b')

subplot (212)

pzmap (i1,'b')

Как видно, процесс имеет экспоненциальный характер.

i1=tf ([0.005 0.0005 0.005 20],[0.005 2.0005 0.205 2])

subplot (211)

step (i1,'b')

subplot (212)

pzmap (i1,'b')

i1=tf ([50 5 50 20],[50 7 50.2 2])

subplot (211)

step (i1,'b')

subplot (212)

pzmap (i1,'b')

Как видно, процесс имеет экспоненциальный характер.

i1=tf ([50 5 50 20],[150 17 150.2 2])

subplot (211)

step (i1,'b')

subplot (212)

pzmap (i1,'b')

i2=impulse (i1)

Как видно, процесс имеет экспоненциальный характер.