- •Предисловие
- •Глава 1. Дифференциальные уравнения в частных производных, их классификация и методы решения. Задача Коши
- •Дифференциальные уравнения в частных производных
- •Дифференциальные уравнения в частных производных второго порядка
- •Понятие об общем решении уравнения в частных производных
- •Начальные и граничные условия
- •Задача Коши для дифференциальных уравнений первого порядка
- •Задача Коши для дифференциальных уравнений второго порядка
- •Классификация уравнений второго порядка в частных производных
- •Численные методы решения дифференциальных уравнений в частных производных
- •Использование метода сеток для решения параболических уравнений в частных производных
- •Использование метода сеток для решения гиперболических уравнений
- •Использование метода сеток для решения эллиптических уравнений
- •Глава 2. Уравнение колебания струны
- •Вывод уравнения колебания струны
- •Методы решения уравнения колебания струны
- •Метод Даламбера (метод бегущих волн) для бесконечной струны
- •Фазовая плоскость
- •Метод продолжений для полубесконечной струны
- •Метод продолжения для конечной струны (начальная и конечная точки жёстко закреплены)
- •Метод Фурье (метод стоячих волн) или метод разделения переменных
- •Упражнения
- •Глава 3. Решение уравнений теплопроводности и Лапласа методом Фурье
- •Вывод уравнения теплопроводности
- •Метод Фурье для конечного стержня
- •Упражнения.
- •Уравнение Лапласа
- •Задача Дирихле для уравнения Лапласа в круге
- •Задача Дирихле для уравнения Лапласа в кольце
- •Уравнение Лапласа в прямоугольнике
- •Приложения
- •Практикум в среде matlab
- •Лабораторная работа № 1
- •Задача Коши для систем обыкновенных дифференциальных уравнений
- •1. Метод Эйлера.
- •2. Решение уравнений p-го порядка и систем обыкновенных дифференциальных уравнений.
- •3. Методы Рунге–Кутты.
- •4. Решение задачи Коши с помощью встроенных функций matlab.
- •Задания для выполнения лабораторной работы
- •Лабораторная работа № 2 Волновые уравнения
- •1. Задача Коши для неограниченной струны.
- •2. Задача Коши для полуограниченной струны. Метод продолжений.
- •3. Создание анимации средствами matlab.
- •4. Фазовая плоскость. Бегущие волны.
- •Индивидуальные задания
- •Рекомендуемая литература
- •Оглавление
Упражнения.
Решить уравнение теплопроводности при заданных начальных и граничных условиях:
Ответ:
.
Ответ:
.
Ответ:
.
Ответ:
.
Ответ:
.
Уравнение Лапласа
При исследовании стационарных процессов различной физической природы (колебания, теплопроводность, диффузия и др.) обычно приходят к уравнениям эллиптического типа. Наиболее распространенным уравнением этого типа является уравнение Лапласа
.
Функция u называется гармонической в области Г, если она непрерывна в этой области вместе со своими производными до 2-го порядка и удовлетворяет уравнению Лапласа.
При изучении свойств гармонических функций были разработаны различные математические методы, оказавшиеся плодотворными и в применении к уравнениям гиперболического и параболического типов. Мы будем искать решение краевых задач для простейших областей методом разделения переменных. Решение краевых задач для уравнения Лапласа может быть найдено методом разделения переменных в случае некоторых простейших областей (круг, прямоугольник, шар, цилиндр и др.). Рассмотрим некоторые из них.
Задача Дирихле для уравнения Лапласа в круге
Найти
функцию u,
удовлетворяющую
уравнению Лапласа внутри круга и
граничному условию
на
границе круга, где
- заданная
функция,
- полярный угол.
Введем
полярную систему координат
с
началом в центре круга:
Уравнение Лапласа в полярных координатах имеет вид
.
Решим уравнение методом разделения переменных, т.е. будем искать частное решение уравнения в виде
.
Подставляя предполагаемую форму решения в уравнение в полярных координатах, получим:
.
Отсюда
.
Получим два обыкновенных дифференциальных уравнения:
Определим знак .
1.
Пусть
,
например
.
Рассмотрим
уравнение
.
Характеристическое уравнение имеет вид
-
это
решение не подходит, так как при изменении
угла
на
величину
однозначная
функция
должна
вернуться к исходному значению
(условие
периодичности).
Отсюда
следует, что
,
т.е.
является периодической функцией угла
с
периодом
.
2.
Пусть
,
тогда
.
-
это
решение подходит при условии
.
Рассмотрим второе уравнение системы:
.
Пусть
,
тогда
.
Получаем:
-
решение уравнения в
общем
случае.
3.
Пусть
,
например
.
Тогда
решение уравнения
:
.
Рассмотрим второе уравнение системы
.
Функцию
будем
искать в виде
.
Тогда уравнение принимает вид
;
;
;
.
Следовательно,
-
решение уравнения, где С и D
- постоянные.
Для
решения внутренней задачи надо положить
,
так как,
если
,
то функция
обращается
в бесконечность
при
и не является гармонической функцией
внутри
круга. Итак, частные решения нашей задачи
найдены:
.
Вид общего решения
.
Удовлетворим краевому условию:
.
Считая,
что
задана
как функция угла
,
возьмем ее разложение в ряд Фурье:
,
где
;
;
.
Будем использовать формулы Эйлера:
;
.
Подставляя выражения для коэффициентов Фурье в решение и меняя порядок суммирования и интегрирования, получим:
Подставляя в это выражение фомулы Эйлера, получаем интегральную формулу, дающую решение задачи
