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

6 сем / Метод_указания_к_расч_заданию-2023 (1)

.pdf
Скачиваний:
0
Добавлен:
26.01.2024
Размер:
8.5 Mб
Скачать

Рис. 4.16. График процесса x(t ) при K = 5K1 = 10 , λ = 0,84 и НУ ( x0 4 , v0 1)

Из графика рис. 4.16 видно, что в случае увеличения K1 в 5 раз, λ в 2 раза и НУ (x0 4 , v0 1) период автоколебаний системы равен 2,1

сек, а амплитуда равна 44. То есть существенно увеличилась амплитуда автоколебаний, примерно в 5 раз, при незначительном увеличении их периода с 2,0 сек до 2,1 сек.

5. Проведем исследование автоколебаний в системе приближенным амплитудно-частотным методом (методом Гольдфарба).

В соответствии с методом Гольдфарба пересечение амплитуднофазовой характеристики линейной части модели Wлч (jω) и инверсной

характеристики эквивалентного комплексного коэффициента усиления нелинейного элемента z(A) свидетельствует о наличии решения уравнения автоколебаний

Wлч(jω) z(A).

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

W12 (p ) W1(p )W2 (p )

32

,

p 2 8 p 16

W30 (p ) 1 pT0 1 2 p .

31

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

 

 

 

Wнэ (A) a (A) jb(A) ,

 

 

где a (A)

2B

(

A2 h 2

A 2 c 2 )

10

 

 

( A2

25

A2 4),

πA 2

 

πA2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b(A)

2B (h c)

 

 

 

30

.

 

 

 

 

 

πA2

 

 

 

 

 

 

 

 

 

 

πA 2

 

 

Инверсная характеристика комплексного коэффициента нелинейного элемента имеет вид:

z (A)

1

 

1

.

Wнэ (A)

a (A) jb(A)

 

 

 

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

W лч (jω) W1(jω)W2 (jω) W3 (jω) 1 K1

K 2

( jωT0

1)

(jω)2 ajω b

 

64 jω 32

.

 

 

 

 

 

 

(jω)2 8 jω 16

 

 

Для построения АФХ линейной части системы и инверсной характеристики эквивалентного комплексного коэффициента усиления нелинейного элемента воспользуемся ППП Mathcad. После задания выражений линейной части модели и инверсной характеристики комплексного коэффициента усиления нелинейного элемента в рабочем поле Mathcad добавим объект «график X-Y», с помощью которого отобразим на комплексной плоскости АФХ линейной части модели и инверсной характеристики комплексного коэффициента усиления нелинейного элемента.

Перед объектом «график X-Y» укажем требуемый диапазон значений и шаг расчетов для частоты ω и амплитуды A .

ω := 0.02, 0.03..1000

А := 5.1, 5.2..60

Результат построения представлен на рис. 4.17. Из него следует, что имеется пересечение двух характеристик – это свидетельствует о наличии автоколебаний.

32

 

 

 

5

Im(Wlc(i ))

10

5

0

Im( z(A))

 

 

 

 

 

 

5

Re(Wlc(i )) Re( z(A))

Рис. 4.17. АФХ Wлч (jω) и [ z (A) ] для п. 5 а) задания

Для определения параметров автоколебаний отсчитаем значения параметров ω и A на соответствующих характеристиках в точке их пересечения. В результате получим (рис. 4.18):

ω := 0.02, 0.03..4.35 А := 5.1, 5.2..51

 

 

 

5

Im(Wlc(i ))

10

5

0

Im( z(A))

 

 

 

 

 

 

5

Re(Wlc(i )) Re( z(A))

Рис. 4.18. Определение параметров автоколебания для п. 5 а) задания

Из рисунка 4.18 видно, что точка пересечения двух характеристик соответствует амплитуде колебаний A 51 и частоте ω 4,35 сек-1, откуда период автоколебаний получается равным 1,44 сек. Поскольку при

33

пересечении характеристика [ z (A) ] выходит из области, ограниченной Wлч (jω) , то автоколебания являются устойчивыми.

Рассмотрим следующий случай задания – п. 5 б): система с параметрами п. 3 при λ=λmin 0,42.

В этом случае изменяется только параметр c 2,9 нелинейного элемента. Найдем эквивалентный комплексный коэффициент усиления нелинейного элемента:

 

 

 

 

Wнэ (A) a (A) jb(A) ,

 

 

 

a (A)

2B

(

A2 h 2

 

A 2 c 2

)

 

10

 

( A 2

 

25

A2 8,41) ,

 

 

πA2

 

πA 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b(A)

2B (h c)

 

 

21

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

πA 2

 

 

πA2

 

 

 

Инверсная характеристика комплексного коэффициента нелиней-

ного элемента имеет вид:

 

1

 

 

 

 

 

1

 

 

 

 

 

 

 

z (A)

 

 

 

 

 

.

 

 

 

 

Wнэ (A)

 

a (A)

jb(A)

 

 

 

 

 

 

 

 

 

 

 

АФХ линейной части модели и инверсной характеристики комплексного коэффициента усиления нелинейного элемента будут иметь

вид (рис. 4.19):

 

ω := 0.02, 0.03..1000

А := 5.1, 5.2..60

 

 

 

5

Im(Wlc(i ))

10

5

0

Im( z(A))

 

 

 

 

 

 

5

 

 

Re(Wlc(i )) Re( z(A))

 

 

Рис. 4.19. АФХ Wлч (jω) и [ z (A) ] для п. 5 б) задания

34

Из рисунка 4.19 видно, что имеется пересечение двух характеристик, что свидетельствует о наличии автоколебаний. Согласно критерию Гольдфарба они являются устойчивыми. Для определения параметров автоколебаний отсчитаем значения параметров ω и A на соответствующих характеристиках в точке их пересечения. В результате получим:

A 51,

ω 4,39 сек-1 ,

откуда период автоколебаний получается равным 1,43 сек.

Рассмотрим последний случай задания – п. 5 в): система с параметрами п. 4 при 5K1 10, min 0,84.

Вэтом случае изменяется параметр c 0,8 нелинейного элемента

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

 

 

 

 

Wнэ (A) a (A) jb(A) ,

 

a (A)

2B

(

A2 h 2

A2

c 2 )

10

(

A2 25

A2 0,64) ,

πA2

 

 

 

 

 

 

 

πA 2

 

 

 

 

 

 

b(A)

2B (h c)

 

42

.

 

 

 

 

 

 

 

 

 

 

 

 

πA2

πA 2

 

Инверсная характеристика комплексного коэффициента нелинейного элемента имеет вид:

z (A)

1

 

1

.

Wнэ (A)

a (A) jb(A)

 

 

 

Определим комплексный коэффициент усиления линейной части системы:

W

лч

(jω)

jω5K1K 2T0 5K1K1

 

320 jω 160

.

 

 

 

(jω)2 ajω b

 

(jω)2 8 jω 16

Построим АФХ линейной части модели и инверсной характеристики комплексного коэффициента усиления нелинейного элемента

(рис. 4.20).

Из рисунка 4.20 следует, что в системе имеются устойчивые автоколебания. Для определения их параметров отсчитаем значения параметров ω и A на соответствующих характеристиках в точке их пересечения. В результате получим, что амплитуда колебаний A 255, частота ω 4,44 сек-1, откуда период автоколебаний получается равным 1,41 сек.

ω := 0.02, 0.03..1000

А := 5.1, 5.2..300

35

 

 

 

20

 

 

 

10

Im(Wlc(i ))

40

20

0

 

 

 

Im( z(A))

 

 

 

 

 

 

10

 

 

 

20

 

 

 

30

 

Re(Wlc(i )) Re( z(A))

 

 

Рис. 4.20. АФХ Wлч (jω)

и [ z (A) ] для п. 5 в) задания

6. Для составления сравнительной таблицы необходимо пересчитать полученные в п. 5 величины амплитуды автоколебаний, поскольку из рис. 3.9 видно, что они соответствуют сигналу y (t ) , а в п. 2 – 4 рас-

сматривается съем данных относительно x(t ) . В результате пересчета

амплитуды автоколебаний посредством формулы Ax

 

A

 

, где ω

1

T

2

 

 

0

 

 

– частота автоколебаний, была составлена сравнительная таблица результатов исследования автоколебаний различными методами (табл. 4.1).

Таблица 4.1

 

 

Исходные параметры

λmin

5K1 , min

Моделирование

Амплитуда

8,4

8,2

44

Период, с

2,1

2

2,1

 

Метод Гольдфарба

Амплитуда

5,8

5,7

28,5

Период, с

1,44

1,43

1,41

 

Сравнивая результаты, приведенные в табл. 4.1, можно отметить, что качественная картина изменения параметров автоколебаний в зависимости от параметров системы имеет одинаковый характер. Различия амплитуд и периодов автоколебаний, полученных двумя методами, составляют приблизительно 30 %. Объяснить различия можно тем, что линейная часть системы не является идеальным фильтром низких частот, что лежит в основе исходного предположения метода Гольдфарба.

36

СПИСОК РЕКОМЕНДУЕМОЙ ЛИТЕРАТУРЫ

Основной

1. Ким, Д.П. Теория автоматического управления. Т.2. Многомерные, нелинейные, оптимальные и адаптивные системы / Д.П. Ким. – 2-е изд., испр. и доп. – М.: Физматлит, 2007.

Дополнительный

2.Теория автоматического управления: Нелинейные системы, управления при случайных воздействиях: учебник / А.В. Нетушил и др. – М.: Высшая школа, 1983.

3.Ягодкина, Т.В. Исследование САУ с использованием прикладного пакета MatLab. Лабораторный практикум: / Т.В. Ягодкина, С.А. Хризолитова, В.М. Беседин и др. – М.: Издательский дом МЭИ, 2007.

37

Производственно-практическое издание

Беседин Валерий Михайлович Державин Отто Михайлович Сидорова Елена Юрьевна

АНАЛИЗ НЕЛИНЕЙНЫХ СИСТЕМ АВТОМАТИЧЕСКОГО УПРАВЛЕНИЯ

Методические указания к расчетному заданию

Редактор Т.А. Феоктистова Компьютерная верстка Ю.В. Макарова

38

Подписано в печать

22.02.22.

Печать цифровая.

Формат 60х84 1/16

Печ. л. 2,5.

Тираж 75 экз.

Изд. № 21у-106

Заказ

Оригинал-макет подготовлен в РИО НИУ «МЭИ». 111250, г. Москва, ул. Красноказарменная, д. 14.

Отпечатано в типографии НИУ «МЭИ». 111250, г. Москва, ул. Красноказарменная, д. 13.

39

ДЛЯ ЗАМЕТОК

40

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