Волновое уравнение

Поделись знанием:
(перенаправлено с «Уравнение колебания струны»)
Перейти к: навигация, поиск

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





Вид уравнения

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

<math>\Delta u=\frac{1}{v^2}\frac{\partial^2 u}{\partial t^2}</math>,

где <math>\Delta</math> — оператор Лапласа, <math>u=u(x,t)</math> — неизвестная функция, <math>t\in \mathbb R</math> — время, <math>x\in \mathbb R^n</math> — пространственная переменная, <math>v</math> — фазовая скорость.

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

<math>\frac{\partial^2 u}{\partial x^2}=\frac{1}{v^2}\frac{\partial^2 u}{\partial t^2}</math>.

Оператор Д’Аламбера

Разность <math>\Delta - \frac{1}{v^2} \frac{\partial^2}{\partial t^2}</math> называется оператором Д’Аламбера и обозначается как <math>\square</math> (разные источники используют разный знак). Таким образом, с использованием оператора Д’Аламбера (даламбертиана) однородное волновое уравнение записывается как:

<math>\square u = 0</math>

Неоднородное уравнение

Допустимо также рассматривать неоднородное волновое уравнение

<math>\frac{\partial^2 u}{\partial t^2}=v^2\Delta u + f</math>,

где <math>f = f(x,t)</math> — некая заданная функция внешнего воздействия (внешней силы).

Стационарным вариантом волнового уравнения является уравнение Лапласа (уравнение Пуассона в неоднородном случае).

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

<math> u(x,t) = v(x) e^{i\omega t}\ </math> или <math> u(x,t) = v(x)\, \mathop{\rm cos}\,(\omega t)\ </math>.

Решение волнового уравнения

Существует аналитическое решение гиперболического уравнения в частных производных. В евклидовом пространстве произвольной размерности оно называется формулой Кирхгофа. Частные случаи: для колебания струны (<math>\mathbb{R}^1</math>) — формула Д’Аламбера, для колебания мембраны (<math>\mathbb{R}^2</math>) — формула Пуассона.

Формула Д'Аламбера

Решение одномерного волнового уравнения (здесь <math> v = a </math> — фазовая скорость)

<math>u_{tt}=a^2 u_{xx} + f(x,t)\quad</math> (функция <math>f(x,t)</math> соответствует вынуждающей внешней силе)

с начальными условиями

<math>u(x,0)=\varphi(x),\quad u_t(x,0)=\psi(x)</math>

имеет вид

<math>u(x,t)=\frac{\varphi(x+at)+\varphi(x-at)}{2}+\frac{1}{2a}\int\limits^{x+at}_{x-at}{\psi(\alpha)d \alpha}+\frac{1}{2a}\int\limits^t_0\int\limits^{x+a(t-\tau)}_{x-a(t-\tau)} f(s, \tau)ds d\tau</math>

Интересно заметить, что решение однородной задачи

<math>u_{tt}=a^2 u_{xx}</math>,

имеющее следующий вид

<math>u(x,t)=\frac{\varphi(x+at)+\varphi(x-at)}{2}+\frac{1}{2a}\int\limits^{x+at}_{x-at}{\psi(\alpha)d \alpha}</math>

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

<math>u(x,t)= f_1(x+at) + f_2(x-at)</math>

где

<math> f_1(x)= \frac{\varphi(x)}{2} + \frac{1}{2a}\int\limits^{x}_{0}{\psi(\alpha)d \alpha} </math>
<math> f_2(x)= \frac{\varphi(x)}{2} + \frac{1}{2a}\int\limits^{0}_{x}{\psi(\alpha)d \alpha} </math>

В таком случае говорят, что решение представлено в виде суммы бегущих волн, а функции <math>f_1(x)</math> и <math>f_2(x)</math> — это профили волн, бегущих, соответственно, влево и вправо. В рассматриваемом случае профили волн со временем не изменяются.

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

Задача на полупрямой

Рассмотрим однородное уравнение колебаний на полупрямой <math>[0; +\infty)</math>

<math>u_{tt} = a^2 u_{xx} </math>

с закрепленным концом:

<math>u(0,t) = 0 </math>

и начальными условиями

<math>u(x,0)=\varphi(x),\qquad u_t(x,0)=\psi(x)</math>

для того, чтобы задача имела решение, необходима согласованность начальных условий и граничного условия, а именно:

<math>\varphi(0) = 0,\qquad \psi(0) = 0</math>

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

<math>\varphi(-x)=-\varphi(x),\qquad \psi(-x)=-\psi(x) \qquad \forall x \in [0, +\infty)</math>

В силу того, что начальные условия <math>\varphi(x), \psi(x)</math> — нечетные функции, логично ожидать, что и решение <math>u(x,t)</math> будет нечетной функцией. В этом можно непосредственно убедиться, рассмотрев решение в виде формулы Д’Аламбера. Поэтому полученное решение u(x, t) будет удовлетворять начальным условиям и граничному условию <math>u(0,t) = 0 </math> (последнее следует из нечетности функции).

Показанный прием широко используется (не только для волнового уравнения) и называется метод отражения. Например, можно рассмотреть волновое уравнение на полупрямой, но с граничным условием второго рода на конце <math>x = 0</math>:

<math>u_x(0,t)=0</math>.

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

Методы решения в ограниченной одномерной области

Метод отражений

Рассмотрим одномерное однородное волновое уравнение на отрезке <math>[0,a]</math>

<math>u_{tt}=a^2 u_{xx} </math>

с однородными граничными условиями первого рода (то есть при закрепленных концах)

<math>u(0,t)=0 \qquad u(a,t)=0</math>

и начальными условиями

<math>u(x,0)=\varphi(x),\quad u_t(x,0)=\psi(x) \qquad \forall x \in [0, a]</math>

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

<math>\varphi(2na + x) = \varphi(x) \qquad \psi(2na + x) = \psi(x) \qquad \forall x \in [0,a] \quad \forall n \in Z</math>
<math>\varphi(2na - x) = - \varphi(x) \qquad \psi(2na - x) = -\psi(x) \qquad \forall x \in [0,a] \quad \forall n \in Z</math>

При рассмотрении неоднородного волнового уравнения:

<math>u_{tt}=a^2 u_{xx} + f(x,t)</math>

используются ровно те же соображения, и функция <math>f(x,t)</math> продолжается таким же образом.

Метод Фурье

Снова рассмотрим одномерное однородное волновое уравнение на отрезке <math>[0,l]</math>

<math>u_{tt}=a^2 u_{xx} </math>

с однородными граничными условиями первого рода

<math>u(0,t)=0 \qquad u(l,t)=0</math>

и начальными условиями

<math>u(x,0)=\varphi(x),\quad u_t(x,0)=\psi(x) \qquad \forall x \in [0, l]</math>

Метод Фурье основывается на представлении решения в виде (бесконечной) линейной комбинации простых решений задачи вида

<math>X(x)T(t)</math>, где обе функции зависит только от одной переменной.

Отсюда другое название метода — метод разделения переменных.

Нетрудно показать, что для того, чтобы функция <math>u(x,t)=X(x)T(t)</math> была решением уравнения колебаний и удовлетворяла граничным условиям, необходимо, чтобы выполнялись условия

<math>X(0) = 0 \qquad X(l) = 0</math>
<math> a^2 X(x) = - \lambda X(x) </math>
<math> T(t) = - \lambda T(t) </math>

Решение задачи Штурма-Лиувилля на <math>X(x)</math> приводит к ответу:

<math>X_n(x) = \sin \left( \frac{\pi n x}{l} \right) \qquad n \in \mathbf{N}</math>

и их собственным значениям <math>\lambda_n = \left(\frac {\pi n a}{l}\right)^2</math>

Соответствующие им функции <math>T</math> выглядят как

<math>T_n(t) = \alpha_n \sin ( \sqrt \lambda_n t ) + \beta \cos ( \sqrt \lambda_n t ). </math>

Таким образом, их линейная комбинация (при условии, что ряд сходится) является решением смешанной задачи

<math>u(x,t) = \sum_{n=1}^{+\infty} X_n(x)T_n(t)

= \sum_{n=1}^{+\infty} \left( \alpha_n \sin ( \sqrt \lambda_n t ) + \beta_n \cos ( \sqrt \lambda_n t ) \right) \sin \frac{\pi n x}{l}. </math>

Разложив функции <math>\varphi(x), \psi(x)</math> в ряд Фурье, можно получить коэффициенты <math>\alpha_n, \beta_n </math>, при которых решение будет обладать такими начальными условиями.

Метод учета волн

Снова рассмотрим одномерное однородное волновое уравнение на отрезке <math>[0,a]</math>

<math>u_{tt}= u_{xx}, </math>

однако на сей раз положим однородные начальные условия

<math>u(x,0) \equiv 0,\quad u_t(x,0) \equiv 0 \qquad \forall x \in [0, a]</math>

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

<math>u(0,t)=\mu(t) \qquad u(a,t)=\nu(t)</math>

Решение записывается в виде

<math>u(x,t)= \sum_{k=0}^{+\infty} \biggl[ \mu(t - x - 2ka) - \mu(t + x - (2k+2)a) \biggr] + \sum_{k=0}^{+\infty} \biggl[ \nu(t + x - (2k+1)a) - \nu(t - x - (2k+1)a) \biggr]</math>

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

<math> \mu(t-x), </math>

которая, добегая за время а до правого конца, отражается и дает вклад

<math> \mu(t+x-2a), </math>

через время а снова отражается и дает вклад

<math> \mu(t-x-2a), </math>

Этот процесс продолжается бесконечно долго, суммируя вклады всех волн и получаем указанное решение. Если нас интересует решение на промежутке <math>[0,T]</math>, то мы можем ограничиться лишь первыми <math> \lceil T / a \rceil</math> слагаемыми.

См. также

Напишите отзыв о статье "Волновое уравнение"

Ссылки