.

Построение приближенного решения нелинейного уравнения методом Ван-дер-Поля

Язык: русский
Формат: реферат
Тип документа: Word Doc
89 710
Скачать документ

Министерство науки и образования Украины

Днепропетровский национальний университет

механико-математический факультет

кафедра дифференциальних уравнений

КУРСОВАЯ РАБОТА

“ПОСТРОЕНИЕ ПРИБЛИЖЕННОГО РЕШЕНИЯ НЕЛИНЕЙНОГО УРАВНЕНИЯ

МЕТОДОМ ВАН-ДЕР-ПОЛЯ”

Допускается к защите Исполнитель

Заведующий кафедрой ДР студентка гр. МЕ-98-2

Поляков М.В. Билан О.Ф.

«____» ________ 2002г. подпись ___________

подпись ___________

Научный руководитель

Профессор Остапенко В.А.

«____» ________ 2002г.

подпись ___________

Рецензент

Доцент Бойцун Л.Г.

«____» ________ 2002г.

подпись ___________

Днепропетровск

2002

Содержание

Содержание ……………………………………………………………………….….2

Реферат ……………………………………………………………………………….3

Annotation …………………………………………………………………………….4

Введение ………………………………………………………………………………5

Метод Ван-Дер-Поля …………………………………………………………7

Метод усреднения Ван-дер-Поля ……………………………………………7

Обоснование метода Ван-дер-Поля

Л. И. Мандельштамом и Н. Д. Папалекси ……………………………………….13

2. Решение уравнения ………………………………………………………………22

Выводы ………………………………………………………………………………29

Список использованной литературы ………………………………………………30

Реферат

Выпускная работа 30 стр., 5 источников.

Выпускная работа «Построение приближеного решения нелинейного уравнения
методом Ван-дер-Поля» посвящена эффективному способу решения нелинейных
задач теории колебаний с одной степенью свободы. Метод Ван-дер-Поля
обладает большой наглядностью и удобен для проведения расчетов.

Работа содержит теоретические выкладки по методу Ван-дер-Поля,
обоснование метода Мандельштамом и Папалекси и построение приближенного
решения уравнения:

.

Работа интересна для специалистов в области прикладной математики,
механики, физики и для студентов старших курсов.

Annotation.

The graduation paper “Approximated solution building of nonlinear
equation by Van-der-Pol’s method” is dedicated to very effective way of
nonlinear problems solution of oscillations theory with one degree of
freedom. Van-der-Pol’s method possesses the great visuality and is
comfortable for calculations.

The work contains theoretical part by Van-der-Pol’s method, the
validation of Mandelshtam and Papalexy method and approximated solution
building of the equation:

.

This work is very interesting for the experts in domain of applied
mathematics, mechanics, physics and for students of senior courses.

Введение.

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

Суть асимптотических методов заключается в том, что при их применении
достигается синтез простоты и точности за счет локализации: в
окрестности некоторого предельного состояния находится упрощенное
решение задачи, которое тем точнее, чем меньше эта окрестность.

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

Эффективность асимптотических методов признана всеми в самых разных
областях прикладной математики.

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

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

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

1. МЕТОД ВАН-ДЕР-ПОЛЯ

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

1.1. Метод усреднения Ван-дер-Поля.

В своих исследованиях Ван-дер-Поль рассматривал, главным образом,
уравнения с малым положительным параметром ? вида

(1)

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

Уравнение (1) условимся называть квазилинейным, а колебания, которые
оно описывает, — квазилинейными колебаниями. Функция f может быть
весьма общего вида и, в частности, даже разрывной.

Уравнение

(1.2)

называется порождающим. Оно описывает гармонические колебания. Общее
решение этого уравнения:

х=acos(?t+?),

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

Для получения приближенного решения уравнения (1) при достаточно малых
значениях параметра ? Ван-дер-Поль предложил особый прием, названный им
методом «медленно меняющихся» коэффициентов, аналогичный одному из
методов, применявшихся еще Лагранжем в небесной механике. Он представил
истинное решение уравнения (1) в виде функции, выражающей гармонические
колебания:

х=acos(?t+?)
(2)

с медленно меняющимися амплитудой а и фазой ?, которые должны находиться
из системы дифференциальных уравнений с разделяющимися переменными:

, (3)

составленными по определенному правилу. Уравнения (3), так называемые
«укороченные уравнения» Ван-дер-Поля, позволяют сравнительно просто
получить приближенное решение исходного уравнения (1). В частности,
задача отыскания периодического решения уравнения (1) сводится к
значительно более простой задаче нахождения состояния равновесия
системы, описываемой «укороченными уравнениями» (3).

Перейдем к составлению «укороченных уравнений» для рассматриваемого
уравнения (1), или эквивалентной ему системы двух уравнений первого
порядка

(4)

Прежде всего, заметим, что при ?=0 уравнение (1) превращается в
дифференциальное уравнение обычного гармонического осциллятора, и тогда
решение системы (4) имеет вид:

(5)

где а и ?— постоянные интегрирования.

Будем отыскивать решение уравнения (4) при достаточно малых значениях
параметра ? в виде выражений (5), но уже считая а и ? не постоянными, а
некоторыми функциями времени. Для этого будем рассматривать выражения
(5) не как решения уравнения (4) при ? = 0, а как формулы замены старых
переменных х и у на новые переменные а и ?.

.

Продифференцировав выражения (5) по t, подставим значения производных в
уравнениях (4). Принимая во внимание формулы (5), получаем систему
уравнений относительно производных, новых переменных а и ?:

(6)

, находим систему уравнений:

(7)

Система дифференциальных уравнений (7) эквивалентна рассматриваемой
исходной системе (4) или, что то же самое, уравнению (1) .

разделены.

Усредняя правые части системы (7) мы получим (8):

(8)

Таким образом, «укорочёнными уравнениями» для системы (7) являются
уравнения (3), где

(9)

разделены.

Интегрируя первое из уравнений этой системы, мы находим закон изменения
амплитуды. Очень часто в прикладных задачах бывает достаточно найти
только зависимость амплитуды от времени. В рассматриваемой теории для
этого достаточно найти решение уравнения первого порядка (в общем случае
нелинейного).

Определение фазы сводится к квадратурам. Наибольший интерес обычно
представляет не сама фаза, а скорость ее изменения в зависимости от
амплитуды. Ответ на этот вопрос дает непосредственно второе уравнение
системы (8).

(которые мы будем называть переменными Ван-дер-Поля) и к замене точных
уравнений (7) укороченной системой (8).

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

(10)

Заметим, что уравнение (10) совпадает с одним из тех уравнений, которое
мы получили бы, если бы рассматривали уравнение (1) как квазилинейное и
разыскивали периодические решения методом Пуанкаре.

Трансцендентное уравнение (10) может совсем не иметь действительных
решений. Это будет означать, что в системе стационарные колебания
невозможны. Уравнение (10) может иметь одно или несколько решений, в
случаях консервативных систем оно удовлетворяется тождественно. В самом
деле, в этом случае функция f зависит только от переменной x, поэтому
уравнение (10) примет вид:

(10а)

, то под знаком интеграла стоит полный дифференциал:

.

,

.

Обоснование метода Ван-дер-Поля

Л. И. Мандельштамом и Н. Д. Папалекси.

Рассмотрим систему стандартного вида

(s=1,2) (1)

Уравнение Ван-дер-Поля также можно привести к системе стандартного вида:

(2)

Сделаем замену

,

(3)

.

:

от t не зависят.

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

, (s=1,2).

Обе системы, приближенная и точная, решаются при начальных условиях

(4)

Для задач Коши (1) и (4), (3) и (4) справедлива следующая теорема:

и L>0:

(5)

(6)

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

Решение задач Коши (1) и (4), (3) и (4) существует и единственно.
Поэтому решение (1) и (4) будем искать методом приближений.

Обозначим

(*)

.

Пусть

(7)

:

.

поэтому

(7a)

находятся те значения t, для которых будет существовать решение (1) –
(4) и оно не выйдет за пределы области G. Это характеризуется так

Из теоремы Пикара следует, что при всех таких t приближенное выражение
сходится к решению задачи Коши:

,

— остаточный интервал.

С учетом возможности такого разбиения

, то последнее выражение перепишется в виде:

,

где с учетом (4)

.

Вычислим

То есть

(8)

Мы можем сказать, что в (8), все, что стоит под знаком суммы

Так как

,

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

Поэтому:

, (9)

(10)

— удовлетворяет условию Липшица, поэтому мы можем воспользоваться этим,
переходя к оценкам

(11)

(12)

, тогда:

(13)

Оценим

(14)

.

, тогда последнее неравенство

(последняя оценка получена с помощью неравенства (11)).

(15)

(16)

Можно увидеть следующую закономерность

(17)

Используя формулу (13), далее получим:

(18)

(19)

Обозначим через

Так как мы пользовались условиями Липшица, нужно убедиться, что
приближения не выходят из области G.

— по теореме Пикара это не выходит за пределы области G, то есть

В силу плотности числовой прямой

(20)

Проверим, вышло ли первое приближение за пределы области G. Пользуясь
оценками (19) и (20), имеем:

Возьмем

,

тогда

Аналогично проверяем второе приближение

Возьмем

, тогда

И если

,

если

, то получим:

(21)

будем выбирать из условия (21), то использование условия Липшица
законно.

с помощью (21) и

Решение уравнения

Рассмотрим уравнение

(1)

Данное уравнение второго порядка описывает колебательное движение. Здесь
? – некоторая действительная постоянная, а ? – малый параметр.

тогда получим систему

(2)

, согласно формулам

(3)

.

(4)

Подставим (4) в (2), учитывая (3).

(5)

,

тогда имеем:

(6)

Система (6) полностью эквивалентна уравнению (1). Соответствующая
системе (6) усредненная система имеет вид

(7)

имеют вид:

то есть

Таким образом имеем

или

(8)

Чтобы найти в явном виде закон изменения амплитуды в зависимости от
времени, необходимо решить первое уравнение системы (8):

:

.

Сделаем замену

,

:

,

,

;

.

Отсюда находим

(9а)

Колебания представятся следующим образом (находим выражение для
приближенного значения x в явном виде)

(9)

.

=0, амплитуда останется равной нулю для любого t, и, следовательно,
получим х=0, то есть тривиальное решение уравнения (1). Это тривиальное
решение, очевидно, соответствует статическому режиму, то есть отсутствию
колебаний в системе.

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

. Это решение соответствует стационарному (установившемуся)
динамическому режиму:

(10)

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

, приводят к уравнению

=0

.

;

соответствует устойчивому предельному циклу.

. Все эти рассуждения следуют из теоремы Мандельштама и Папалекси.

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

, (s=1,2) .

,

(s=1,2)

:

непрерывны.

для системы (2) имеют вид:

.

непрерывно дифференцируемы по t, а следовательно удовлетворяют условию
Липшица.

.

определяется уравнением (9а)).

Выводы

В рамках теории Ван-дер-Поля нельзя уточнить полученные решения.

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

Список использованной литературы

Ю.А. Митропольский Метод усреднения в нелинейной механике, «Наукова
думка» Киев — 1971г.

Н.Н. Моисеев Асимптотические методы нелинейной механики, М.: Наука,
1981г.

А. Найфэ Методы возмущений. Издательство «Мир», Москва—1976г.

А. Найфэ Введение в методы возмущений. «Мир», 1984г.

Андронов А.А., Витт А.А., Хайкин С.Э. Теория колебаний. Физматгиз,
М.,1959г.

PAGE

PAGE 3

Нашли опечатку? Выделите и нажмите CTRL+Enter

Похожие документы
Обсуждение

Ответить

Курсовые, Дипломы, Рефераты на заказ в кратчайшие сроки
Заказать реферат!
UkrReferat.com. Всі права захищені. 2000-2020