Добавил:

student_tipo Опубликованный материал нарушает ваши авторские права? Сообщите нам.

Вуз:

Балаковский институт техники, технологии и управления

Предмет:

Идентификация и диагностика систем

Файл:

лекции / Идентификация объектов управления - УП - Семёнов-Артамонов-Брюхачев - 2003 - 215

.pdf

Скачиваний:

297

Добавлен:

17.02.2014

Размер:

1.61 Mб

Скачать

☆

<<< < Предыдущая 1 23 / 223 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 > Следующая >>>

Знание собственных значений и векторов матрицы системы позволяет осуществлять линейные преобразования (2.6) в пространстве состояний, придавая различный вид этой матрице. В задачах управления наиболее часто используются следующие канонические виды матрицы А [30].

1. Диагональная каноническая форма

		λ	0 ..	0
		1	λ2 ..	0
Aд =		0	λ2 ..	0	,
Aд =	.. .. ..			..	,
		0	0 ..
		0	0 ..	λn
где λi - собственные различные значения матрицы А.
Диагональной канонической		форме		соответствует
структурная схема рис. 2.3.
U(t)

(2.18)

b1	b2	bn-1	bn
x1	x2	xn-1	xn
1	1	1	1
p	p	p	p
λ1	λ2	λn−1	λn
c1	c2	cn-1	cn

y(t)

Рис. 2.3 2. Каноническая форма управляемости

Aу = −

1	..	0
0	..	0		(2.19)
..	..	..	. ,	(2.19)
..	..	..
− d1 ..		− dn−1
− d1 ..		− dn−1

где di – коэффициенты характеристического уравнения матрицы (2.17).

Структурная схема для канонической формы управляемости показана на рис.2.4.

Y(t)

U(t)

1 p

		+
			dn-1
				xn-1
				xn-1
1					1
	p					p

1 p

cn	cn-1	c2	c1
+	+		+

Рис. 2.4.

3. Каноническая форма наблюдаемости

		0	0 ..
Aн =		1	0 ..
Aн =	.. .. ..
	0		... 1

− b	0
	0
− b1			(2.20)
..		. .	(2.20)
..
− bn−1
− bn−1

Этой канонической форме соответствует структура рис. 2.5

U(t)

	b2
x1	+
1
p

bn-1

	bn
xn-1	+

an-1

xn	y(t)
1
p

Рис. 2.5.

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

		λ	0	..	0
		1	λ2
Aд =		0	λ2	..	0	= T	−1	AT ,	(2.21)
Aд =	..		..	..	..	= T		AT ,	(2.21)
		0	0	..
		0	0	..	λn

где λ1 , λ2 ,....λn - собственные значения матрицы А, а ее нормированные собственные векторы являются столбцами матрицы Т. Условием нормирования является равенство единице модуля собственных векторов.

Собственные векторы и собственные значения матрицы А однозначно определяют динамические свойства системы, задаваемые уравнением (2.10).

Пример 2.2. Рассмотрим приведение многомерного объекта, задаваемого уравнением

dx1		= a x + a x						2		+b u			1

dt	11			1	12					11

dx2		= a	21	x	+ a	22	x		2	+b	22	u		2	(2.22)

dt				1

y= c1 x1 +c2 x2

кканонической диагональной форме. Матрицы объекта равны

− 0,5	3		1	0	; C = (3 2).
A =	− 2,5	; B =		2	; C = (3 2).
− 0,25	− 2,5	0		2

Вычислим собственные значения матрицы А из условия (2.16)
	−0,5 −λ	3	= 0	(2.23)

	−0,25	−2,5 −λ
Раскрывая полученный определитель, получим характеристическое
уравнение объекта
λ2 +3λ + 2 = 0.				(2.24)

Корни этого уравнения или собственные значения матрицы А равны

λ1 = −1; λ2 = −2

Найдем собственные векторы матрицы, подставляя в (2.23) вычисленные собственные значения. Для первого собственного значения имеем:

(−0,5 −λ1 )t11 +3t21				= 0			= 0.	(2.25)
−0,25t	11	+(−2,5	−λ	1	)t	21

При λ1 = −1 получаем следующую систему уравнений для вычисления первого собственного вектора

0,5t11 +3t21 = 0 . −0,25t11 −1,5t21 = 0

Откуда t11 = −6t21 .

Конкретные значения первого собственного вектора определяются условием нормировки

t112 +t212 =1.

(2.26)

Подставляя сюда, решение уравнений t11 = −6t21 получим

t11 =	6	; t21 = −	1
	37		37

Аналогично найдем и второй собственный вектор t12 = − 25 ;t21 = 15 .

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

		6	−	2
T =		37	−	5	.
	−	1		1
	−	37		5
		37		5
Найдем обратную матрицу T−1
		37		37
		4		2
T−1 =		4		2	.
		5	3	5
		5	3	5
		4		2

Не трудно убедиться, что

T × T−1		1	0
T × T−1	= I =		1	.
	0		1

Новый вектор координат q задается линейным преобразованием
				x = Tq	.	(2.27)
Подставляя его в уравнения объекта (2.22), получим
		T	dq	= ATq + Bu .		(2.28)

			dt
Умножая обе части уравнения на T −1 слева будем иметь
I	dq	= T−1ATq + T−1Bu .				(2.29)

	dt
Матрица A д = T −1 AT , в соответствии с (2.21), будет						иметь

диагональный вид, где в главной диагонали стоят ее собственные значения. Действительно

				37		37						6	−	2
							− 0,5		3			37	−	5	−1		0
A							− 0,5		3						−1		0
A		= 4				2	×			×					=			.
	д		−	5	−	3 5		− 0,25	− 2,5		−	1		1		0	− 2
			−	5	−	3 5		− 0,25	− 2,5		−	37		5		0	− 2
				4		2						37		5

Вычислим новую матрицу управления Bд = T−1B
			37	37							37
								1	0				37
B		= 4		2				1	0
B	д	= 4		2			×		= 4
	д		5 3 5								5
			5 3 5					0 2			5	3 5
			4	2						−	4	3 5
			4	2							4
и матрицу измерения
					6			−	2
					37				5		16	−	4
Сд = (3 2)×					37				5	=		−
				−	1				1		37		5
				−		37			5		37		5
						37			5
Тогда в новых координатах q1,							q2	уравнения объекта примут вид

dq1

= λ q

+ b

д11

+ b

д12

;

1 1

dq2

= λ2 q2

+ bд21u1 + bд22u2 ; .

(2.30)

y = cd1q1 + cd 2 q2 .

Проведем аналогичные расчеты, используя систему MATLAB, используя сначала операторы матричных преобразований, а затем - оператор canon, который осуществляет линейные преобразования моделей систем заданных в пространстве состояний и поддерживает две канонические формы модальную (диагональную) и присоединенную (форму наблюдаемости)

Программа расчета

A=[-.5 3;-.25 -2.5] B=[1 0;0 2]

C=[3 2]

D=[0 0]

[T Ad]=eig(A) % Вычисление собственных векторов и собственных значений матрицы А To=inv(T) % Вычисление обратной матрицы Т

I=T*To % Умножение матрицы Т на Т-1

Ad=To*A*T % Вычисление коэффициентов канонической диагональной модели

Bd=To*B

Cd=C*T

sys=ss(A,B,C,D) % Модель исходной системы в пространстве состояний sd=ss(Ad,Bd,Cd,D) % Модель приведенной системы в пространстве состояний [sysd,Td]=canon(sys,'modal') % Получение канонической диагональной (модальной) формы

[sysn,Tn]=canon(sys,'companion') % Получение канонической формы наблюдаемости

Результаты расчетов с использованием матричных преобразований

T =

0.9864 -0.8944 -0.1644 0.4472

Ad =
-1	0
0	-2
To =
1.5207		3.0414
0.5590		3.3541
I =
1.0000		-0.0000
-0.0000		1.0000
Ad =
-1.0000		0
0.0000		-2.0000
Bd =
1.5207		6.0828
0.5590		6.7082
Cd =
2.6304		-1.7889

Результаты расчетов в пространстве состояний

a =		Исходная модель
a =	x1	x2
	x1	x2
x1	-0.5	3
x2	-0.25	-2.5
b =	u1	u2
	u1	u2
x1	1	0
x2	0	2
c =	x1	x2
	x1	x2
y1	3	2
d =	u1	u2
	u1	u2
y1	0	0

Модель, преобразованная к диагональной форме a =

	x1	x2
x1	-1	0
x2	0	-2
b =	u1	u2
	u1	u2
x1	1.5207	6.0828

c =	x2	0.55902	6.7082
c =		x1	x2
		x1	x2
d =	y1	2.6304	-1.7889
d =		u1	u2
		u1	u2
	y1	0	0
Матрица перехода в новую систему координат
Td =
1.5207		3.0414
0.5590		3.3541
Модель, преобразованная к форме наблюдаемости
a =		x1	x2
		x1	x2
	x1	0	-2
b =	x2	1	-3
b =		u1	u2
		u1	u2
	x1	1	-4
c =	x2	0	-8
c =		x1	x2
		x1	x2
d =	y1	3	-2
d =		u1	u2
		u1	u2
	y1	0	0
Матрица перехода в новую систему координат
Tn =
1	-2
0	-4

2.3 Структурированные модели.

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

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

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

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

2.Все звенья, из которых состоит система, является детектирующими. В детектирующем звене выходная величина зависит только от входной. Если выходная величина звена оказывает влияние на входную, то звено называется недетектирующим.

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

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

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

dnY

dn−1Y

+L+a

+a Y =b

dmU

dm−1U

+L+b U,

(2.31)

0 dtn

1 dtn−1

n−1 dt n

0 dtm

1 dtm−1

где a0 ,a1 ,K, an ;b0 ,b1 ,K,bn - постоянные коэффициенты; n - порядок системы.

Для реальных физически реализуемых систем управления m < n .

Подвергая (2.31) преобразованию Лапласа при нулевых начальных условиях получим алгебраическое уравнение, связывающее изображения по Лапласу от входной X(p) и выходной Y(p) величины объекта

(a0 pn +a1pn−1 +K+an)Y(p)=(b0 pm +b1pm−1 +K+bm) X(p),

где p – оператор Лапласа

Последнее уравнение можно представить в виде:

Y (p)	=	b0 pm + p1 Pm−1 +K+ bm
			.
X ( p)		a0 pn + a1 pn−1 +K+ an

(2.32)

(2.33)

Это отношение называется передаточной функцией объекта и обозначается символом W(p).

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

Зная передаточную функцию системы или звена можно легко получить дифференциальное уравнение в форме (2.31), справедливо также и обратное утверждение.

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

Пусть имеется многомерный объект управления со структурной схемой рис. 2.1 б. По аналогии с одномерными системами (2.32) можно записать:

Q( p)y( p) = R( p)u( p) + S( p)f ( p) ,					(2.34)
где Q(p)-квадратная матрица операторных коэффициентов размера n × n
q	( p); q	( p);....q	( p)
11	12	1n
q21	( p); q22	( p);....q2n		( p)	,
Q( p) =	............................				,
	............................
qn1	( p); qn2 ( p);....qnn ( p)

R(p)- прямоугольная матрица операторных коэффициентов размера n × k

<<< < Предыдущая 1 23 / 223 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 > Следующая >>>

Соседние файлы в папке лекции

#
17.02.20141.61 Mб297Идентификация объектов управления - УП - Семёнов-Артамонов-Брюхачев - 2003 - 215.pdf
#
17.02.20141.23 Mб170ИДС.doc
#
17.02.20141.6 Mб115Лекции по ИДС.doc
#
17.02.2014167.42 Кб128ОСНОВНЫЕ СВЕДЕНИЯ ОБ ИДЕНТИФИКАЦИИ И МОДЕЛИРОВАНИИ.doc
#
17.02.2014853.5 Кб94Случайные сигналы.doc
#
17.02.2014100.35 Кб82Тесты Идентификация и диагностика 1.doc