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

книги из ГПНТБ / Конторович М.И. Нелинейные колебания в радиотехнике. (Автоколебательные системы)

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

в правой

части

(7)

приведет

к малым погрешностям

(порядка

|i) для

всех

времен

Аналогичные со­

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

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

В дополнение полезно сказать, что во многих зада­ чах, связанных с процессом установления колебаний, длительность процесса зависит от величины [i и тем больше, чем меньше р,. Поэтому необходимо не огра­

ничиваться малыми [it, а рассматривать

большие отрез­

ки времени. Метод последовательных

приближений

в его первоначальной форме такой возможности не дает.

Сказанное выше по поводу метода ММА при над­ лежащих условиях может быть строго доказано. Ви­ димо, наиболее общая формулировка задачи и строгое доказательство того, что решение уравнений первого приближения в методе ММА (укороченных уравнений) для времен порядка L/[i могут быть сделаны сколь угод­ но близкими к решениям точных уравнений, если [i выбрать достаточно малым, имеется в книге Н. Н. Бого­ любова и Ю. А. Митропольского «Асимптотические ме­ тоды в теории нелинейных колебаний».

В некоторых специальных случаях (при наличии устойчивых стационарных состояний) укороченные урав­

нения

дают решения,

мало

отличающиеся

от

решения

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

на

бесконечном

интервале

времени.

Доказательство

этого

положения,

принадлежащее

Н. Н. Боголюбову, приведено в указанной выше

моногра­

фии (см. также приложение

1).

 

 

 

 

 

 

 

4.2.1. Случай двух близких частот

 

 

 

 

В

предыдущих

рассуждениях

предполагалось,

что

в правой части (7) § 4.2

первый член

не содержит

яв­

но

t. Однако если

среди частот т будут близкие друг

к

другу или среди

частот

vr

найдется

близкая

к одной

из частот cofe (в том смысле,

что соответствующая

раз­

ность

будет иметь

порядок малости

и,), то условие

(5)

§ 4.2 не будет выполняться.

В этом

случае

надлежит

соответствующие

медленно

осциллирующие

члены

110

4.2.1

изъять из Fft(i' и отнести их к первому слагаемому Ми. Вид уравнений (4) § 4.2 в этом случае остается преж­ ним, с той лишь разницей, что теперь Mh явно будет зависеть не только от Л, но и от /. Однако легко видеть, что предыдущие рассуждения останутся в силе и уко­ роченные уравнения по-прежнему сохранят вид (8)

§4.2.

4.3.Уравнения со свободным членом в правой части

Рассмотрим теперь уравнение, имеющее вид

L(y) = anyW

+ ап-гУ^ + ... + а0у =

 

= рР(У,

У', .... У(п~1]\ 0 + / ( 0

(1)

и отличающееся от (1) § 4.1 только наличием f(t) в пра­ вой части. Эта функция в дальнейшем предполагается периодической с периодом Т и может быть представ­ лена в виде ряда Фурье

 

00

 

 

 

S = — 0 0

 

 

где <о = 2я/7\

 

 

 

Покажем теперь, что уравнение

(1) сводится к урав­

нению (1) § 4.1, и, следовательно, вся изложенная

мето­

дика может быть к нему

применена. Нам придется рас­

смотреть два случая.

 

 

 

4.3.1. Нерезонансный

случай

 

В этом случае ни одна

из разностей | СОЙ.—sco| не рав­

на нулю для любых СОЙ

из числа

определяемых

уравне­

нием (4) § 4.1 и любых | s | = 0 , 1, .. . Рассмотрим

вспо­

могательное уравнение

 

 

 

L(y0)=f(t).

 

(1)

Будем искать его периодическое решение с периодом Т. Полагая

s=—oo

4 . 3 . 1 .

I l l

после подстановки в (1) получаем

 

 

 

 

00

 

4

-

dsL(jsm)ei;

= 4 -

£

^

 

 

s=—оо

 

S = — 0 0

 

где под L(p)

подразумевается anpn

+ an-zpn'2 + - • . + а0.

Отсюда

непосредственно

находим

 

 

 

 

de = Ce/L(jas).

 

(2)

Формула

 

(2) однозначно

определяет

коэффициенты

ds, так как ее знаменатель нигде в нуль не обращается.

Теперь положим у = уо + £,, тогда

из (1)

следует:

an|<n) + an-2|(B -2 >+...+"aoS = |iO(g,

I', ....

l ^ ; t),

(3) причем Ф — функция, образующаяся из F после под­ становки у = уо + Ь,. Полученное уравнение (3) имеет та­ кой вид, как и (1) § 4.1.

4.3.2. Резонансный случай

Здесь одна из комбинаций |со/гso)| равна нулю или имеет порядок малости ц. Теперь уже (1) § 4.3 не имеет периодического решения, так как один из коэффициен­ тов (2) п. 4.3.1 обращается в бесконечность (или велик) В этом случае уравнение (1) § 4.3, вообще говоря, не будет иметь решений, остающихся ограниченными при стремлении \\ к 0. Учитывая это обстоятельство, обычно считают, что «резонансный» член в (1) § 4.3 сам имеет малую амплитуду (порядка малости р,), что во многих случаях вызывается требованием существа рассматри­ ваемой физической задачи.

Пусть, например, член с номером г является резо­ нансным. Тогда напишем

Д 0

=

~ И S

 

c ^

t + ^™t}=M)

+

i-Cre'rnt.

 

 

s=

— oo

 

 

 

 

Относя

теперь

в

уравнении

('1) § 4.3 резонансный член

к F,

получаем

 

 

 

 

 

 

 

 

any(1t)

+

я » .2У{п

~2) + - + а0у

=

 

 

 

F(y,

у ' , . . . ;

0 + 4г е/ Г А "]+Л(0.

112

4 . 3 . 2 .

Таким образом, мы вновь пришли к уравнению ви­ да (1) §4 . 1, рассмотренному раньше.

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

4.4.Обобщение предыдущих результатов.

Получение второго приближения

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

ние и-го порядка вида (1) § 4.1 и было показано, что его можно привести к системе п уравнений первого по­ рядка с малым параметром в правой части. В дальней­ шем будем эту систему называть стандартной. Далее составлялась система укороченных уравнений, дающих решение задачи в первом приближении. Теперь, уже исходя из заданной стандартной системы, постараемся более точно сформулировать условия, при которых мо­ гут быть получены решения первого приближения, и по­ кажем, как могут быть получены более точные решения.

Пусть имеется стандартная система вида

^.=^[Ф1 (Л, 0 +

Ф,(А

')].

О)

где |л]>0— малый параметр

и Л,

Ф п

Ф2 я-мерные

ком-

плекснозначные векторы-столбцы

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

соот­

ветственно равными А\ Ф{к),

Ф

 

/ г = 1 , 2,

п.

Ищется непрерывное решение системы (1), обращаю­ щееся при t = t0 в заданную вектор-функцию A(ta). Не уменьшая общности рассуждений, можем в дальнейшем положить /0 = 0.

Будем считать, что в рассматриваемом векторном пространстве существует надлежащим образом опреде­ ленная «-мерная область D, содержащая A (to) в каче­ стве своей внутренней точки, в пределах которой функ­ ции CPi и Ф2 удовлетворяют следующим условиям при всех t^to:

 

\d>i{A, t)\<M-

 

2 (Л,

t)\<M-

(2)

 

 

{- дФ2*> (A,

t)

 

^ф[к)(А,

/) < / И ;

J

длт

dt

(3)

3 4.4.

- 1 2

113

где k=\, 2, . . . , п; М — некоторое положительное число, не зависящее от Л и от ц, или функция от ц, остающаяся ограниченной при сколь угодно малых р. При вычисле­

нии интегралов (3) Л

считается постоянным

(интегри­

рование

ведется лишь

по

явно входящему

Прямые

скобки

в этих формулах

(и далее), отнесенные к век­

тору или матрице, означают норму, определяемую как

сумму модулей

отдельных

составляющих

этих

величин;

в случае

скаляра — модуль заключенной

в

скобки

ве­

личины.

Далее

предполагается, что Ф 1

удовлетворяет

условию

Липшица

\Ф\(А

+ r|, t)— Ф4 (Л,

* ) | < / ( | г | | ,

^>

А+ц^О

при всех t^O и

К — некоторая

константа,

не

зависящая

от р.

 

 

 

 

 

 

 

Область D может

быть

определена различным обра­

зом, например, так:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(4)

где

— некоторые

положительные числа

(не

зависят

ни от t,

ни от д.).

 

 

 

 

 

 

 

Очевидно, что при таком определении термин «об­ ласть означает лишь множество величин Л</;> (или векторов Л), удовлетворяющих ( 4 ) . Если хотя бы для одной из величин AW неравенство (4) превращается в равенство, мы будем говорить, что вектор Л лежит на границе области D.

Условимся в дальнейшем множество функций, удов­ летворяющих условию (3), обозначать как М0 и писать, например, так: Ф20.

В соответствии с тем, что говорилось в предыдущих параграфах, естественно составить укороченные уравне­

ния, соответствующие системе (1), путем

отбрасывания

в правой части слагаемых класса

М0, т. е.

написать

dfl/d* = (!(!>! (В,

t),

(5)

где вектор Л заменен на В.

Основой для дальнейшего будет следующее утвер­ ждение, доказательство которого приводится в прило­

жении 1.

 

 

 

 

 

Пусть вектор В,

удовлетворяющий

(5) и

заданным

начальным

условиям

при всех 0^,t^L/ix,

где

L — неко­

торое положительное

число, лежит в

области

D, такой,

1 Это значит, что каждое Л<*>(<)—Л<*>(0)

не выходит

в комплекс­

ной плоскости

за пределы

круга (включая точки

окружности) радиу­

сом

 

 

 

 

 

114

 

 

 

 

§4 . 4 .

что |5<*>(0—Я( Л , (0)|<#<*>—Д, k=\, 2, . . . . n, и Л — не­ которое сколь угодно малое, но отличное от нуля поло­

жительное число (область

D является

«расширенной»

по отношению к Dt). Тогда

для вектора

ц=А—В, опре­

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

имеет место

следующая оценка по

порядку величины:

 

h K Q i ,

 

(6)

где С — некоторая константа,

пропорциональная

М, т. е.

оцененная по

норме разность

между

решениями

точных

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

Следует

обратить внимание

на

то,

что

М,

входящее

в (3), может зависеть от «размеров»

выбранной

обла­

сти D и, в частности, во многих случаях M<Mi\B(t)

—В (0)|max,

где

Mi — абсолютная константа,

не завися­

щая ни от

(я,

ни от t,

ни от

«размеров»

D;

\B(t)—•

—В (0)| max — наибольшее

значение,

которого

достигает

величина при всех рассматриваемых

значениях / и

B^D.

Тогда оценка приобретает

вид

 

 

 

 

 

 

 

 

| г | | < О г | Л - Л ( 0 ) | т а х .

 

 

 

(7)

При некоторых условиях оценка (7) оказывается значи­ тельно менее грубой, чем (6).

Следует обратить внимание на то, что сказанное от­ носится лишь к конечным интервалам времени t, хотя и сколь угодно большим при достаточно малых ц. Если рассматривать фиксированное и все 0 = £ ^ < о о , то соот­ ветствующие теоремы будут формулироваться иначе и требовать дополнительных условий. Сейчас мы на этом

останавливаться не будем

и

рассмотрим

этот

случай

в приложении

1.

 

 

 

 

 

 

 

Отметим теперь еще некоторые обстоятельства, от­

носящиеся к тому случаю, когда рассматривается

конеч­

ный интервал t длиной L/ц.

Очевидно, что в этом

случае

нет необходимости

требовать,

чтобы

условия

(2)

и (3)

выполнялись

при

любых t,

а

лишь

для

тех,

которые

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

Однако этот случай не является единственным. Пусть, например, имеется функция (векторная) F(A, t) такая,

§4.4.

8*

115

что на всем интервале 0 ^ / < L / L I она сама и все ее пер­ вые частные производные по компонентам Л ограничены, тогда

$Р'Р(А,

t)dt

 

^\\F\mJ<*\F\naL<v,M,

о

 

|

о

где l^lmax наибольшее значение нормы вектора F. Аналогично получаем для этого случая и второе нера­ венство (3). Таким образом, при интервале 0 ^ / < L / L I функции, ограниченные вместе со своими первыми част­ ными производными, умноженные на малый параметр ц, также принадлежат к классу М0. Отсюда, в частности, следует, что входящие в правую часть слагаемые, имею­ щие порядок малости р 2 (и при выполнении других ука­ занных условий), при составлении укороченных уравне­ ний также могут быть отброшены.

4.4.1. Построение уравнений второго приближения

Для получения уравнений второго приближения по­ ложим A=Ai + ^a, где А и а — векторы-столбцы а ( 0 ) = 0 . Подставляя в (1) §4.4, получаем

 

— + ц ™ = j , [Ф, (Л, +

va) +

Ф2 (Л, +

w)].

 

Подчиним А уравнению первого приближения, т. е.

положим

 

 

 

 

 

 

 

 

 

 

 

 

dAildt

=

\iOi{Al)

 

 

 

 

и, следовательно,

 

 

 

 

 

 

 

 

 

^ - = Ф ,

(Л, +

ца) -

Ф, (Л,) +

Ф2 (Л, +

т )

= .

 

 

^

^ - а

+ ФЛЛО +

^ С Л , , ца),

 

0)

причем

здесь

 

 

 

 

 

 

 

 

 

dQx\dA

— квадратная

матрица

\\дФ(1)/дА<к> ||;

4=1,

2

п.

 

 

 

 

 

 

 

I/,

Ф3 (Л„ pa) =

j -

1 [Ф, (Л, + И

-

Ф, (Л,) _

а ] +

 

 

 

 

 

 

 

 

 

 

(2)

116

4 . 4 . 1 .

Введем теперь

t

о

тогда

^ I ~й~Аа

(3)

Если при определении а4 ограничиться величинами ну­ левого порядка малости, что при вычислении Л дает погрешность на всем интервале 0 ^ / < L / | x , меньшую, чем C\i2, то последнее слагаемое в (3) можно упростить следующим образом: величина, стоящая в (2) в квадрат­ ной скобке, имеет порядок ц.2 и, следовательно, даже после деления на ц. согласно изложенному ранее при­ надлежит к классу MQ. Кроме того, если Фг имеет огра­ ниченные частные производные по компонентам Л, по­ следнее слагаемое (2) можно представить в виде

Фа ( Л + ца) Фг (А,)

_дфг

[j.

дА

где i|? ограниченная функция при всех Л ; + дл, находя­ щихся в пределах D, и при рассматриваемых t. Тогда согласно сказанному ранее цдр принадлежит к классу М0. Таким образом, укороченное уравнение для а можно написать так:

~diH > [ $ - + & ] [ « . + J « . h .dt. o

(4)

1

5

 

решение которого ищется при условии ai(0)=0.

 

Приведенные выше оценки,

конечно, имеют

силу

лишь при условии ограниченности всех частных произ­ водных от (Di и Ф2 до второго порядка включительно в D и для всех рассматриваемых значений t. Таким

образом, можем

написать

0f Ф2 [A, (r), t) dt 1

 

Л =

Л, + |» a, +

 

4 . 4 . 1 .

 

J

117

 

 

(з'десь уже при интегрировании At(t) рассматривается как функция от t, а не как постоянная). При этом по­ грешность не превосходит величину, имеющую порядок малости и2 .

Если оставить в стороне вопрос о вычислительных трудностях, можно к уравнению (1) вновь применить аналогичную процедуру и получить более высокую сте­ пень точности, причем на основании сформулированной выше основной теоремы можно будет утверждать, что погрешность результата при аналогичных условиях бу­ дет величиной порядка малости ц.3. Отметим еще обстоя­ тельство, полезное для дальнейшего. Если рассматривать конечные интервалы t, не зависящие от ц, то величина а определяется с ошибкой порядка ц2 , а не и, как это предполагалось раньше. Действительно, члены, которые

мы отбросили в (3)

при переходе

к

( 4 ) , имели

вид

[i2ty(t),

где ty(t)—ограниченная

в рассматриваемом

ин­

тервале

функция. Если

0 ^ / < Г ,

где

Т

не зависит от р.,

мы можем, не меняя результатов, относящихся к этому интервалу, продолжить его до L/u. и выбрать при этом 1|з(0=0 при t>T. Тогда

и,

следовательно, величина М,

входящая

в неравенства

(3)

§ 4.1, будет иметь вид

|л| , ф|тах7'

и

согласно

(6)

|r)| ^ C | i = Cip2 , где Ci — константа, не

зависящая от

и.

Заметим еще, что сама величина а в этом случае будет иметь порядок малости ц.

4.5. О применении метода ММА к уравнениям с медленно меняющимися коэффициентами

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

Пусть дано уравнение

(/С> + а п _ , ( / 0 - < > +

. . .

+a0y=f(t),

(1)

где коэффициенты as(t)

являются

медленными функциями от

t

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

являют­

ся малыми величинами.

 

118

§4 . 5 .

Подобно тому, как это делалось раньше в Начале настоящем главы, умножим обе части рассматриваемого уравнения на е^? 1 и проинтегрируем по t в пределах от 0 до со. Можем написать

оо

 

 

d

d

^asy(^e~

rtdt

= — (/(«->) (0) а, (0) -

j V s

- •) -jj (a,e~ *') dt =

о

 

б

 

 

 

{/(•' - ' ) (0) as (0) + (/f»-« (0)

d

 

-

- j j -

(a.*- P')f=o +

 

 

и

 

 

и, продолжая

по аналогии, найдем

 

 

00

 

5

 

 

^a.yi')e-PUt

= JJjrf'-4(0) ^ г т "

( - I ) " +

о

 

ft=l

 

 

 

 

00

 

 

о

и, следовательно, после преобразования уравнение приобретает вид

ооп

j W ) £ (-1 )8 iK(a-e'vi)

dt+

Оs=0

п

S

00

Учитывая, что

-Zf(be-")

= -e-"{p:-D)

Ь,

где b — произвольная функция от t; D — оператор дифференцирова­ ния, можем написать

-jjjrVe-*)=-D[e-rt(p-D)

b] =

e-'4p{p-D)~

— D(p — D)\ Ь=е-Р'(р>—

D)2b

и, продолжая, найдем

 

 

' ~ ш ( 6 е - р , ) = e~pt

( _ 1 ) 8 ( р ~ D)'b'

§4 . 5 .

119

 

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