Исследование устойчивости разомкнутой и замкнутой системы. Необходимое и достаточное условие устойчивости

Страницы работы

Содержание работы

1. Исследование устойчивости разомкнутой и замкнутой системы.

1) Устойчивость разомкнутой системы определяется корнями D(p)  ( Q(p) ):

p = 0                          – на границе устойчивости

p = -1/T0 = -0,2 < 0

Необходимое и достаточное условие устойчивости состоит в том, чтобы вещественные части всех корней характеристического уравнения были бы отрицательными.

2) Устойчивость замкнутой системы определим по знакам характеристического уравнения:

 - для системы 2 порядка положительность коэффициентов является необходимым и достаточным условием устойчивости.

Замкнутая система устойчива.

2. Построение границ областей устойчивости в плоскости параметров .

Выделение областей устойчивости в плоскости двух параметров:

Обозначим: А – T0, В – k1

Листинг 2.1. Выполнение построения  областей устойчивости в плоскости параметров системы в Matlab’е.

>> q_string = '[y, 1 + x*120, 120]';

>> range_x = [-0.1:0.01:0.25];

>> range_y = [-0.1:0.01:1.2];

>> dcom( q_string, range_x, range_y );

Рис 2.1. Области устойчивости в плоскости двух параметров системы

Система устойчива только в области параметров D(0).

( область D(1) – T0 < 0, область D(2) – k1 < -1/120 )

3. Построение линий равной степени устойчивости в плоскости тех же параметров.

     

 

Листинг 3.1. Файл с исходными данными (data.m), для выполнения функции Rtanalti в Matlab’е.

function [R, Q, area] = data( a )

%DATA - задание исходных данных

%~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

%

% Файл создается пользователем

%

% a(a1, a2)  - вектор варьируемых параметров системы

% R, Q  - числитель и знаменатель передаточной функции

% Задание области вариации параметров a1,a2. Формат: [минимум, шаг, максимум]

area = [ 0 0.1 6; 0 0.01 1 ];

% Задание коэффициентов числителя и знаменателя п.ф., зависящих от параметров

% Записываются по убыванию степеней

R = [120*a(1) 120];                                 % Числитель передаточной функции

Q = [a(2) 1+120*a(1) 120];    % Характеристический полином

a(1) и a(2) – варьируемые параметры (kи T0 оответственно ).

Листинг 3.2. Выполнение построения  линии равной степени устойчивости в плоскости тех же параметров в Matlab’е.

>> rtanalti('data')

Рис 3.1. Линии равной степени устойчивости

Рис 3.2. Линии равной степени колебательности

4. Построение переходного процесса для единичного ступенчатого воздействия на входе при нулевых начальных условиях для оптимальных значений параметров. Определение показателей качества переходного процесса.

Исследование чувствительности системы к изменению этих параметров.

Амплитуда

t,c

Рис.4.1. Переходная характеристика на выходе инерционного пропорционально дифференцирующего звена при ступенчатом входном воздействии (для первой точки – t1).

Текущие значения параметров системы:

a1 = 0.010141

a2 = 0.010125

Значение степени устойчивости: 97.924376

Значение степени колебательности: 0.000000

Величина перерегулирования: 0 %

Время переходного процесса: 0.025 с

Амплитуда

t,c

Рис.4.2. Переходная характеристика на выходе инерционного пропорционально дифференцирующего звена при ступенчатом входном воздействии (для второй точки – t2).

Текущие значения параметров системы:

a1 = 0.008446

a2 = 0.004392

Значение степени устойчивости: 70.408383

Значение степени колебательности: 0.000000

Величина перерегулирования: 0 %

Время переходного процесса: 0.03 с

Амплитуда

 t,c

Рис.4.3. Переходная характеристика на выходе инерционного пропорционально дифференцирующего звена при ступенчатом входном воздействии (для третьей точки – t3).

Текущие значения параметров системы:

a1 = 0.008636

a2 = 0.013716

Значение степени устойчивости: 74.231975

Значение степени колебательности: 0.766607

Величина перерегулирования: 4 %

Время переходного процесса: 0.025 с

Амплитуда

 t,c

Рис.4.4. Переходная характеристика на выходе инерционного пропорционально дифференцирующего звена при ступенчатом входном воздействии (для четвертой точки – t4).

Текущие значения параметров системы:

a1 = 0.012056

a2 = 0.013716

Значение степени устойчивости: 89.189164

Значение степени колебательности: 0.315948

Величина перерегулирования: 1 %

Время переходного процесса: 0.025 с

Амплитуда

 t,c

Рис.4.5. Переходная характеристика на выходе инерционного пропорционально дифференцирующего звена при ступенчатом входном воздействии (для пятой точки – t5).

Похожие материалы

Информация о работе