Асимптотические методы решения дифференциальных уравнений

Асимптотические методы решения дифференциальных уравнений

Содержание:

По этой ссылке вы найдёте полный курс лекций по математике:

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

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

Следствие теоремы 2. Пусть при (t,x) выполнены условия теоремы 2, и при решение задачи (1) проходит в области D; t0 £ [tl9t2]- Тогда решение x(t, р) задачи (1) при t разлагается по формуле Тейлора по степеням параметра ц до цт включительно: Асимптотические методы решения дифференциальных уравнений Здесь х(Ь9ц) и v^t) — n-мерные вектор-функции, v0(t) = 0) есть решение системы (1) при \i = 0, оно считается известным.

Чтобы найти надо подставить разложение (12) в систему (1) и начальные условия, и разложить правые части по степеням /х до /хт включительно. Далее надо приравнять коэффициенты при одинаковых степенях /х. Получается для vx. 9vm система дифференциальных уравнений с начальными условиями. Последовательно решая уравнения системы и пользуясь начальными условиями, находим . vm(£). Пример 2. Найти разложение решения задачи по степеням параметра /х до /х2 включительно.

Решение примера. Правая часть уравнения в области х > 0 име- ет производные любого порядка по ж,/х. Условия теоремы 2 выполнены для любого т, пока решение задачи (13) с \l = О проходит в области х > 0. При = 0 задача (13) принимает вид dx/dt = t/x, 1, и имеет решение ж(*) = t, оно проходит в области х > 0 при t > 0. Поэтому v0(*) = t (t > 0).

Разложение х = t + /iv1 +/i2v2 + о(/л2) подставляем в уравнение и начальные условия (13), члены порядка o(/i2) не пишем. Разлагаем дробь в (14) по степеням /х, члены с цк, к > 2, не пишем. Подставляем это в (14) и приравниваем коэффициенты слева и справа при одинаковых степенях параметра /х: при v\ = -^-2t\ М0 = 4; (16> при Ц = + «2(l) = i (17) Здесь начальные условия получены из (15). Все дифференциаль- ные уравнения для v<9. vm всегда линейные. Из (16) получаем Vj = — у.

Возможно вам будут полезны данные страницы:

Подставляя это в (17), находим v2 = + fj. Итак, Так как условия теоремы 2 выполнены для любого m ^ 2, то следующий член разложения имеет вид p3v3(t) и, не находя vy в (18) вместо o(/i2) можно написать 0(fi3). Задачи для упражнений: [ | 3« | Отыскание периодических решений. Нижеследующие лемма 2 и теорема 3 дают условия существования периодических решений соответственно для линейной системы с периодической правой частью и для нелинейной системы, близкой к линейной, и указывают методы отыскания таких решений.

Лемма:

Пусть при 0 вектор-функция x(t) — решение уравнения х’ = f(t, х), где вектор-функция f и все dfjdxj непрерывный f(t+p, х) = f(t, х). Если х(р) = ж(0), то решение x(t) продолжается на интервал (-оо, оо) с периодом р. Доказательство. Так как Асимптотические методы решения дифференциальных уравнений то продолженная с периодом р функция x(t) 6 С1.

Она всюду удовлетворяет данному уравнению, ибо для любого k € Z имеем Лемма 2. Если для всех собственных значений матрицы А имеем то система х1 = Ах + f(t) для каждой непрерывной функции f(t) с периодом р имеет (и только одно) решение с периодом р. Условие (19) называется условием отсутствия резонанса. Доказательство. Пусть v(£) — частное решение данной системы с v(0) = 0. В силу теоремы 5 § 9 и следствия 1 § 15 общее решение имеет вид х = etAb + v(t), где Ь — произвольный вектор из Rn.

Чтобы это решение имело период р, по лемме 1 надо, чтобы х(р) = ж(0). То есть Это — линейная алгебраическая система относительно неизвестных координат вектора Ь. Для существования единственного решения достаточно, чтобы det (ерА — 1 • Е) Ф 0, то есть чтобы матрица не имела собственных значений, равных 1. Если АР. АП — собственные значения матрицы А> то согласно замечанию в имеет собственные значения Для А = а + /3i имеем = е?» (cos р/3 + tsinp/5). Это число равно 1 только в случае а = 0, рр = 2*кку к = 0,±1,±2.

Поэтому при условии (19) имеем Теорема 3. Пусть функции f(t), g(t9 х, р) непрерывны при имеют период , где т ^ 1. Пусть выполнено условие (19) и решение x°(t) с периодом р уравнения х Ах + f(t) содержится в области D. Тогда при всех достаточно малых |/*| система имеет решение периода р по t, стремящееся Такое решение единственно и принадлежит классу Ст по ц. -7- Доказательство.

Пусть х(Ь\Ъ,ц) — решение системы (20) с начальным условием ж(0; ц) = Ь. По лемме 1 оно будет иметь период р, если Докажем, что при малых \l существует Ь € R», удовлетворяющее уравнению (21). Функция ж(р;Ь, ц) Е Ст по b, \i в силу теоремы 2. При /х = 0 уравнение (20) линейное, как в лемме 2, уравнение (21) принимает вид и имеет единственное решение Ь. Далее, якобиан левой части равенства (21) по координатам Ъ<9. Ьп вектора b при /х = 0 совпадает с детерминантом (22), значит, не равен нулю.

Тогда по теореме о неявных функциях уравнение (21) при достаточно малых имеет решение Ь = Ь(ц), стремящееся к Ь° при такое решение единственно и Ъ(ц) 6 Ст.

Пусть дано уравнение с постоянными коэффициентами а,- и непрерывными функциями /, д периода р по t и д € Ст по у, /1, а корни А характеристического уравнения удовлетворяют условию (19). Тогда для отыскания решения периода р не нужно переходить от уравнения (24) к системе, можно сразу отыскать решение в виде (12), где теперь v^t) — скалярные функции с периодом р. 1 Пример 3. Найти с точностью о(ц2) периодическое решение 1 | уравнения [ Решение примера.

Здесь р = 2т, А2 + 3 = О, Л = ±iV3 Ф 2*ki/p = ki (к Е Z), условие (19) выполнено. Ищем периодическое решение в виде х = v0 + /avx + p2v2 + . Подставляя в уравнение (25) и приравнивая коэффициенты при одинаковых степенях /*, получаем систему уравнений . Надо найти решения t;0, Vj, v2 с периодом 2т. Для каждого из этих уравнений надо найти лишь частное решение (методой неопределенных коэффициентов), так как по теореме 3 при выполнении условия (19) решение с периодом р единственно.

Замена х = а?0 + ру дает . Так как F(x°) = 0, то по формуле Тейлора Остаточный член г Е С (ибо другие члены в равенстве принадлежат Cm+1), г = ц2д(у, ц). Получаем систему вида (20) Если собственные значения матрицы А удовлетворяют условию (19) (нет резонанса), то по теореме 3 система (27) при достаточно малых |/*| имеет решение с периодом р. I Пример 4. Рассмотрим уравнение Решение примера. При \i = 0 положения равновесия х< = 1 и х2 = -1.

Периодическое решение

Найдем периодическое решение, близкое к х = 1, Замена х = 1 + /ху Дает Здесь р = 2х, А = -1 ± i Ф 2*ki/p (к € Z), условие (19) выполнено. Поэтому при малых ц уравнение (29) имеет решение периода 2тг и вида у = v0($) + fivx(t) +..где все имеют период 2т. Подставляя это в (29), получаем, как в примере 3, Асимптотические методы решения дифференциальных уравнений Отсюда находим . Следовательно, Уравнение (28) при малых ц имеет и другое решение с периодом 2х.

Оно близко к неустойчивому положению равновесия х2 = -1 и отыскивается аналогичным способом. Можно доказать, что оно неустойчиво. ч Задачи для упражнений: [12], § 18, № 1079-1083. | 5« | Естественно возникает вопрос, в каких случаях разложения по степеням параметра /х, полученные в следствиях теорем 2 и 3, можно продолжить до бесконечного ряда Тейлора, сходящегося к искомому решению при малых /1.

Пуанкаре об аналитической зависимости решения от параметра, см. [13], гл. 1, §6, теорема 1.3 и [2], гл.6, §2, теорема 6.2.1′ и §3, п. 1. I О методах исследования устойчивости периодических решений, получаемых методом малого параметра, см. [2], гл. 7, § 3 и [33], гл. 3, § 10-15. В частности, при условиях теоремы 3 асимптотическая устойчивость периодического решения при достаточно малых обеспечена, если для матрицы А все собственные значения А.

Имеют Re At. а неустойчивость — если есть хотя бы одно Re At. > 0. Поэтому в примере 4 при малых /1 решение (30) асимптотически устойчиво, а периодическое решение, близкое к х = -1 (для него А = — 1 ± л/3), — неустойчиво. Метод отыскания периодических решений при резонансе, то есть когда условие (19) не выполнено, изложен в [13), гл. 2, §8, пункты 2 и 3; примеры там же; в [2], гл. 5, §3, п. 2 и в [33], гл. 2, §6, 7. Об отыскании периодических решений автономной системы х’ = Ах -f pf(x,p) в случае, когда при р = 0 периодическое решение известно, см. [13], гл.2, §8, пункт 4; [2], гл.5, §3, п. 3 и [33], гл.2, §11-13.

Метод малого параметра применялся к широкому кругу задач, в частности, в [33], главы 4-8. Методы последовательных приближений для уравнений с малым параметром разработаны в [24]. Существенно отличным от предыдущих является случай, когда малый параметр является множителем при производной, например, fix = /(*, у), у =д(х, у). Здесь нет непрерывной зависимости от /1 при р 0, и решения имеют другие свойства, см., например, [13], гл.4, §6 и [15], гл.10, §3,4. Известно много работ, в которых подробно исследуются такие случаи.

Присылайте задания в любое время дня и ночи в ➔

Официальный сайт Брильёновой Натальи Валерьевны преподавателя кафедры информатики и электроники Екатеринбургского государственного института.

Все авторские права на размещённые материалы сохранены за правообладателями этих материалов. Любое коммерческое и/или иное использование кроме предварительного ознакомления материалов сайта natalibrilenova.ru запрещено. Публикация и распространение размещённых материалов не преследует за собой коммерческой и/или любой другой выгоды.

Сайт предназначен для облегчения образовательного путешествия студентам очникам и заочникам по вопросам обучения . Наталья Брильёнова не предлагает и не оказывает товары и услуги.

Курсовая работа: Асимптотика решений дифференциальных уравнений

по дисциплине “Эффективные алгоритмы исследования моделей естествознания”

на тему: «Асимптотические решения дифференциальных уравнений по малому параметру. Регулярные возмущения»

Применения регулярного возмущения

1. Асимптотическое поведение решений дифференциальных уравнений с малым параметром

1.1 Асимптотическое поведение решений системы

2. Регулярные возмущения

2.1 Асимптотические методы

2.2 Регулярные возмущения решений задачи Коши для обыкновенных дифференциальных уравнений

2.3 Существование решении возмущенной задачи

Невозможно представить себе современную науку без широкого применения математического моделирования. Сущность этой методологии состоит в замене исходного объекта его «образом» — математической моделью — и дальнейшем изучении модели с помощью реализуемых на компьютерах вычислительно-логических алгоритмов. Этот «третий метод» познания, конструирования, проектирования сочетает в себе многие достоинства как теории, так и эксперимента. Работа не с самим объектом (явлением, процессом), а с его моделью дает возможность безболезненно, относительно быстро и без существенных затрат исследовать его свойства и поведение в любых мыслимых ситуациях (преимущества теории). В то же время вычислительные (компьютерные, симуляционные, имитационные) эксперименты с моделями объектов позволяют, опираясь на мощь современных вычислительных методов и технических инструментов информатики, подробно и глубоко изучать объекты в достаточной полноте, недоступной чисто теоретическим подходам (преимущества эксперимента).

Неудивительно, что методология математического моделирования бурно развивается, охватывая все новые сферы — от разработки технических систем и управления ими до анализа сложнейших экономических и социальных процессов.

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

Применения регулярного возмущения

Выходные пучки лазеров часто имеют квазирегулярную модуляцию волнового фронта (ВФ). В газовых лазерах с движущейся активной средой такую модуляцию могут вызывать неоднородности, возникающие под действием периодической сопловой решетки [1], под влиянием страт и доменов в газовом разряде [2], в результате наложения ударных волн [3,4], а также под действием ряда других физических факторов. Модуляция ВФ выходных лазерных пучков в литературе чаще всего рассматривается как фактор, влияющий, прежде всего на расходимость излучения.

Гораздо меньше внимания уделяется анализу метаморфоз структуры ВФ, условиям появления и взаимосвязи каустических и фазовых дислокационных образований в лазерных пучках. Такого рода образования регистрируются в излучении лазеров с самыми разными оптическими резонаторами [5,6]. В настоящей работе рассматриваются качественные изменения амплитудно-фазовой структуры лазерных пучков, первоначально обладающих плавной регулярной модуляцией ВФ.

Общее представление о характере рассматриваемых процессов можно получить на примере известной задачи [7] о распространении безграничной волны, фаза которой в начальной плоскости меняется по гармоническому закону. Амплитуда такой волны имеет следующий вид:

где m — параметр, характеризующий глубину фазовой модуляции; х — поперечная координата; а — период модуляции. На расстоянии z от начальной плоскости поле можно представить в виде суперпозиции плоских волн [7]:

где — функция Бесселя порядка — волновое число. Это поле является частным случаем самовоспроизводящихся полей, свойства которых нашли применение в лазерной технике [8,9]

Используя для расчета характеристик поля его разложение по плоским волнам (2), а также лучевой метод из работы [10], можно установить основные особенности трансформации первоначального распределения амплитуды и фазы. Расчеты показывают, что даже при малой глубине модуляции фазы и равномерном распределении интенсивности в начальной плоскости дифракционные эффекты приводят к значительному пространственному перераспределению интенсивности.

Перераспределение наиболее заметно вблизи плоскостей Эти плоскости располагаются между плоскостями, в которых, согласно эффекту Тальбо, воспроизводится первоначальное равномерное распределение интенсивности. Так, при m = 0.1 контраст картины распределения интенсивности , а при m = 0.5 контраст К = 2.82. С превышением определенной критической глубины модуляции в структуре волны происходят качественные изменения..

На рис.1 приведены распределения амплитуды А и фазы Ф на расстояниях при разных первоначальных глубинах модуляции фазы. Видно, что при превышении критической глубины модуляции появляются линии с нулевыми амплитудами. В распределении фазы им соответствуют КД, обусловленные скачкообразным изменением фазы на π. КД располагаются симметрично относительно осей клювообразных каустик. Клювы каустик, находящиеся сначала вблизи плоскостей , с дальнейшим увеличением глубины фазовой модуляции приближаются к плоскостям воспроизведения первоначальной структуры. При этом растет и число КД.

Их расположение по отношению к образующим каустик соответствует рассчитанной на основе интеграла Перси фазовой структуре поля, приведенной в работе [11].

Рис. 1. Распределение амплитуды А (1,2) и фазы Ф (3,4) по поперечной координате х для безграничной волны на расстоянии стрелками указано положение клювов каустик.

Продольная структура распределения интенсивности излучения показана на рис. 2 для т = 1.2. Из него видно, что фазовая модуляция вызывает формирование каналов, вытянутых вдоль направления распространения, в которых интенсивность излучения существенно превышает среднюю. Оси этих областей совпадают с осями симметрии клювообразных каустик.

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

Здесь функция ) совпадает с функцией при замене поперечной координаты х поперечной координатой у С —константа.

Структура эквифазных линий в начальной плоскости на рис.3,а построена с помощью формулы (3) для С = 0.2 и m = 2. Ход линий свидетельствует о наличии плавных регулярных возмущений волнового фронта. На рис.3,6 изображена структура эквифазных линий на расстояниях

ВД располагаются в точках пересечения эквифазных линий. Они образуют своеобразные квадруполи каждый из которых состоит из четырех ВД. Две из них имеют положительный знак (являются «правыми»), две — отрицательный знак (являются «левыми»). Квадруполи окружают оси каустик.

В отличие от КД каждая из которых строго говоря формируется в определенной плоскости z = const, ВД характеризуются определенной продольной длиной. Как и КД дислокации винтового типа возникают лишь при превышении глубиной первоначальной модуляции волнового фронта некоторого критического значения. Если обозначить через разность между максимальной и минимальной фазами в начальной плоскости (при модуляции по одной координате совпадает с ni) то ВД будут возникать когда >

Все вышеперечисленные эффекты были проанализированы применительно к пространственно-ограниченному пучку с гауссовым профилем распределения интенсивности. В основу расчета была положена формула (2), в которой суперпозиция плоских волн была заменена системой распространяющихся под углом друг к другу raусовых мод свободного пространства [12]. Горловины мод располагались в начальной плоскости

Расчеты показали, что переход к более точной модели гауссова пучка с периодической модуляцией ВФ не вносит существенных качественных изменений в данные о преобразовании амплитудно-фазового распределения, по крайней мере на расстояниях сопоставимых с характерной длиной . Как и в случае безграничной волны дислокации ВФ начинают формироваться в ближней зоне, когда глубина модуляции фазы превышает . Сказанное иллюстрирует рис. 4, который является аналогом рис. 1 для гауссова пучка. Отношение радиуса пучка в горловине к периоду модуляции, а равно пяти.

Из сравнения рис. 1 и 4 видно, что имеющиеся в них различия проявляются в дифракционном «замывании» части дислокаций. Различия усиливаются с ростом координаты z по мере того, как ухудшается периодическая воспроизводимость первоначальной структуры поля. Это видно, в частности, из рис. 5, на котором изображено продольное распределение интенсивности для параметров

Распределения, приведенные на рис.2 и 5, близки лишь в ближней зоне, для которой характерны узкие зоны, где концентрируется энергия светового потока. В дальней зоне дифракции перекрытие гауссовых угловых компонент излучения ослабевает, и структура излучения кардинальным образом отличается от структуры безграничной волны: излучение представляет собой «веер» пучков, интенсивность которых убывает с увеличением угла наклона.

Фазовая модуляция гауссова пучка по двум поперечным координатам, если ее глубина превышает указанную выше критическую глубину, приводит к появлению на волновом фронте БД. Как и в безграничной волне, эти БД обладают определенной продольной длиной, увеличивающейся с ростом глубины модуляции. Это свойство БД значительно облегчает их экспериментальное обнаружение. В дальней зоне дифракции вследствие изменения фазы в начальной плоскости по двум координатам будут формироваться два веера пучков, располагающихся во взаимно перпендикулярных плоскостях.

Заметим в заключение, что результаты выполненного анализа могут быть частично перенесены и на случай нерегулярной плавной модуляции ВФ, если длина рассматриваемой пространственной области сопоставима с величиной аЦк, где ап — характерный размер нерегулярных возмущений ВФ. В частности, это относится к образованию в световом поле каналов с повышенной интенсивностью и к появлению дислокаций волнового фронта при превышении фазовыми возмущениями определенного значения.

Таким образом, плавные возмущения ВФ играют важную роль в трансформации амплитудно-фазового профиля излучения и в формировании каустических и дислокационных образований. Появление каустик и дислокаций волнового фронта носит пороговый характер и непосредственно связано с глубиной первоначальной модуляции фазы. Для практики важным является то, что появление указанных образований в лазерном пучке сопряжено с формированием узких каналов, в которых интенсивность излучения значительно превышает среднюю.

Работа выполнена при финансовой поддержке государственной научно-технической программы «Физика квантовых и волновых процессов» (проект 1.61) и физического учебно-научного центра «Фундаментальная оптика и спектроскопия».

1. Асимптотическое поведение решений дифференциальных уравнений с малым параметром

Многие колебательные системы описываются дифференциальными уравнениями с малым параметром при производных:

или, в векторной форме

где — малый положительный параметр, — неизвестные функции времени t, характеризующие данную систему.

В работах ( х ) — ( 5 ) находится асимптотика решений системы (1.1) в случае, когда при каждом z любое решение системы «быстрых движений» **

при приближается либо к устойчивому положению равновесия, либо к устойчивому предельному циклу.

Но возможны случаи, когда система «быстрых движений» (1.2) может не иметь асимптотически устойчивых положений равновесия и изолированных предельных циклов. Такова, например, гамильтонова система. Целью настоящей работы и является изучение этих случаев. Так, в § 2 с точностью до величин порядка О (г) находится решение системы (1.1), для которой соответствующая система «быстрых движений» гамильтонова и к = 2, т. е. находится решение системы

Асимптотические формулы для решения этой системы находятся для области, где траектории соответствующей гамильтоновой системы «быстрых движений» при каждом векторе z замкнуты (в случае невырожденного центра в рассматриваемую область включается и сам центр). Метод исследования системы (1.3) таков: сначала рассматривается система «быстрых движений» (1.4), а затем система (1.3) после соответствующей замены переменных усредняется вдоль решений (1.4). Оказывается, что уравнение с малым параметром и. при старшей производной и с пропущенной в основном члене Q (п — 1)-й производной, исследованное В.М. Волосовым (при п — 2 — в работе ( 12Г ), при F

О — в ‘работах ( 8 ) — ( п )) методом конечных разностей, является частным случаем системы (1.3). Поэтому результаты работ ( 8 ) — ( 12 ) (эти результаты сформулированы в § 3 настоящей работы) следуют из результатов § 2.

Метод построения решения уравнения (1.5) при п = 2 с любой наперед заданной точностью в случае, когда известно общее решение (в форме разложения в тригонометрический ряд Фурье) соответствующего невозмущенного уравнения был дан в работе Ю.А. Митрополъским.

Задача исследования системы (1.3) с точки зрения работ ( 3 ) — ( 4 ) и вывода из нее известных результатов В.М. Волосова [работы ( 8 ) — ( 12 )] относительно уравнения (1.5) была поставлена Л.С. Понтрягиным в его докладе на семинаре В.И. Смирнова в Ленинграде в середине апреля 1957 г.

Выражаю глубокую благодарность Л.С. Понтрягину за ценные указания, советы и постоянное внимание к настоящей работе.

1.1 Асимптотическое поведение решений системы

Система (1.3) в векторной форме имеет вид:

глк, в быстром времени

При е = 0 система (2.1′) переходит в гамильтонову систему

являющуюся системой «быстрых движений» для системы (2.1). 1. Изучение системы (2.2). Пусть функции

определены и непрерывны вместе со всеми своими первыми частными производными в некоторой области G эвклидова пространства E ₂ + i переменных х, у, zi . zi . Как известно, система (2.2) имеет первый интеграл

и (2.3) представляет собой семейство всех фазовых траекторий системы(2.2) на кажтгой плоскости z = constобласти G .

Возьмем некоторую точку (х, у, z ) из G, не являющуюся положением равновесия системы (2.2). По известной теореме существования и единственности решений системы обыкновенных дифференциальных уравнений, через эту точку пройдет только одна фазовая траектория системы

(2.2). Уравнение этой траектории запишется в виде:

Докажем следующее утверждение.

Пусть траектория (2.4) замкнута и целиком лежит внутри области G . Тогда в пространстве E ₂ + i существует некоторая окрестность G этой траектории (2.4) такая, что

1)фазовые траектории системы (2.2), проходящие через точки G , замкнуты и целиком лежат в G;

2)уравнение (2.3) при каждой паре (/г, z ) определяет одну и только одну фазовую траекторию системы (2.2), расположенную в G;

3)на каждой фазовой траектории (2.3) системы (2.2), лежащей в G, можно выбрать по одной точке , гладко зависящей от

В самом деле, в силу известных свойств гамильтоновой системы, в пространстве E ₂ + i существует некоторая окрестность G траектории (2.4) ( Gd G), в которой выполняется условие 1). Выделим из G ту окрестность траектории (2.4), в которой выполняются и условия 2), 3). Для этого возьмем поверхность, пересекающую каждую плоскость z = const области G

о о по нормали в точке (х, у, z ) к фазовой траектории системы (2.2), проходящей через эту точку. Уравнение этой поверхности имеет вид:

Следовательно, точка (х, у, z , h ) эвклидова пространства £»₂ +z переменных х, у, z , h удовлетворяет системе

Левые части системы (2.5) определены и непрерывны вместе со всеми своими частными производными в области Г: (#, у, z ) £ G, —ос 0 такое, что при любых eg (0, е₀ ], t 6 [*<>> L ] решение <ф ( t , е), h ( t , е), z ( t , г)> системы (2.15) с начальными условиями

и решение <ф ( t , e), h ( t ), z ( t )> усредненной системы (2.17) с теми же начальными условиями

связаны следующим образом: точка <h ( t , e), z ( t , г)> остается в некоторой и выполняются соотношения:

окрестность решения)). А так как, по (2.13),

и так как точка <h (£, е), z ( t , е)> остается в G_hp CZG_h , то на отрезке [ t_Q , L ] при любом 8 g (0, е₀ ] решение <х (t , е),?/ (£, е), z (£, г)> системы (2.1) остается в G, причем, по свойству 3,

и потому соотношения (2.39), (2.40) доказывают первую часть теоремы 1. Докажем вторую часть теоремы 1. По формуле конечных приращений, из (2.41) получаем:

Возникает вопрос, как ведут себя решения системы (2.1) во всей указанной окрестности Go (включая и положения равновесия g (z), z> системы (2.3)). На этот вопрос отвечают теорема 1 и нижеследующие теоремы 2 и 3.

ТЕОРЕМА 2. Пусть в окрестности Go выполнены условия теоремы 1, касающиеся гладкости правых частей системы (2.1). Тогда найдется число 8° у> О, такое, что при любом г £ (0, е°] (е° и Н

сохраняет вид системы (2.1), но дает условия (2.54). Следовательно, в силу.

Это решение на конечном промежутке времени [ t ₀ , L ] составляет некоторое замкнутое ограниченное множество F_Q CZG ₀ и поэтому найдется ро > 0 такое, что G₀₀ С G ₀ ( G_{Q 0} — р₀ -окрестность F ₀ ).

Следовательно, по формуле Тейлора, примененной к функциям

(формула Тейлора применима в G₀₀ относительно х, у, так как прямолинейный отрезок, соединяющий любые две точки (я, у, z ) и (0, 0, z ) из Goo, содержится в Goo, поскольку каждое сечение области G₀₀ плоскостью z = const представляет собой круг с центром в точке (0, 0, z ), по определению Goo).

Функция О 2 (х, у, е), в силу указанной в условиях теоремы гладкости правых частей системы (2.1), является однородной квадратичной относительно х, у, е с ограниченными в Gooкоэффициентами, и поэтому

С другой стороны, по формуле Тейлора, в силу (2.54) имеем в G₀₀

то соотношение (2.61), в силу (2.57), дает на [£₀ , t (е) ]:

Но, по (2.56) — (2.58) и (2.63),

откуда следует, что на отрезке

Но так как, в силу

т. е. окончательно, по (2.64), (2.67),

2. Регулярные возмущения.

2.1 Асимптотические методы

Пусть задано банахово пространство и отображение .

Определение . Будем ряд называть асимптотическим рядом для функции , если для любого найдутся числа и такие, что

при (2.1)

Пример 1 . Если функция имеет производные всех порядков в точке , то справедливо формула Тейлора

(2.2)

Ряд Тейлора может расходиться на любом отрезке , но он будет асимптотическим рядом для функции . Действительно,

(2.3)

Пример 2 . Рассмотрим функцию

Интегрируя по частям, получаем

Ряд расходится при любом , но является асимптотическим для функции , так как

Замечание . Асимптотический ряд может быть полезен при вычислении значений функции при малых или больших значениях параметра .

Рассмотрим функцию примера 2. Вычисляя интеграл численно, получаем при

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

0.0015633, -0.0004366, 0.0001633, -0.0000766, 0.0000433, -0.0000287, 0.0000217,

-0.000186, 0.0000177, -0.0000186, 0.00002133, -0.0000266, 0.0000357, -0.0000515,

0.0000793, -0.0001299, 0.0002257, -0.0004145, 0.0008020

Наилучшее приближение дает девятая частичная сумма.

На рис. 1 изображен графически характер приближения частичных сумм к значению. На горизонтали оси откладывается номер , по вертикали частичная сумма .

Пусть банаховы пространства и при задано семейство операторов . Рассмотрим при уравнение . Допустим, что это уравнение при каждом имеет единственное решение . Уравнение будем называть вырожденным . Допустим, что вырожденное уравнение имеет единственное решение . Будем говорить, что вырождение регулярное, если

при (2.4)

Если (18.4) не выполняется, то говорят, что вырождение сингулярное.

Распространена еще и такая терминология: Уравнение называют уравнением возмущений для уравнения . Если условие (2.4) выполнено, то говорят о регулярных возмущениях. В противном случае речь идет о сингулярных возмущениях. Сам термин «теория возмущений» возник в рамках небесной механики. В следующем параграфе будет исследована задача о регулярных возмущениях для обыкновенных дифференциальных уравнений.

2.2 Регулярные возмущения решений задачи Коши для обыкновенных дифференциальных уравнений

Рассмотрим задачу Коши

(2.2.1)

Функция непрерывна по переменной и бесконечно дифференцируемая по переменным и при , , .

Предполагается, что вырожденная задача

(2.2.2)

имеет единственное решение при , причем .

(2.2.3)

и воспользовавшись тем, что функция удовлетворяет уравнению (2.2.2) запишем систему уравнений для функции в виде

(2.2.4)

(2.2.5)

(2.2.6)

Будем искать решение задачи Коши (2.1.4) в виде формального ряда по степеням малого параметра

(2.2.7)

Для определения неизвестных функций получаем рекуррентную систему задач Коши для линейных уравнений (уравнений в вариациях)

(2.2.8)

Уравнение (2.2.8) называют уравнением в вариациях.

Вычислим две первых функции

(2.2.9)

Подставляя разложения (2.2.7) и (2.2.8) в уравнения (2.2.4),получаем рекуррентную систему уравнений

(2.2.10)

Все уравнения (2.2.4) имеют одинаковую структуру

, (2.1.11)

Столбцы фундаментальной матрицы образуют фундаментальную систему решений. При помощи формулы Коши получим решение в виде

(2.2.12)

Линейный оператор

(2.2.13)

Покажем, что ряд (2.2.3) асимптотический для решения . Положим

(2.2.14)

Применяя формулу Тейлора, получаем

(2.2.15)

где функции те же, что и в формуле (19.8), а

(2.2.16)

Подставляя представление (2.2.14) в уравнение (2.2.4), воспользовавшись представлением (2.2.15) и формулами (2.2.8), получаем уравнение для функции .

(2.2.17)

(2.2.18)

Из формулы (2.2.6) получаем

и формула (2.2.18) может быть записана в виде

(2.2.19)

Так как вторые производные функции ограничены, то функция удовлетворяет условию Липшица и

(2.2.20)

Вспоминая определение оператора , получаем функциональное уравнение

(2.2.21)

Используя принцип сжатых отображений, покажем, что уравнение (2.1.21) при имеет единственное решение, и справедливо неравенство . Тем самым будет доказано, что ряд является асимптотическим рядом для функции , являющейся решением задачи Коши (2.2.1).

Пусть . Так как частные производные равномерно непрерывны, то из (2.2.17)- (2.2.20) получаем оценки

при . Таким образом, шар радиуса отображается в себя при.

Используя (2.2.20), получаем

Используя равномерную непрерывность частных производных, получаем

Уменьшая, если нужно, получаем, что при оператор является оператором сжатия. Следовательно,

и ряд асимптотический для решения задачи Коши (2.1.1).

2.3 Существование решении возмущенной задачи

Результаты, полученные обладают той особенностью, что справедливость асимптотического представления гарантируется на некотором сегменте [0,T], определяемом свойствами правой части (2.3.1), одновременно с существованием и единственностью как невозмущенного, так и возмущенного уравнений.

Можно ставить вопрос иначе. Допустим, что решение невозмущенной задачи (2.3.2) существует, единственно и принадлежит некоторой области G пространства переменных y(t,μ) при, 0≤t≤T.Величину T в данном случае можно, например, установить непосредственно из явного вида y(t). Будет ли при достаточно малых μ решение задачи (2.3.1) также существовать на всем [0,T] и подчиниться формуле (1.3)? Ответ на этот вопрос дает следующая

Теорема 1 .2. Пусть в области

непрерывны и равномерно ограничены:

Пусть решение y ( t ) задачи (2.3.2) существует, единственно на [0, T ] и принадлежит . Тогда при каждом достаточно малом μ решение y ( t , μ) задачи (2.3.1) также существует, единственно на [0, T ] принадлежит G , и имеет место равномерный относительно предельный переход

(2.3.8)

Доказательство. Перейдем в (2.3.1) к новой неизвестной функции . Имеем

Перейдем к эквивалентному интегральному уравнению

(2.3.9)

где причем . Здесь и в

дальнейшем бесконечно малые при μ →0 величины будем обозначать

ω(μ), ω₁ (μ) и т. д. Применим к уравнению (2.3.9) метод последовательных приближений и докажем, что ▲(t,μ) существует на сегменте [0,Т] и .Это очевидно, равносильно утверждению теоремы 1.2.

Построим последовательные приближения обычным образом

Предварительно заметим, что так как y=y(t) принадлежит G для кривая , где при достаточно малом μ. также принадлежит G для

Положим Тогда

(2.3.10)

В равномерной сходимости последовательности ( k ) ▲ к решению ▲(t,μ) уравнения (2.3.9) можно убедиться совершенно, может в пределе при k→∞ появиться равенство. Поэтому , что равносильно (2.3.8).

Замечание. Теорема доказана для скалярного случая, но аналогичное утверждение справедливо и для случая, когда y — вектор.

2.3.2 Теорема 2.3.2 остается справедливой, если имеет место возмущение не только в уравнении, во и в начальных условиях, т. е. имеет вид

1. Понтрягин Л.С. Асимптотическое поведение решений системы дифференциальных уравнений с малым параметром при высших производных, Известия Ак. паук СССР, серия метем, 21(1957), 605—626.

2. Мищенко Е.Ф., Понтрягии Л.С. Периодические решения систем дифференциальных уравнении, близкие к разрывным, Доклады Ак. наук СССР, 102, № 5 (1955), 889-891.

3. Мищенко Е.Ф., Асимптотическое вычисление периодических решении систем дифференциальных уравнении, содержащих малые параметры при производные. Известия Ак. наук СССР, серия матем., 21 (1957), 627—654.

4. Мищенко Е., Понтрягии Л.С. Вывод некоторых асимптотических оценок для решений дифференциальных уравнений с малым параметром пр» производных, Известия Ак. наук СССР, серия матем., 23(1959), 643—660.

5. Тихонов А. И-, Системы дифференциальных уравнений, содержащие малые параметры при производных, Матем. еборн., 31(73): 3 (1952), 574—586.

6. Боголюбов Н.И., Митропольский Ю.А., Асимптотические методы в теории нелинейных колебаний, Москва, 1955.

7. Митропольскнй Ю.А., Нестационарные процессы в нелинейных колебательных системах, Изд. АН УССР, 1955.

Асимптотические методы для линейных обыкновенных дифференциальных уравнений, Федорюк М.В.

Асимптотические методы для линейных обыкновенных дифференциальных уравнений, Федорюк М.В.

В книге содержатся асимптотические методы решения линейных обыкновенных дифференциальных уравнений. Рассмотрен ряд важных физических приложений к задачам квантовой механики, распространения волн и др. Для математиков, физиков, инженеров, а также для студентов и аспирантов университетов и инженерно-физических вузов.

ПРЕДИСЛОВИЕ.
В настоящей книге изложены основные результаты асимптотической теории обыкновенных линейных дифференциальных уравнений и систем, относящиеся к поведению решений с малыми параметрами при старших производных и к поведению решений при больших значениях аргумента. Литература по этим вопросам велика и разрозненна, но методы доказательств довольно однотипны, так что этот материал хорошо укладывается в монографию справочного типа. Мы ограничились только однородными уравнениями. Асимптотику решений неоднородного уравнения можно получить из асимптотики фундаментальной системы решений, применяя методы асимптотических оценок интегралов.

Оглавление.
Глава I.Аналитическая теория дифференциальных уравнений.
Глава II.Уравнения второго порядка на вещественной оси.
Глава III.Уравнения второго порядка в комплексной плоскости.
Глава IV.Уравнения второго порядка с точками поворота.
Глава V.Уравнения и системы n-го порядка.
Литература.

Бесплатно скачать электронную книгу в удобном формате, смотреть и читать:
Скачать книгу Асимптотические методы для линейных обыкновенных дифференциальных уравнений, Федорюк М.В. — fileskachat.com, быстрое и бесплатное скачивание.

Скачать pdf
Ниже можно купить эту книгу по лучшей цене со скидкой с доставкой по всей России. Купить эту книгу