🌐 Матфизика

Как решать уравнения в частных производных: методы и разбор

Q: Чем уравнение в частных производных отличается от обычного дифференциального?

В ОДУ одна независимая переменная, решение — кривая, и производные берутся по ней одной. В УЧП переменных несколько, производные частные, а решение — поверхность или поле, например температура в стержне для каждой точки и каждого момента. Обычные дифференциальные уравнения при этом никуда не уходят: временной множитель $T(t)$ в методе Фурье ищется именно как ОДУ.

Q: Что делать, если уравнение неоднородное — есть правая часть?

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

Q: Когда метод характеристик лучше разделения переменных?

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

Q: Обязательно ли запоминать собственные значения?

Полезно, но вывод важнее: на контрольной их легко восстановить. Подставьте $X = \sin(kx)$ в условия $X(0) = X(L) = 0$: второй конец требует $\sin(kL) = 0$, откуда $kL = n\pi$ и $\lambda_n = (n\pi/L)^2$. Для стержня с утеплёнными концами (условия Неймана) появятся косинусы и лишнее значение $\lambda_0 = 0$ — проверьте сами.

Классификация по дискриминанту B²−4AC, выбор метода и полный разбор колебания струны: ряд с коэффициентами через интегралы и чек-лист проверки краевых условий.

Сложность: Продвинутый 2 схемы действий

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

Логика такая: классификация подсказывает физику задачи, физика — тип условий, тип условий — метод. Струна колеблется (гиперболика), тепло растекается (параболика), поле установилось (эллиптика). Перепутать нельзя: к эллиптическому уравнению начальные условия не ставят вовсе, а гиперболическому нужны сразу два.

Шаг 1. Классифицируйте уравнение#

Уравнение второго порядка $A u_{xx} + B u_{x y} + C u_{y y} + \dots = 0$ (здесь $x$ и $y$ — независимые переменные, старшие производные — только вторые) классифицируется одним числом — дискриминантом $D = B^{2} - 4 A C$ . Считать его проще, чем дискриминант квадратного трёхчлена: коэффициенты стоят прямо перед производными.

D = B^{2} - 4 A C : D > 0 — гиперболическое, D = 0 — параболическое, D < 0 — эллиптическое

классификация уравнения второго порядка — первый шаг любого решения УЧП

Тип	Эталон	Физика	Метод по умолчанию
гиперболическое, $D > 0$	волновое $u_{tt} = a^{2} u_{xx}$	колебания, волны	характеристики (задача Коши), разделение переменных (отрезок)
параболическое, $D = 0$	теплопроводность $u_{t} = a^{2} u_{xx}$	диффузия, остывание	разделение переменных, преобразование Фурье
эллиптическое, $D < 0$	Лапласа $u_{xx} + u_{y y} = 0$	стационарное поле, прогиб мембраны	разделение переменных, только краевые условия

три типа уравнений матфизики и их эталоны — спутать типы значит решать чужую задачу

Проверим эталоны. Волновое $u_{tt} - a^{2} u_{xx} = 0$ : перед $u_{xx}$ стоит $- a^{2}$ , то есть $A = - a^{2}$ , $B = 0$ , $C = 1$ , и $D = 4 a^{2} > 0$ — гиперболическое (подробный урок). Теплопроводность $u_{t} - a^{2} u_{xx} = 0$ : второй порядок только у $u_{xx}$ , $A = - a^{2}$ , $C = 0$ , $D = 0$ — параболическое (урок). Уравнение Лапласа $u_{xx} + u_{y y} = 0$ : $A = C = 1$ , $D = - 4 < 0$ — эллиптическое. Тип может зависеть от точки: у уравнения Трикоми $u_{xx} + y u_{y y} = 0$ при $y > 0$ он эллиптический, при $y < 0$ — гиперболический, и такие уравнения честно называются смешанными.

Шаг 2. Выберите метод под область и условия#

Классификация говорит, чего ждать от решения; метод выбирают область и условия. Инструментов три:

Разделение переменных — конечная область (отрезок, прямоугольник) с однородными краевыми условиями: термин и урок метода.
Метод характеристик — гиперболические уравнения и задача Коши на всей прямой: урок метода характеристик; его вершина — формула Даламбера.
Преобразование Фурье — бесконечная или полубесконечная область: урок о преобразовании Фурье для УЧП.

Сколько условий ставить? Эллиптическому — только краевые, по всей границе. Параболическому — одно начальное плюс краевые. Гиперболическому — два начальных (отклонение и скорость) плюс краевые. Это не прихоть преподавателя: без второго начального условия у волнового уравнения останутся неопределёнными коэффициенты при синусах времени, и задача будет иметь бесконечно много решений. Постановка — половина ответа; термин краевая задача систематизирует варианты.

Шаг 3. Разделение переменных: задача о струне целиком#

Струна длины $L = 1$ м закреплена на концах, скорость всех точек в начальный момент нулевая, отклонение задано: $f (x) = sin (π x)$ . Полная постановка:

u_{tt} = a^{2} u_{xx}, u (0, t) = u (1, t) = 0, u (x, 0) = f (x), u_{t} (x, 0) = 0, 0 < x < 1

волновое уравнение с закреплёнными концами: два начальных и два краевых условия; a² = T/ρ — отношение натяжения струны к её линейной плотности

Ищем решение произведением $u (x, t) = X (x) T (t)$ — свой множитель для координаты, свой для времени. Допущение смелое, но метод Фурье на нём стоит: суперпозиция потом склеит любое начальное условие.
Подставляем: $X T^{''} = a^{2} X^{''} T$ . Делим на $a^{2} X T$ : слева функция только от $t$ , справа — только от $x$ . Равенство при любых $x$ и $t$ возможно, лишь когда обе части — одна и та же константа $- λ$ .
Краевые условия переезжают на множитель: $u (0, t) = X (0) T (t) = 0$ при любом $t$ даёт $X (0) = 0$ , аналогично $X (1) = 0$ . Получили задачу Штурма–Лиувилля $X^{''} + λ X = 0$ , $X (0) = X (1) = 0$ — систематически она разобрана в уроке о собственных функциях.
Разбираем знаки $λ$ : при $λ < 0$ решением будут гиперболические функции, и оба краевых условия дают тривиальное $X \equiv 0$ ; случай $λ = 0$ тоже вырождается в ноль. Нетривиальные решения живут только при $λ_{n} = n^{2} π^{2}$ : $X_{n} (x) = sin (nπ x)$ , $n = 1, 2, 3, \dots$
Временной множитель: $T_{n}^{''} + a^{2} n^{2} π^{2} T_{n} = 0$ — гармонические колебания с частотой $ω_{n} = nπ a$ : $T_{n} = b_{n} cos (nπ a t) + c_{n} sin (nπ a t)$ . Это обычное линейное ОДУ второго порядка, метод здесь.
Собираем ряд: каждое слагаемое $u_{n} = T_{n} sin (nπ x)$ решает уравнение и краевые условия, а суперпозиция решает тоже: $u (x, t) = \sum_{n = 1}^{\infty} (b_{n} cos (nπ a t) + c_{n} sin (nπ a t)) sin (nπ x)$ .
Начальное отклонение: $u (x, 0) = \sum b_{n} sin (nπ x) = f (x)$ — разложение $f$ в ряд Фурье по синусам (урок): $b_{n} = 2 \int_{0}^{1} f (s) sin (nπ s) d s$ .
Начальная скорость: $u_{t} (x, 0) = \sum c_{n} nπ a sin (nπ x) = 0$ , и раз ряд по синусам нулевой — каждый $c_{n} = 0$ . Нулевые скорости гасят все синусы времени одним махом.

b_{n} = 2 \int_{0}^{1} f (s) sin (nπ s) d s, c_{n} = \frac{2}{nπ a} \int_{0}^{1} ψ (s) sin (nπ s) d s

коэффициенты через начальные условия: f — начальное отклонение, ψ — начальная скорость; здесь ψ = 0, и все c_n = 0

Считаем коэффициенты для $f (x) = sin (π x)$ — и здесь сработает ортогональность синусов: $\int_{0}^{1} sin (nπ s) sin (mπ s) d s = 0$ при $n \neq = m$ .

$b_{1} = 2 \int_{0}^{1} sin^{2} (π s) d s = 2 \cdot \frac{1}{2} = 1$ : среднее значение $sin^{2}$ на целом числе полупериодов — половина.
$b_{2} = 2 \int_{0}^{1} sin (π s) sin (2 π s) d s$ . Произведение синусов превращаем в разность косинусов (шпаргалка по тригонометрии): $sin α sin β = \frac{1}{2} (cos (α - β) - cos (α + β))$ , и интеграл распадается на $\frac{1}{2} \int_{0}^{1} (cos (π s) - cos (3 π s)) d s = 0$ — оба косинуса дают ноль на целом числе полупериодов. Численное интегрирование подтверждает: $- 3 \cdot 1 0^{- 17}$ , ноль.
$b_{3} = 2 \int_{0}^{1} sin (π s) sin (3 π s) d s = 0$ — та же ортогональность. И все остальные $b_{n}$ тоже нули: $f (x) = sin (π x)$ — уже готовая первая мода, раскладывать нечего.

u (x, t) = cos (π a t) sin (π x)

ответ для f(x) = sin(πx): стоячая волна — основная мода; все точки качаются синфазно, узлы сидят на закреплённых концах, период равен 2/a

Проверка подстановкой — секунда дела, а уверенность на всю задачу: $u_{tt} = - (π a)^{2} cos (π a t) sin (π x)$ , $u_{xx} = - π^{2} cos (π a t) sin (π x)$ , и уравнение $u_{tt} = a^{2} u_{xx}$ выполняется тождественно. На концах синус даёт ноль при любом $t$ ; при $t = 0$ косинус равен единице, то есть $u = f$ ; производная по времени в нуле пропорциональна $sin 0 = 0$ . Бонус-контроль: максимум отклонения не превышает $max ∣ f ∣ = 1$ — для нулевой начальной скорости это гарантировано формулой Даламбера.

Тренировочный профиль: $f (x) = x (1 - x)$ . Двойное интегрирование по частям (приёмы — здесь) даёт $b_{n} = \frac{4 ( 1 - ( - 1 ) ^{n} )}{n ^{3} π ^{3}}$ : у нечётных $b_{1} = \frac{8}{π ^{3}} \approx 0, 2580$ , $b_{3} = \frac{8}{27 π ^{3}} \approx 0, 0096$ , у чётных — нули. Провал на высоких гармониках ( $1/ n^{3}$ ) — наглядный ответ, почему параболическая форма звучит глухо, а чистый $sin (π x)$ — как камертон: у него вообще одна гармоника.

Ряд Фурье: сумма синусов

Добавляй гармоники!

Серые волны — гармоники $\frac{4}{π n} sin (n x)$ для нечётных $n$ . Фиолетовая — их сумма. Добавляйте слагаемые: сумма всё больше похожа на прямоугольную волну. В этом смысл ряда Фурье: любую периодическую функцию можно собрать из синусов.

Число гармоник3Анимация фазы

Каждая серая волна — «нота», а их сумма — «аккорд». Именно так аудиокодеки (MP3) и JPEG-сжатие раскладывают сигнал на частоты.

Двиньте ползунок до n=12 — углы квадрата всё ещё чуть «звенят». Это явление Гиббса!

Шаг 4. Чек-лист проверки решения#

Начальные условия: подставьте $t = 0$ — ряд обязан свестись к $f (x)$ внутри отрезка, а $u_{t} (x, 0)$ — к начальной скорости. Для одного слагаемого это мгновенно, для ряда — проверьте пару коэффициентов интегралом.
Краевые условия: $x = 0$ и $x = L$ обязаны давать ноль при любом $t$ , причём каждый член ряда по отдельности. Не даёт — вы выбрали не ту систему собственных функций (косинусы вместо синусов, например).
Само уравнение: подставьте общий член ряда и продифференцируйте. Для временного множителя всё сводится к $T_{n}^{''} = - a^{2} λ_{n} T_{n}$ — знаки обязаны сойтись.
Физический контроль: гиперболическое уравнение бережёт энергию — максимум отклонения не должен расти со временем; параболическое сглаживает контрасты; эллиптическое достигает максимума на границе, а не внутри области.

Типичные ошибки

Первая — потерять второе начальное условие: без

u_{t} (x, 0)

коэффициенты

c_{n}

повиснут в воздухе, и решение не будет единственным. Вторая — принять

λ < 0

: гиперболический косинус не пересекает ноль дважды, и краевые условия выкинут решение. Третья — забыть делитель в

c_{n}

: из

u_{t} (x, 0) = \sum c_{n} nπ a sin (nπ x)

следует

c_{n} = \frac{2}{nπ a} \int_{0}^{1} ψ sin (nπ x) d x

, а не

2 \int_{0}^{1} ψ sin

. Четвёртая — знак в дискриминанте: у волнового уравнения

A = - a^{2}

отрицателен, поэтому

D = 4 a^{2} > 0

; потерянный минус превращает гиперболику в псевдопараболику.

Проверь себя+15 XP

К какому типу относится уравнение $u_{y y} - 2 u_{xx} = 0$ ?

Проверь себя+15 XP

Струна, $u_{t} (x, 0) = 0$ , начальное отклонение $f (x) = sin (3 π x)$ . Каково решение?

Частые вопросы

Чем уравнение в частных производных отличается от обычного дифференциального?

В ОДУ одна независимая переменная, решение — кривая, и производные берутся по ней одной. В УЧП переменных несколько, производные частные, а решение — поверхность или поле, например температура в стержне для каждой точки и каждого момента. Обычные дифференциальные уравнения при этом никуда не уходят: временной множитель $T (t)$ в методе Фурье ищется именно как ОДУ.

Что делать, если уравнение неоднородное — есть правая часть?

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

Когда метод характеристик лучше разделения переменных?

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

Обязательно ли запоминать собственные значения?

Полезно, но вывод важнее: на контрольной их легко восстановить. Подставьте $X = sin (k x)$ в условия $X (0) = X (L) = 0$ : второй конец требует $sin (k L) = 0$ , откуда $k L = nπ$ и $λ_{n} = (nπ / L)^{2}$ . Для стержня с утеплёнными концами (условия Неймана) появятся косинусы и лишнее значение $λ_{0} = 0$ — проверьте сами.

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