Влияние распределения корней на вид переходного процесса передаточных функций

ЛАБОРАТОРНАЯ РАБОТА №6

Цель работы: Изучить влияние распределения корней на вид переходного процесса передаточных функций  в программе MATLAB .

Задание:

По заданному преподавателем варианту (см. таблицу 1)  найти корни передаточных функций и построить графики переходных функции.

Таблица 1.

№ п/п	Функция
	W=1/(5s^3+25s^2+12s+2)
	W=1/(12s^3+30s^2+12s+1)
	W=2/(10s^3+20s^2+10s+1)
	W=3/(5s^3+30s^2+7s+1)
	W=1/(7s^3+20s^2+2s+1)
	W=4/(12s^3+10s^2+12s+2)
	W=2/(5s^3+20s^2+7s+1)
	W=3/(12s^3+20s^2+5s+3)
	W=2/(10s^3+30s^2+10s+1)
	W=1/(20s^3+20s^2+12s+1)
	W=2/(15s^3+30s^2+10s+1)
	W=1/(30s^3+20s^2+12s+2)
	W=1/(15s^3+10s^2+25s+1)
	W=1/(5s^3+25s^2+10s+2)
	W=2/(30s^3+10s^2+25s+2)
	W=2/(5s^3+10s^2+7s+1)
	W=3/(12s^3+15s^2+15s+3)
	W=2/(10s^3+30s^2+10s+1)
	W=1/(20s^3+20s^2+20s+1)
	W=2/(15s^3+30s^2+10s+1)
	W=1/(30s^3+20s^2+30s+2)
	W=1/(15s^3+5s^2+10s+1)

Общие положения

Повторим, что переходной функцией называется реакция объекта, описанного передаточной функцией W(s), на единичное ступенчатое воздействие. Для её получения в программе Matlab имеется специальная функция step(w(s)).

Установлено, что переходная функция будет сходящейся (амплитуда убывает со временем), если корни передаточной функции будут находиться в левой части комплексной плоскости, то есть все они отрицательны. На рисунке 1 это a1, a2+b, a2-b.

Если хотя бы один из них положителен, то переходный процесс будет расходящимся. На рисунке 1 это a3, a4+с, a4-с. Если корень находится на мнимой оси, то переходный процесс будет гармоническим с постоянной амплитудой.

Исследуем следующую передаточную функцию:

Момент перехода процесса из сходящегося в расходящийся для передаточной функции третьего порядка W(s) можно определить из следующего условия:

k1*k2 = k0*k3

Если  k1*k2 > k0*k3, то переходный процесс сходящийся, а если k1*k2 < k0*k3, то переходный процесс расходящийся.

Рассмотрим пример:

Напишем программу для вычисления корней и построения переходной функции.

num=[1];                  % числитель передаточной функции

den=[20 4 5 1];       % знаменатель передаточной функции

w=tf(num,den);       % определение передаточной функции

pole(w)                     % вычисление корней передаточной функции

step(w,150,'k'),grid on % построение графика переходной функции

Корни передаточной функции напечатаны в окне команд:

ans =

0.0000 + 0.5000i

0.0000 - 0.5000i

-0.2000         

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

Изменим коэффициент при первой производной с 5 на 7.

num=[1];                  % числитель передаточной функции

den=[20 4 7 1];       % знаменатель передаточной функции

w=tf(num,den);       % определение передаточной функции

pole(w)                     % вычисление корней передаточной функции

step(w,150,'k'),grid on % построение графика переходной функции

Все корни, в этом случае, оказываются отрицательными,

-0.0269 + 0.5843i

-0.0269 - 0.5843i

-0.1461   

Переходный процесс, показанный на рисунке 3, в этом случае, сходящийся, то есть амплитуда колебаний со временем   уменьшается.

Содержание отчета:

3. График переходной  функции по заданию и расчет корней.

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

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

4.  Выводы к каждому графику.

Дата добавления: 2023-10-06