|
|
|
2. |
Второй метод |
Предположим, что выполняются условия |
|
|
»<1-г\и\, |
|
|
|
(7-54) |
где г ^ 1 , р — произвольно малая |
положительная |
вели |
чина. |
|
|
|
|
|
Уравнение (7-44) позволяет записать: |
|
|
|Хп-^т |-^7 (/l ||Xn+m—i |“h ••• |
|fm |j Хп |, |
(7-55) |
откуда |
п+пг—1 |
|
|
|
|
|
|
|
|
|
|Мг+т| = |
^Г V |
|Х,-|. |
|
(7-56) |
|
i=n |
|
г раз |
|
|
Величина переменной |
системы |
в |
меньше, чем |
сумма последовательности абсолютных значений, если
^ < 1 —r|/i|.
З а м е ч а н и е . Условия (7-47) и (7-54) не являются достаточ ными, обеспечивающими асимптотическую устойчивость. Необходимо
дополнить их условием асимптотической устойчивости, |
например |
условием (7-45), или же, что то же самое, положить |
1, как это |
было сделано выше. |
|
Если принять r=m, то устанавливается эквивалентность усло вий (7-46) и (7-54). Наибольшая область начальных условий опре
деляет такую область, что длительность переходного |
процесса аТ |
с точностью Р удовлетворяет неравенству |
|
й о да‘ > 1- |
(7‘57> |
б) СИСТЕМА НЕЛИНЕЙНЫХ РАЗНОСТНЫХ УРАВНЕНИЙ
Можно также установить закон постепенного убыва ния решений системы двух разностных уравнений перво го порядка и отыскать область начальных условий, для которой этот закон удовлетворяется. С. Мира [Л. 7-8] рассматривает систему
■*(« + 1)= аил;(/г) + а12г/(«); J |
(7.58^ |
У(« + 1 )= а а1л:(я) + а22у(п), J |
|
где величины а.ц являются функциями х(п), у(п) и я. Предположим, что для конечного я для данных на
чальных условий |
существует единственное |
решение |
х(п) = у (п ) -— 0 и |
что оно асимптотически |
устойчиво, |
а критический случай исключен. |
|
f
Наложим на коэффициенты а^- условия
i |
I -Н! |
I |
ь > |
1 |
|
|
|
! |
(7-59) |
I ^12 |
I “Ь | ^"22 | “-С £ > |
k^>- \ J |
и отыщем неравенства, связывающие переменные х(п),
*/(«)•
Уравнения (7-58) позволяют записать:
\х(п+1) |< |ац| |х{п) |+ |а12|\у{п) |; |
|
|г/(п+1) |^ |a2i| |х(п) j -f ja22| |^(н) |, |
(7-60) |
откуда
\х{п+\) |+ \у{п+\) |<f|an| + |a2i|J|Jc(n) |+
+[ai2|+ |«22131у(п) | |
(7-61) |
и, если принимать во внимание соотношения |
(7-59), |
I * (* + 1) i + 1У (п + 1) 1< |
(7-62) |
Соотношения (7-59) позволяют найти область началь ных условий в дискретной фазовой плоскости х(п), у(п), гарантирующую требование, чтобы изменение перемен ных происходило бы по (7-62).
Таким образом, мы получаем фиксированную после довательность кривых, ограничивающих область, для ко
торой значения |
не |
являются явными функциями п. |
|
в) |
СИСТЕМА МАТРИЧНЫХ УРАВНЕНИИ |
Ф. Лоран {Л. 7-9] обобщил вышеописанный критерий на случай нелинейных систем, описываемых разностным матричным уравнением. Однако, сославшись на соответ
ствующие публикации, опустим изложение |
этих работ. |
1. Конечно-разностное уравнение первого |
Для уравнения первого порядка |
порядка |
|
xn+i+ fixn^0 |
(7-63) |
условия |
|
p < l-r | fi| ; /->1 |
(7-64) |
требуют, чтобы удовлетворялось неравенство |
|
Эти неравенства можно было бы очень легко дока зать прямым способом.
2. Конечно-разностное уравнение второго порядка
Пример 7-5. Рассматриваем нелинейное уравнение |
|
^n+2+ /l^n+l+ /2^n=0. |
(7-65) |
Условие |
|
Р < 1 - '1 М - '21Ы |
(7-66) |
накладывает на а ограничение
— га > 1
100 г >
Неравенство (7-66) может быть изображено в плоскости коэф фициентов (рис. 7-5). На рис. 7-5 производится сравнение областей устойчивости для линейного и нелинейного уравнений.
Области устойчиЬости
линейной системеI Область устойчибости ' начальны х услобий
Области -1
асимптотичес кой устойчивости нелинейной системы
Iх on-ej<~ jfxa+/j+l XocjJ
Рис. 7-5.
Условия (7-54) запишутся:
|
РК I — г |f, |
(7-67) |
|
р/< 1—г |f2 |
|
: } |
Соответствующие области асимптотической устойчивости изобра жены на рис. 7-5.
Пример 7-6 [Л. 7-7]. Обратимся снова к уже исследованной не линейной системе (рис. 7-6). Соотношения между входными и вы ходными сигналами различных элементов
—Wn=kiXn+i', |
(7-68) |
|
. |
Tidy |
ktw. |
(7-69) |
|
У |
dt |
|
|
e = s — y. |
|
(7-70)' |
|
х = Ф (е ) = f(e )e . |
(7-71) |
Рис. 7-6. |
f
Примем k0= fe,£2, D — e г/г'< 1 , s = 0.
Уравнение в конечных разностях, описывающее систему, будет:
Bn+2[( 1—D) Kof (бп+i) —D— 1]бп+1-\-Dbu=0. |
(7-72) |
Достаточное условие асимптотической устойчивости, как это сле |
дует из гл. 6, определяется неравенством |
|
|
|
|
2D |
|
|
г |
|
|
2 |
|
(7-73) |
|
1 _ о |
< |
k ° f (s) < |
j _ о |
' |
|
|
Для линейной системы kBf (е) = const достаточное и необходимое |
условие будет: |
|
|
|
14-D |
|
|
|
|
|
|
|
|
|
(7-74) |
|
*of ( « ) < 2 Т- ± 0 . |
|
Соотношение (7-66) преобразуется для этого частного случая |
гЮ + rD - f г — 1 |
|
|
|
|
rW + rD + г + 1 |
r( 1 - D ) |
< |
kof («) < |
■ |
|
r ( l - D ) |
(7-75) |
Оно иллюстрируется на рис. 7-7 |
(а — область устойчивости не |
линейной системы; |
б — область с длительностью переходных процес |
|
|
|
сов |
аТе; |
в — область |
устойчивости |
|
|
|
линейной системы). Двойной штри |
|
|
|
ховкой обозначается область, в кото |
|
|
|
рой для системы, описываемой (7-72), |
|
|
|
обеспечивается |
время |
переходного |
|
|
|
процесса, |
равное аТ, при Р, если |
|
|
|
|
|
|
|
100га> |
(7-76) |
|
|
|
|
Пример 7-7 [Л. 7-8]. Рассмотрим |
|
|
|
систему |
|
|
|
|
|
|
|
%п+1 “ |
6 , 2 уп, |
|
(7-77) |
|
|
|
|
|
|
|
|
|
Рис. 7-7. |
|
|
|
У п + 1 = 0,1х„ (1 + 0,2х"). |
|
|
|
|
|
|
|
|
|
Примем |
|
|
|
|
|
|
|
|
|
I * п + 1 |
I + Уп + |
|
|
|
х (и) |
|
У (п) 1 |
(7-78) |
I 1 |
I |
< |
|
' |
|
|
|
|
|
|
Область переменных, для которой удовлетворяются эти соотно |
шения, определяется как |
|
|
|
|
|
|
|
|
|
0,1[1 + 0 ,2 х 2( л ) ] < ^ |
’ |
(7-79) |
|
I х (п) |
|< |
Vm |
|
|
(7-80) |
|
—ту- |
|
Первое |
выражение |
этой |
системы накладывает следующее |
условие: |
|
|
|
|
|
\ y ( n ) \ < ~ V W . |
(7-81) |
3 |
а м е ч а н ii е. |
Только |
что рассмотренные |
методы хотя и н |
ограничены, но и не полностью решают поставленную задачу. При каждом конкретном использовании их в соответствии с изложенными работами следует искать область начальных условий, наиболее соот ветствующую решаемой задаче.
Интересно отметить, что вышеприведенные условия всегда более жестки, чем условия асимптотической устойчивости. С другой сторо ны, чем меньше требуемая длительность переходного процесса аТ, тем меньше должна быть область, внутри которой должны нахо диться коэффициенты нелинейного уравнения, что легко объяснимо..
7-4. ПРЕДЕЛЬНЫЕ ЦИКЛЫ КОЛЕБАНИЙ НЕЛИНЕЙНЫХ ИМПУЛЬСНЫХ СИСТЕМ
За последние годы отмечаются большие успехи в обла сти исследований нелинейных импульсных систем, при чем работы ведутся главным образом в направлении исследования устойчивости.
Вопросы, связанные с рассмотрением предельных ко лебаний, носят совершенно иной характер; ниже будет сделана попытка сделать их краткий обзор без подроб ного рассмотрения доказательств предложенных крите риев, (для более подробного знакомства с ними даются ссылки на литературу). Вопреки тому, что в настоящее’ время многие авторы заинтересовались задачей исследо вания периодических или почти периодических колеба ний в системах с различными нелинейностями, в боль шинстве случаев в нелинейной импульсной системе су ществует непериодическое неустойчивое колебание (либо; колебание с бесконечным периодом).
Теоретическое исследование таких колебаний из-за их: непериодичности весьма затруднено, и, как нам известно,, этому вопросу посвящено небольшое число публика ций.
а) ПЕРИОДИЧЕСКИЕ И КВАЗИПЕРИОДИЧЕСКИЕ РЕШЕНИЯ СИСТЕМ КОНЕЧНО-РАЗНОСТНЫХ УРАВНЕНИЕ
А. Халанай (Л. 7-10] получил некоторые результаты: по установлению существования периодических и квазипериодических решений систем конечно-разностных урав нений. Мы остановимся на основных идеях этих работ.
г
Рассмотрим |
квазилинейную |
систему |
следующего |
вида: |
Х(га+1)=А(га)Х(га)+Р[Х(п), га]; |
|
(7-82) |
|
|
X и |
F — векторы; |
А — квадратная матрица, |
причем А |
и F — периодические функции га с периодом N. |
|
|
Предположим, что линейная система |
|
|
|
|
|
|
Х(га+1) = А(га)Х(га) |
|
|
|
(7-83) |
не имеет никакого |
другого периодического |
решения |
с периодом N, кроме нулевого решения. |
|
|
|
|
Пусть уравнение |
|
|
|
|
|
|
Х ( п ) = |
N—1 |
|
гаг], |
|
(7-84) |
|
G(ra, rai)F[X(/n), |
|
где |
|
|
m =zО |
|
|
|
|
|
G(n, гаг) = Х(га, 0);[I—X (N—0)h 1Y (гаг, n—1) + |
|
|
|
для m ^ n — 1; |
|
+ Y (гаг, ra— 1) |
|
|
|
(7-85) |
|
|
|
|
|
|
|
|
G{n, гаг) = |
X(«, 0)[I—X(N, 0 )h 1Y(m, « — 1) |
(7-86) |
для гаг^га; здесь |
I — единичная |
матрица; |
Y(гаг, ra)— та |
кая |
функция, |
что Y(гаг, |
ra) = I, |
Y(гаг— 1, |
ra) — |
= Y(гаг, n)A(m). |
|
|
|
|
|
|
|
Легко показать, что решение задачи о существовании периодических решений системы (7-82) эквивалентно решению задачи о существовании периодических реше ний для уравнения (7-84).
Пусть
F(0, га) — Gj(ra); F(X, ra)-F(0, га) = ==pF1i[X(n)]. (7-87)
Система (7-82) запишется в этом случае как
X ( n + l ) = A ( n ) X ( n ) + G l(n) +
+ JiF1[X(n), |
га], |
(7-88) |
а уравнение (7-84) |
|
|
|
N~ I |
N—1 |
|
|
X (га) = Yi G (п> т ) G, (гаг) |
Р* 2 |
G [га, |
гаг] F! [X (гаг), гаг]. |
т=.0 |
т=0 |
|
|
(7-89)
Если
га—1
Ф («) = Ц G (га, гаг) (гаг), т=0
где Ф (п) — периодическое решение системы
X { n + l ) = A ( n )X ( n )+ G l(n), |
(7-90) |
то уравнение (7-84) приводится к виду
N -1
Х(я) = ф (я) + |Л£ Q (я, т) F, [X (яг), т\. (7-91) т=0
Решение строится в виде постепенно сходящейся по
следовательности |
|
Х0(п.) = ф(пу, |
|
N -1 |
(7-92) |
X j(я) = Ф (я) -{- fA £ G (я, |
т\. |
т=0 |
|
Теорема 1. Рассмотрим систему
X (п+ 1) = A (n) X (п) + Gi (п) + pFijX (п) , п],
где A, Gi, Ft — периодические функции с периодом N. Если система
Х (п + 1 )= А (п )Х («)
не имеет никаких других периодических решений с пе
риодом N, кроме нулевого, то |
для достаточно малого |
f pi| возможно существование |
единственного периоди |
ческого решения с периодом N, которое можно получить |
методом последовательных приближений. |
Халанай также предложил следующую теорему. |
Теорема 2. Если в системе, |
описанной уравнением |
(7-82), функция F(X, п] непрерывна относительно X и если
|Р[Х, Я Ц < Н | Х | ) ; ^ < ^ Г .
где М — константа, определяемая только лишь системой (7-88), то в системе (7-82) возможно периодическое ре шение с периодом N.
В этой теореме предполагается, что в системе (7-82) не существует никакого другого периодического решения с периодом N, кроме нулевого.
Теорема 3. Рассмотрим систему |
|
Х (« + 1) =А(п)Х(п) +F(n), |
(7-93) |
матрицы которой А и А-1, вектор F — квазипериодические функции.
Если нулевое решение системы
Х («+ 1) = А (п )Х (п )
является равномерно асимптотически устойчивым, то си стема (7-93) допускает квазипериодическое решение, равномерно асимптотически устойчивое и удовлетворяю щее ограничению вида
|Х(я) | |
sup|F(n) |. |
(7-94) |
Теорема 4. Исследуем систему
Х(п+1)=Хо[л, Х(п)] + еХ1[пХ(«), е], |
(7-95) |
где Х0 и Х( — периодические функции п с периодом N. Предположим, что для е= 0 в системе возможно пе риодическое решение Х0(я), которое является равномер но асимптотически устойчивым. Тогда для достаточно малого |е| система, возможно, имеет периодическое ре
шен:-! з, которое при е— й) стремится к Х0(я). Теорема 5. Рассмотрим систему
Х(п + 1) — А(п)Х(п) + eX i[п, Х(я), п]; е>0 (7-96)
и предположим, что А, А-1 и Xt — квазипериодические функции относительно п, и что система
Х (я + 1) = А (п)Х(п) |
(7-97) |
допускает все квазипериодические решения.
Пусть Х(п, 0) — матрица решений системы (7-97) при Х(0, 0) = Е и
Z[п, Х(п), е]=
= Х-1(л+1, OJXJn, |
Х(л, 0), |
z(n), |
г]; |
|
(7-98) |
|
|
N |
|
|
|
|
|
|
|
X Z (П, |
Z, е ) |
|
|
|
|
z 0(2, е) = иш |
---------. |
|
|
(7-99) |
где Z0 — решение уравнения |
Z0(г0, |
0) = |
0, |
такое, что |
сооственные величины |
|
„ |
дЪ0 ( г , 0) |
1 |
имеют |
матрицы Н = — |
— - |
|
отрицательные части. |
Тогда |
система |
|
02 |
| г= г„ |
(7-96) |
допускает |
квазипериодическое решение, которое для г— Я) стремит ся к решению Х(п, 0)г0:
Х(п) = Х(«, O)z(n).
Уравнение (7-96) запишется тогда так:
z(n+ 1) = z(n ) +eZ[n, z(n), е].
З а м е ч а н и е . Доказательство этих теорем очень сложно и гро моздко. Кроме того, отметим, что эти теоремы позволяют судить лишь о существовании периодических и квазипериодических решений систем конечгго-разностных уравнений; периодическое решение не может быть получено только лишь приближенным образом.
б) ОПРЕДЕЛЕНИЕ И УСТОЙЧИВОСТЬ СИММЕТРИЧНЫХ ПЕРИОДИЧЕСКИХ КОЛЕБАНИИ
Я. 3. Цыпкин [Л. 7-11] предложил метод для опреде ления симметричных периодических решений нелинейно го разностного уравнения и исследовал их устойчивость.
1. Определение колебаний
Для того чтобы получить периодические решения основной частоты для уравнения с конечными разностя ми, необходимо произвести трудоемкие вычисления. По этому Цыпкин предлагает приближенный метод, пригод ный, когда частота колебаний по крайней мере в 3 раза ниже, чем частота импульсного элемента.
Рассмотрим нелинейное уравнение с конечными раз ностями вида
т |
т ' |
i)\ |
|
(П - f - Т) = |
^ (Зг- р ( п |
( I - 1 0 0 ) |
1=П i=0
у[п)=<Цх(п)1 (7-101)
где Ф(х) — симметричная нелинейная функция и m>m'. Отыщем симметричные периодические решения х(п)
с частотой co = it/lV (где N — половина периода колеба ний, N — целое).
С физической точки зрения можно доказать, что эти решения могут быть представлены разложением в ряд Фурье
м, |
|
х (п) = 0,5 2 ] а’ь exp j / -|r nj • |
(7-102) |
Значения а\—a^exp (jcpk) — искомые комплексные амплитуды. Первый индекс указывает, что суммирование производится лишь по нечетным значениям, и
N — 1, если N четное;
Ni — N, если N нечетное.
Если х(п) — периодическая функция, то у(п) будет периодической и сможет быть записана в виде полинома Фурье
у(п) = Ф [х («)] = Ф 10,5 |
a'ftexp |
~. knn |
|
|
} - |
^k = N t
N1
= 0 ,5 ^ |
|
, knn |
(7-103) |
b'kexp^j ~w |
k -N , |
|
|
|
Коэффициенты этого полинома могут быть определе^ |
ны по формуле гармонического анализа Бесселя |
|
|
N—1 |
|
(7-104) |
Ь\ = Ькехр ЦЩ = |
Ф [х (/)] ехр [ - } ~ , |
|
<=о |
|
|
и если N — нечетное для k= ±:N, то |
|
|
N- |
1 |
|
b'±N= b±N^ -JT |
ф Iх (0] ехр [—//it]. |
(7-105) |
i=0
Подставляя (7-102) и (7-103) в уравнение (7-101) и приравнивая коэффициенты при одинаковых степенях, получаем систему алгебраических уравнений, опреде ляющих желаемые комплексные амплитуды:
k = s l’ 3’ - |
’ |
(?-106) |
или |
|
|
т |
|
|
Q (/ ^ ) = S e' exp[ /^ ?] : |
I |
[(7-10 |
,=* |
|
i= 0 |
* |
|