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

книги / Теория автоматического управления техническими системами

..pdf
Скачиваний:
15
Добавлен:
12.11.2023
Размер:
15.06 Mб
Скачать

Рис. 5.18. Выделение областей устой­ чивости, £>[1], D[2], ДЗ]—области, соответствующие различным значе­ ниямкорией

на один

больше, так

как при этомпроисходит движение в сторону штри­

ховки.

Следовательно,

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

все три

корня лежат слева от мнимой оси. Здесь

существеннытолько дей­

ствительны значения к, принадлежащие области

устойчивости. Они опре­

деляются отрезком оси

U, лежащим внутри области

D[3]. Условиюустой­

чивости

рассматриваемой системыотвечаю значения

0<£<1.

Контрольные вопросы

1.Дайте строгое определение устойчивости динамической системы (САР) по Ляпунову. Что такое критерий устойчиво­ сти?

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

3.В чем состоит смысл принципа аргумента?

4.Сформулируйте и проиллюстрируйте на примерах частот­ ные критерии Михайлова и Найквиста, используя АФХ.

5.Сформулируйте и проиллюстрируйте частотный критерий устойчивости Найквиста с использованием ЛЧХ.

6.Что такое запасы устойчивости по модулю и фазе? Что они характеризуют?

7.Что такое полоса пропускания системы?

8.Как построить области устойчивости САР в плоскости одного параметра? Что дает проектировщику САР выделение областей устойчивости?

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

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

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

частотному методу анализа, который базируется на связи ве­

щественной частотной характеристики САР, вычисленной при

разомкнутой цепи обратной связи, с переходной функцией си­

стемы. Это позволяет использовать логарифмические характе­

ристики разомкнутой

САР, необходимые для анализа как

устойчивости по амплитудно-фазовому критерию, так и качест­

ва регулирования.

САР в установившемся режиме при

Для оценки точности

действии детерминированных сигналов и возмущений исполь­

зуют метод коэффициентов ошибок.

6.1. Методы анализа качества САР

Качество САР или качество регулирования являются обоб­ щенными показателями системы, характеризующими: во-пер­ вых, переходный процесс; во-вторых, статическую точность си­ стемы при некоторых типовых воздействиях; в-третьих, точность

соответствующие движению с постоянными скоростью и уско­

рением, а также гармонические воздействия и др.).

Для оценки стохастических САР применяют вероятностные

методы, определяя

динамическую точность САР по значению

среднеквадратической

ошибки или среднеквадратического от­

клонения регулируемой переменной. Вопросы анализа и синте­

при медленно изменяющихся сигналах [19, 12].

Для оценки качества детерминированных САР выбирают ти­

повые (тестовые) воздействия, являющиеся наиболее неблаго­ приятными или характерными для данной САР (например, сту­ пенчатые функции: единичные; дельта-функции; воздействия,

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

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

тальным или расчетным кривым переходного процесса в САР; косвенно, т. е. по каким-либо динамическим характеристикам или параметрам системы. Оценки, полученные непосредственно по кривой переходного процесса САР, называют прямыми (пер­ вичными) показателями качества, а оценки, определяемые дру­ гими способами —косвенными показателями.

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

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

Основные методы определения переходных процессов и ана­ лиза качества линейных САР: 1) частотный; 2) корневого годо­ графа (распределения и траектории нулей и полюсов переда­ точной функции системы); 3) логарифмического корневого го­ дографа; 4) интегральных оценок. Рассмотрим их.

1. Частотный метод. Основан на рассмотрении преобразова­

ния Лапласа X(s) для регулируемой величины при чисто мни­

мых значениях аргумента s=jcо, а также на связи, существую­

щей между частотными характеристиками замкнутой

(разомк­

нутой) системы

и переходным

процессом. Одно из

основных

различий между прямым методом анализа переходных процес­

сов, основанным на преобразовании Лапласа, и частотным за­

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

с вычислением

корней характеристического уравнения системы,

а второй (как

и частотный

метод анализа устойчивости) —

графоаналитическим, не требующим вычисления корней. .При

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

ференциальных уравнений всей системы в целом или отдель­ ных ее элементов. Этот метод позволяет: а) проводить полный анализ динамики, а также решать многие вопросы синтеза корректирующих устройств или регулятора САР; б) учитывать особенность САР, заключающуюся в том, что анализ систем в разомкнутом состоянии обычно проще, чем в замкнутом; в) осу­ ществлять анализ устойчивости, качества, а также переходных процессов и систем любого порядка (как одно-, так и многокон­ турных, содержащих не только сосредоточенные, но и распреде­ ленные параметры); г) решать вопросы анализа и синтеза си­ стем при непрерывно изменяющихся воздействиях.

2. Метод корневого годографа. Основан на связи между расположением нулей и полюсов передаточной функции системы в замкнутом и разомкнутом состоянии и на изучении их пере­ мещения на плоскости s при изменении параметров системы. Если в процессе проектирования САР были получены характе­ ристики переходного процесса, не соответствующие предъявляе­ мым техническим требованиям, то изменением положения кор­ ней характеристического уравнения можно изменить показате­ ли качества. Метод корневого годографа позволяет проанализи­ ровать, как меняются корни уравнения при изменении от -—<»до +оо линейно-входящего параметра системы, и показы­ вает, как нужно изменить эти корни для получения требуемых характеристик.

3. Метод логарифмического корневого годографа. Основан на анализе свойств замкнутой системы по логарифмическим комплексным частотным характеристикам разомкнутой систе­ мы, т. е. характеристикам, построенным для комплексных зна­ чений аргумента s=o+/cù в выражении для передаточной функ­ ции \P(s).

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

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

6.2. Частотный метод определения переходных функций линейных непрерывных САР

Первичные показатели качества. Поведение системы в пере­

ходном процессе, вызванном

типовым воздействием, стремя­

щимся с течением времени к постоянному установившемуся значению, можно охарактеризовать при помощи некоторых па­ раметров, называемых первичными показателями качества САР. В частности, при ступенчатом воздействии и нулевых начальных условиях переходный режим системы характеризуется следую­ щими показателями качества (рис. 6.1):

Рис. 6.1. Определение показателей качества САР

а) установившимся значением лгуСт=х(°°) (т. е. статическим отклонением регулируемой переменной x{t), которое определяет статическую ошибку (точность) системы

8уст==ХВх—Хуст, б) временем регулирования ТР(т. е. временем переходного

процесса Тая), которое характеризует быстродействие системы.

Время Гр определим как наименьшее значение интервала вре­ мени tp—to, отсчитываемого от начала приложения воздействия /0до момента tp, после которого

Ix{t)-Xycrl^A»

где А—заданная малая постоянная величина (значение Аобыч­

но назначают равным 0,05 хуст);

 

в) положительным перерегулированием хтат (т. е. макси­

мальным динамическим

отклонением) регулируемой

величины:

. gi=*maxr-*ycT.jQQ^;

 

 

 

*уст

N величины x{t) в течение

времени

г) числом колебаний

переходного процесса, т. е. в

интервале tp—to.

 

Три показателя: хусг,

Гр,

а (или N) —определяют качество

переходных процессов минимально-фазовых систем, т. е. систем с минимально-фазовыми или «левыми» передаточными нулями

и полюсами.

 

 

 

 

 

 

 

Показатели качества выбирают в зависимости от техниче­

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

необходимым качеством, если она удовлетворяет заданным ус­

ловиям качества, а переходный про­

'//////////,

 

цесс не выходит из области допусти­

 

мых значений

(на рис. 6.2

граница

 

lL~

области допустимых значений обозна­

 

чена штриховкой).

метод

опре­

 

77777

Графоаналитический

 

деления переходного

процесса,

рас­

1

1

 

смотренный далее, может быть реали­

 

зован на ЭВМ с подсистемой

визу­

 

 

t

ализации или с использованием обыч­

Рис. 6.2. Область

ных графических средств. В первом

допу­

случае данный

метод является

осно­

стимых

значений

пере­

вой для математического

и

про­

ходной функции

граммного обеспечения ЭВМ.

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

возможно решение различных задач анализа и синтеза систем.

Математической основой для частотного метода определения и анализа переходных процессов является преобразование

Фурье. Оно позволяет получить на основании дифференциаль­

ных уравнений с учетом начальных условий и приложенных к системе воздействий (или с использованием экспериментальных

данных) некоторые функции, называемые обобщенными частот­

ными характеристиками САР.

координаты

Можно показать, что изменение обобщенной

или регулируемой переменной x(t) для широкого

класса воз­

действий g(0 как при ненулевых, так и при нулевых началь­

ных условиях определяют выражениями [19]:

 

00

(6.1)

JC(t)-*„(*)+ !-$ Rr(©) cos(ùtd(ù ;

je(t)=xn(t)—- ^ Sr (cd) sin ©tda,

(6.2)

d

 

где xn(t) —нерегулярная часть функции x{t), которая при Н-оо является расходящейся, или периодической, функцией, либо по­ стоянной величиной; R,(a) и Sr(©) —обобщенная вещественная и мнимая частотные характеристики процесса, получающиеся, если преобразование Лапласа для функции X-(s) при s=/© представить в виде

*r(/®)=*r(®)+/S,(®).

Здесь Xr (s) —регулярная часть функции X(s), т. е. когда все полюсы ее расположены в левой полуплоскости и limj:r(9=0.

Формулы (6.1) и (6.2) можно переписать в виде

&x(t)=x(t)—xn(t)= ^Rr(a>)cost<ùd<ù\

(6.3)

о

 

Дx(t)=x(t)—xn(t) = - д2 Jsr(©)sin*©d©.

(6.4)

 

 

о

 

регу­

Левы части уравнений (6.3) и(6.4) представляю собой отклонения

лируемой

переменной

от нерегулярной составляющей xn(t) переходного

процесса.

Вслучае,

если функция Л(s) не содержит особенностей

во

всей

правой полуплоскости ина мнимой оси,

 

 

ИшsX(s)=limх (f)=0

 

 

j-*0

<-+о

 

(6.3) и

ифункция x(t) имеет только регулярнуючасть. Поэтому выражения

(6.4) сводят к виду

 

 

 

 

Rr (ш) cos ttodiВ,

 

(6.5)

или

о

 

 

 

2_

 

 

 

x{t).

sin t(ûd(ù, t>0.

 

(6.6)

л

 

Если Rr{<ù) совпадает с вещественной частью#(©) выражения X(/со), а 5г((о) —с мнимой частьюS(<o), то обобщенны частотные характеристики определяют следующим выражениями:

R(©)=P(a)i>»(w)-Q(©)Q„(ш)+Рв(©); 5(и)=Р(со)Q,(©)-fQ(со)рв(хо)+Qh(и),

где Р(ы) и Q(<o)—вещественная и мнимая частотные характеристики си­ стемысоответственно; />,(©) и Q«(со) —вещественная и мнимая частотные характеристики воздействия; Рп(<о) и Qn(со)—вещественная и мнимая частотные характеристики начальных условий.

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

а все остальные полюсы.расположеныв левой s-полуплоскостн, а также в

случае, когда

такой

полюс

отсутствует, формулы(6.5) и (6.6)

имеют вид

2

F R (©)

sin t(ùd<ù> t>0>

(6.7)

у (0=~ )

 

«ли

2

(*5 (®)

 

 

 

cos t<ùd<ù, t>0,

(6.8)

x(t)=R (0) + —^

 

 

где /? (0)—значение

о

 

 

 

R (ю) при <0=0.

 

При единичном ступенчатомвоздействии g(t)=l(t) инулевых началь­ ных условиях функциями R((ù) и S(at), входящим в формулы(6.7) и (6.8), являются собственны вещественная />(<в) и мнимая Q(o>) частотные харак­ теристики замкнутой системы. Типовые вещественная Р(а) и мнимая Q(<o) частотные характеристики САР с единичной обратной связьюпоказанына

Рис. 6.3. Типовые вещественная (а) и мнимая (б) частотные характеристики САР с единичной обратной связью

рис. 6.3. Функция x(t) в выражениях (6.7) и (6.8) представляет собой пере­ ходнуюфункциюh(t) системы. Поэтому формулы(6.7) и(6.8) примут вид

Л(О=4"^о

 

s,n t(ùd<ùt

 

(6.9)

или

2

СО(со)

 

 

 

в.

(6.10)

h{t)=P (0) + —у

m- cos

 

 

о

 

процесса

Таким образом, первымшагомпри вычислении переходного

по формулам

(6.7)

и

(6.8) является

определение частотных характеристик

Р(ю) или Q(со)

системы.

 

 

Определение вещественной и мнимой частотных характери­ стик замкнутой системы по ЛЧХ разомкнутой системы. Взаим­ ную связь между АФЧХ замкнутой и разомкнутой САР опреде­

ляют выражением

 

ф(/ш)

W(Jcù)

(6.11)

1+W{j<ùY

Как было показано ранее, АФЧХ разомкнутой системы №(/ю) •может быть представлена соответствующими амплитудной H((ù) и фазовой 0(©) частотными характеристиками следующим образом:

Подставляя выражение (6.12) в (6.11), получим

(6.12)

,g

^ И cos 9(о))+/Я(со) sin 8 (ю)__р(

. Q( .

' 1+//(со) cos 0 (со) н- у//(со) sia 0 (to)

'

/+74:V *'

'

Отделяя в правой части выражения (6.13) вещественную часть

P(ûï) от мнимой Q(ш), найдем

 

Р(0)):

Нг (ш) + Н(ш) cos9 (m) .

(6.14)

Нг (©) +2Н(со) cos 0 (<в) + Г

^

Н (ш) sin Q(со)

(6Л5>

Н2(и) +2Н (.си) cos 0 ((о) + Г

Геометрическое место точек Р(©)=const=РСсогласно (6.14) определяют уравнением

Нг (и) +Н (a)cos 8 (а) _р

(W) +Ш(to) cos 0 (ш) + 1~

или

(6.16)

Я(o)cos0(cû)(1—2РС) =Я2(©) (Р0-1)+Ре.

Точно так же геометрическое место точек в соответствии с вы­ ражением (6.15)

Q(<û) = const=QC,

 

определяют уравнением

 

Н(to)sin 9 (<а)___ _

 

Я* (©)+2Н(to) cos 0 (со) +1

 

ИЛИ

(6.17)

Я(©)sin 0(со)—2<2еЯ(и)cos0(ш= Qc+QCH2(q).

Номограммы можно построить: в соответствии с выражени­ ем (6.16)—для определения вещественной частотной характе­ ристики замкнутой системы [18, 20]; в соответствии с выраже­ нием (6.17)—для мнимой. Эти номограммы строят в коорди­ натах фаза —амплитуда, т. е. в следующих единицах отсчета: град —по оси абсцисс, дБ —по оси ординат. Параметры Р и Qномограмм—безразмерные. На номограммы наносят точки,, соответствующие значениям амплитудной характеристики L{(ù)] разомкнутой системы (в дБ) и фазовой характеристики 0(g)) (в град) для выбранной частоты ©<, и определяют веществен ную P((ù) и мнимую Q(©) частотные характеристики, соответст­ вующие частоте со,-.

Определение вещественной Р(©) и мнимой Q(co) частотных характеристик замкнутой системы по логарифмическим ампли­ тудной А(ш) и фазовой ср(со) частотным характеристикам замк­ нутой системы. При вычислении переходных процессов в САР в случае, когда возмущающее воздействие может быть приложе­ но в любой точке системы, а также для системы с неединичной обратной связью, вещественную и мнимую частотные характе­ ристики P((ù) и Q{(ù) следует определять по частотным харак­ теристикам замкнутой системы:

Р(<а)=Л(о)со5ф(ш);

(6.18)

Q(©) =A(©)sinqp(o).

.(6.19).

Р(а>) пользуются номограммой, вычисленной по формуле (6.18), при определении мнимой частотной характеристики Q(co) —

номограммой, вычисленной по формуле (6.19). Эти номограммы

построены в координатах фаза —амплитуда замкнутой системы Фи LA соответственно. Параметры Р и Qномограмм —без­ размерные.

Имея зависимость LA(<ù) от ф(ш), построенную на кальке в масштабе номограммы, и наложив ее на номограмму, можно

найти вещественную и мнимую частотные характеристики замк­ нутой системы.

Пересечение кривой LA(ш) =/[ф(о))] с кривой номограммы Реили Qe при заданной частоте ш=©, означает, что при со=

=ю, вещественная частотная характеристика Р(ю) имеет зна­

чение Р{, а мнимая—значение Q,. Положительным фазовым

углам соответствуют положительные значения Q(<a), а отрица­

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

 

6.3. Определение переходных процессов методом

трапецеидальных частотных

характеристик

Формулы (6.9) и (6.10) применимы

для вычисления пере­

ходных процессов (при нулевых начальных условиях), стремя­

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

Рассмотрим задачу вычисления переходной функции (пере­ ходной характеристики) САР, т. е. задачу нахождения переход­ ного процесса как реакции системы на единичное ступенчатое воздействие.

Первый шаг при вычислении переходной характеристики по формуле (6.9) или (6.10) заключается в нахождении веществен­ ной P(<ù) или мнимой Q((ü) характеристики системы. Различные варианты задания и способы определения функций Р(ш) и

Q(со) были рассмотрены ранее.

 

Следующим шагом является вычисление интеграла в выра­

жениях

(6.9) и (6.10). Для приближенного нахождения пере­

ходной

характеристики САР применяют

графоаналитический

метод, заключающийся в аппроксимации

функций Р(ю) или

Q(©), заданных в виде графиков, трапецеидальными частотны­

Чтобы облегчить процесс построения функций Р(а>) и ф(ш)

по ЛЧХ LA((ù) и ф(со), применяют специальные номограммы.

Для нахождения вещественной частотной характеристики

ми характеристиками. Метод может быть реализован либо на ЭВМс графическим дисплеем в диалоговом режиме, либо с ис­ пользованием обычных графических средств.

Типовая трапецеидальная частотная вещественная характе­ ристика. Эта характеристика (рис. 6.4) определяется следующи­ ми параметрами: высотой г0; интервалом равномерного пропу­ скания частот ©d; интервалом пропускания частот ©п. Отноше-

дальная частотная характери­

Рис. 6.5. Аппроксимация

вещественной

частотной

стика

характеристики САР тра­

пецеидальными

частот­

 

ными характеристиками

ние ©d/©n=x характеризует наклон боковой стороны типовой

трапецеидальной частотной характеристики.

 

Если го=1, ©„=1, то трапецеидальную частотную характе­

ристику называют единичной. Коэффициент наклона этой харак­

теристики может быть любым (от 0 до 1). Единичной вещест­

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

переходная

функция МО-

 

 

Вычисленные значения h в функции времени t для различных

х приведены в таблице /^-функций [см.: 2, 10, 19]. Зная коэф­

фициент наклона х, по этой таблице можно

определить hHи

гтабл. Для перехода от ^-функции процесса

(на выходе САР),

соответствующего единичной частотной

характеристике, к за­

данной следует значение функции МО

умножить на г0, а зна­

чение аргумента *табл>найденное по таблице, разделить на ©„.

Так, например, если r0=5, ©d=15, ©„=20, то x=©d/©„= =0,75; по таблице Мфункций для х=0,75 находим значения

^табл» МО и составляем табл. 6.1.

 

Таблица 6.1

 

 

 

 

Пример определения *дейсто. и h(t)__________

*табл

МО

*табл

*(О-МОг.

'действ" ^

 

 

ши

 

0,0

0,0000

0,00

0,000

1,0

0,5344

0,05

2,672

2,0

0,9383

0,10

4,691

Далееопределяем /действ и h(t)

 

Аппроксимация вещественной частотной характеристики тра­

пецеидальными характеристиками. Кривая Р(©) может

быть

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

трапе-