Операторная форма записи дифференциального уравнения это

Передаточная функция. В ТАУ часто используют операторную форму записи дифференциальных уравнений

В ТАУ часто используют операторную форму записи дифференциальных уравнений. При этом вводится понятие дифференциального оператораp = d/dt так, что, dy/dt = py, а p n = d n /dt n . Это лишь другое обозначение операции дифференцирования. Обратная дифференцированию операция интегрирования записывается как 1/p. В операторной форме исходное дифференциальное уравнение записывается как алгебраическое:

Не надо путать эту форму записи с операционным исчислением хотя бы потому, что здесь используются непосредственно функции времени y(t), u(t) (оригиналы), а не их изображения Y(p), U(p), получаемые из оригиналов по формуле преобразования Лапласа. Вместе с тем при нулевых начальных условиях с точностью до обозначений записи действительно очень похожи. Это сходство лежит в природе дифференциальных уравнений. Поэтому некоторые правила операционного исчисления применимы к операторной форме записи уравнения динамики. Так оператор p можно рассматривать в качестве сомножителя без права перестановки, то есть py yp. Его можно выносить за скобки и т.п.

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

Дифференциальный оператор W(p) называют передаточной функцией. Она определяет отношение выходной величины звена к входной в каждый момент времени: W(p) = y(t)/u(t), поэтому ее еще называют динамическим коэффициентом усиления. В установившемся режиме d/dt = 0, то есть p = 0, поэтому передаточная функция превращается в коэффициент передачи звена K = bm/an.

Знаменатель передаточной функции D(p) = aop n + a1p n — 1 + a2p n — 2 + . + an называют характеристическим полиномом. Его корни, то есть значения p, при которых знаменатель D(p) обращается в ноль, а W(p) стремится к бесконечности, называются полюсами передаточной функции.

Числитель K(p) = bop m + b1p m — 1 + . + bm называют операторным коэффициентом передачи. Его корни, при которых K(p) = 0 и W(p) = 0, называются нулями передаточной функции.

Звено САУ с известной передаточной функцией называется динамическим звеном. Оно изображается прямоугольником, внутри которого записывается выражение передаточной функции. То есть это обычное функциональное звено, функция которого задана математической зависимостью выходной величины от входной в динамическом режиме. Для звена с двумя входами и одним выходом должны быть записаны две передаточные функции по каждому из входов. Передаточная функция является основной характеристикой звена в динамическом режиме, из которой можно получить все остальные характеристики. Она определяется только параметрами системы и не зависит от входных и выходных величин. Например, одним из динамических звеньев является интегратор. Его передаточная функция Wи(p) = 1/p. Схема САУ, составленная из динамических звеньев, называется структурной.

Дата добавления: 2015-12-26 ; просмотров: 857 ; ЗАКАЗАТЬ НАПИСАНИЕ РАБОТЫ

Формы записи дифференциальных уравнений

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

. (2.9)

В этом уравнении — выходная переменная (управляемая (регулируемая) величина) САУ, — входная переменная САУ. Правая часть уравнения (3.1) записана относительно управляющего воздействия , однако используются формы записи уравнения относительно задающего воздействия , возмущения или нескольких входных воздействий.

Применяется также операторная форма записи уравнения (2.9):

. (2.10)

В этом уравнении через « » обозначен оператор дифференцирования .

Заметим, что по сложившейся традиции символ « » используется также в преобразованиях Лапласа и Карсона-Хевисайда, но является комплексным числом .

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

(2.11)

или в операторной форме

. (2.12)

Первая стандартная символическая форма записи уравнения (2.11) имеет следующий вид:

, (2.13)

где ; ; ; .

Форма (2.13) представляет собой операторно-структурное описание системы, т.е. в виде операторов звеньев, составляющих структурную схему системы (далее эти понятия разъясняются), и связей между ними. В этой форме — постоянные времени звена, измеряемые в секундах; — передаточный коэффициент звена.

Из изложенного выше следует, что уравнение (2.9) в этой форме перепишется в следующем виде:

, (2.14)

где ; ; .

Во второй стандартной форме записи дифференциального уравнения используется передаточная функция системы, которая для рассматриваемого примера (2.11) имеет вид

.

Передаточная функция САУ, поведение которой во времени описывается уравнением (2.9), имеет следующий вид :

.

В формуле (2.15) через и обозначены изображения (по Лапласу) выходной и входной переменных САУ при нулевых начальных условиях и равенстве нулю внешних возмущений, а через и — полиномы относительно комплексной переменной .

Вторая стандартная форма записи дифференциального уравнения имеет следующий вид:

или . (2.16)

В (2.16) и являются полиномами (символическими) относительно оператора .

Из сравнения первой и второй стандартных форм записи дифференциальных уравнений следует, что с математической точки зрения различие между этими формами весьма несущественно и состоит лишь в различном представлении коэффициентов уравнений. В ТАУ принято называть уравнения вида (2.9) — (2.14), (2.16) уравнениями типа «вход-выход».

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

Координаты состояния х1, х2 , . ,хn
u
y

Рис.2.2 — Состояние системы

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

Множество возможных положений этого вектора образует векторное пространство , называемое пространством состояний системы. В переменных состояния САУ описывается векторно-матричным уравнением

, (2.17)

где — квадратная матрица коэффициентов (ее называют также собственной параметрической матрицей системы); — входная матрица (матрица управления) системы; — выходная матрица системы;

— вектор переменных состояния — внутренних координат системы;

— вектор входных переменных (управляющих и возмущающих);

— вектор наблюдаемых или выходных переменных; размерности матриц , , , соответственно, ( ), ( ), ( ).

Процессы в САУ в свободном движении (без внешних воздействий) согласно уравнению (2.17) описываются векторно-матричным уравнением с характеристическим уравнением , где — единичная матрица, или в развернутом виде системой дифференциальных уравнений

с характеристическим уравнением

. (2.18)

Эти уравнения при определенных начальных условиях дают возможность изучить процессы в системе путем их решения численными методами с использованием ЭВМ.

Разработаны различные способы перехода от уравнений типа «вход-выход» к уравнениям состояния вида (2.17) и наоборот. Один из наиболее распространенных способов состоит в следующем. Пусть САУ описывается уравнением (2.9). Введем обозначения

, , . , ,

.

С помощью этих обозначений преобразуем уравнение (3.1) к следующему виду:

, (2.19)

где ; ;

; .

В нашем примере и являются скалярными величинами. В общем случае (2.17) — это, соответственно, вектор наблюдаемых или выходных переменных и вектор входных переменных (управляющих и возмущающих), поэтому в (2.19) матрицы и выродились в вектор-столбец и вектор-строку соответственно.

Система уравнений (2.19) представляет собой описание линейной непрерывной системы в пространстве состояний . Уравнения (2.19) с матрицей называют уравнениями в форме Фробениуса.

Если , то

; .

Форма уравнений (2.19) с подобными матрицами и называется в ТАУ канонической формой фазовой переменной.

Задание 1

1.1. По дифференциальному уравнению системы:

Для каждого типового звена 1 – 12 (таблицы 2.1) в соответствии с его параметрами вывести дифференциальное уравнение, операторное уравнение, и выражение передаточной функции.

1.2Математическое описание типового звена системы автоматического регулирования записать в трех формах записи дифференциальных уравнений.

Первая стандартная символическая форма операторно-структурное описание системы, т.е. в виде операторов звеньев.

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

Третья стандартная форма записи дифференциального уравнения — переменные состояния.

Таблица 2.1 – Исходные коэффициенты

№ п.п.Наименование звенаа 0а 1а 2b0b1Примечания
Безынерционное (пропорциональное)к
Инерционное 1-го порядка (апериодическое)Тk
Инерционное 2-го по- рядка (апериодическое)Т2 2Т1kТ1³2Т2
Инерционное 2-го по- рядка (колебательное)Т2 2Т1kТ1 T

Задание 2

2.1Для каждого звена (таблицы 2.2) по его передаточной функции записать дифференциальное уравнение.

2.2 Математическое описание типового звена системы автоматического регулирования записать в трех формах записи дифференциальных уравнений.

ВарПередаточная функцияЗначения параметров передаточной функции
а0=1; а1=5; а2 =1,2; а3 =0,9; а4=0,5;в0=1;в1=3; в2=0,8;в3=0,3
а0=1; а1=5; а2 =1,2; а3 =0,9; в0=1;в1=3; в2=0,8
а0=1; а1=5; а2 =1,2; в0=1; в1=3;
а0=1; а1=5; а2 =1,2; а3 =0,9; а4=0,5; в0=10
а0=1; а1=5; а2 =1,2; а3 =0,9; в0=10
а0=1; а1=5; а2 =1,2; а3=0,9;а4=0,5;в0=1;в1=3; в2=0,8;в3=0,3
Т0=2; Т1=4; Т2=1,1;Т3=0,9
Т0=2; Т1=4; Т2=1,1;Т3=1,1;Т4=,9
К= 10;Т1=4; Т2=1,1;Т3=0,9
К= 10; Т2=1,1;Т3=0,9 Т4=0,9
Т0=0,7; Т1=3;Т2=1,2;Т3=0,9;Т4=0,8;Т5=0,5
К=10 Т0=0,7; Т1=3;Т2=1,2;Т3=0,9;Т4=0,8;Т5=0,5;

Задание №3

3.1 Для заданной схемы необходимо составить операторное уравнение для каждого элемента схемы САУ.

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

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

3.4 Сформировать математическое описание систем автоматического регулирования в виде третьей стандартной формы записи дифференциального уравнения — В переменных состояния САУ описываемых векторно-матричным уравнением.

Схема, показанная на рисунке 2.2, представляет собой САР температуры в помещении. Объектом регулирования (ОР) в дан­ной системе является помещение, для которого регулируемая ве­личина — температура внутри помещения Ө, регулирующее (уп­равляющее) воздействие — температура воздуха ӨК, поступающего из калорифера, возмущающее воздействие — изменения внешних факторов f(в общем случае изменение температуры атмосферного воздуха, его влажности, скорости ветра). При исследовании сис­темы в качестве основного возмущения следует рассматривать из­менение температуры окружающего воздуха.

Воспринимающим органом — ВО (датчиком, чувствительным элементом) в данной САР является терморезистор RД, включен­ный в мостовую схему, обеспечивающую с помощью резистора RОзадание необходимого значения температуры в помещении и выполняющую также функции сравнивающего органа — СО (эле­мента сравнения). Усиление сигнала разбалансаΔU(сигнала рас­согласования) измерительной мостовой схемы обеспечивается посредством усилителя. Усиленный сигнал Uобеспечивает вра­щение двухфазного исполнительного двигателя, который изменя­ет перемещение клапана (заслонки) на трубопроводе подачи парав калорифер, чем достигается изменение температуры воздуха на входе калорифера — регулирующего воздействия на объектерегулирования.

1 — помещение; 2 — теплообменник (калорифер), 3 — измерительная мостовая схема; 4 — двухфазный ис­полнительный двигатель, 5 — дифференциальный магнитный усилитель; 6 — клапан (заслонка)

Рис. 2.2. Схема САР температуры

Динамические свойства объекта регулирования и элементов системы описываются следующими уравнениями:

где То, Т2, Т3, Т4 — постоянные времени, с; Ө — значение температуры воздуха в помещении, °С, Ө к — значение температуры воздуха на выходе калорифера, °С; к, к1, к2, к3, к4— коэффициенты передачи; f— возмущающее воздействие на объекте регулирования; Uд —падение напряжения на термодатчике, В; ΔU— напряжение на выходе мостовой схемы (сигнал рассогласования), В; μ. — линейное перемещение клапана, см; U0 — задающий сигнал, В.

Значения параметров элементов САР по вариантам даны в таб­лице 2.3.

Заданное значение температуры в помещении Ө = 20 °С.

Значения параметров элементов САР

ВариантТ0, сТ2, скк1, В/ 0 Ск4к2, см/(В*с)f,. 0 СК3, °С/см
0,060,20,20,002-11
0,070,250,30,001
0,080,30,250,0018-8
0,090,350,20,002
0,100,40,20,002-5
0,500,180,250,003
0,0550,190,40,0035
0,060,170,40,0025-15
0,060,250,20,0016
0,080,40,150,0014-18

Примечание. Для всех вариантов постоянные времени Т3 = 20 с, Т4=0,5 с.

Схема САР, приведенная на рисунке 2.3, обеспечивает стаби­лизацию угловой скорости электродвигателя постоянного тока который совместно с рабочим механизмом является объектом ре­гулирования. Регулируемая величина объекта — угловая скорость двигателя ω, регулирующее воздействие — напряжение Uг,пода­ваемое от генератора на якорь двигателя. Возмущающее воздейст­вие на объекте регулирования — момент сопротивления Мс, соз­даваемый рабочим механизмом. Угловая скорость двигателя ωконтролируется тахогенератором, сигнал которого Uтг, пропор­циональный скорости, сравнивается с задающим сигналом U3. Сигнал рассогласования ΔU = U3— UTг усиливается магнитным усилителем и воздействует на обмотку возбуждения генератора, выполняющего функции исполнительного органа (элемента).

Динамические свойства объекта регулирования и элементов САР описываются следующими уравнениями:

гдеТд, Ту, Tv — постоянные времени, с; Кд, Км, Ктг, Ку, Кг — коэффициенты передачи соответствующих элементов систем

1 — задающий потенциометр; 2 — магнитный усилитель; 3 — генератор; 4 — двигатель; 5 — тахогенератор; 6 — рабочий механизм

Рис. 2.3. Схема САР угловой скорости электродвигателя

Значения параметров элементов САР

Вари­антТу, сКуКгТг, сКд, рад/ с*ВТд, сКм рад/ с*Н* мМс, Н*мКгг, В*с/ рад
0,0204,02,00,101,00,50,021,0
0,0155,01,80,120,950,600,030,9
0,0184,51,70,150,850,700,04
0,0226,01,50,200,80,800,050,7
0,0205,81,60,161,50,650,060,6
0,0254,22,00,251,40,750,070,5
0,0203,52,70,221,30,800,080,4
0,0286,22,10,301,20,750,020,5
0,0186,52,30,161,00,500,0130,6
0,0147,02,50,201,250,800,0150,7

Значения параметров объекта регулирования и элементов сис­темы для различных вариантов указаны в таблице 2.4. Заданное значение угловой скорости ω = 40 рад/с.

На рисунке 2.4 изображена схема САР давления Р в ресивере (воз­духосборнике) 1, который является в данной системе объектом регу­лирования. Давление в ресивере регулируется посредством изменения количества воздуха Q, зависящего от положения заслонки 2, т.е. от ее линейного перемещения Х3, которое можно рассматривать как регу­лирующее воздействие на входе объекта регулирования. Внешним возмущением, вызывающим отклонение регулируемой величины — давления Р, является изменение расхода сжатого воздуха Qc.

Рис 2.4 Схема САР давления Р в ресивере

Давление в данной системе контролируется с помощью сильфонного датчика 3, выход­ная величина которого — пере­мещение Хс сильфона 5 одно­значно зависит от разности сил ΔF= F0— Fp, где Fp— сила, соз­даваемая давлением Р, F0— си­ла натяжения пружины 6, кото­рое можно изменять винтом 7.

Перемещение сильфона Хсс помощью потенциометрического преобразователя 4 преобразуется в электрический сигнал — напряжение U, которое усиливается электронным усилителем 8. Выходной сигнал усилителя Uyуправляет электромагнитным при­водом 9, связанным с заслонкой 2,

В данной САР сильфонный датчик выполняет функции вос­принимающего, задающего и сравнивающего органов. Как вос­принимающий орган он контролирует давление Р, преобразуя его в силу Fp. Задание требуемого давления в ресивере обеспечивается посредством силы F0. Как сравнивающий орган сильфон обеспе­чивает сравнение величин F0 и Fp, в результате чего, как отмеча­лось ранее, получается ΔF= F0 — Fp — сигнал рассогласования.

Динамические свойства объекта регулирования и элементов САР описываются следующей системой уравнений:

заслонкой

Физическая сущность переменных, входящих в уравнения, от­ражена выше в описании схемы САР. Параметры T0, T1, T2, T3 и К0, Кq, Кв, Кc, Кп, Ку, К3 — соответственно постоянные времени и ко­эффициенты передачи. Их размерности и значения по вариантам даны в таблице 2.5. Требуемое значение давления Р = 500 кПа.

Значения параметров элементов САР

Вари­антТ0,сКо КПа/ммТ1,сТ2,сКс мм/НКв Н/кПаКQ, Кпа*с/м 3ΔQC, м3/сКп В/ммКуТ3К3 Мм/В
1,30,20,0452,50,50,10,20,01
0, 250 ,042,50,5-0, 20,20,01
0,63,50,340,0222,50,50,30,20,01
4,80,250,0352,50,5-0,150,20,01
0,74,50,30,042,50,50,120, 90,01
0,83,50,180, 0252 ,50,5-0,20 ,20,01
0,44,40,250,032,50,50,110,20,01
0,655,50,20,022,50,5-0,120,20,01
0, 70, 40 ,0252,50,50,140,20,01
0,550,250,0352,50,5-0,140,20,01

На электрических станциях при производстве электроэнергии предъявляют определенные требования к стабильности частоты f генерируемой ЭДС. Частота f однозначно определяется угловой скоростью ω рабочего колеса гидротурбины. В связи с этим гид­ротурбины на электростанциях оснащают САР угловой скорости. На рисунке 2.5 показана схема одного из вариантов такой САР.

В данной системе объектом регулирования является гидротур­бина 1, регулируемой величиной — угловая скорость ω .Она при постоянном расходе воды изменяется в зависимости от нагрузки на валу турбины, т. е. от мощности Р, которая потребляется от ге­нератора 2 (с увеличением мощности угловая скорость снижается, с уменьшением — возрастает). Таким образом, мощность Р явля­ется внешним возмущающим воздействием на объекте регулиро­вания. Для регулирования угловой скорости предусмотрена за­слонка 3, с помощью которой изменяется расход воды через тур­бину. Он однозначно зависит от вертикального перемещения X заслонки. Следовательно, перемещение заслонки X можно рас­сматривать как регулирующее воздействие объекта регулирова­ния. Угловая скорость ω контролируется посредством тахогенератора 4, ЭДС Е которого сравнивается с задающим напряжением U0. Сигнал рассогласования Δ U через усилитель 5 управляет по­средством электродвигателя 6 и редуктора 7 заслонкой 3.

Рис. 2.5 Схема САР угловой скорости рабочего колеса гидротурбины

Динамические свойства элементов САР описываются следую­щей системой уравнений:

Формы записи линейных неоднородных дифференциальных уравнений

4.4. Формы записи линейных неоднородных дифференциальных уравнений

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

Рассмотрим операторную форму записи на примере. Пусть неоднородное линейное дифференциальное уравнение, описывающее динамику САУ, имеет вид

. (4.27)

Операцию дифференцирования обозначим оператором p,

.

Используя эти обозначения, (4.27) перепишется как

,

. (4.28)

— назовем собственным оператором;

— назовем операторами воздействия.

В итоге получим компактную форму записи в операторном виде

. (4.29)

Теперь получим форму записи неоднородного линейного дифференциального уравнения n-го порядка в преобразованиях Лапласа. Это уравнение в стандартном виде представлено как

. (4.30)

Используя свойства преобразования Лапласа, получим

.

, (4.31)

можем записать сразу решение дифференциального уравнения (4.30) в преобразованиях Лапласа

или ,

где изображение решения дифференциального уравнения или выходной величины y(t) равно произведению дробно-рациональной функции W(s) и изображения по Лапласу входной величины x(t).

уравнение ;

решение уравнения в изображениях

;

при .

Тогда решение исходного уравнение в преобразованиях Лапласа

,

или, используя таблицы обратного преобразования Лапласа, получим

4.5. Виды передаточных функций САУ

Отношение оператора воздействия к собственному оператору называется передаточной функцией в операторной форме.

Передаточная функция относительно входной величины x(t)

. (4.32)

Передаточная функция относительно возмущающего воздействия f(t)

(4.33)

Используя передаточные функции (4.32) и (4.33), уравнение (4.27) можно переписать

(4.34)

Уравнения (4.32)-(4.33) называют уравнениями в символической или операторной форме записи.

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

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

Найдем передаточные функции в форме изображений Лапласа для звена, описанного уравнением (4.27)

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

, (4.35)

где , , .

; (4.36)

W1(s) называют передаточной функцией по управляющему воздействию,

W2(s) – передаточной функцией по возмущающему воздействию.

Сравнивая передаточные функции (4.32) и (4.33) с (4.36), заметим, что они полностью совпадают при замене аргументов p на s. Это совпадение возможно только для стационарных систем. Если , , зависят от времени, то. динамические системы называются нестационарными.

Используя передаточные функции (4.36) и уравнение (4.35),можно записать решение исходного уравнения в преобразованиях Лапласа

(4.36)

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

Следует отметить, что линейные дифференциальные уравнения с постоянными коэффициентами не выше второго порядка, описывающие

поведение звеньев и САУ, записываются в стандартной форме. В этом случае все члены, содержащие выходную переменную и ее производные, записывают в левой части уравнения, все остальные члены — в правой, а коэффициент при выходной переменной делают равным единице.

Уравнение (4.27) в стандартной форме принимает вид:

,

где — имеют размерность времени и их называют постоянными времени;

и — передаточные коэффициенты.

Стандартную форму остальных звеньев рассмотрим ниже.

4.6. Характеристики САУ и их элементов

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

4.6.1. Частотные характеристики

Частотные характеристики с физической точки зрения являются

функциями параметров реакции САУ и ее элементов на гармонические воздействия. Для определения частотных характеристик рассмотрим решение уравнения (4.30) при

Решение уравнения (4.30) имеет вид

,

где ya(t) — общее решение однородного уравнения;

(t)- частное решение неоднородного уравнения.

Составляющая уa(t) определяет свободные движения системы или переходной процесс, который со временем затухает, если система устойчива. Составляющая yb(t) определяет вынужденное движение системы под воздействием x(t) . Найдем yb(t) при заданном гармоническом воздействии x(t). Согласно формуле Эйлера, представим

, (4.38)

где , . (4.39)

Найдем отдельно каждое из этих решений.

Для этого определим производные от x1(t) с учетом того, что .

Уравнение (4.30) после его преобразования по Лапласу и подстановки в правую часть выражений (4.40) примет вид

.

На основе свойств неоднородного линейного дифференциального уравнения, частное решение будем искать в виде

, (4.41)

где оператор не зависит от времени. Тогда, подставляя (4.41) в (4.40), получим

.

Откуда

. (4.42)

Функцию , которую можно получить на основе передаточной функции W(s) при подстановке в нее , называют частотной передаточной функцией. Частотная передаточная функция является комплексной функцией от действительной переменной , где — частота синусоиды входного сигнала x(t).

Функцию представим как

, (4.43)

, (4.44)

, (4.45)

Тогда частное решение неоднородного уравнения (4.30) при x1(t)

. (4.46)

Аналогично получаем решение при x2(t)

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

Частотную характеристику как комплексную функцию можно отобразить на комплексной плоскости в виде годографа (рис.4.3).

При заданной будет изображаться точкой С. Модуль вектора . Угол. Кривая, которую описывает конец вектopа ОС при изменении от 0 до называют амплитудно-фазовой час­тотной характеристикой (АФЧХ) (рис.4.3),

U()- вещественная частотная функция;

V() — мнимая частотная функция;

A()=/W(j)/- амплитудно-частотная функция;

()=arg/W(jw)/- фазо частотная функция.

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

Кроме перечисленных частотных характеристик, используют еще и логарифмические частотные характеристики. Логарифмическая амплитудно-частотная функция определяется как . График зависимости логарифмической амплитудно-частотной функции L() от логарифма частоты называют логарифмической амплитудной частотной характеристикой (ЛАЧХ).

При построении ЛАЧХ по оси абсцисс откладывают частоту в логарифмическом масштабе. Для построения ЛАЧХ по оси ординат откладывают в дБ величину 20lg A() .Для построения ЛФХ по оси ординат откладывают значение .

Шкалу частот по оси абсцисс разбивают в масштабе lg(), что позволяет на всем диапазоне частот выделить интервалы частот на октаву или на декаду.

Особенность lg масштаба состоит в том, что отрезок, изображающий октаву или декаду, имеет одну и ту же длину, не зависящую от частоты и равную соответственно lg 2 или I (pис.4.4):

Рис.4.4 • Логарифмический масштаб

где по оси абсцисс пишут не значение 1g(), а значение , по оси ординат откладывают L() =20 lg А() Логарифмической фазо-частотной характеристикoй (ЛФЧХ) называют график зависимости фазовой частотной функции от логарифма

частоты (при построении графиков по оси абсцисс пишут не значение, а ).

Единица измерения L()- децибел (дБ), единица измерения — декада, которая определяется как интервал, на котором частотa изменяется в 10 раз.

Ось ординат при построении ЛАЧХ строится через произвольную точку (обычно при =I).

Зависимость (4.42) можно представить как

,

(4.49)

Откуда возможно получить аналитические зависимости для определения частотных характеристик на основе (4.43) и (4.49):

— вещественная частотная функция;

— мнимая частотная функция;

— амплитудно-частотная функция;

— фазовая частотная функция.

5.2. Временные характеристики

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

Переходной функцией системы (звена) называют функцию, описывающую изменение выходной величины системы (звена) при подаче на ее вход единичного ступенчатого воздействия при нулевых начальных условиях. Переходную функцию обычно обозначают h(t). Иначе, переходная функция h(t) есть функция, описывающая реакцию системы (звена) на единичное ступенчатое воздействие при нулевых начальных условиях.

Аналитически единичное ступенчатое воздействие можно описать так называемой единичной функцией .

График переходной функции — кривая зависимости функции h(t) от времени t — называют переходной характеристикой.

Импульсной переходной или весовой функцией (функцией веса) системы (звена) называют функцию, описывающую реакцию системы (звена) на единичное импульсное воздействие при нулевых начальных условиях; обозначают эту функцию через (t). График импульсной переходной функции называют импульсной переходной характеристикой.

Переходную и импульсную переходную характеристику называют временными характеристиками.

При определении весовой функции было использовано понятие единичного импульса. Физически единичный импульс можно представить как очень узкий импульс, площадь которого равна единице. Математически он описывается функцией , которую называют дельта функцией. Дельта-функция и единичная функция относятся к так называемым обобщенным функциям. Теория обобщенных функций — сравнительно новый раздел функционального анализа. Отметим, что в рамках теории обобщенных функций любые встречающиеся в приложении функции обладают производными любого порядка. В частности, существует производная и от единичной функции, которая равна дельта — функции, [1(t)]1= . Обладает производными любого порядка и сама дельта-функция.

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

дельта-функция есть функция, которая обладает следующими свойствами —

(4.50)

(4.51)

Производные от дельта-функции можно определить по следующим соотношением:

(4.52)

(4.53)

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

Найдем изображение Лапласа от дельта-функции и ее производных. При этом преобразование Лaпласа будем трактовать как предельное соотношение

(4.54)

Используя соотношения (4.50)-(4,53), нетрудно получить

, . . (4.55)

Рассмотрим линейное неоднородное дифференциальное уравнение n- го порядка с постоянными коэффициентами в общем виде

(4.56)

В изображении Лапласа уравнение принимает вид

, (4.57)

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

Легко проверить, используя формулу (4.55) ,что уравнение (4.57) справедливо и в тех случаях, когда или .

В соответствии с определением весовой функции при переменная. Так как , то (4.57) можно переписать как

(4.58)

Таким образом передаточная функция равна изображению Лапласа от весовой функции и соответственно

(4.59)

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

Установим связь между весовой и переходной функциями.

Так как , то уравнение (4.57) при ,принимает вид:

(4.60)

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

(4.61)

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

(4.62)

Эта формула, как и уравнение (4.57), справедлива только при нулевых начальных условиях.

4.7. Элементарные звенья САУ и их характеристики

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

. (4.63)

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

Из курса алгебры на основании теоремы Безу известно, что полином произвольного порядка можно разложить на простые множители вида

, . (4.64)

Поэтому передаточную функцию (4.63) можно представить, как произведение простых множителей вида (4.64) и простых дробей вида

, , . (4.65)

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

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

Так как , , то для модуля и аргумента комплексного числа имеем

, .

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

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

Частотные и временные функции этого типового эвена имеют вид:

, , ,

, , , .

Ha рис. 4.5 представлены некоторые из характеристик пропорционального звена: амплитудно-фазовая частотная характеристика (4.5 а) — это точка К на действительной оси; фазовая частотная

20 lgK K

K U w t

Рис.4.5 Характеристики пропорционального звена

характеристика (или АФЧХ) совпадает с положительной осью частот; логарифмическая амплитудная частотная характеристика (рис. 4.56) параллельна оси частот и проходит на уровне. Переходная характеристика (рис.4.5в) параллельна оси времени и проходит на уровне .

Интегрирующее звено. Интегрирующим называют звено, которое описывается уравнением или передаточной функцией . Частотная передаточная функция .

Остальные частотные и временные функции имеют вид:

, , , ,

, , .

АФЧХ (рис.4.6а) интегрирующего звена совпадает с отрицательной мнимой полуосью. ЛФЧХ (рис.4.66) параллельна оси частот и проходит на уровне : сдвиг фазы не зависит от частоты и равен .

ЛАЧХ (рис.4.6б) — наклонная прямая, проходящая через точку с координатами и . Как видно из уравнения при увеличении частоты на I декаду ордината , уменьшается на 20 дБ. Поэтому наклон ЛАЧХ равен -20 дБ/дек (читается: минус двадцать децибел на декаду).

Переходная характеристика представляет собой прямую, проходящую через начало координат с угловым коэффициентом наклона, равным k. (рис.4.6в).

а) б) в)

jV U L(w) (w) h(t)

Рис 4.6 Характеристики интегрирующего звена

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

Частотные и временные функции этого звена имеют вид

, , , , ,

, , .

jV а) L(w) (w) б)

Рис.4.7 Характеристики дифференцирующего звена

АФЧХ (рис 4.7а) совпадает с положительной мнимой полуосью. ЛФЧХ (рис 4.7б) параллельна оси частот и проходит на уровне , то есть сдвиг фазы не зависит от частоты и равен /2.

ЛАЧХ есть прямая линия, проходящая через точку с координатами =1, и имеющая наклон 20 дБ/дек (читается: плюс двадцать

децибел на декаду): увеличивается на 20 дБ при увеличении частоты на одну декаду.

Апериодическое звено. Апериодическим эвеном первого порядка называют звено, которое описывается уравнением

(4.66)

или передаточной функцией

. (4.67)

Это звено также называют инерционным звеном первого порядка. Апериодическое звено в отличие от вышерассмотренных звеньев характеризуется двумя параметрами: постоянной времени T и передаточным коэффициентом k.

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

. (4.68)

Умножив числитель и знаменатель на комплексно-сопряженное знаменателю выражение, получим

, . (4.69)

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

Так как модуль числителя частотной передаточной функции (4.68) равен k, а модуль знаменателя ,то

(4.70)

Аргумент числителя равен нулю, а аргумент знаменателя . Поэтому

(4.71)

Решив дифференциальное уравнение (4.66) при и нулевом начальном условии , получим переходную характеристику . Весовая функция или импульсная переходная характеристика

.

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

ЛАЧХ представлена на рис 4.8б. На практике обычно ограничиваются построением так называемой асимптотической ЛАЧХ (ломаная линия на том же рис 4.86). В критических случаях, когда небольшая погрешность может повлиять на выводы о состоянии исследуемой системы, рассматривают точную ЛАЧХ. Впрочем, точную ЛАЧХ можно легко построить по асимптотической ЛАЧХ, если воспользоваться следующей зависимостью ( L — разность между асимптотической и точной ЛАЧХ):

T = 0,10 0,25 0,40 0,50 1,0 2,0 2,5 4,0 10,0

L= 0,04 0,25 0,62 0,96 3,0 0,96 0,62 0,25 0,04

Частоту , при которой пересекаются асимптоты, называют сопрягающей частотой. Точная и асимптотическая ЛАЧХ

Рио.4.8 Характеристики апериодического звена

наиболее сильно отличаются при сопрягающей частоте; отклонение при этой частоте примерно 3 дБ.

Уравнение асимптотической ЛАЧХ имеет вид:

Оно получается из уравнения (4.71), если в нем под корнем при пренебречь первым слагаемым, а при вторым слагаемым.

Согласно полученному уравнению, асимптотическую ЛАЧХ можно строить следующим образом: на уровне частоты провести прямую, параллельно оси частот, а далее через точку с координатами и — прямую под наклоном — -20 дБ/дек.

По АФЧХ или ЛАЧХ легко определить параметры Т и k аперио­дического звена (рис.4.86).

ЛФЧХ изображена на рис. 4.86. Эта характеристика асимптотически стремится к нулю при и к при . При фазо — частотная функция принимает значение —, то есть . ЛФЧХ всех апериодических звеньев имеют одинаковую форму и могут быть получены на основе одной характеристики параллельным сдвигом вдоль оси частот влево или вправо в зависимости от значения постоянной времени T. Поэтому для построения ЛФЧХ апериодического звена можно воспользоваться шаблоном, представленном на рис.4.8г.

Переходная характеристика апериодического звена (рис.4.8в) представляет собой экспоненциальную кривую, по которой можно определить параметры этого звена: передаточный коэффициент k определяется по установившемуся значению ; постоянная времени T равна значению t, соответствующему точке пересечения касательной, построенной на переходной характеристике в начале координат, с ее асимптотой (рис 4.8в).

Форсирующее звено. Форсирующим звеном или форсирующим звеном первого порядка называют звено, которое описывается уравнением

,

или передаточной функцией

.

Это звено, как и апериодическое, характеризуется двумя параметрами: постоянной времени T и передаточным коэффициентом k.

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

.

Остальные частотные и временные функции имеют вид:

, , , ,

, , .

АФЧХ есть прямая, параллельная мнимой оси и пересекающая действительную ось в точке U=k.(рис. 4.9а). Как и в случае апериодического звена, на практике ограничиваются построением асимптотической ЛАЧХ. Частоту , соответствующую точке излома этой характеристики, называют сопрягающей частотой. Асимптотическая ЛАЧХ при параллельна оси частот и пересекает ось ординат при , а при имеет наклон +20дБ/дек.

ЛФЧХ форсирующего звена можно получить зеркальным отображением относительно оси частот ЛФЧХ апериодического звена и для ее построения можно воспользоваться тем же шаблоном и номограммой, которые используются для построения последней.

Колебательное, консервативное и апериодическое второго порядка звенья. Звено, которое можно описать уравнением

(4.72)

или в другой форме

(4.73)

где, , .

Передаточная функция этого звена

(4.74)

Это звено является колебательным, если ;-консервативным, если ;- апериодическим звеном второго порядка, если . Коэффициент называют коэффициентом демпфирования.

Колебательное звено . Частотная передаточная функция этого звена

.

Умножив числитель и знаменатель на комплексно-сопряженное выражение, получим вещественную и мнимую частотные функции колебательного звена:

,

Фазовая частотная функция, как это видно из АФЧХ (рис 4.10б), изменяется монотонно от 0 до — и выражается формулой

(4.75)

ЛФЧХ (рис.410б) при асимптотически стремится к оси частот, а при к прямой . Ее можно построить с помощью шаблона. Но для этого необходимо иметь набор шаблонов, соответствующих различным значениям коэффициента демпфирования.

Амплитудная частотная функция

и логарифмическая амплитудно-частотная функция

.

Уравнение асимптотической ЛФЧX имеет вид

(4.75)

где — сопрягающая частота. Асимптотическая ЛАЧХ (рис.4.106) при параллельна оси частот, а при имеет наклон — -40 дБ/дек.

Рис. 4.10 .Характеристики колебательного звена

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

,

, ,

.

.

По переходной характеристике (рис.4.10в) можно определить параметры колебательного звена следующим образом.


источники:

http://megaobuchalka.ru/5/12593.html

http://pandia.ru/text/78/455/43428.php