книги из ГПНТБ / Конторович М.И. Нелинейные колебания в радиотехнике. (Автоколебательные системы)
.pdfв правой |
части |
(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, |
ни от д.). |
|
|
|
|
|
|
|
Очевидно, что при таком определении термин «об ласть D» означает лишь множество величин Л</;> (или векторов Л), удовлетворяющих ( 4 ) . Если хотя бы для одной из величин AW неравенство (4) превращается в равенство, мы будем говорить, что вектор Л лежит на границе области D.
Условимся в дальнейшем множество функций, удов летворяющих условию (3), обозначать как М0 и писать, например, так: Ф2^М0.
В соответствии с тем, что говорилось в предыдущих параграфах, естественно составить укороченные уравне
ния, соответствующие системе (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 |
|