6 сем / Метод_указания_к_расч_заданию-2023 (1)
.pdfРис. 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, 2λ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 , 2λ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