Теория систем автоматического управления. В.А. Бесекерский, Е.П. Попов, 1975. Теория систем автоматического управления. В.А. Бесекерский, Е.П.. 3 Примеры непрерывных автоматических систем
 Скачать 25.93 Mb.
|
|
u >u 1 ) скорость р х имеет тот те знак те. отклонение возрастает по величине) либо когда скорости р х имеет противоположный знак, но мала (| v | ). Система с переменной структурой. Как уже указывалось вначале книги (§ 2.3), системы с переменной структурой содержат в себе специальное переключающее устройство для изменения структуры регулятора, которое срабатывает в зависимости от размеров и знаков входных величин. Примеры переключающих устройств приведены схематически на рис. 16.27, где КЭ — ключевой элемент, БИС — блок изменения структуры. Уравнение принято [42] записывать в виде (16.70) Функция Ψ может строиться по-разному. Например риса) Для случая, указанного на рис. 2.9 и 2.10 (16.72) Под символами и β могут также иметься ввиду различные выражения в простейшем случае постоянные (16.73) в другом случае (16.74) и любые другие, в том числе и нелинейные. Основная же характерная нелинейность здесь состоит в самом факте автоматического переключения в зависимости от состояния входных величин. ГЛАВА ТОЧНЫЕ МЕТОДЫ ИССЛЕДОВАНИЯ УСТОЙЧИВОСТИ И АВТОКОЛЕБАНИЙ § 17.1. Фазовые траектории и метод точечных преобразований Понятие о фазовом пространстве, о фазовых траекториях и их типах было уже дано выше. В данном параграфе на примерах построения фазовых траекторий для простейших систем второго порядка будут проиллюстрированы некоторые важные особенности процессов в нелинейных системах автоматического регулирования. Пример 1. Возьмем систему автоматического регулирования с объектом без самовыравнивания и с приводом регулирующего органа, имеющим постоянную скорость. Уравнение регулируемого объекта без самовыравнивания будет (17.1) Для регулятора без массы и демпфера с жесткой обратной связью, те. при ζ ξ ξ η σ ϕ δη = − = − = , , , получим (17.2) где ζ ξ σ ϕ η , , , , — относительные изменения регулируемой величины, смещений чувствительного элемента, регулирующего органа, элемента обратной связи и управляющего золотника риса коэффициент. Привод регулирующего органа пусть имеет постоянную скорость в двух вариантах 1) с мгновенным переключением (рис. 16.22, ж) при переходе управляющего элемента (золотника, струйной трубки) через нейтральное положение ( σ = 0); 2) с зоной нечувствительности рис. 16.22, з) вследствие наличия перекрытия золотника или струйной трубки. В первом случае уравнение привода регулирующего органа будет (17.3) а во втором (17.4) Возьмем фазовую плоскость (х, у, приняв (17.5) Из уравнений (17.1), (17.2) и (17.5) имеем (17.6) Следовательно, переключения привода в первом варианте ( σ = 0) будут иметь место при (17.7) что соответствует прямой АВ (риса) на фазовой плоскости, причем согласно (17.16) значениям σ >0 соответствует часть плоскости слева от прямой АВ, а σ < 0 — справа. На основании первого из соотношений (17.6) с учетом (17.3) при σ < 0 получаем (17-8) а из (17.5) (17.9) откуда находим уравнения фазовых траекторий (17.10) или, после интегрирования, Это есть семейство парабол, показанное на риса справа от линии АВ (они симметричны относительно оси х. Так как (17.8) и (17.9) являются проекциями скорости V изображающей точки М на оси хи у, то имеем v у < О, а знак v x совпадает со знаком у. В соответствии с этим на риса укажем стрелочками направление движения изображающей точки М по фазовым траекториям. Аналогичным путем легко строятся параболы слева от прямой АВ. В результате, как видно из общего расположения фазовых траекторий (риса, получается устойчивая система с затухающим колебательным переходным процессом. Но число колебаний будет конечным. В самом деле, здесь имеется особый отрезок Св который вливаются все фазовые траектории. Чтобы выявить поведение системы на этом отрезке, вспомним, что для него согласно (17.7) и (17.5) Следовательно, попав на отрезок С, изображающая точка не может с него уйти, и система будет апериодически приближаться к установившемуся состоянию, те. изображающая точка будет сползать по отрезку С к началу координат 0. Таким образом, имевший место вначале колебательный переходный процесс после конечного числа колебаний вырождается в этот так называемый скользящий процесс. Крайние точки особого отрезка С определяются, очевидно, как точки, в которых прямая АВ касается одной из парабол соответственно правого и левого семейств. Поэтому, подставив значения dx dy изв выражение (17.10), найдем точку С По найденной картине расположения фазовых траекторий можно качественно представить себе кривую переходного процесса ϕ (t) при любых начальных условиях. Начальными условиями определяется начальное положение изображающей точки Ми тем самым — определенная фазовая траектория, иллюстрирующая протекание процесса. Она показывает (риса) максимальное отклонение регулируемой величины фтаХ1 максимальную скорость (р ϕ ) mах ; а также все последующие отклонения, число колебаний и т. п. Рассмотрим теперь туже систему, нос учетом зоны нечувствительности. В этом случае переключениям привода (при σ = -b и σ = +b) на фазовой плоскости соответствуют согласно (17.6) две наклонные прямые (рис. 17.1, б Между этими прямыми | σ | < b, правее их σ σ | < b из (17.4), (17.6) и (17.5) получаем откуда (при упрямые, параллельные оси х в полосе АВ на рис. 17.1, б. При | σ | > b получим прежние параболы. В результате снова система оказывается устойчивой и имеет колебательный переходный процесс, но вместо особой точки О получаем особый отрезок ух, те. установившееся состояние определяется неоднозначно. Это соответствует тому, 'что регулятор может находиться в равновесии в любом месте внутри зоны нечувствительности. Здесь точно также возможен скользящий процесс, как x в случае риса. В данном примере система оказывается устойчивой при любых значениях параметров и при любых начальных условиях. Однако здесь для получения системы второго порядка была проведена грубая идеализация уравнений регулятора (пренебрежение массой и демпфированием. Пример 2 . Допустим, что требуется стабилизировать угловое положение некоторого тела, когда сопротивлением среды его вращению можно пренебречь. Уравнение объекта будет (17.11) где J— момент инерции тела, ϕ — угол поворота тела, w — его угловая скорость, М — управляющий момент со стороны исполнительного органа системы стабилизации. Уравнение регулятора (системы стабилизации) запишем в виде (17.12) где М — постоянная положительная величина, Ф ( ϕ ,w) — нелинейный закон регулирования, осуществляемый при помощи логического устройств потому же простейшему принципу, что и на рис. 16.26, стой лишь разницей, что по углу ϕ фазовая плоскость ограничена значениями + π итак как это составляет один полный оборот тела (рис. 17.2). Изобразим процесс регулирования на фазовой плоскости. Уравнение всей системы согласно (17.11) и (17.12) будет (17.13), где обозначено причем с имеет физический смысл величины углового ускорения, сообщаемого данному телу постоянным моментом М 1 Умножив почленно уравнение (17.13) на выражение получим дифференциальное уравнение фазовой траектории (17.14) Это уравнение легко интегрируется внутри участков, на которых Ф = const. В результате для каждого отдельно взятого участка уравнение фазовой траектории будет (17.15) где ни н — значения фи сов начальной точке данного участка. Зададим начальные условия процесса Для данной начальной точки процесса (см. рис. 17.2) имеем Ф = 0. Поэтому на первом участке процесса согласно (17.15) уравнение фазовой траектории будет Этот участок движения с постоянной скоростью заканчивается в точке 1 (рис. 17.2), где происходит включение исполнительного органа (Ф = -1). Следовательно, для второго участка процесса (после точки 1) из (17.15) получим уравнение фазовой траектории (17.16) так как в начальной точке 1 этого участка н = b, w н = w 0 . Фазовая траектория (17.16) — парабола, ось которой совпадает с координатной осью ϕ . Это соответствует равнозамедленному движению ( c dt dw − = ) Изображая параболу графически, доводим ее до границы ϕ = π (участок 1—2 на рис. 17.2), причем в точке 2 согласно (17.16) (17.17) Это значение переносим в точку 2' (для вращающегося тела π ϕ ± = это одна и та же точка. Здесь происходит выключение исполнительного органа (Ф = 0). Поэтому дальнейшее движение согласно (17.15) пойдет с постоянной скоростью до точки 3 (рис. 17.2). Таким образом, в рассмотренной начальной части процесса регулирования тело совершило один полный оборот, нов конце этого оборота скорость вращения его стала меньше начальной. В точке 3 снова включается исполнительный орган (Ф = - 1), в результате чего фазовая траектория будет (17.18) так как в точке 3 b = ϕ , w н = w 2 . Допустим, что соответствующая уравнению (17.18) парабола 3 — 4 не доходит до границы π ϕ = . Это означает, что тело больше не совершит полного оборота, а начнет (сточки А) возвращаться в сторону нулевого положения. В точке 4 (рис. 17.2) имеем скорость w 4 =- b 1 . Следовательно, из (17.18) угловая координата ее будет где w 2 определяется по формуле (17.17). Дальше (4 — 5) процесс пойдет-с постоянной скоростью (так как Ф = 0), после чего тело войдет в установившийся автоколебательный режим, определяемый предельным циклом (5—6—7—8). Уравнение параболы 7—8 согласно (17.15) будет Отсюда амплитуда угловых автоколебаний а ф, как значение ф прибудет, а амплитуда колебаний скорости Она равна зоне нечувствительности измерителя угловой скорости b 1 , в то время как амплитуда угловых колебаний (17.19) несколько больше зоны нечувствительности измерителя угла b. Период автоколебаний t а можно вычислить как сумму времен где t хол и t раб — времена участков (6—7) + (8—5) и (5—6) + (7—8) соответственно. По законам равномерного и равнозамедленного движений соответственно получаем Итак, установившийся режим стабилизации в данной системе является автоколебательным. Однако уравнение системы (17.13) справедливо только для идеальной системы стабилизации. Всякое реально имеющееся запаздывание в работе усилительно- преобразовательного и исполнительного устройств приведет к увеличению амплитуд автоколебаний по сравнению с полученными здесь значениями. Решение задачи с учетом постоянных времени системы управления будет дано в следующей главе. Пример 3 . Уравнения системы автоматического регулирования курса водяной торпеды в упрощенном варианте имеют вид линейная часть (16.40) и (16.41), те) и нелинейное звено (возьмем сначала один случай) (17.21) Покажем, что здесь равновесное установившееся состояние системы с постоянным значением ψ = 0 неустойчиво, но будет иметь место устойчивый автоколебательный процесс. Возьмем фазовую плоскость (х, ус координатами х = - ψ , у = р (угол отклонения и угловая скорость отклонения оси торпеды от заданного курса. Уравнения (17.20) и (17.21) перепишутся в виде (17.22) Из сравнения этих уравнений с упрощенными уравнениями системы регулирования температуры в конце § 16.1 видна их полная аналогия. Поэтому здесь, также как ив случае рис. 16.15, установившийся процесс движения торпеды будет автоколебательным, причем картина фазовых траекторий будет иметь вид, показанный на риса При этом кривая АВ предельного цикла, соответствующая автоколебательному процессу, определяется из уравнения (16.31) с таким значением произвольной постоянной C 1 , чтобы выполнялось условие (17.23) так как только в этом случае и получится замкнутая кривая предельного цикла АВD (риса. Определив таким образом С найдем амплитуду автоколебаний а как значение х при у = 0, те. согласно (16.31) Значения же (17.23) дают амплитуду q колебаний скорости у. Можно все это определять и графически прямо по чертежу (риса. Период автоколебаний остается неизвестным. Введем теперь в характеристику нелинейного звена (рулевой машинки) зону нечувствительности, как показано на рис. 17.3, б, в. Так, на том участке характеристики рис. 17.3, 6), где = δ 0, из (17.22) следует, что что соответствует наклонным прямым внутри полосы ЕЕ на фазовой плоскости (рис. 17.3, б. Аналогичная полоса НGG 1 Н 1 будет ив нижней части плоскости. Все остальное заполняется такими же кривыми, как на риса. В результате с увеличением зоны нечувствительности размеры предельного цикла, а значит, и амплитуда автоколебаний уменьшаются. При b 1 = 0 предельный цикл вырождается в точку О. При дальнейшем увеличении зоны нечувствительности характеристика нелинейного звена и картина фазовых траекторий принимают вид, показанный на рис. 17.3, в. Здесь автоколебания отсутствуют и становится устойчивым установившийся процесс с постоянным значением ψ . Ранее неустойчивый особый отрезок F 1 G теперь стал устойчивым. Дальнейшее увеличение зоны нечувствительности приводит к расширению отрезка F 1 G, тек увеличению установившейся ошибки системы из-за слишком широкого участка равновесия. Пример 4 . Рассмотрим вибрационный регулятор напряжения, уравнения которого были составлены в § 16.2, а именно (17.24) причем уравнение нелинейного звена (регулирующего органа) (17.25) В качестве ординаты фазовой плоскости здесь удобнее взять не скорость отклонения регулируемой величины dt dU , как делалось раньше, а вторую переменную 2 I ∆ . Итак, примем для этой задачи (17.26) Тогда уравнения (17.24) преобразуются к виду (17.27) (17.28) где согласно (17.25), (17.26) и (17.28) имеем (17.29) следовательно, первое из этих условий имеет место ниже прямой (риса второе — выше нее. В первом случае переключение реле происходит при у = i 1 ; те. на прямой С (риса во втором случае — при у, те. на прямой Е. Чертеж сделан в предположении, что k 1 k 2 r 1 В результате получаем, что выше линии ЕFСD будет (17.30) а ниже линии ЕFСD (17.31) Рассмотрим сначала верхнюю область. Для нее, деля (17.28) нас учетом (17.30) получим уравнение фазовых траекторий (17.32) которое можно представить в виде и проинтегрировать, применив вспомогательную подстановку где z — новая переменная вместо у. В результате найдем следующее уравнение фазовых траекторий (при Т Т (17.33) (17.34) где С — произвольная постоянная, при γ = 1 решение будет иметь другой вида прибудет аи эти решения не будут исследоваться. Чтобы представить себе всю совокупность фазовых траекторий, можно провести на фазовой плоскости прямую (17.35) и ко всем ординатам этой прямой добавлять придавая С произвольные значения (каждому значению Сбудет соответствовать определенная фазовая траектория. Это будут параболы степени γ с осью (17.37) и с единым началом в точке Н (рис. 17.4), имеющей координаты На рис. 17.4 показаны все ветви этих парабол, лежащие выше линии ЕFСD (так как только там справедливы данные выкладки. Направления стрелок на полученных фазовых траекториях определяются тем, что проекция скорости изображающей точки справа от прямой (17.37) согласно (17.27) будет отрицательна, а слева — положительна проекция же dt dy y = υ согласно (17.28) выше прямой ух будет отрицательна, а ниже — положительна (во всех точках прямой у k 2 x касательные к фазовым траекториям горизонтальны. Аналогично строятся и все фазовые траектории ниже линии ЕFСD, так каких дифференциальное уравнение отличается от (17.32) только заменой +r 1 на —r 1 согласно (17.31). В результате на рис. 17.4 видим, что все фазовые траектории, исходящие из особого отрезка ОС, расходятся, а все траектории, идущие от краев чертежа, сходятся. Как те, таки другие асимптотически приближаются к установившемуся предельному циклу, обозначенному на чертеже жирной замкнутой кривой (чичевицеобразной). Это соответствует тому, что установившийся процесс в системе является автоколебательным, причем размеры предельного цикла аи а представляют собой амплитуды автоколебаний соответственно регулируемого напряжения U ∆ и тока в обмотке электромагнита реле Определить фазовую траекторию, образующую этот предельный цикл, можно как такую кривую (17.33), у которой (17.38) чем определяется значение произвольной постоянной С. Значение x (17.38) для этой кривой и дает искомую амплитуду а. Амплитуда же а определяется как ордината пересечения кривой предельного цикла с прямой ух (ибо, как было показано ранее, в точках этой прямой касательные к фазовым траекториям горизонтальны. Из чертежа (рис. 17.4) видно, что предельный цикл лежит левее точки L и охватывает точку С. Поэтому х с < ах, те. амплитуда автоколебаний регулируемого напряжения заключена в интервале где аи определяются формулами (17.34). Амплитуда же a I2 будет немного больше Пример 5. Рассмотрим следящую систему с сухим трением в управляемом объекте, для которой уравнения были написаны в § 16.3. Уравнение регулируемого объекта (16.52) как нелинейного звена при отсутствий линейного трения (с = 0) имеет вид (17.39) При написании уравнения линейной части системы (16.53) пренебрежем постоянными времени (чтобы иметь возможность рассматривать уравнение всей системы как уравнение второго порядка, а именно (17.40) Подставив это в уравнения объекта (17.39) и обозначив получим уравнение всей следящей системы в целом (17.41) (17.42) За координаты фазовой плоскости примем, как обычно, х = β , у = р. Условие у 0 их ikiiсiiсi1 , при котором согласно (17.42) будет рте. система будет в равновесии, изображается на фазовой плоскости отрезком АВ (рис. 17.5). Вне этого отрезка согласно (17.41) необходимо отдельно рассмотреть два случая 0 ≥ = β p y и 0 ≤ = β p y , те. верхнюю и нижнюю половины фазовой плоскости. При 0 ≤ y из (17.41) имеем Это уравнение совпадает с уравнением (16.23), но со сдвигом на величину х Следовательно, ниже оси а надо нанести такие же кривые, как на рис. 16.9, или как на рис. 16.11, б, но со сдвигом начала координат в точку А, что и сделано на риса и б соответственно. Аналогичные кривые наносятся и выше оси x, но только со сдвигом начала координат в точку В (рис. 17.5), так как согласно (17.41) при y > 0 имеем уравнение В обоих случаях (риса и б) система устойчива, причем в первом случае переходный процесс состоит из конечного числа затухающих колебаний управляемого объекта, а во втором случае имеем апериодическое движение. Положение равновесия объекта определяется неоднозначно, объект может остановиться в любой точке особого отрезка АВ (рис. 17.5), как это было уже ранее при наличии зоны нечувствительности (см. пример 1). Особый отрезок АВ определяется соотношением | М вр | с я <с, где с — абсолютное значение момента сухого трения при движении управляемого объекта. Заметим, что произведенное здесь упрощение уравнений системы хотя и позволило решить их точно, но это решение, дающее в результате устойчивость системы при любых числовых значениях параметров системы, неполно отражает действительную картину явлений в данной нелинейной системе. Метод точечного преобразования. Как видно было из примеров, фазовая траектория обычно складывается из отдельных кусков, представляющих решение уравнений системы по участкам. Пусть (риса) граничными линиями между кусками фазовых траекторий являются ось х, линия FG и линия LN. Возьмем начальное положение изображающей точки М где-нибудь на полуоси Ох. Один этап движения системы состоит в переходе изображающей точки на линию РО, ограничивающую этот этап, в некоторое положение М (риса. Следующий этап переводит изображающую точку в положение М на полуоси ОН, затем в положение М 3 на кривой Ми, наконец, в положение М на исходной полуоси Ох. Каждому положению М (хна полуоси Ох соответствует определенное положение точки М (хна кривой FG. Это называется точечным преобразованием полупрямой Ох в кривую FG. Для краткости ему присваивают название, например преобразование Дальше (риса) идет точечное преобразование кривой FG в полупрямую ОН, названное Е затем — точечное преобразование S - полупрямой ОН в кривую LN и преобразование Е- кривой LN в исходную полуось Ох. Все это в целом (или, как говорят, преобразование ЕЕ) называется точечным преобразованием полупрямой Ох самой в себя. Это преобразование в данном случае записывается в виде определенной зависимости где через x 4 их обозначены абсциссы точек M 4 и М (риса. Если при любом х оказывается x 4 < х, тов системе будет затухающий процесса если x 4 > х — расходящийся процесс. Если же возможно равенство ,х 4 =х 0 , тона фазовой плоскости получится предельный цикл, который, как известно, может изображать либо устойчивый автоколебательный процесс, либо границу устойчивости системы в малом, либо может соответствовать особому случаю бифуркации (см. ниже. В тех случаях, когда общая картина фазовых траекторий разделяется на две симметричные части, достаточно исследовать только половину всего точечного преобразования. В рассматриваемом случае верхняя полуплоскость симметрична нижней относительно начала координат. Поэтому достаточно рассмотреть только первую половину преобразования (Е, те. точечное преобразование Ох в полупрямую ОН, и выразить его в виде зависимости (17.43) причем условие наличия предельного цикла на фазовой плоскости будет |x 2 | = х при х. Пусть, например, зависимость (17.43) имеет вид кривой, показанной на рис. 17.6, б. Проведем на этом графике еще прямую изначала координат под углом 45° к координатным осям. Если она пересечет кривую, тов точке пересечения получим | x 2 хо Чтобы определить, какому типу предельного цикла это соответствует, надо взять на оси абсцисс начальную точку х сначала слева, а затем справа от точки пересечения и проследить ход точечного преобразования, как показано стрелками на рис. 17.6, б. В данном случае процесс сходится с обеих сторон к точке пересечения. Следовательно, здесь будет устойчивый предельный цикл, соответствующий автоколебательному процессу в системе. При этом абсцисса точки пересечения (рис. 17.6, б) дает амплитуду автоколебаний. Можно поступить и иначе. Допустим, преобразование S + от точки М к точке выполняется достаточно просто, но оказывается, что по полученной определенной точке М находить соответствующую точку М труднее, чем по заданному положению М 2 определять М. Тогда будем подходить к кривой FG с двух сторон, задавая одновременно точку М (x 0 ) на полуоси Охи точку М ( - хна полуоси ОН и находя соответствующие точки Ми М (рисе. В результате получим точечные преобразования полупрямых Охи ОН в кривую FG, выраженные некоторыми зависимостями Изобразив это в виде двух кривых (рис. 17.6, г) анализируем их тем же способом, как и кривую с прямой на рис. 17.6, б. Такие графики (17.6, 6 и г) называют диаграммами точечного преобразования. Они соответствуют в данном случае устойчивому предельному циклу, те. наличию установившегося автоколебательного процесса в системе. Другие возможные типы диаграмму точечного преобразования показаны на рис. 17.7. При этом риса соответствует неустойчивому предельному циклу он ограничивает область начальных условий хау, при которых система оказывается устойчивой относительно установившегося состояния с постоянным значением регулируемой величины (х = 0). При начальных же условиях х, у, выходящих за контур этого предельного цикла, система неустойчива (система устойчива в малом и неустойчива в большом. Рис. 17.7, в соответствует двум предельным циклам, из которых меньший неустойчива другой больший) устойчив. Следовательно, при начальных условиях х, y 0 , расположенных внутри первого предельного цикла, система устойчива, как ив предыдущем случае, а при всяких других начальных условиях она стремится к установившемуся автоколебательному процессу, который определяется вторым предельным циклом. Этот случай может выродиться в случай, изображенный на рис. 17.7, б, когда оба предельных цикла сливаются в один полу устойчивый. Подобные особые случаи называются бифуркационными. Наконец, на рис. 17.7, гид изображены случаи, когда на диаграмме точечного преобразования кривая | х | = хане пересекается с прямой, проведенной под углом 45° к осям. Это означает соответственно устойчивость (г) и неустойчивость (д) системы при любых начальных условиях (до которых справедливы исследуемые уравнения системы. Заметим, что изложенное выше является лишь качественным рассмотрением, так как в нем отсутствует время t. Остается неизвестным течение процессов во времени, период частота) автоколебаний. Чтобы получить полное решение задачи, нужно к данному рассмотрению добавить еще решение дифференциальных уравнений на отдельных участках во времени (как в методе припасовывания, см. § 16.1). Поэтому в методе точечного преобразования вводится соответствующий временной параметр (здесь это не рассматривается. Для систем выше второго порядка вместо фазовой плоскости придется иметь дело с фазовым пространством и с точечным преобразованием не линий, а поверхностей. Там появляются новые особенности процессов. Однако ввиду чрезвычайной сложности таких построений рассматривать их не будем. Метод изоклин. Выше были рассмотрены такие примеры нелинейных систем второго порядка, для которых фазовые траектории легко находятся интегрированием уравнений по участкам. В тех случаях, когда интегрирование затруднено, ход фазовых траекторий, хотя бы качественно, можно проследить с помощью так называемого метода изоклин (без интегрирования уравнений. Количественно этот способ имеет сравнительно низкую точность. Применение его пока ограничено системами второго порядка. Изоклиной называется такая линия, во всех точках пересечения которой с фазовыми траекториями последние наклонены под одними тем же определенным углом коси абсцисс х. Так, если известно дифференциальное уравнение фазовых траекторий (17.44) то для получения изоклины нужно положить Уравнение изоклины, следовательно, будет (17.45) где с обозначает определенный тангенс угла наклона фазовых траекторий. Каждому заданному значению с соответствует своя изоклина. Например, часто встречается нелинейное уравнение которое можно записать в виде тогда дифференциальное уравнение фазовых траекторий будет выглядеть так а уравнение изоклин Задавая различные значения с (при заданном К, для каждого из них строим поэтому уравнению кривую на фазовой плоскости — изоклину (сплошные кривые на рис. 17.8). Затем на каждой кривой наносим стрелочки под углами ас коси абсцисс (на рис. 17.8 указаны значения с для каждой кривой. Из этих стрелочек и составляются искомые фазовые траектории некоторые из них изображены на рис. 17.8 пунктиром. В данном случае получается устойчивый предельный цикл, что соответствует автоколебаниям в системе. |