Концепция решения линейных уравнений в частных производных
В математике фундаментальное решение для линейного оператора в частных производных L представляет собой формулировку на языке теории распределения старой идеи функции Грина (хотя, в отличие от функций Грина, фундаментальные решения не учитывают граничные условия).
В терминах дельта-функции Дирака δ ( x ) фундаментальное решение F является решением неоднородного уравнения
LF знак равно δ ( Икс ) .
Здесь F априори только предполагается распределением .
Эта концепция уже давно используется для лапласиана в двух и трех измерениях. Он был исследован для всех измерений лапласиана Марселем Риссом .
Существование фундаментального решения для любого оператора с постоянными коэффициентами — важнейший случай, напрямую связанный с возможностью использования свертки для решения произвольной правой части — было показано Бернаром Мальгранжем и Леоном Эренпрейсом . В контексте функционального анализа фундаментальные решения обычно разрабатываются с помощью альтернативы Фредгольма и исследуются в теории Фредгольма .
Пример
Рассмотрим следующее дифференциальное уравнение Lf = sin( x ) с
![{\displaystyle L={\frac {d^{2}}{dx^{2}}}.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Фундаментальные решения могут быть получены путем решения LF = δ ( x ) в явном виде:
![{\displaystyle {\frac {d^{2}}{dx^{2}}}F (x)=\delta (x)\,.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Поскольку для функции единичного шага (также известной как функция Хевисайда ) H имеем
![{\displaystyle {\frac {d}{dx}}H (x)=\delta (x)\,,}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle {\frac {d}{dx}}F(x)=H(x)+C\,.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
CC = −1/2После интегрирования и выбора новой константы интегрирования равной нулю, имеем![{\displaystyle {\frac {dF}{dx}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle F(x)=xH(x)-{\frac {1}{2}}x={\frac {1}{2}}|x|~.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Мотивация
Как только фундаментальное решение найдено, найти решение исходного уравнения несложно путем свертки фундаментального решения и желаемой правой части.
Фундаментальные решения играют важную роль также при численном решении уравнений в частных производных методом граничных элементов .
Приложение к примеру
Рассмотрим оператор L и дифференциальное уравнение, упомянутое в примере:
![{\displaystyle {\frac {d^{2}}{dx^{2}}}f(x)=\sin(x)\,.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Решение исходного уравнения можно найти путем свертки (обозначена звездочкой) правой части с фундаментальным решением : ![{\ displaystyle f (x)}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\ displaystyle \ грех (х)}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\textstyle F(x)={\frac {1}{2}}|x|}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle f(x)=(F*\sin)(x):=\int _{-\infty }^{\infty }{\frac {1}{2}}|xy|\sin(y) \,ди\,.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Это показывает, что необходимо соблюдать некоторую осторожность при работе с функциями, которые не имеют достаточной регулярности (например, компактный носитель, интегрируемость L 1 ), поскольку мы знаем, что искомое решение есть f ( x ) = −sin( x ) , в то время как приведенное выше интеграл расходится для всех x . Однако оба выражения для f равны как распределения.
Пример, который более наглядно работает
![{\displaystyle {\frac {d^{2}}{dx^{2}}}f(x)=I(x)\,,}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Iхарактеристическая (показательная) функция[0,1]IF ( x ) = | х |/2![{\displaystyle (I*F)(x)={\begin{cases}{\frac {1}{2}}x^{2}-{\frac {1}{2}}x+{\frac {1 }{4}},&0\leq x\leq 1\\|{\frac {1}{2}}x-{\frac {1}{4}}|,&{\text{иначе}}\end {случаи}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
IДоказательство того, что свертка является решением
Обозначим свертку функций F и g как F ∗ g . Скажем, мы пытаемся найти решение Lf = g ( x ) . Мы хотим доказать, что F ∗ g является решением предыдущего уравнения, т. е. мы хотим доказать, что L ( F ∗ g ) = g . При применении дифференциального оператора L к свертке известно, что
![{\displaystyle L(F*g)=(LF)*g\,,}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
LЕсли F является фундаментальным решением, правая часть уравнения сводится к
![{\displaystyle \delta *g~.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Но поскольку дельта-функция является единичным элементом для свертки, это просто g ( x ) . Подводя итоги,
![{\ displaystyle L (F * g) = (LF) * g = \ delta (x) * g (x) = \ int _ {- \ infty } ^ {\ infty } \ delta (xy) g (y) \ ,dy=g(x)\,.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Следовательно, если F является фундаментальным решением, свертка F ∗ g является одним решением Lf = g ( x ) . Это не значит, что это единственное решение. Можно найти несколько решений для разных начальных условий.
Фундаментальные решения некоторых уравнений в частных производных
С помощью преобразования Фурье можно получить следующее:
Уравнение Лапласа
Для уравнения Лапласа
![{\displaystyle [-\Delta ]\Phi (\mathbf {x},\mathbf {x} ') = \delta (\mathbf {x} -\mathbf {x} ')}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \Phi _{\textrm {2D}}(\mathbf {x} ,\mathbf {x} ')=- {\frac {1}{2\pi }}\ln |\mathbf {x} - \mathbf {x} '|,\qquad \Phi _ {\textrm {3D}}(\mathbf {x} ,\mathbf {x} ')={\frac {1}{4\pi |\mathbf {x } -\mathbf {x} '|}}~.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Экранированное уравнение Пуассона
Для экранированного уравнения Пуассона
![{\displaystyle [-\Delta +k^{2}]\Phi (\mathbf {x},\mathbf {x} ')=\delta (\mathbf {x} -\mathbf {x} '),\quad k\in \mathbb {R},}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \Phi _{\textrm {2D}}(\mathbf {x} ,\mathbf {x} ')={\frac {1}{2\pi }}K_{0}(k|\mathbf { x} -\mathbf {x} '|),\qquad \Phi _{\textrm {3D}}(\mathbf {x} ,\mathbf {x} ')={\frac {\exp(-k|\ mathbf {x} -\mathbf {x} '|)}{4\pi |\mathbf {x} -\mathbf {x} '|}},}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
модифицированная функция Бесселя![{\displaystyle K_{0}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
В более высоких измерениях фундаментальное решение экранированного уравнения Пуассона дается потенциалом Бесселя .
Бигармоническое уравнение
Для бигармонического уравнения
![{\displaystyle [-\Delta ^{2}]\Phi (\mathbf {x},\mathbf {x} ')=\delta (\mathbf {x} -\mathbf {x} ')}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \Phi _{\textrm {2D}}(\mathbf {x},\mathbf {x} ')=- {\frac {|\mathbf {x} -\mathbf {x} '|^{2 }}{8\pi }}\ln |\mathbf {x} -\mathbf {x} '|,\qquad \Phi _{\textrm {3D}}(\mathbf {x} ,\mathbf {x} ' )={\frac {|\mathbf {x} -\mathbf {x} '|}{8\pi }}~.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Обработка сигнала
В обработке сигналов аналог фундаментального решения дифференциального уравнения называется импульсной характеристикой фильтра.
Смотрите также
Рекомендации