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

книги из ГПНТБ / Видаль П. Нелинейные импульсные системы

.pdf
Скачиваний:
6
Добавлен:
24.10.2023
Размер:
10.96 Mб
Скачать

 

 

 

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 и

что оно асимптотически

устойчиво,

а критический случай исключен.

 

231

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)

требуют, чтобы удовлетворялось неравенство

 

232

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

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.

233

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)

 

—ту-

 

234

Первое

выражение

этой

системы накладывает следующее

условие:

 

 

 

 

 

\ y ( n ) \ < ~ V W .

(7-81)

3

а м е ч а н ii е.

Только

что рассмотренные

методы хотя и н

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

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

7-4. ПРЕДЕЛЬНЫЕ ЦИКЛЫ КОЛЕБАНИЙ НЕЛИНЕЙНЫХ ИМПУЛЬСНЫХ СИСТЕМ

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

Вопросы, связанные с рассмотрением предельных ко­ лебаний, носят совершенно иной характер; ниже будет сделана попытка сделать их краткий обзор без подроб­ ного рассмотрения доказательств предложенных крите­ риев, (для более подробного знакомства с ними даются ссылки на литературу). Вопреки тому, что в настоящее’ время многие авторы заинтересовались задачей исследо­ вания периодических или почти периодических колеба­ ний в системах с различными нелинейностями, в боль­ шинстве случаев в нелинейной импульсной системе су­ ществует непериодическое неустойчивое колебание (либо; колебание с бесконечным периодом).

Теоретическое исследование таких колебаний из-за их: непериодичности весьма затруднено, и, как нам известно,, этому вопросу посвящено небольшое число публика­ ций.

а) ПЕРИОДИЧЕСКИЕ И КВАЗИПЕРИОДИЧЕСКИЕ РЕШЕНИЯ СИСТЕМ КОНЕЧНО-РАЗНОСТНЫХ УРАВНЕНИЕ

А. Халанай (Л. 7-10] получил некоторые результаты: по установлению существования периодических и квазипериодических решений систем конечно-разностных урав­ нений. Мы остановимся на основных идеях этих работ.

235-

г

Рассмотрим

квазилинейную

систему

следующего

вида:

Х(га+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

236

где Ф (п) — периодическое решение системы

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)

2 3 7

матрицы которой А и А-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)

допускает

238

квазипериодическое решение, которое для г— Я) стремит­ ся к решению Х(п, 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)

239

Значения а\—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

*

 

2 4 0

Соседние файлы в папке книги из ГПНТБ