BSc: DifferentialAndPartialDifferentialEquationsADV

Уравнения в частных производных

Квалификация выпускника: бакалавр

Направление подготовки: 09.03.01 - “Информатика и вычислительная техника”

Направленность (профиль) образовательной программы: Математические основы ИИ

Программу разработал(а): В.А.Гордин

1. Краткая характеристика дисциплины

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

2. Перечень планируемых результатов обучения

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

Пререквизиты (Предварительные знания у слушателей)

Матанализ, линейная алгебра (в объеме книги И.М.Гельфанда «Лекции по линейной алгебре» кроме главы о тензорах), начала функционального анализа, курс обыкновенных дифференциальных уравнений, методы вычислений..

Общая характеристика результата обучения по дисциплине

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

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

3. Структура и содержание дисциплины

№ п/п	Наименование раздела дисциплины	Содержание дисциплины по темам
Тема 1.	Введение. Вывод уравнений неразрывности, диффузии, теплопроводности. Метод характеристик для урчп 1 порядка.	Понятие плотности сплошной среды. Перенос плотности потоком с заданной скоростью. Изменение плотности импульса под действием распределенных сил. Уравнения движения. Уравнение Эйлера – Хопфа. Характеристики. Градиентная катастрофа. Законы Фика и Фурье. Уравнения диффузии и теплопроводности. Начальные и граничные условия. Интенсивность движения купеческих караванов. Уравнение Блэка – Шоулза – Мертона для определения справедливой цены европейских опционов..
Тема 2.	Задача Коши. Существование, единственность, корректность. Теорема Коши – Ковалевской, которая не гарантирует корректность.	Уравнения, разрешенные (и не разрешенные) относительно старшей производной по времени. Уравнения и системы типа Коши – Ковалевской. Сколько нужно начальных условий при решении во всем пространстве. Аналитические начальные условия. Теорема Коши – Ковалевской для аналитических решений задачи Коши «в малом». Понятие корректности для задачи Коши. Преобразование Фурье для урчп с постоянными коэффициентами. Некорректность задачи Коши для уравнения теплопроводности в обратную сторону по времени. Условия на бесконечности по переменной x. Теорема Тихонова. Пример Адамара – задача Коши для уравнения Лапласа. Формула Пуассона для решения краевой задачи.
Тема 3.	Краевые и смешанные задачи. Метод разделения переменных для уравнения теплопроводности.	. .
Тема 4.	Уравнение струны. Формула Даламбера. Метод разделения переменных для смешанной краевой задачи.	Вывод уравнения струны и граничных условий. Задача на всей прямой. Многомерное уравнение лакуны и диффузия волн – зависит от четности размерности пространства. Задача на отрезке. Разделение переменных. Колебания амплитуд. Законы сохранения. Отражение и неотражение от границ.
Тема 5.	Свойства собственных функций.	Ряд Фурье. Проекция на первые собственные функции. Неравенство Бесселя и равенство Парсеваля. Сходимость разложения функции в метриках $L^{2}$ и в $C$. Скорость сходимости и убывания коэффициентов. Явление Уилбрахама - Гиббса. Ортогонализация многочленов в Failed to parse (syntax error): {\displaystyle L^2\(-1, 1\)} . Многочлены Лежандра. Формула Родрига. Задача Штурма - Лиувилля. Приведение к самосопряженному виду. Собственные числа и функции. Асимптотика собственных чисел при больших номерах. Простота спектра для оператора 2 порядка. Чередование нулей собственных функций. Разделение переменных для параболической и гиперболической задач на отрезке. Разностная аппроксимация задачи Штурма - Лиувилля. Компактные конечно-разностные схемы. Уравнение Бесселя нулевого порядка. Ограниченное в нуле решение и его разложение в ряд Тейлора. Отрицательность спектра при условии ограниченности на бесконечности. Оператор Лапласа на плоскости в полярных координатах. Функции Бесселя – собственные для задачи Дирихле в круге. Асимптотика ограниченного решения в начале координат. Уравнение колебаний стержня. Варианты граничных условий. Принцип наименьшего действия. Амплитуда вероятностей. Уравнение Шрёдингера.
Тема 6.	Простейшие экологические модели.	Логистическое уравнение. Устойчивая и неустойчивая стационарная точки. Возможные границы отлова. Жесткая и мягкая модели. Опасность оптимизации в жесткой модели. Гибкие планы отлова. Модель Лотки – Вольтера. Стационарные точки и исследование их устойчивости. Первый интеграл системы (два способа построения). Сравнение теории с экспериментальными данными. Ограничения модели. Задача о двух видах, конкурирующих за общий ресурс.
Тема 7.	Конечно-разностные уравнения и системы.	Последовательность Фибоначчи. Конечно-разностные уравнения и системы. Пространство решений линейного конечно-разностного уравнения n-го порядка n-мерно. Общее решение для уравнения с постоянными коэффициентами. Случаи простых и кратных корней характеристического уравнения. Матрица Лесли и предельное распределение популяции по возрастам. Игра с конечной суммой (блуждание частицы). Марковские цепи. Пример нелинейного конечно-разностного уравнения. Метод Герона и метод Ньютона. Метод Ньютона и сверхсжатие для некратных корней.
Тема 8.	Устойчивость и неустойчивость стационарных точек систем.	Устойчивость положения равновесия при t→+∞ для дифференциальных и разностных систем с постоянными коэффициентами. Теория Ляпунова – исследование устойчивости стационарных точек нелинейных систем (без док.). Примеры, когда спектральный метод бессилен. Метод функции Ляпунова.
Тема 9.	Введение в разностные схемы для решения задачи Коши.	Разностные схемы для решения задачи Коши. Схема Эйлера, Эйлера с пересчетом, схема центральных разностей. Схемы Рунге – Кутты. Метод экстраполяции Ричардсона для повышения точности схемы.
Тема 10.	Семейства траекторий и решение типа бегущей волны.	Уравнение в вариациях. Гомотопия. Метод стрельбы (пристрелки) для решения краевой задачи. Изменение фазового объема в окрестности траектории. Дивергенция векторного поля. Теорема Лиувилля, спектр матрицы системы и разбегание траекторий. Вывод уравнения неразрывности. Характеристики. Уравнение переноса и решение типа бегущей волны. Неоднородное уравнение переноса и изменение решения вдоль характеристики. Применение фундаментальной системы решений к интегрированию неоднородных систем – метод Лагранжа вариации постоянных.
Тема 11.	Специальные типы систем ОДУ.	Теория Флоке. Матрица монодромии и мультипликаторы. Анализ устойчивости для решений уравнений с периодическими коэффициентами. Гамильтоновы системы. Сохранение гамильтониана. Сохранение фазового объема. Примеры. Принцип наименьшего действия.
Тема 12.	Сингулярные возмущения и релаксационные колебания.	Системы вида $\varepsilon d_{t}{\overrightarrow {X}}^{\varepsilon }={\vec {g}}({\vec {X}}^{\varepsilon },{\vec {Y}}^{\varepsilon },\varepsilon )$, $d_{t}{\overrightarrow {Y}}^{\varepsilon }={\vec {f}}({\vec {X}}^{\varepsilon },{\vec {Y}}^{\varepsilon },\varepsilon )$, где $0<\varepsilon \ll 1$.Поверхность вырождения: устойчивые и неустойчивые участки. Фазовый портрет. Чередование быстрых и медленных движений. Периодические режимы. Уравнение Ван дер Поля.
Тема 13.	Конечно-разностные уравнения и системы.	Формула Тейлора. Интерполяционный полином Лагранжа. Численная производная двухточечная и четырехточечная формулы для численной производной. Численное интегрирование. Оценки погрешностей численной производной и интеграла. Определение неподвижной точки. Теорема о неподвижной точке. Доказательство, примеры. Сжимающее отображение. Определение. Теорема Банаха о неподвижной точке. Доказательство, примеры. Метод Ньютона для решения нелинейных уравнений. Условия сходимости метода. Примеры.
Тема 14.	Системы с несколькими первыми интегралами.	Центральная сила. Примеры. Секториальная скорость и доказательство второго закона Кеплера. Эффективная потенциальная энергия. Сохранение энергии радиального движения. Определение фазы. Апоцентр и перицентр. Условие периодичности орбиты. Первый и третий законы Кеплера, Бертрана и Кенига – без док. Сохранение импульса для замкнутой системы. Движение центра масс. Сохранение момента импульса замкнутой системы. Случай сохранения проекции момента импульса в некоторых незамкнутых системах. Сохранение энергии в задаче N-тел. Полная интегрируемость в задаче 2 тел.
Тема 15.	Обобщенные функции.	Основные функции (варианты пространства). Пример Коши и разбиение единицы.Топология в пространстве основных функций. Неметризуемость этой топологии.Умножение обобщенной функции на гладкую.Дельта-функция.Обобщенные функции типа функции и другие.Носитель и сингулярный носитель. Примеры.Слабое дифференцирование обобщенных функций. Примеры.Решение ОДУ в пространстве обобщенных функций. Сравнение с классическим решением. Преобразование Фурье обобщенных функций.
Тема 16.	Простейшие задачи вариационного исчисления.	Задача согласования информации о координате и скорости. Задача Дидоны. Цепная линия. Катеноид. Брахистохрона. Геодезические. Принцип Ферма и рефракция. Непрерывные и гладкие функционалы. Первая вариация. Необходимое условие экстремума гладкого функционала – равенство нулю первой вариации. Примеры недостаточности этого условия. Вывод уравнения Эйлера. Граничные условия трансверсальности. Условный экстремум и множители Лагранжа. Вторая вариация. Задачи со старшими производными, несколькими функциями и в частных производных. Принцип наименьшего действия. Примеры.

4. Методические и оценочные материалы

Формы контроля: Контроль знаний студентов включает формы текущего и итогового контроля. Текущий контроль осуществляется в виде коротких контрольных работ в начале многих занятий и контрольной работы в середине триместра. Кроме того будет выдано несколько домашних работ, которые должны выполняться студентом в течение одной недели. Если она сделана в течение второй недели, оценка за нее делится пополам. После второй недели 10-балльная оценка за сданную работу - нулевая. Итоговый контроль осуществляется в виде двух экзаменов, один из которых теоретический, а второй – обсуждение письменной работы и решение задач на компьютере. Веса обоих экзаменов равные. Итоговая оценка по 100-балльной шкале формируется по формуле ${\text{О}}{\text{итог}}=0.3\times {\text{О}}{\text{предв}}+0.2\times {\text{О}}{\text{контр}}+0.5\times {\text{О}}{\text{экз}}$. Округление происходит только для итоговой отметки. Кроме того на протяжении курса студентам выдаются домашние задания, где решение требует комбинированного подхода: аналитические соображения + численная компьютерная реализация. Задачи, как правило, содержат индивидуальный параметр Y или параметры.

Задания для практических занятий:

1.Доказать, что в некоторой точке на интервале (-1;1) производная функции $f(x)=(1-x^{2})^{100}$ равна 1.

2.Как повлияет на пропускную способность трубы со стационарным ламинарным течением жидкости а) удвоение перепада давления на концах трубы; б) удвоение диаметра трубы?

3.Приведите определение метрического пространства, линейного нормированного пространства, евклидова пространства. Каковы соотношения между этими понятиями?

4.Лестница прислоняется к трубе, как показано пунктиром на рисунке. Радиус трубы равен 1. Нижний конец лестницы находится в точке x=1+k/10. Определить координаты точки касания лестницы и трубы. Здесь k – номер студента в списке.

5.Из плоскости вырезан круг радиуса 1 с центром в начале координат. Точка A=<1+Y/1000,0> на этой плоскости соединена натянутой нитью с точкой B=<0,1+Y/100>. Определить длину нити. При каком Y нить будет иметь форму отрезка?

6.Методом разделения переменных решить уравнение${\frac {dx}{dt}}=Ax^{\alpha }$ при $A=\pm Y$,α=0,1/2,1,2. Здесь Y – номер студента в списке. Под словом решить понимается описание общего решения при произвольных начальных данных. Для каждого из вариантов построить несколько траекторий.

7.Известно, что в момент $t_{1}=7$ первая популяция имела численность 2Y, а вторая Y. Первая популяция возрастает со временем согласно уравнению ${\frac {dx}{dt}}=Yx$, а вторая ${\frac {dz}{dt}}=2Yz$. В какой момент численности обеих популяций будут (или были) равными?

8.Построить график решения уравнения ${\frac {dx}{dt}}=Yx^{2}$ с начальным условием Y при $t_{1}=0$. В какие моменты решение будет в два раза больше и в два раза меньше начального значения?

9.Рассмотрим три дифференциальных уравнения первого порядка вида: ${\frac {dz}{dt}}=a|z|^{b}$, где a=2Y,_ b=1-Y/100,_ 1,_ 1+Y/100 с начальным условием при t=0: z(0)=Y. Методом разделения переменных найти решение в максимально возможных пределах в обе стороны по t. Существенен ли знак модуля в уравнении ? Построить графики решений.

10.Для уравнения фон Берталанфи с α=1,_ β=Y определить время удвоения объема при различных начальных данных.

11.Для уравнения Гомперца финальная масса вдвое больше начальной. Временной масштаб λ=Y. Определить константу r.

12.Привести к диагональному виду оператор с матрицей A=

13.Решить систему дифференциальных уравнений ${\frac {d{\vec {X}}}{dt}}=A{\vec {X}}$ с упомянутой выше матрицей А и начальным условием X ⃗(0)=<0,0,0,1>. Построить графики компонент решения на отрезке [0,3].

14.То же задание для системы с матрицей B=A-2E, E – единичная матрица.

15.Ответить еще раз на вопросы 8-10, используя метод Рунге - Кутты. Сравнить полученные графики с аналитическими. Исследовать зависимость погрешности численного решения от шага схемы Р-К и времени интегрирования.

16.В шар объемом V вписан цилиндр. При каком радиусе объем цилиндра максимален? При каком радиусе максимальна площадь его поверхности ?

17.Чтобы удержать груз на канате, перекинутом через балку, нужна сила ≥10кг, а чтобы начать подтягивать на свою сторону ≥10+ Y кг. Определить вес груза. Определить коэффициент трения, если угол обхвата $\omega =10Y$ градусов.

18.В следующих задачах начальные данные <x,y> для системы двух дифференциальных уравнений первого порядка пробегают единичную окружность: $x^{2}(0)+y^{2}(0)=1$. Требуется описать (и нарисовать кривые – геометрическое место точек) множество решений в следующие моменты времени t=-1, 1, 2. Для ориентировки: если А – нулевая матрица, то все три искомые кривые совпадают с единичной окружностью, а если А – единичная матрица, то это окружности радиуса e^t. Для каждой задачи указать, имеется ли у системы первый интеграл?

19.Для уравнений химической кинетики при $\nu =\mu =1,\quad k1=1,\quad k2=2$ на фазовой области описать множество начальных условий, для которых реакция полностью заканчивается за время Y.

20.Пружинный маятник с трением описывается уравнением ${\ddot {x}}+k{\dot {x}}+\omega ^{2}x=0$. Пусть $\omega =1,\quad k=2\pm {\frac {Y}{100}}$. Построить графики решения при нескольких различных начальных данных. Построить фазовые портреты.

21.Длина физического маятника без трения Y см. Ускорение свободного падения g=9,8 м/сек^2. Определить период малых колебаний. Численно определить амплитуду колебаний при которых период вдвое и втрое больше. Для уравнения идеального маятника ${\frac {d^{2}x}{dt^{2}}}+Y\sin {x}=0$ энергетическим методом построить траектории. Интеграл для периода колебаний вычислить методом Симпсона. Построить график зависимости периода колебаний от его амплитуды.

22.Построить методом Рунге – Кутты траектории (несколько - с разными начальными данными) для уравнения маятника с трением ${\frac {d^{2}x}{dt^{2}}}+{\frac {dx}{dt}}+Y\sin {x}=0$.

23.Тот же вопрос для ${\frac {d^{2}x}{dt^{2}}}+11{\frac {dx}{dt}}+Y\sin {x}=0$

24.Для уравнения ${\frac {d^{2}x}{dt^{2}}}+{\frac {dx}{dt}}+Yx=\sin {\omega t}$ подобрать частоту $\omega$ так, чтобы амплитуда вынужденных колебаний была максимальна. Построить графики и траектории для решения (с нулевыми начальными условиями) для этой частоты, а также для ${\frac {\omega }{2}}$ и $2\omega$.

25.Построить траектории для уравнения ${\frac {d^{2}x}{dt^{2}}}+{\frac {dx}{dt}}+Yx^{3}=\sin {\omega t}$

26.Пусть функция $\Psi$ - ступенчатая, попеременно на отрезках равной длины принимающая значения 1 и -1, $\Psi (s)={\text{sign}}\sin(s)$.

27.Для уравнения ${\frac {d^{2}x}{dt^{2}}}+{\frac {dx}{dt}}+Yx=\Psi (\omega t)$подобрать частоту $\omega$ так, чтобы амплитуда вынужденных колебаний была максимальна. Построить графики и траектории для решения (с нулевыми начальными условиями) для этой частоты, а также для ${\frac {\omega }{2}}$ и $2\omega$.

28.Построить траектории для уравнения ${\frac {d^{2}x}{dt^{2}}}+{\frac {dx}{dt}}+Yx^{3}=\Psi (\omega t)$.

29.Определить численно зависимость периода и сдвига фаз между компонентами решения от амплитуды периодических решений системы Лотки – Вольтера при $\alpha =\beta =\gamma =1,\quad \delta =Y$. Построить графики решений для нескольких вариантов начальных данных.

30.Для уравнений химической кинетики построить графики зависимости решения от времени при $\nu =1,\quad \mu =2,\quad k1=10,\quad k2=Y,\quad n1(0)=5,\quad n2(0)=7$.

31.Для многочлена $P(x)=x^{3}+\alpha x^{2}+Yx$ определить значения $\alpha _{*}$, при которых имеется вырожденная стационарная точка. Как различаются многочлены со значениями $\alpha$ по разные стороны от $\alpha _{*}$? Построить графики для примеров.

32.Параметрическим называется резонанс вследствие изменения коэффициентов уравнения. Для уравнения второго порядка $x''+x'+Y(1+a\sin {\omega t})x=0$ численным экспериментом определить, при каких значениях параметров $a,\omega$ теряется устойчивость состояния покоя ? Указать границы областей параметров, где ПР наблюдается. Для нескольких вариантов параметров построить графики решения.

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

34.Для системы уравнений Лотки - Вольтерры с параметрами из задачи 31 определить погрешность (изменение первого интеграла в конечный момент времени по отношению к начальному) в зависимости от периода и шага разностной схемы Рунге - Кутты 4-го порядка с постоянным шагом. Время интегрирования – 100 периодов.

35.Использовать экстраполяционный метод Ричардсона для уменьшения этой погрешности. В каком диапазоне шагов метод неэффективен и почему?

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

37.Пусть $f(x,y)=\sin {x}+\cos {(x+Yy)}$. С помощью леммы Морса исследовать поведение функции в окрестности стационарных (критических) точек. Нарисовать изолинии f.

38.Пусть$f(x,y,z)=\sin {x}+\cos {(x+Yy)}+\sin {(x+Yy+Y^{2}z)}$ Исследовать поведение функции в окрестности стационарных точек.

39.Исследовать системы ${\dot {x}}=y\pm Yy^{2},\quad {\dot {y}}=-x\pm x^{2}$ методом разделения переменных. Устойчивы ли стационарные точки этих систем? Сопоставить с исследованием устойчивости методом Ляпунова.

40.Для уравнения пружинного маятника с трением ${\ddot {x}}+{\text{sign}}({\dot {x}})+Yx=0$ с начальным условием <Y,Y> определить число колебаний до остановки.

41. Методом Ньютона исследовать уравнение $\sin {z}=Yz,\quad z\in \mathbb {C}$

42.Рассмотрим линейный дифференциальный оператор $(x^{2}+3x-Y){\frac {d^{2}}{dx^{2}}}+(2x-7){\frac {d}{dx}}$. Докажите, что он оставляет инвариантным подпространство многочленов степени не выше 10. Вычислить матрицу этого оператора в базисе, составленном из этих 11 мономов. Вычислить спектр этого оператора. То же для многочленов степени не выше 12.

43.Для дифф. уравнения Бесселя степени 2: $r^{-1}{\frac {d}{dr}}r{\frac {d}{dr}}u+(Y-4r^{-2})u=0$ построить решение, ограниченное в нуле. При малых r строить разложением в ряд Ньютона, а потом при некотором $r=\varepsilon$ использовать полученные $u(\varepsilon ),\quad u'(\varepsilon )$в качестве начальных данных для метода Рунге - Кутты, каковым интегрировать до 10/\sqrt Y. Оценить зависимость погрешности от числа членов ряда Ньютона, выбора $\varepsilon$ и шага схемы. Применить метод Ричардсона для повышения точности.

44.Для сетки {-Y,0,1,2,3} построить многочлен степени 5, который во всех узлах обращается в нуль, а первая производная которого в левой точке равна 1. То же для правой точки. Построить графики.

45.На единичной окружности $x\in [0,2\pi )$задана равномерная сетка из Y+10 точек. Значения сеточной функции равны значениям функции sin(x). Вычислить в точках сетки первую и вторую производные по компактной схеме и сравнить с истинным результатом. Построить графики.

46.На маятник сбоку дует ветер. Поэтому уравнение для его колебаний принимает вид: ${\ddot {x}}+A\sin {x}+B\cos {x}=0$. Нужно a) Объяснить физический смысл В; b) Построить фазовый портрет и проинтегрировать уравнение при А=1, В=Y. c)При этих же значениях параметров определить зависимость периода от амплитуды и сравнить со случаем В=0.

47.Для уравнения маятника с трением ${\frac {d^{2}x}{dt^{2}}}+{\frac {dx}{dt}}+Yx=0$ рассмотрим решения с начальными данными <1,0> и <0,1>. На фазовой плоскости построить соответствующие параллелограммы при t=-1, 1, 3. Построить график зависимости вронскиана от времени.

48. То же для начальных данных в круге $(x-Y)^{2}+{\dot {x}}^{2}<{\frac {Y}{2}}$. Приложить распечатку программы.

49.Для уравнения ${\frac {d^{2}x}{dt^{2}}}+{\frac {dx}{dt}}+Y(x+\varepsilon \sin {x})=0$, зависящего от параметра$\varepsilon$ рассмотрим решения с начальными данными <1,0> и <0,1>. На фазовой плоскости построить для |t|< 1 решения при $\varepsilon =0,0.01,0.02$. Построить разности решений для 0,01 и 0, для 0,02 и 0. Построить решения для уравнения в вариациях, полагая $\varepsilon \rightarrow 0$. Для тех же значений $\varepsilon ==0,_{0},01,_{0},02$. вычислить его решение. Сравнить оба метода. Оценить скорость нарастания погрешности метода уравнения в вариациях со временем.

50. Для уравнения второго порядка $x''+x'+Y(1+a\sin {\omega t})x=0$ численным экспериментом определить, в зависимости от значений параметров $a,\quad \omega$ след и определитель матрицы монодромии. Определите мультипликатор с наибольшим модулем. Постройте кривые, на которых он равен 1, - они являются границами устойчивости нулевого решения (и параметрического резонанса). Сопоставить с результатами задачи, полученными в 53.

51.Для уравнения ${\ddot {x}}+Yx=f(t)$ методом вариации постоянных получить общее решение.

52.Шарик с нулевой начальной скоростью под действием силы тяжести без трения движется под действием силы тяжести по желобу, имеющему форму параболы, причем z(0)=1, z(Y)=0. Требуется с помощью численных экспериментов определить, какая из парабол обеспечивает наименьшее время для достижения конечной точки.

53. Рассмотрим гамильтонову систему с двумя степенями свободы с гамильтонианом

$H(x_{1},x_{2},p_{1},p_{2})={\frac {1}{2}}c({\vec {x}})|p|^{2}-{\frac {1}{2}}c^{-1}({\vec {x}})$ Выписать систему Гамильтона. Траектории этой гамильтоновой системы (бихарактеристики), точнее их проекции на плоскость <x_1,x_2>, описывают движение лучей в среде с переменной скоростью с. Докажите, что при с=const лучи – прямые. Если рассматривается точечный источник лучей (скажем, из начала координат), то начальные данные образуют двумерное подпространство <0,0,p_1,p_2>. Предположим, что среда «слоистая»: $c=f(x2)={\begin{cases}1&{\text{if }}|x2|>1\\1+Y\sin {\pi x_{2}}&{\text{if }}|x_{2}|<1\end{cases}}$. Вычислить геометрию лучей. Определить критический угол $\alpha =\left|\arctan {\frac {p_{1}(0)}{p_{2}(0)}}\right|$, при котором лучи не покидают волновод $|x_{2}|<1$.

Контрольные вопросы для подготовки к промежуточной аттестации:

Для оценки качества освоения дисциплины можно использовать задачи, приведенные в задачнике Филиппов А.Ф. Сборник задач по дифференциальным уравнениям. 2008.

Несколько тысяч задач имеется в тексте книг:

1. В.И.Арнольд: Обыкновенные дифференциальные уравнения. М., ``Наука, 1984, 2002.

2. В.А.Гордин: Дифференциальные уравнения. Какие явления они описывают и как их решать. М. Издательский Дом ВШЭ. 2016.

3. Гордин В.А. Математика, компьютер, прогноз погоды и другие сценарии математической физики. М., ФИЗМАТЛИТ. 2010, 2013.

4. Гордин В.А. Прикладная математика. Искусство и ремесло вычислений. Готовится к изданию в М. Издательский Дом ВШЭ. 2024.

Перечень учебно-методического обеспечения дисциплины

Основная литература

Тема 1.

1. В.И.Арнольд: Обыкновенные дифференциальные уравнения. М., ``Наука, 1984, 2002.
   
2. Гельфонд А.О. Исчисление конечных разностей. М., Наука, 1967, УРСС 2012.
   
3. В.А.Гордин: Дифференциальные уравнения. Какие явления они описывают и как их решать. Изд. Дом ВШЭ, М., 2016.
   
4. Петровский И.Г. Лекции по теории обыкновенных дифференциальных уравнений. М., Наука - Физматлит, 1970, Изд. МГУ, 1984. М.: ФИЗМАТЛИТ, 2009, ЛитРес, 2022.

Тема 2.

1. В.И.Арнольд: Обыкновенные дифференциальные уравнения. М., ``Наука, 1984, 2002.
   
2. В.А.Гордин: Дифференциальные уравнения. Какие явления они описывают и как их решать. Изд. Дом ВШЭ, М., 2016.
   
3. Петровский И.Г. Лекции по теории обыкновенных дифференциальных уравнений. М., Наука - Физматлит, 1970, Изд. МГУ, 1984. М.: ФИЗМАТЛИТ, 2009, ЛитРес, 2022.
   

Тема 3.

1. В.А.Гордин: Дифференциальные уравнения. Какие явления они описывают и как их решать. Изд. Дом ВШЭ, М., 2016.
   
2. Петровский И.Г. Лекции по теории обыкновенных дифференциальных уравнений. М., Наука - Физматлит, 1970, Изд. МГУ, 1984. М.: ФИЗМАТЛИТ, 2009, ЛитРес, 2022.
   
3. Г.М.Фихтенгольц: Курс дифференциального и интегрального исчисления. т. 1. М., Физматгиз, 1963, Лань 2023.
   
4. К.Чен, П.Джиблин, А.Ирвинг: MATLAB в математических исследованиях. М., ``Мир, 2001.
   

Тема 4.

1. В.И.Арнольд: Обыкновенные дифференциальные уравнения. М., ``Наука, 1984, 2002.
   
2. В.А.Гордин: Дифференциальные уравнения. Какие явления они описывают и как их решать. Изд. Дом ВШЭ, М., 2016.
   
3. Оболенский А. Ю. Лекции по качественной теории дифференциальных уравнений — М.; Ижевск, 2006.
   
4. Петровский И.Г. Лекции по теории обыкновенных дифференциальных уравнений. М., Наука - Физматлит, 1970, Изд. МГУ, 1984. М.: ФИЗМАТЛИТ, 2009, ЛитРес, 2022.
   

Тема 5.

1. В.И.Арнольд: Обыкновенные дифференциальные уравнения. М., ``Наука, 1984, 2002.
   
2. В.А.Гордин: Дифференциальные уравнения. Какие явления они описывают и как их решать. Изд. Дом ВШЭ, М., 2016.
   
3. В.А.Гордин: Математика, компьютер, прогноз погоды и другие сценарии математической физики. М., ФИЗМАТЛИТ, 2010, 2013.
   
4. Петровский И.Г. Лекции по теории обыкновенных дифференциальных уравнений. М., Наука - Физматлит, 1970, Изд. МГУ, 1984. М.: ФИЗМАТЛИТ, 2009, ЛитРес, 2022.
   
5. Н.М.Матвеев. Методы интегрирования обыкновенных дифференциальных уравнений. М.: Высшая школа, 1967.
   

Тема 6.

1. В.И. Арнольд, «Жесткие и мягкие математические модели», М., МНЦМО, 2000.
   
2. В.Вольтерра: Математическая теория борьбы за существование.} М., ``Наука", 1976.
   
3. В.А.Гордин: Дифференциальные уравнения. Какие явления они описывают и как их решать. Изд. Дом ВШЭ, М., 2016.
   
4. А.А.Самарский, А.П.Михайлов: Математическое моделирование. Физматлит, М., 2002.
   
5. Ю.М.Свирежев, Д.О.Логофет: Устойчивость биологических сообществ. ``Наука, М., 1978.
   

Тема 7.

1. Гельфонд А.О. Исчисление конечных разностей. М., Наука, 1967, УРСС 2012.
   
2. В.А.Гордин: Как это посчитать? М., МЦНМО, 2005.
   
3. В.А.Гордин: Дифференциальные уравнения. Какие явления они описывают и как их решать. Изд. Дом ВШЭ, М., 2016.
   

Тема 8.

1. Гельфонд А.О. Исчисление конечных разностей. М., Наука, 1967, , УРСС 2012.
   
2. В.А.Гордин: Дифференциальные уравнения. Какие явления они описывают и как их решать. Изд. Дом ВШЭ, М., 2016.
   
3. А.А.Андронов, А.А.Витт, С.Э.Хайкин: Теория колебаний, М., Физматгиз, 1959, ``Наука, 1981.
   
4. Зайцев В.Ф., Полянин A.Д. Справочник по нелинейным обыкновенным дифференциальным уравнениям. М., Факториал, 1997.
   
5. Камке Э. Справочник по обыкновенным дифференциальным уравнениям. М., Наука, 1976; Спб. ``Лань, 2003.
   

Тема 9.

1. Н.С.Бахвалов, Н.П.Жидков, Г.М.Кобельков: Численные методы. ``Наука, М., 1987.
   
2. В.А.Гордин: Дифференциальные уравнения. Какие явления они описывают и как их решать. Изд. Дом ВШЭ, М., 2016.
   
3. Рябенький В.С. Введение в вычислительную математику. М., Наука, 1994, Физматлит, 2008.
   

Тема 10.

1. Н.М.Матвеев. Методы интегрирования обыкновенных дифференциальных уравнений. М.: Высшая школа, 1967.
   
2. В.А.Гордин: Дифференциальные уравнения. Какие явления они описывают и как их решать. Изд. Дом ВШЭ, М., 2016.
   
3. В.А.Гордин: Математика, компьютер, прогноз погоды и другие сценарии математической физики. М., ФИЗМАТЛИТ 2010, 2013.
   

Тема 11.

1. В.И.Арнольд: Обыкновенные дифференциальные уравнения. М., ``Наука, 1984, 2002.
   
2. В.А.Гордин: Дифференциальные уравнения. Какие явления они описывают и как их решать. Изд. Дом ВШЭ, М., 2016.
   

Тема 12.

1. Айнс Э.Л. Обыкновенные дифференциальные уравнения. Харьков, Государственное научно-техническое издательство Украины, 1939, М. «Факториал Пресс», 2005.
   
2. Федорюк М.В. Асимптотические методы для линейных обыкновенных дифференциальных уравнений. М., Наука, 1983.
   

Тема 13.

1. Е.Ф.Мищенко, Н.Х.Розов Дифференциальные уравнения с малым параметром и релаксационные колебания. М., Наука, 1975.
   

Тема 14.

1. В.И.Арнольд. Обыкновенные дифференциальные уравнения. М., ``Наука, 1984, 2002.
   
2. В.А.Гордин. Дифференциальные уравнения. Какие явления они описывают и как их решать. Изд. Дом ВШЭ, М., 2016.
   
3. Л.Д.Ландау, Е.М.Лифшиц. Механика, ``Наука, М., 1965, 1988.
   

Тема 15.

1. М.С.Агранович. Обобщенные функции. М.: МЦНМО, 2008.
   
2. В.А.Гордин. Математика, компьютер, прогноз погоды и другие сценарии математической физики. М.: Физматлит, 2010; 2013.
   
3. Прикладная математика. Искусство и ремесло вычислений. Готовится к изданию в М. Издательский Дом ВШЭ. 2024.
   
4. Л.Шварц. Математические методы для физических наук. М.: Мир, 1965.
   
5. Г.Е.Шилов. Математический анализ. Второй спецкурс. М.: Наука, 1965.
   

Тема 16.

1. Буслаев В.С. Вариационное исчисление. Л.: Изд. ЛГУ, 1980.
   
2. И.М.Гельфанд, С.В.Фомин Лекции по вариационному исчислению. М.: Физматгиз, 1961
   
3. В.А.Гордин. Математика, компьютер, прогноз погоды и другие сценарии математической физики. М.: Физматлит, 2010; 2013.
   

Дополнительная литература

Тема 1.

1. Пайтген Х.-О., Рихтер П.Х. Красота фракталов. Образы комплексных динамических систем. М., Мир, 1993.
   

Тема 2.

1. Абрамовиц М., Стиган И. Справочник по специальным функциям. М., Наука, 1979.
   
2. Шилов Г.Е. Математический анализ. Специальный курс. М., Физматгиз, 1960.
   

Тема 3.

1. Камке Э. Справочник по обыкновенным дифференциальным уравнениям. М., Наука, 1976; Спб. ``Лань, 2003.
   
2. Зайцев В.Ф., Полянин A.Д. Справочник по нелинейным обыкновенным дифференциальным уравнениям. М., Факториал, 1997.
   

Тема 5.

1. Айнс Э.Л. Обыкновенные дифференциальные уравнения. Харьков, Государственное научно-техническое издательство Украины, 1939, М. «Факториал Пресс», 2005.
   
2. А.А.Андронов, А.А.Витт, С.Э.Хайкин: Теория колебаний, М., Физматгиз, 1959, ``Наука, 1981.
   
3. В.И.Арнольд, Ю.С.Ильяшенко: Обыкновенные дифференциальные уравнения. В сб. Динамические системы. т.1. Современные проблемы математики. Фундаментальные направления, М., ВИНИТИ, 1985.
   
4. Зайцев В.Ф., Полянин A.Д. Справочник по нелинейным обыкновенным дифференциальным уравнениям. М., Факториал, 1997.
   
5. Камке Э. Справочник по обыкновенным дифференциальным уравнениям. М., Наука, 1976; Спб. ``Лань, 2003.
   
6. Л.Д.Ландау, Е.М.Лифшиц. Механика, ``Наука, М., 1965, 1988.
   

Тема 6.

1. Братусь А.С., Новожилов А.С., Платонов А.П. Динамические системы и модели биологии. М., ФИЗМАТЛИТ, 2010.
   
2. А.С.Братусь, С.В.Дрожжин, Т.С.Якушкина. Математические модели эволюции и динамики репликаторных систем. М., Ленанд, 2022.
   

Тема 7.

1. Бабенко К.И. Основы численного анализа. М., Наука, 1986, 2002.
   
2. Пайтген Х.-О., Рихтер П.Х. Красота фракталов. Образы комплексных динамических систем. М., Мир, 1993.
   

Тема 8.

1. Пайтген Х.-О., Рихтер П.Х. Красота фракталов. Образы комплексных динамических систем. М., Мир, 1993.
   

Тема 9.

1. Камке Э. Справочник по обыкновенным дифференциальным уравнениям. М., Наука, 1976; Спб. ``Лань, 2003.
   
2. Р.П.Федоренко: Лекции по вычислительной физике. М., МФТИ, 1994, 2008.
   

Тема 10.

1. Братусь А.С., Новожилов А.С., Платонов А.П. Динамические системы и модели биологии. М., ФИЗМАТЛИТ, 2010.
   
2. И.М.Гельфанд. Некоторые задачи теории квазилинейных уравнений. Успехи математических наук, 1959. Т.14. № 2, С.87-158.
   

Тема 11.

1. Айнс Э.Л. Обыкновенные дифференциальные уравнения. Харьков, Государственное научно-техническое издательство Украины, 1939, М. «Факториал Пресс», 2005.
   
2. В.И.Арнольд: Математические методы классической механики. ``Наука, М., 1989.
3. Л.Д.Ландау, Е.М.Лифшиц: Механика, ``Наука, М., 1965, 1988.
   

Тема 12.

1. Камке Э. Справочник по обыкновенным дифференциальным уравнениям. М., Наука, 1976; СПб.: "Лань, 2003.
   

Тема 14.

1. А.А.Андронов, А.А.Витт, С.Э.Хайкин: Теория колебаний, М., Физматгиз, 1959, ``Наука, 1981.
   
2. В.И.Арнольд: Математические методы классической механики. ``Наука, М., 1989.
   

Тема 16.

1. Блисс Г.А. Лекции по вариационному исчислению. М.: Иностранная литература, 1950.
   
2. Лаврентьев М.А., Люстерник Л.А. Курс вариационного исчисления. М.--Л.: ГИТТЛ, 1950.
   

Методические указания для обучающихся по освоению дисциплины

Вид учебных занятий/деятельности	Деятельность обучающегося
Лекция	Написание конспекта лекций: кратко, схематично, последовательно фиксировать основные положения лекции, выводы, формулировки, обобщения; помечать важные мысли, выделять ключевые слова, термины. Обозначить вопросы, термины или другой материал, который вызывает трудности, пометить и попытаться найти ответ в рекомендуемой литературе. Если самостоятельно не удается разобраться в материале, необходимо сформулировать вопрос и задать преподавателю на консультации, во время семинарского (практического) занятия.
Практическое (семинарское) занятие	При подготовке к семинарскому (практическому) занятию необходимо проработать материалы лекций, основной и дополнительной литературы по заданной теме. На основании обработанной информации постараться сформировать собственное мнение по выносимой на обсуждение тематике. Обосновать его аргументами, сформировать список источников, подкрепляющих его. Во время семинарского (практического) занятия активно участвовать в обсуждении вопросов, высказывать аргументированную точку зрения на проблемные вопросы. Приводить примеры из источниковой базы и научной и/или исследовательской литературы.
Устный/письменный опрос	Отвечать, максимально полно, логично и структурировано, на поставленный вопрос. Основная цель – показать всю глубину знаний по конкретной теме или ее части.
Подготовка к промежуточной аттестации	При подготовке к промежуточной аттестации необходимо проработать вопросы по темам, которые рекомендуются для самостоятельной подготовки. При возникновении затруднений с ответами следует ориентироваться на конспекты лекций, семинаров, рекомендуемую литературу, материалы электронных и информационных справочных ресурсов, статей. Если тема вызывает затруднение, четко сформулировать проблемный вопрос и задать его преподавателю.
Самостоятельная работа	Самостоятельная работа состоит из следующих частей: 1) чтение учебной, справочной, научной литературы; 2) повторение материала лекций; 3) составление планов устных выступлений; 4) подготовка видеопрезентации. При чтении учебной литературы нужно разграничивать для себя материал на отдельные проблемы, концепции, идеи. Учебную литературу можно найти в электронных библиотечных системах, на которые подписан АНО Университет Иннополис.
Контрольная работа	При подготовке к контрольной работе необходимо проработать материалы лекций, семинаров, основной и дополнительной литературы по заданной теме.
Разработка отдельных частей кода	Разработать часть кода, исходя из поставленной задачи и рекомендаций преподавателя. При выполнении работы рекомендуется обращаться к материалам лекций и семинарских (практических) занятий. Если возникают затруднения, необходимо проконсультироваться с преподавателем.
Выполнение домашних заданий и групповых проектов	Для выполнения домашних заданий и групповых проектов необходимо получить формулировку задания от преподавателя и убедиться в понимании задания. При выполнение домашних заданий и групповых проектов необходимо проработать материалы лекций, основной и дополнительной литературы по заданной теме.

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

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

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