Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

МОТС ИТМО

.pdf
Скачиваний:
77
Добавлен:
14.04.2015
Размер:
11.04 Mб
Скачать

0

0

5

0

 

0

1

U

 

; S

4

 

; R

 

.

0

0

 

2

10

7

12.2.Определить, разрешимы ли

 

матричные

уравнения вида

MS RM U относительно матрицы

M ,

а в случае разрешимости

решить первым способом при:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a)

U

 

0

0

 

;

S

5

0

 

; R

 

0

 

1

;

 

 

 

 

 

 

 

 

 

 

 

 

 

0

0

 

 

 

4

2

 

10

7

 

 

 

 

 

0

0

 

2

10

 

 

0

10

 

 

б)

U

 

 

 

;

S

0

 

 

; R

 

 

 

 

7

 

;

 

 

 

0

0

 

 

5

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

в)

U

 

0

0

 

;

S

0

1

 

; R

 

10

0

 

;

 

 

 

 

 

 

 

 

 

7

 

 

 

 

 

 

1

0

 

4

5

 

 

 

 

1

 

 

 

 

2

1

 

 

 

0

1

 

3

0

 

 

г)

U

 

 

 

 

;

S

 

 

 

; R

0

 

 

;

 

 

 

4

2

 

10

7

 

 

2

 

 

 

 

2

4

 

 

 

0

1

 

 

 

4

 

0,707

д)

U

 

 

 

; S

 

 

 

; R

 

 

 

 

 

3

.

 

 

5

10

 

 

8

9

 

0,707

 

 

12.3. Решить

 

матричное

уравнение

 

MS RM U , используя

кронекеровское произведение матриц и приведение уравнения к векторно – матричному виду для матриц S ,R и U из примеров а) - д)

задачи 12.2.

12.4. Используя матричное уравнение (12.12), найти матрицы модального преобразования M , приводящие к диагональному виду следующие матрицы.

 

1

2

2

 

 

 

1 1

2

 

 

3 1

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a)

1

2 1 ;

 

 

б)

 

4 4

 

2

 

;

в)

1 2

3

;

 

2

0 3

 

 

 

 

1 2

 

3

 

 

 

2 5

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5

2 5

 

 

6 3

8

 

 

1 2

3

 

г)

2

1

8

;

д)

 

3 2

3

;

е)

7

1

5

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

3

2

 

 

1 0

 

4

 

 

2 3

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

5 6

 

 

 

2

3

 

2

 

3 1

0

 

ж)

4

3 1

;

 

з)

 

4 6

 

4

; и)

 

4 2

3 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

2 2

 

 

 

1 4

 

7

 

 

4 2

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

12.5. Путем перестановки произвольных столбцов матрицы преобразования M убедиться, что в диагональных матрицах, подобных матрицам а) и) задачи 12.4, происходит согласованная с перестановкой столбцов перестановка собственных чисел на диагонали.

229

12.6. Найти матрицу преобразования подобия M для следующих пар подобных матриц:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a)

5

1

38

81

 

;

б)

17

6

 

14

60

 

 

,

 

 

 

 

,

3

 

 

 

9

1 16

34

 

 

45

16

13

 

3

2

1 24

11

22

 

 

3

2 5 6 20

34

в)

 

2

2

2 , 20

8

20

; г)

2

6 10 , 6 32

51

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

6

5 12

6

10

 

 

1

2 3

4 20

32

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6

6

15

 

 

37

20

4

 

 

 

 

 

 

д)

 

1

5 5

,

 

34

17

4

 

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

2

 

2

 

119

70

11

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

6

15

 

13

70 119

 

 

 

 

 

 

е)

 

1

3 5

,

 

 

4

19 34

 

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

2 4

 

 

4

20 35

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6

10 19

 

4 41

4

26

7

 

 

 

 

1

6 8

 

 

3 14

13

91

18

 

 

ж)

 

 

4 6

 

 

,

4

25

 

.

 

 

 

 

1

 

 

2 40

8

 

 

 

 

 

1 1

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

0 0

0

 

0

1

 

 

 

12.7. Используя формулу

 

(12.24)

для

матрицы

M ,

решить

матричное уравнение (12.6)

частного

вида

(12.10),

где

матрица

Г диагональная матрица, для следующих композиций спектра

Г

матрицы Г, и значений матрицы A, B, H :

 

 

 

 

 

 

 

 

 

 

 

 

 

a) Г {2, 5},

A

 

0

 

1

 

 

0

1

1 ;

 

 

 

 

 

 

 

 

 

, B

 

, H

 

 

 

 

 

 

10

 

7

 

 

1

 

 

 

 

 

 

 

 

 

б) Г { 3, 4},

A

 

7

 

1 , B

2 ,

H 3

5 ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

10

 

0

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

1

2

 

1

 

 

 

 

 

 

в) Г { 1, 2, 6}, A

2

 

 

1

2 , B 2 , H 4

5

6 ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

1

 

1

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

1

1

 

 

1

 

 

 

 

 

 

 

 

 

 

г) Г { 2, 4, 6},

A 2

 

4

1 , B

1 , H 1

 

1

1 ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

1

4

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

5

 

3

 

 

1

0

 

 

 

 

1

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

д) Г { 1, 3, 5},

A 1

2

 

3 , B 0

2 , H

1

 

;

 

 

 

 

 

 

 

3 15

 

12

 

 

3

4

 

 

0 1 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

230

 

 

 

4

 

 

 

 

 

 

 

 

 

 

4

2

 

1

0

 

1

1

2 .

е) Г { 3, 7, 11},

A

2

2

1

, B 1

1 , H

 

 

 

 

 

 

 

 

 

 

 

 

 

4

4 2

0

1

2

1

1

 

 

 

 

 

 

 

 

 

 

 

 

Решение вариантов задач

Задача12.1. Определить, разрешимо ли матричное уравнение

MS RM U

относительно матрицы

M,

если матрицы S, R и U

принимают значения

 

 

 

 

 

0

0

 

5

0

 

0

1

U

 

; S

 

 

; R

 

.

0

0

 

4

2

10

7

Решение. Как указывалось в раздела, уравнение (12.13) при нулевой правой части имеет нетривиальное решение M 0 в случае, если спектры матрицы S и –R совпадают, т.е. {S} { R}. Таким

образом, первым шагом решения задачи является установление факта совпадения спектров матриц S и –R . Характеристические уравнения матриц S и –R принимают вид

5

0

 

2

7 10 0,

det( I S ) det

 

 

 

4

 

 

2

 

 

 

 

1

2 7 10 0.

det( I R) det

 

 

10

7

 

 

 

Из приведенных характеристических уравнений матриц S и –R видно, что они идентичны, а, следовательно, спектры матриц

совпадают.

Таким

образом,

уравнение

 

MS RM 0 имеет

нетривиальное решение. Найдем матрицу M , используя четвертый

(поэлементный) способ.

Уравнения MS RM 0

в поэлементной

форме запишется:

 

 

 

 

 

 

 

m11

m12 5

0

0

1 m11

m12

 

0

0

 

 

 

 

 

m22

 

 

.

m21

m22 4

2

10

7 m21

 

0

0

Перемножив матрицы и сложив их поэлементно, получим матричное соотношение

5m11 4m12 m21

2m12

m22

 

0

0

 

10m11 4m22

 

 

 

 

.

2m21

10m12 5m22

0

0

Если две матрицы равны, то равны их элементы с одинаковыми индексами, в силу чего для элементов матриц M получим систему скалярных соотношений

5m11 4m12 m21 0, 2m12 m22 0; 2m21 10m11 4m22 0, 10m12 5m22 0.

231

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

5m11 4m12 m21 0; 2m12 m22 0.

Система уравнений относительно элементов матрицы M оказалось недоопределенной, так как содержит два уравнения относительно четырех неизвестных. Зададим два недостающих условия, положив, к примеру, m11 1, m12 1. Матрица M в результате

принимает вид

1

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

M

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Матрица M 1 , обратная M , равняется M 1

 

2

 

1

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

1

 

 

 

Покажем, что матрицы S

 

и

 

 

R ,

 

связанные

друг

с другом

матричным уравнением (12.13),

 

при известных матрицах

M и M 1

могут быть найдены друг из друга с помощью отношений подобия:

 

 

S M 1RM, R MSM 1 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Действительно,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

M 1RM

2 1

 

0

1

1

 

 

1

 

 

 

2

1 1

 

 

2

 

 

5 0

 

S

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

1

 

 

 

 

 

 

 

 

1 1 10

 

7 1

 

2

 

 

3

 

 

4

 

4 2

 

и, наоборот,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

MSM 1

1

 

1

 

5

0

 

2

1

 

 

0

1

R .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 2

4

2 1 1

10 7

 

 

 

 

 

 

В заключение заметим, что произвольность назначения двух

элементов матрицы

M не сказалось на

выполнении

соотношения

подобия, т.е. на связи матрицы

 

S и R . Положим теперь

m11 3

и

m 5, в этом случае матрицы M и M 1 принимают вид

 

 

 

 

12

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

5

, M 1

2 11

 

1 11

 

 

 

 

 

 

 

 

 

 

 

 

 

 

M

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

5

10

 

 

 

 

1 11

 

3 55

 

 

 

 

 

 

 

 

 

 

 

 

 

Используя соотношения подобия матриц S

 

и R ,

для последних

получаем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

5 0

 

 

M 1RM

2 11 1 11

 

1 3

5

 

 

 

S ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 11

3 55 10

 

7 5

10

 

 

4 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3 5

 

 

5 0

 

 

2 11 1 11

 

0 1

 

 

 

MSM 1

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

R .

 

 

 

5

10

2 1 11

 

3 55

10

 

7

 

 

232

13*. МОДЕЛИ ТРАЕКТОРНОЙ ПАРАМЕТРИЧЕСКОЙ ЧУВСТВИТЕЛЬНОСТИ ДИНАМИЧЕСКИХ ОБЪЕКТОВ И

СИСТЕМ

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

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

13.1.Модели траекторной параметрической чувствительности непрерывных динамических объектов и систем

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

В связи с тем, что наблюдение дополнительного движения осуществляется с помощью дополнительной динамической системы с фиксированными параметрами, именуемой моделью траекторной чувствительности (МТЧ), аппарат дает возможность разработчику при формировании объекта управления, представляющего собой функциональное объединение физического (технологического) процесса, регулирующих органов и устройств измерения компонентов вектора состояния и выхода, сравнивать конфигурации ОУ на предмет оценки потенциальной стабильности показателей качества проектируемой системы в условиях неопределенности параметров. Анализ управляемости составной системы «номинальный ОУ – МТЧ» по выходу модели траекторной чувствительности с помощью аппарата матриц управляемости по состоянию и выходу МТЧ позволяет ранжировать параметры по степени достижимости стабильности показателей качества систем с использованием возможностей

233

неадаптивных алгоритмов управления, рационально распределять ресурсы управления, решать задачу «оптимального номинала» агрегатов ОУ. Применительно к спроектированной системе аппарат ФТЧ позволяет как на траекторном, так и на структурном уровне оценивать эффект введения в состав системы устройств управления (УУ) в условиях параметрической неопределенности, проводить сравнения альтернативных вариантов УУ. Применение аппарата функций траекторной чувствительности к дискретным динамическим системам дает возможность как траекторно, так и структурно оценивать влияние таких чисто «дискретных» параметров, как интервал дискретности и запаздывания вывода из микроЭВМ вычисленного управления.

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

вектором состояния x Rn , вектором выхода y Rm , а также вектором

q квазистационарных параметров (q(t) 0) ,

который

претерпевает

вариацию q относительно вектора q0

номинальных значений

параметров

так, что

q q q, q RP .

Чтобы

обеспечить

 

 

0

 

 

прозрачность

трактовки

результатов,

будем

использовать

безразмерную форму представления элементов q j вектора параметров

q ( j 1, p) .

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

x t,q q0

q x t x t,q0 , q ,

(13.1)

y t,q q0

q y t x t,q0 , q ,

(13.2)

 

 

 

где x t x t, q0 ; y t y t,q0 . В выражениях (13.1), (13.2)

x t и

y t представляют собой номинальные траектории непрерывной

динамических объекта или системы соответственно по состоянию и

выходу, x t,q0 , q

и y t,q0 , q дополнительные движения

объекта или системы по состоянию и выходу, определяемые вариацией

q , а также номинальным значением q0 вектора параметров.

Будем

полагать справедливыми две гипотезы: первая – о малости

 

 

 

q

 

 

 

 

нормы

 

 

 

 

вариации q вектора

 

параметров, вторая – о непрерывной

дифференцируемости по вектору параметров q в точке

q q0

траекторий x t, q и y t, q в каждый момент времени. Тогда (13.1) и

(13.2) принимают вид

 

 

 

 

 

 

 

 

 

 

 

 

x t, q x t

x t, q

 

q q

q O2 q ,

(13.3)

 

 

 

 

q

 

x

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

234

 

y t, q y t

 

y t, q

 

 

q q

q Oy2 q ,

(13.4)

 

 

q

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

в которых выполняются соотношения

 

 

 

 

 

 

O2 q

 

 

 

 

 

 

 

 

 

 

 

 

O2 q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

lim

 

 

x

 

0; lim

 

 

 

 

 

y

 

 

 

 

 

 

0 .

(13.5)

 

 

 

 

 

q

 

 

 

 

 

 

 

 

 

 

q

 

 

 

 

 

 

q

0

 

 

 

 

 

 

 

q

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Если воспользоваться (13.3)–(13.5), то для дополнительных

движений x t,q0 , q и

y t,q0 , q

параметрически

возмущенной

системы можно записать:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x t,q0 , q t q ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(13.6)

y t, q0 , q t q .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(13.7)

Матрицы Якоби вида

t и

t именуются

матрицами

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

 

 

 

x t, q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

|q q

 

 

 

 

 

t row j t

q

 

; j 1, p

,

(13.8)

 

 

 

 

0

 

 

 

 

 

 

 

j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y t, q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

|q q

 

 

 

 

 

 

row j t

q

 

; j 1, p

,

(13.9)

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

j

 

 

 

 

 

 

 

где

j (t)

и j (t)

являются

функциями

траекторной

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

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

максимальное M t и минимальное m t сингулярные числа матрицы t , задают значение нормы максимальной и минимальной полуосей элипсоидных покрытий дополнительных движений (13.6) и (13.7), порожденных сферой q 1, а элементы правого сингулярного базиса матрицы t задают сочетания вариаций параметров, порождающие максимальную и минимальную полуоси этого покрытия.

235

Конструирование модели траекторной чувствительности проиллюстрируем на примере линейного непрерывного ОУ, матричные компоненты модельного представления которого зависят от

вектора параметров q .

x t,q q x t,q q u t ; x 0,q x 0 ; y t,q C q x q t,q , (13.10)

где x Rn , u Rr , y Rm q,t . Продифференцируем выражение (13.10) по j -му компоненту q j вектора параметров q в точке q q0 . Поменяем порядок дифференцирования по времени t и параметру q j в

левой части первого уравнения (13.10) так, что получим цепочку равенств

 

 

x t, q

 

 

 

 

 

 

 

 

 

dx

 

t, q

 

 

 

 

 

 

 

 

d

x

 

t, q

 

 

 

 

j t ,(13.11)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q j

 

 

 

 

 

 

 

 

 

 

 

 

 

q j

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

q j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q q0

 

 

 

 

 

 

 

 

 

q q0

 

 

dt

 

 

 

q q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

а также введем обозначения

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

A q

 

 

 

; B

 

 

 

 

B q

 

 

 

 

 

; C

 

C q

 

 

 

,

 

(13.12)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q j

 

 

 

 

 

 

 

 

 

 

 

q j

 

 

 

 

 

 

 

 

 

 

q j

 

 

 

 

 

q j

 

 

 

 

 

 

 

 

 

q j

 

 

 

 

 

 

 

 

 

 

 

q j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q q0

 

 

 

 

 

 

 

 

 

 

q q0

 

 

 

 

 

 

q q0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A q

 

q q

A; B q

 

q q

 

 

B; C q

 

q q

 

C ,

 

 

 

 

 

 

 

 

 

(13.13)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x t, q

 

q q

 

x t ;

y t, q

 

q q

 

 

y t .

 

 

 

 

 

 

 

 

 

 

 

 

 

(13.14)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Теперь для

 

j -й модели траекторной чувствительности получим

представление

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j t A j t Aq

x t Bq j u t ; j

C j t Cq j

x t ,

(13.15)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

МТЧ (13.15) будет генерировать функции траекторной

чувствительности j

t

 

по состоянию

 

и

j

t

по выходу,

если ее

дополнить моделью номинального ОУ (см. рисунок 13.1), полученной из (13.10) при q q0 :

x t Ax t Bu t ; x 0 ; y t Cx t .

(13.16)

Нетрудно видеть из (13.10) и (13.16), что динамическая модель дополнительных движений (13.6) и (13.7) с точностью до

мультипликативной составляющей q j j 1, p по выходам j j иj t совпадает системой из p МТЧ (13.15).

Анализ возможности сведения дополнительных движений к нулю хотя бы в асимптотике сводится к анализу управляемости МТЧ вида (13.15). Для этих целей сконструируем агрегированную систему с

составным вектором xj col x, j размерности dim x 2n , которая

236

объединением (13.15) и (13.16), получает векторно-матричное представление

xj (t) Aj xj (t) Bju(t); x j (0) col{x(0),0}

 

 

 

 

 

 

 

 

 

(13.17)

x(t) Cxj xj (t); y(t) C j xj (t); j (t) C j xj (t); j (t) C j xj (t)

(13.18)

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

0

 

 

 

 

 

 

 

 

 

 

B

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Aj

 

A

A

;

 

Bj

 

B

 

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(13.19)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q j

 

 

 

 

 

 

 

 

 

 

 

q j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C

x j

I

n n

0

n n

; C

 

 

C

0

m n

; C

 

0

n n

I

n n

; C

C C

.(13.20)

 

 

 

 

 

 

 

 

j

 

 

 

 

 

 

 

 

 

j

 

 

 

 

 

 

 

 

 

j

 

q j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x t

 

 

 

 

 

 

 

 

 

 

 

 

 

y t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

x t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

B

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

s

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Aq j

 

 

 

 

 

 

Cq

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j t

 

 

 

 

 

 

j t

 

 

 

 

 

 

j t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

B

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q j

 

 

 

 

 

 

 

 

 

 

 

 

s

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рисун к 2.4

 

 

 

 

 

 

 

 

 

 

Рисунок 13.1 Модель траекторной

чувствительности, дополненная

 

 

 

 

 

 

 

g t

 

моделью номинального ОУ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x t

 

 

 

 

 

 

 

 

 

 

 

 

y t

 

 

 

y t, q j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

x t

C

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

G

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Если провести агрегирование номинального ОУ (13.16) и всех p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

s

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

МТЧ (13.15) путем введения вектора x col{x,

 

; j

1, p

} размерности

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j y t, q j

 

dim x ( p 1)n, то векторно-матричное представление такой системы

получает вид

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F

 

 

 

 

 

 

 

 

 

 

 

q j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x(t) Ax(t) Bu(t); x(0) col{x(0),F

(0) 0; j

1, p

},

 

 

 

 

 

q j

j

 

 

 

 

Cq

j

 

 

 

 

 

 

 

 

 

 

 

 

 

x(t) C

x(t); y(t) Cx(t); (t) C x(t); (t) C x(t) ,

x

 

 

 

t

 

 

j

t

 

 

 

 

j

 

 

 

 

 

 

 

 

 

 

 

где

 

 

Gq j

1

 

 

 

 

 

C

 

 

j t

 

 

 

s

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

 

 

 

 

 

 

 

 

 

0n np

 

 

 

 

 

 

 

 

A

 

 

 

 

 

 

 

 

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

col{Aqj ; j 1, p} diag{Ajj

A; j

1, p}

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F

237

(13.21)

(13.22)

(13.23)

Рисунок 2.5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(13.24)

 

 

 

 

 

 

 

 

 

 

 

 

 

B col{B, Baj ; j 1, p}; Cx row In n

0n np ,

 

 

 

 

 

 

 

 

,

 

 

 

 

(13.25)

C row C

 

0m pn ; C 0np n

Inp np

 

 

 

 

 

col {C

 

 

 

 

 

C; j

 

} ,

 

C

; j

1, p

}

 

diag{C

 

1, p

(13.26)

 

 

qj

 

 

 

 

 

 

 

jj

 

 

 

 

 

 

(t) col{ j (t); j

1, p

}; (t) col{ j (t); j

1, p

}.

(13.27)

Нетрудно видеть, что с ростом числа варьируемых параметров

заметно растет размерность

dim x ( p 1)n агрегированной системы

(13.21), (13.22), что может породить проблемы вычислительной устойчивости. В этой связи аддитивная природа дополнительных движений по состоянию (13.06) и выходу (13.07) позволяет p раз воспользоваться агрегированной системой (13.17), (13.18) размерности

dim xj 2n для всех j 1, p .

Для оценки достижимости нулевой траекторной чувствительности к вариациям параметра q j ( j 1, p) , а также ранжирования параметров

по возможным затратам ресурсов управления для достижения нечувствительности траектории проектируемой системы к этим вариациям проведем анализ управляемости системы (13.17), (13.18) по вектору состояния j МТУ и ее выходу j . Первая задача решается на

тройке матриц (C j , Aj , Bj ) , а вторая – на тройке матриц (C j , Aj , Bj ) . Для этих целей сформулируем следующее утверждение.

Утверждение 13.1 (У13.1). Если тройка матриц (C j , Aj , Bj )

полностью управляема для всех

j

1, p

в том смысле,

что матрица

управляемости

 

 

 

 

 

 

 

 

 

 

 

 

 

C B

 

 

 

 

 

C A2n 1B

, j

 

 

(13.28)

W

 

C

A

B

 

1, p

y j

j

j

j

j

 

j

j

j

j

 

имеет ранг,

равный

n (rang Wy j

n) ,

то в системе управления,

полученной агрегированием параметрически возмущенного ОУ (13.10) и УУ, содержащего в своем составе номинальный ОУ (13.16) и реализующего закон управления по вектору дополнительного движения x(t,q0 , q j ) , достижима в асимптотике траекторная

нечувствительность вектора состояния x(t)

к

вариациям всех

 

 

 

 

 

 

компонентов q j ( j

1, p

)

вектора параметров

q

относительно их

номинальных значений в смысле выполнения условия

 

 

 

 

lim x(t, q0 , q j ) 0

j

1, p

 

 

(13.29)

t

 

 

 

 

 

с наперед заданным темпом.

 

238

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]