Авсиевич_КР

		WЖ(s) = WР(s)*WПКУ(s),
где WЖ(s) – желаемая передаточная функция; WР(s) – располагаемая передаточная
функция; WПКУ(s) – передаточная функция корректирующего устройства.
Перейдем к ЛАХ
	LЖ(s) = LР(s) + LПКУ(s)
	LПКУ(s) = LЖ(s) – LР(s)
Для построения ЛАХ корректирующего устройства необходимо взять разность
ЛАХ желаемой и располагаемой систем (рис. 17.).
100	I	II		III
80
60
40
20
0						Lp(s)
						Lж(s)
-200,1	1		10	100	1000
-40		1	T
-60
-80				2
-100
		Рис. 17. Графики располагаемой и желаемой ЛАХ.

Производим вычитание:

а). Участок графика до частоты I – пройдет с наклоном (-20) – (-20) = 0 дб/дек. б). Участок графика от I до II – с наклоном (-60) – (-20) = -40 дб/дек.

в). Участок графика от II до 1 – с наклоном (-20) – (-20) = 0 дб/дек.

г). Участок графика от 1 до Т – с наклоном (-20) – (-40) = 20 дб/дек. е). Участок графика от Т до III – с наклоном (-20) – (-60) = 40 дб/дек. ж). Участок графика от III до 2 – с наклоном (-80) – (-60) = -20 дб/дек.

з). Участок графика от 2 – с наклоном (-80) – (-80) = 0 дб/дек.

Получаем график переходного процесса (рис. 18.).

31

80
70
60
50
40
30
20
10
0
0,1	1	10	100	1000

Рис. 18. График корректирующей ЛАХ с параметрами

I = 0.5 рад/с	TI = 2 c
II = 3.16 рад/с	TII = 0.32 c
1 = 5.13 рад/с	T1 = 0.195 c
T = 20 рад/с	TT = 0.05 c
III = 100 рад/с	TIII = 0.01 c
2 = 200 рад/с	T2 = 0.005 c

= 0,009 с

Определяем передаточную функцию корректирующего устройства:

		(S)	T II S 1 2 * T S 1 * T S 1 * T S 1
	WПКУ		1	T	2
			T I S 1 2	* T III	S 1 3

Задание 2. Исследование качества переходных процессов САУ

Для САУ в соответствии с вариантом задания построить переходную и импульсную характеристику замкнутой системы и определить следующие параметры качества:

1)величину перерегулирования;

2)время переходного процесса;

3)статическую ошибку;

4)время регулирования;

5)коэффициент ошибки по положению;

6)коэффициент ошибки по скорости;

7)коэффициент ошибки по ускорению;

8)интегральные оценки качества.

Теоретические сведения

32

L (t) W ( p) .

Графическое представление переходных и импульсных функций называют временными характеристиками. Временные характеристики представляют процессы, происходящие в динамическом и статическом режимах. Переходной функцией h(t) называют функцию, описывающую сигнал на выходе при условии, что на вход подано единичное ступенчатое воздействие при нулевых начальных условиях. График переходной функции, представляющий собой зависимость функции h(t) от времени t, называют переходной характеристикой. В том случае, если амплитуда единичного ступенчатого воздействия отлична от единицы, получают разновидность переходной характеристики, которая называется кривой разгона.

Импульсной дикцией или весовой функцией (t) называют функцию, описывающую реакцию на единичное импульсное воздействие при нулевых начальных условиях. График зависимости функции (t) от времени называют импульсной переходной (импульсной характеристикой).

Аналитическое определение переходных функций и характеристик основано на следующих положениях. Если задана передаточная функция системы или составной части W(p) и известен входной сигнал x(t), то выходной сигнал y(t) определяется следующим соотношением:

L y(t) L x(t) W ( p) .

Таким образом, изображение выходного сигнала L y(t) представляет собой произведение передаточной функции на изображение входного сигнала L x(t) . Сигнал y(t) в явном виде получился после перехода от изображения L y(t) к оригиналу y(t). Для большинства случаев линейных систем и составных элементов разработаны таблицы, позволяющие производить переход от изображений к оригиналу и обратно.

Так как изображение единичного ступенчатого воздействия равно 1 p , то

изображение переходной функции определяется соотношением:

L h(t) W ( p) . p

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

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

Таким образом, передаточная функция является изображением импульсной

	функции.
	Так как	(t)	d1(t)	,	то между импульсной и переходной функциями существует
			dt
						dh(t)
	следующая зависимость:						(t)
						dt

Методы построения переходной функции

Из выше сказанного известно, что изображение переходной функции, имеет вид

H ( p) W ( p) 1p .

33

Оригинал переходной функции может быть получен использованием точных и приближенных методов. Будем рассматривать только точные методы, связанные с применением обратного преобразования Лапласа.

Метод построения переходных функций с использованием таблиц преобразования Лапласа

Построение переходной функции с использованием таблиц преобразования Лапласа начинается с представления изображения переходной функции в виде произведения передаточных функций типовых звеньев

B( p)

n

B ( p)

H ( p)

i

.

pD( p)

pD

( p)

i 1

i

Затем полученное выражение преобразуется в сумму передаточных функций с

неопределенными коэффициентами

A

m

Bi

l

Ck p Dk

H ( p)

.

(7)

p

p

i

p 2 b p c

i 1

k 1

k

k

В этом выражении A, Bi, Ck, Dk – неопределенные коэффициенты, i – вещественные корни уравнения D(p)=0, а выражения p2+bkp+ck – соответствуют комплексно – сопряженным корням характеристического уравнения. Выражение (7) необходимо привести к общему знаменателю и числитель полученного выражения приравнять числителю изображения исходной переходной функции B(р). Приравнивая члены при одинаковых степенях оператора p в левой и правой частях, получим систему уравнений относительно неопределенных коэффициентов. Вычислив значения этих коэффициентов, обратимся к таблицам преобразования Лапласа и определим оригиналы слагаемых формулы (7). Отметим, что оригиналы, соответствующие первым двум членам суммы (7), определяются непосредственно из таблиц преобразований Лапласа (Приложение 1). Для определения оригиналов для третьего члена этой суммы может потребовать некоторых преобразований.

Метод построения переходных функций с использованием формулы обратного преобразования Лапласа

Оригинал переходной функции может быть получен использованием формулы обратного преобразования Лапласа, так называемой формулы разложения, которая в общем случае имеет вид

n		1		n	B( p )	p t	(ni 1)
h(t)			lim ( p pi ) i		i	e i	,
	(ni				pi D( pi )
i 1		1)! p pi

где pi – корни уравнения pD(p)=0, ni- кратность корней. В случае только простых корней, когда среди них имеется m вещественных корней и l пар комплексно – сопряженных корней, формула разложения принимает вид

B(0)

m

B( p )

l

h(t)

i

2 Ak e k t cos(k t k ) .

D(0)

.

i 1

p

D( p

)

k 1

i

i

В этом выражении k

и

k

–

вещественная и мнимая части комплексно –

сопряженных корней, а

.

d

D( p) . Амплитуда и фаза колебательных составляющих

D( p)

dp

определяются следующим образом:

34

		B( pi )			U k		jVk ,
			.
		pi D( pi )
							arctg	Vk	.
A	U 2		V 2		;	k
k		k	k				U k

При вычислениях по приведенным выше формулам, в первую подставляется только один из пары комплексно – сопряженных корней. При вычислении фазового сдвига необходимо учитывать квадрант, в котором находится вектор Akej k.

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

	Wз	( p)					k kc (T1 p 1)							,
			d	0	p 4	d	1	p3 d	2	p2 d	3	p d	4

где все коэффициенты заданы или вычислены ранее. Корни характеристического уравнения равны:

1 0,49 ; 2 4,608 ;

3,4 4,538 j41,37 j .

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

H ( p) Wз ( p)

1

k kc (T1 p 1)

p

d

0

( p

)( p

2

)( p 2

bp c)

1

,

A

B1

B2

C3 p D3

p

p

p

2

p 2

bp c

1

где b=2 , c= 2+ 2.

Приведем полученное выражение к общему знаменателю и приравняем числитель этого выражения к числителю исходного изображения переходной функции. Приравняв члены при одинаковых степенях оператора p в правой и левой частях, получим систему линейных уравнений относительно неопределенных коэффициентов. Из этой системы сразу определяется А=1, после чего решается система 4-го порядка.

B1 B2 B3 1

13,684B1 9,566B2 5,098C3 D3 14,174

1773,892B1 1736,517B2 2,258C3 5,098D3 1780,5977981,379B1 848,714B2 2,258D3 1807,108

Решение этой системы дает:

B1 0,1318 ; B2 0,8789 ;

C3 0,01076 ; D3 4,02 .

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

C3 p D3

C3 p D3

p D3 / C3

p

D3 / C3

C3

C3

C3

p 2 bp c

( p )2 2

( p )2 2

( p )2 2

( p )2 2

p

C3 (D3 / C3 )

.

C3

( p )2 2

( p )2 2

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

35

h(t) 1 0,131e 0,49t 0,8789e 4,608t 0,0107e 4,538t cos(41,37t) 0,098e 4,538t sin(41,37t) .

2) Построение переходной функции с использованием формулы разложения. Вначале определим составляющие процесса, соответствующие вещественным

корням.

	B( p1 )	0,1319	;	B( p2 )	0,8788 .
	.			.
	p1 D( p1 )			p2 D( p2 )

Для колебательной составляющей получим

B( p3 )	0,0147 j0,0472 ,	A3 0,0494 .
.

p3 D( p3 )

Так как вектор этой составляющей находится во втором квадранте, то

3 arctg		0,0472	1,8725 .
	0,0147
Следовательно
h(t) 1 0,1319e 0,49t	0,0878e 4,608t 0,0988e 4,538t cos(41,37t 1,8725) .
Переходные процессы, полученные различными				способами, совпадают с

точностью до арифметических вычислений. Кривая переходной функции показана на рис.15. Переходный процесс практически монотонный. Колебательная составляющая фактически никак себя не проявляет ввиду крайне малой амплитуды. Перерегулирование отсутствует: =0. Время регулирования, определенное при =0,05, приближенно равно 2,2 с., что для системы автоматического регулирования большинстве случаев является вполне приемлемым.

Рис.15. Переходная функция системы регулирования

Показатели качества САУ

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

36

Рис.16. Переходная функция и показатели качества

Используются следующие прямые показатели качества:

1)величина перерегулирования ,

xmax xуст 100% ;

xуст

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

xуст x( ) lim pX ( p) ;

p 0

2) время переходного процесса или время регулирования tp – наименьшее значение времени, после которого имеет место неравенство

x(t) x( ) ,

где - заданная величина, обычно лежащая в пределах =0,02÷0,05;

3) статическая ошибка сm – величина отклонения установившегося значения

регулируемой величины x( ) от требуемого значения N

ст N x( )

или ст lim pE( p) , где E(p) – изображение ошибки;

p 0

4) время регулирования tр – промежуток времени, по истечении которого регулируемая величина первый раз достигает установившегося значения.

Для определения качества системы могут использоваться и другие показатели, соответствующие решаемой задаче, например, число колебаний регулируемой величины за время регулирования, частота и период колебаний и т.д.

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

Коэффициенты ошибок

Точность САУ в установившемся режиме, при относительно медленно изменяющихся воздействиях, может быть оценена с помощью коэффициентов ошибок. Изображение ошибки определяется выражением

E( p) Ф ( p)G( p) ,

где Ф ( p) - передаточная функция по ошибке.

37

Разложим передаточную функцию системы по ошибке в степенной ряд в окрестности точки p=0. Отметим, что при p 0, t и именно поэтому мы говорим о точности в установившемся режиме.

Обозначим: Ki Ci

i!

Ф ( p) C0 C1 p C2!2 p2 ... Ci!i pi ...

и получим

Ф ( p) K			0	K p K			p2 ... K		pi ...,	(8)
					1	2		i
E( p) (K	0	K p K				p2 ... K		pi ...)G( p) .
				1	2		i

Учитывая, что оператор p, умноженный на изображение самой величины, является символом дифференцирования, можно для оригиналов записать

	.		g(2) (t) ... K		g(i) (t) ... .
(t) K	g(t) K g(t) K					(9)
0	1	2		i
Выражение (9) определяет зависимость		ошибки		регулирования от		различных

составляющих входного воздействия, коэффициенты Ki получили название коэффициентов ошибок:

-K0 - коэффициент ошибки по положению;

-K1- коэффициент ошибки по скорости;

-K2 – коэффициент ошибки по ускорению и т.д. Из (8) следует, что

	1		i 1
Ki ( p)		Ф ( p) Km .
	i
	p		m 0

Численные значения коэффициентов ошибок определяются из этого выражения при p 0.

		1			i 1
Ki	( p) lim			Ф ( p) Km .
		p	i
	p 0				m 0

Очевидно, что К0=Ф (0).

Входное воздействие можно представить в виде степенного ряда

g(t) g0 g1t g22 t2 ... gmm! tm ,

где g0 – постоянная величина, характеризующая начальное значение, g1=const – скорость изменения входного воздействия, g2 =const – ускорение и т.д.

Тогда

G( p)	g0		g1		g2	...	gm	.
					2 p3		m! pm 1
	p p2

Пусть передаточная функция разомкнутой системы имеет вид

W ( p) B( p) , p A( p)

где - порядок астатизма системы. Для передаточной функции замкнутой системы по ошибке получим

Ф ( p)		1
		W ( p)
1

Изображение ошибки запишется в виде

	p A( p)	g		0
E( p)
			p
	p A( p) B( p)

		p A( p)						.
p A( p) B( p)
	g			g	2				g	m
	1						...					.
		2				3					m 1
	p			2 p					m! p

38

Отсюда следует, что если порядок астатизма больше порядка старшей производной воздействия, т.е. >m, то ошибка в установившемся режиме будет равна нулю. Если =m, то установившаяся ошибка будет равна постоянной величине, называемой статической ошибкой. И если <m, то при t и . В отношении коэффициентов ошибок последнее выражение позволяет сделать следующие выводы.

1). Если система статическая, т.е. =0, то существуют все составляющие ошибки и все коэффициенты ошибок не равны нулю, т.к. К0 = Ф (0) 0.

2).Система с астатизмом 1-го порядка, =1, не имеет ошибки по положению и

К0=0.

3).Система с астатизмом 2-го порядка, =2, не имеет ошибок по положению и по скорости и К0 =0, К1=0.

Этот список можно продолжить. Таким образом, повышение порядка астатизма повышает точность системы в установившемся режиме. Но повышение порядка астатизма снижает запасы устойчивости, т.к. введение интегрирующих звеньев увеличивает фазовое запаздывание (снижает частоту ). Поэтому на практике порядок астатизма выше второго не применяют, а чаще всего ограничиваются астатизмом первого порядка, используя для повышения точности другие способы.

Интегральные оценки качества

Интегральные оценки характеризуют качество протекания переходных процессов. Наибольшее распространение получили две интегральные оценки

J0	x(t) x( ) 2 dt 2 (t)dt ,	(10)
0	0

J1 x(t)

0

				dx
2		2				2				2
x( )			(		)		dt				(t)
				dt				0

	2	d (t) 2
					dt .	(11)
			dt

Интеграл (10) определяет площадь под кривой квадрата динамической ошибки. Чем меньше этот интеграл, тем быстрее затухает переходный процесс и, следовательно, интеграл J0 служит мерой быстродействия системы. В ряде случаев система, удовлетворяющая условию минимума J0, имеет значительную колебательность переходного процесса. Для уменьшения колебательности можно попробовать воспользоваться оценкой J1 (11). Представим этот интеграл в виде

J1 (t)

0

	d (t)		2
				dt 2 (t)
	dt			0

d (t)				d (t)		2
							2
	dt		(t)			dt		(t)		.
dt									0
		0		dt

Последний член в полученном выражении является постоянной величиной и, если считать, что при t ошибка (t) 0, то он равен 2(0). Минимальное значение интеграл J1 будет иметь, если подынтегральное выражение будет равно нулю, т.е.

(t) d (t) 0 . dt

Решение этого дифференциального уравнения будет

(t) (0)e t .

При подаче на вход системы единичного ступенчатого воздействия начальное значение ошибки (0)=1, и можно рекомендовать следующую методику выбора величины постоянной времени .

39

1) выберем из каких либо соображений время регулирования tp и величину , по

уровню которой выбирается это время, т.е. (tp ) e

t p

;

2)

определим

логарифм

натуральный

от

полученного

выражения

t p

ln , и получим

t p

.

ln

Недостатками интегральных оценок являются невозможность получения прямых показателей качества и высокая сложность вычислительных процедур. Достоинство – это возможность выразить интегральные оценки как функции параметров системы и, воспользовавшись известными методами поиска экстремума, определить значения этих параметров, дающие минимум избранной оценке. Именно это и послужило развитию аналитических методов синтеза систем автоматического управления, основанных на минимизации квадратичных интегральных оценок.

Задание 3. Моделирование САУ

В третьем задании необходимо на базе лабораторных работ выполнить следующие:

1.В соответствии с вариантом в задании 1 смоделировать в пакете MathLab Simulink заданную САУ.

2.Найти частотные характеристики и устойчивость САУ в соответствии с вариантом из таблицей 1.

3.Построить график переходной и импульсной переходной функции.

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

Теоретические сведения

Исследование линейных динамических моделей с использованием ППП

Control System Toolbox системы MatLab 6.x

Для выполнения лабораторной работы используется пакет прикладных программ (ППП) Control System Toolbox 5 [3] системы инженерных расчетов MatLab. ППП

предназначен для работы с линейными или LTI-моделями (Linear Time Invariant Models) систем управления.

В Control System Toolbox имеется тип данных, определяющих динамическую систему в виде комплексной передаточной функции.

Синтаксис команды, создающий LTI-систему c одним входом и одним выходом, в виде передаточной функции:

TF([bm, …, b1, b0], [an, …, a1, a0]),

bm, …, b1 - значения коэффициентов числителя полинома, an, …, a1 - значения коэффициентов знаменателя полинома.

Для выполнения работы могут применяться команды, приведенные в таблице 3.

40