Многомерные и многосвязные системы (86216)

Посмотреть архив целиком















Контрольная работа

«Многомерные и многосвязные системы»




Задание


Для многомерной системы, заданной матрицами А, В, С, получить:

1. Передаточную функцию ;

2. Частотную передаточную функцию ;

3. Годограф;

4. Импульсную характеристику ;

5. Переходную характеристику ;

6. ЛАЧХ ;

7. ФЧХ .

Составить структурную схему системы.




Дано:


;

;

.


Решение:


1. Передаточная функция

Рассматриваем линейную систему с постоянными параметрами:


,

.


Преобразуем по Лапласу матричные уравнения:


; (1)

, (2)


где


; ;


лапласовы преобразования координат состояния , выходных и входных сигналов.

Преобразуем уравнение (1):




Выносим за скобки:



где

единичная матрица.

Умножаем слева на обратную матрицу:



Откуда получаем:


.


Подставляем в уравнение (2):



Получаем:



Выражение называют передаточной функцией системы.

Находим её:



Находим обратную матрицу:



Подставляем:


.


2. Частотная передаточная функция

Для получения частотной передаточной функции производим замену в передаточной функции :


,


получаем:


.


Выделим действительную и мнимую части:


,


для этого умножим числитель и знаменатель на комплексно – сопряжённый знаменатель:


;

;

;

.


3. Годограф

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

Изменяя частоту, производим расчёт действительной и мнимой частей частотной передаточной функции.

Результат расчёта записываем в таблицу 1.


Таблица 1. Расчёт годографа

0

2,8750000

0,0000000

10

-0,0512719

0,4570747

200

-0,00018

0,020008

1

2,7230769

0,9846154

20

-0,0163435

0,2074170

300

-0,000078

0,013336

2

1,9500000

1,9000000

30

-0,0075500

0,1355448

400

-0,000044

0,010001

3

0,8344828

1,9862069

40

-0,0043030

0,1009350

500

-0,000028

0,008001

4

0,2250000

1,5500000

50

-0,0027705

0,0804792

600

-0,000019

0,006667

5

0,0130624

1,1611030

60

-0,0019302

0,0669441

700

-0,000014

0,005715

6

-0,0500000

0,9000000

70

-0,0014209

0,0573176

800

-0,000019

0,005000

7

-0,0645030

0,7269777

80

-0,0010893

0,0501171

900

-0,000009

0,004445

8

-0,0634615

0,6076923

90

-0,0008614

0,0445267

1000

-0,000007

0,004000

9

-0,0578113

0,5216604

100

-0,0006982

0,0400600

2000

-0,000002

0,002000


Можно построить график на комплексной плоскости – рис. 1.



Рис. 1. Годограф


4. Импульсная характеристика

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


.


Найдём полюса передаточной функции:



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

Разложим передаточную функцию на простые дроби:


.


Используя табличные значения, находим:


,

.


Таким образом, получаем:


.


Изменяя время от нуля до 5 секунд, производим расчёт по формуле, результаты заносим в таблицу 2.


Таблица 2. Импульсная характеристика

0

0,5

1

1,5

2

2,5

3

3,5

4

4,5

5

-4

11,28

62,69

100,8

-167,1

-1236

-2395

2097

23854

54578

-15944


Строим график импульсной характеристики – рис. 2.


Рис. 2. Импульсная характеристика


5. Переходная характеристика

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


.


Найдём полюса передаточной функции:


; .


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

Разложим передаточную функцию, делённую на р, на простые дроби:


.


Приводим к общему знаменателю:


.


Приравниваем коэффициенты при равных степенях р:


,

,

.


Откуда находим:


,

,

.


Используя табличные значения, находим:


,

,

.


Таким образом, получаем:


.


Изменяя время от нуля до 5 секунд, производим расчёт по формуле, результаты заносим в таблицу 3.


Таблица 3. Переходная характеристика

0

0,5

1

1,5

2

2,5

3

3,5

4

4,5

5

0

0,654

17,59

62,52

69,32

-243

-1209

-1744

3830

24151

42653


Строим график переходной характеристики – рис. 3.


Рис. 3. Переходная характеристика


6. ЛАЧХ

Для получения ЛАЧХ найдём модуль частотной передаточной функции:


.


далее находим 20 десятичных логарифмов от найденного модуля:


.


Это и есть выражение для ЛАЧХ.

Расчёт значений ЛАЧХ ведём в логарифмическом масштабе. Результаты записываем в таблицу 4. Размерность ЛАЧХ – децибелы (дБ).



Таблица 4. ЛАЧХ

-1

0,1

9,17406

0,1

1,25893

9,20891

1,2

15,8489

-11,426

-0,9

0,12589

9,17482

0,2

1,58489

9,08243

1,3

19,9526

-13,614

-0,8

0,15849

9,17601

0,3

1,99526

8,70564

1,4

25,1189

-15,738

-0,7

0,19953

9,17788

0,4

2,51189

7,83066

1,5

31,6228

-17,818

-0,6

0,25119

9,18077

0,5

3,16228

6,23375

1,6

39,8107

-19,869

-0,5

0,31623

9,18519

0,6

3,98107

3,94960

1,7

50,1187

-21,902

-0,4

0,39811

9,19182

0,7

5,01187

1,26946

1,8

63,0957

-23,923

-0,3

0,50119

9,20135

0,8

6,30957

-1,5050

1,9

79,4328

-25,936

-0,2

0,63096

9,21400

0,9

7,94328

-4,1982

2

100

-27,944

-0,1

0,79433

9,22792

1

10

-6,7459

2,1

125,893

-29,950

0

1

9,23483

1,1

12,5893

-9,1470

2,2

158,489

-31,953


Строим график ЛАЧХ – рис. 4.


Рис. 4. ЛАЧХ


7. ФЧХ

ФЧХ – угол поворота вектора на комплексной плоскости в зависимости от частоты:


.


Расчёт значений ФЧХ ведём в логарифмическом масштабе. Результаты записываем в таблицу 5. Размерность ФЧХ – радианы (рад).


Таблица 5. ФЧХ

-1

0,1

0,03263

0,1

1,25893

0,44997

1,2

15,8489

1,66382

-0,9

0,12589

0,04110

0,2

1,58489

0,58831

1,3

19,9526

1,64958

-0,8

0,15849

0,05177

0,3

1,99526

0,77030

1,4

25,1189

1,63592

-0,7

0,19953

0,06524

0,4

2,51189

0,99225

1,5

31,6228

1,62384

-0,6

0,25119

0,08227

0,5

3,16228

1,22480

1,6

39,8107

1,61359

-0,5

0,31623

0,10383

0,6

3,98107

1,42316

1,7

50,1187

1,60513

-0,4

0,39811

0,13123

0,7

5,01187

1,56064

1,8

63,0957

1,59824

-0,3

0,50119

0,16622

0,8

6,30957

1,63913

1,9

79,4328

1,59268

-0,2

0,63096

0,21126

0,9

7,94328

1,67427

2

100

1,58822

-0,1

0,79433

0,26981

1

10

1,68250

2,1

125,893

1,58466

0

1

0,34696

1,1

12,5893

1,67633

2,2

158,489

1,58182


Строим график ФЧХ – рис. 5.


Рис. 5. ФЧХ



8. Структурная схема системы

Записываем матричные уравнения системы:


;

.


Подставляем исходные данные:


;

.


Производим умножение матриц:


,

,

.


Получили систему уравнений, на основе которой строим структурную схему – рис. 6.


Рис. 6. Структурная схема системы



Часть 2:

Осуществить синтез замкнутой системы с собственными числами

{–1; –4;  5j}.

Построить наблюдатель полного порядка.


Дано:

,

,

.


Решение:


1. Синтез замкнутой системы

Рассматриваем линейную систему с постоянными параметрами:


,

.


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


,


где

входной командный сигнал,

К – матрица коэффициентов обратной связи.

После замыкания эта система имеет структуру, изображённую на рис. 7.


Рис. 7. Структура исходной системы


Движение системы описывается линейным дифференциальным уравнением:


.


Таким образом, динамические свойства системы полностью определяются матрицей А – ВК, её характеристическими числами.

Характеристический многочлен исходной системы равен:


.


Спектр характеристических чисел (корни характеристического многочлена):

.

Желаемый характеристический многочлен замкнутой системы по условию имеет 4 собственных числа, но наша исходная система имеет третий порядок, поэтому одно из собственных чисел необходимо убрать, убираем собственное число (–1), тогда:


.


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


.


Приравниваем коэффициенты при равных степенях многочленов и :

,

,

,

.

Решая полученную систему уравнений, получаем:

,

,

.

Искомое управление принимает вид:


.


Структура синтезированной системы представлена на рис. 8.

Она построена по уравнениям:


,

,

,

,

.


Рис. 8. Структура синтезированной системы


2. Построение наблюдателя полного порядка

Система



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

Найдём структуру асимптотического наблюдателя, для чего определим ошибку восстановления и найдём модель её изменения:


.


Затем потребуем, чтобы при всех и .

Это равенство возможно при:


,

.


Таким образом, структура асимптотического наблюдателя полного порядка определяется моделью вида:


.


На рис. 9 изображена структура системы и её наблюдателя.


Рис. 9. Структура системы с наблюдателем


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

Пусть ошибка восстановления , тогда

.

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

Пусть матрица


,


тогда матрица


.

Полюса наблюдателя определяются уравнением:


.


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

{– 4;  5j},

то расположим полюса наблюдателя в точках:

.

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


,


что будет иметь место тогда, когда:


,

,

.


Решая полученную систему уравнений, получаем:


;

;

.


Находим матрицу:



Модель асимптотического наблюдателя системы принимает вид:


,

,

,

.


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

Она построена по уравнениям:


,

,

,

,

,

,

.


Случайные файлы

Файл
124942.rtf
150880.rtf
138836.rtf
177257.rtf
22481-1.rtf




Чтобы не видеть здесь видео-рекламу достаточно стать зарегистрированным пользователем.
Чтобы не видеть никакую рекламу на сайте, нужно стать VIP-пользователем.
Это можно сделать совершенно бесплатно. Читайте подробности тут.