Линейные уравнения система производная решения

Линейные дифференциальные уравнения в частных производных первого порядка

Линейные однородные уравнения в частных производных первого порядка

Пусть X ₁ , X ₂ , . X_n – заданные функции переменных x ₁ , x ₂ , . x_n .

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

необходимо решить систему обыкновенных дифференциальных уравнений (уравнение характеристик):
:
Далее нужно представить решение в виде:
φ ₁( x ₁ , x ₂ , . x_n ) = C ₁ ,
φ ₂( x ₁ , x ₂ , . x_n ) = C ₂ ,
.
φ_{n- 1} ( x ₁ , x ₂ , . x_n ) = C_{n- 1} ,
где C_k – постоянные.
После чего сразу получаем общее решение:
,
где F – произвольная функция от n – 1 аргументов.

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

Линейные неоднородные уравнения в частных производных первого порядка

Пусть X ₁ , X ₂ , . X_{n+ 1} – заданные функции от переменных x ₁ , x ₂ , . x_n и z .

Чтобы решить линейное неоднородное уравнение в частных производных первого порядка:
,
необходимо решить уравнение характеристик:
.
Решение этой системы нужно представить в следующем виде:
φ ₁( x ₁ , x ₂ , . x_n , z ) = C ₁ ,
φ ₂( x ₁ , x ₂ , . x_n , z ) = C ₂ ,
.
φ_n ( x ₁ , x ₂ , . x_n , z ) = C_n .
После чего сразу получаем общий интеграл в неявном виде:

где F – произвольная функция. Также общий интеграл можно представить в различных вариантах, например:
φ ₁ = F ( φ ₂ , φ ₃ , . φ_n ) ,
φ ₂ = F ( φ ₁ , φ ₃ , . φ_n ) ,
и т. д.

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

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

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

Это линейное однородное уравнение в частных производных первого порядка. Составляем уравнение характеристик:

Это уравнение характеристик содержит три уравнения:
;
;
.
Нам нужно выбрать и решить любые два из них. Тогда третье будет выполнено автоматически.

Выбираем и решаем первое уравнение:

Здесь переменные уже разделены, интегрируем:

Интегралы табличные,

Потенцируем:

Отсюда

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

Или:

Это линейное уравнение. Решаем с помощью интегрирующего множителя. Умножим на x -1 и преобразуем:

Интегрируем:

Подставим полученное ранее выражение C₁ = x y 2 :

Итак, мы нашли два интеграла уравнения характеристик:

Общее решение исходного уравнения в частных производных имеет вид:

где F — произвольная функция от двух аргументов F(φ₁, φ₂) . Найдем ее вид из граничного условия
при .

Рассматриваем решение на границе.
Положим x y = –1 :

Отсюда

На границе
.

Итак, мы нашли, что на границе функция F имеет вид:
F ( φ ₁ , φ ₂ ) = φ ₁ φ ₂ .
Такой же вид она имеет и во всей области
Подставляя
;
,
получаем частное решение исходного уравнения в частных производных с заданным граничным условием:

Общее решение:

где F — произвольная функция от двух аргументов F ( φ ₁ , φ ₂ ) .

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

Найти поверхность, удовлетворяющую данному уравнению
,
и проходящую через данную окружность x + y + z = 0 , x 2 + y 2 + z 2 = a 2 .

Это линейное неоднородное уравнение в частных производных первого порядка. Составляем уравнение характеристик:

Оно содержит три уравнения:
;
;
.
Нам нужно выбрать и решить любые два из них. Тогда третье удовлетворится автоматически. Выбираем первое и второе уравнения.

Решаем уравнение:

Умножаем на 2 z и интегрируем:

Интегралы табличные,

Потенцируем:

Отсюда
x = C ₁ y

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

Или:

Замечаем, что , тогда

Это линейное уравнение. Решаем с помощью интегрирующего множителя. Разделим на y 2 и преобразуем:

Интегрируем:

Подставим полученное ранее выражение и преобразуем:

Итак, мы нашли два интеграла уравнения характеристик:

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

Общий интеграл исходного уравнения в частных производных имеет вид:
F ( φ ₁ , φ ₂) = 0
Но, поскольку F — произвольная функция от двух аргументов, то общий интеграл можно записать также в виде:
φ ₁ = F ( φ ₂) ,
где F — произвольная функция от одного аргумента.

Найдем вид этой функции, рассматривая решение на границе.
На границе, x 2 + y 2 + z 2 = a 2 , .
Из уравнения x + y + z = 0 , z = – ( x + y ) . Подставим в x 2 + y 2 + z 2 = a 2 и преобразуем:
x 2 + y 2 + ( x + y ) 2 = a 2
x 2 + y 2 + x 2 + 2 xy + y 2 = a 2
2 x 2 + 2 xy + 2 y 2 = a 2
Разделив на y 2 , имеем

Итак, мы нашли, что на границе:

.
Подставим в выражение общего интеграла:
φ ₁ = F ( φ ₂)
.
Сделаем подстановку
:
.

Итак, мы нашли, что на границе функция F имеет вид:
.
Такой же вид она имеет и во всей области, тогда
.
Подставляем выражения для φ₁ и φ₂ :

.
Умножим на a 2 y 2 .

Автор: Олег Одинцов . Опубликовано: 23-09-2014

Системы уравнений с частными производными. Характеристики

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

Рассмотрим систему двух квазилинейных уравнений относительно искомых функций :

(2.60)

Коэффициенты этой системы переменные и зависят от х, t, U, V. Введем следующие обозначения: U — искомый вектор; F — вектор правой части; А, В — матрицы коэффициентов:

Запишем систему уравнений (2.60) в векторном виде:

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

Мы не будем повторять сказанное ранее для одного уравнения, а остановимся на одном частном случае системы (2.60), важном для приложений. Речь идет о системах гиперболического типа. Введем матрицу С где α, β — некоторые числа. Тогда определитель этой матрицы

(2.61)

является квадратичной формой относительно α и β,т.е.

(2.62)

где коэффициенты q1, q2, q3, легко выразить через элементы матриц А, В, раскрывая определитель (2.61).

Система уравнений с частными производными первого порядка (2.60) называется гиперболической,если квадратичная форма (2.62) разлагается на вещественные линейные множители:

причем векторы неколлинеарны. Эти векторы в каждой точке плоскости (х,t) образуют два направления, которые называются характеристическими. Линия, касательная к которой в каждой точке имеет характеристическое направление, называется характеристикой. Через каждую точку проходят две характеристики, соответствующие двум характеристическим направлениям. Таким образом, всю плоскость (х, t) можно покрыть двумя семействами характеристик (рис. 2.18).

Рис. 2.18. Характеристики

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

Наряду с гиперболическими системами существуют также параболические (с одним семейством характеристик) и эллиптические (действительных характеристик нет) системы.

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

Не приводя подробных выкладок и опуская сами формулы, изложим идею метода характеристик. Рассмотрим задачу Коши. Пусть при t = 0 заданы начальные значения функций . Выбираем любой отрезок [а,b] на оси х и разбиваем его на части точками (рис. 2.19). В данном случае принято n= 4.

Рис. 2.19. К решению задачи Коши методом характеристик

Из точки А0 проводим характеристику первого семейства, из А1 — второго. Находим точку пересечения В0. Используя некоторые соотношения (характеристические) вдоль отрезков характеристик А0В0 и А1В0,заменяющие исходные уравнения, вычисляем искомые функции в точке В0. Аналогично находим решение в других точках слоя В. При этом в отличие от метода сеток этот слой не является прямолинейным отрезком t= const, а определяется точками пересечения характеристик.

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

При решении краевой задачи используют значения искомых функций на границах. В этом случае расчетная область изменяется: она прилегает к границе х = const, на которой заданы значения функций U(x), V(x). При этом вблизи границы используют характеристики одного семейства, выходящие из границы и попадающие в расчетную область. Если граничные условия задают при двух значениях х, то алгоритм метода характеристик значительно усложняется.

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

Для устранения этого недостатка разработаны так называемые сеточно-характеристические методы. Их идея состоит в том, что сетка фиксируется заранее, а характеристики проводятся «назад» из узлов (j+ 1)-ого слоя до пересечения с j-ым слоем. Значения U, Vв точках пересечения вычисляются путем интерполяции по ранее найденному решению в узлах j-ого слоя.

Геометрическая интерпретация сеточно-характеристического метода показана на рис. 2.20. Здесь точками отмечены заранее выбранные узлы; штриховые линии — отрезки характеристик. Значения функций в точках пересечения А и В находятся интерполированием решения в узлах и Эти значения используют для определения решения в расчетном узле (i, j+1).

Рис. 2.20. Геометрическая интерпретация сеточно-характеристического метода

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

Дифференциальным уравнением с частными производными называется уравнение вида
(1)

связывающее независимые переменные x1, х2, … , хn искомую функцию и = и(х1, х2,…, хn) и ее частные производные (наличие хотя бы одной производной обязательно). Здесь ki,k2,… ,кn — неотрицательные целые числа, такие, что к1 + к2 + … + кп = т.

Порядком дифференциального уравнения называется наивысший порядок входящие в уравнение частных производных. Так, если х, у — независимые переменные, и = и(х, у) — искомая функция, то

— дифференциальное уравнение 1-го порядка;

— дифференциальные уравнения 2-го порядка.

Для упрощения записи пользуются также следующими обозначениями:

Пусть имеем дифференциальное уравнение с частными производными (1) порядка т. Обозначим через С m (D) множество функций, непрерывных в области D вместе со всеми производными до порядка m включительно.

Определение:

Решением дифференциального уравнения (1) в некоторой области D изменения независимых переменных x1, x2…xn,. называется всякая функция и = и(х1, х2,…, xп) ∈ С m (D) такая, что подстановка этой функции и ее производных в уравнение (1) обращает последнее в тождество по x1, x2, …., хп в области D.

Пример:

Найти решение и = и(х,у) уравнения

Равенство (2) означает, что искомая функция и не зависит опт х, но может быть любой функцией от у,

u = φ(y). (3)

Таким образом, решение (3) уравнения (2) содержит одну произвольную функцию. Это — общее решение уравнения (2).

Приме:

Найти решение u = u(z, у) уравнения

Положим = о. Тогда уравнение (4) примет вид = 0. Его общим решением будет произвольная функция v = w(у). Поскольку v= приходим к уравнению = w(у). Интегрируя по у (считая х параметром), получим

где g(x) — произвольная функция. Так как w(у) — произвольная функция, то и интеграл от нее также является произвольной функцией; обозначим его через f(у). В результате получим решение уравнения (4) в виде

u(x, y) = f(y) + g(x) (5)

произвольные дифференцируемые функции).

Решение (5) уравнения с частными производными 2-го порядка (4) содержит уже две произвольные функции. Его называют общим решением уравнения (4), так как всякое другое решение уравнения (4) может быть получено из (5) подходящим выбором функций f и g.

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

Пример:

Множество действительных решений уравнения

исчерпывается функцией u(x, y) = const, а уравнение

вовсе не имеет действительных решений.

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

Линейные дифференциальные уравнения с частными производными. Свойства их решений

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

Пример:

— линейное уравнение; уравнения

Линейное дифференциальное уравнение 2-го порядка для функции двух независимых переменных х, у в общем случае имеет вид
(1)

где А(х, у), В(х, у), …, с(х,у), f(x,y) — функции переменных х, у, заданные в некоторой области D плоскости хОу. Если f(x,y) ≡ 0 в D, то уравнение (1) называется однородным, в противном случае — неоднородным.

Обозначив левую часть уравнения (1) через L[u], запишем (1) в виде

L[u] = f(x, у). (2)

Соответствующее однородное уравнение запишется так:

L[u] = 0. (3)

Здесь L — линейный дифференциальный оператор, определенный на линейном пространстве C 2 (D) функций и = и(х, у).

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

Теорема:

Если и(х, у) есть решение линейного однородного уравнения (3), то си(х, у), где с — любая постоянная, есть также решение уравнения (3).

Теорема:

Если и₁(х, у) и и₂(х, у) — решения линейного однородного уравнения (3), то сумма и₁(х, у) + и₂(x, у) есть также решение этого уравнения.

Следствие:

Если каждая из функций и₁(х, у) и и₂(х, у), u k(x, у) является решением уравнения (3), то линейная комбинация

где c1, c2 …, сk — произвольные постоянные, также является решением этого уравнения.

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

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

Пример:

имеет общее решение k = φ(х), так что решениями его будут, например, функции 1,х,…, х n ,… . В соответствии с этим в линейных задачах для уравнений с частными производными нам придется иметь дело не только с линейными комбинациями конечного числа решений, но и с рядами , членами которых являются произведения постоянных Сп на частные решения иn(х, у) дифференциального уравнения.

Возможны случаи, когда функция и(х, у; λ) при всех значениях параметра λ из некоторого интервала (λо, λ1), конечного или бесконечного, является решением уравнения (3). В этом случае говорят, что решения уравнения зависят от непрерывно меняющегося параметра λ. Если теперь взять функцию С(λ) такую, что первые и вторые производные интеграла

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

L[u] = f (4)

справедливы следующие предложения.

Теорема:

Если и(х, у) есть решение линейного неоднородного уравнения (4), a v(x, у) — решение соответствующего однородного уравнения (3), то сумма и + v есть решение неоднородного уравнения (4).

Теорема:

Принцип суперпозиции. Если и1(х, у) —решение уравнения L[u] = f1, a u2(x,y) — решение уравнения L[u] = f2, то и1 + u2 — решение уравнения L[u] = f1 + f2.

Классификация линейных дифференциальных уравнений второго порядка с двумя независимыми переменными

Определение:

Линейное дифференциальное уравнение второго порядка

в некоторой области Q на плоскости хОу называется

1) гиперболическим в Ω, если

2) параболическим в Ω, если

3) эллиптическим в Ω, если

Пользуясь этим определением, легко проверить, что уравнения

— гиперболические при всех х и у, уравнение

— параболическое при всех х и у, а уравнение

— эллиптическое при всех х и у. Уравнение

— эллиптическое при у > 0, параболическое на линии у = 0 и гиперболическое в полуплоскости у

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

Уравнение гиперболического типа (∆ > 0) преобразуется к вшу

(два канонических вида уравнений гиперболического типа).

Уравнение параболического типа (∆ ≡ 0) преобразуется к виду

(канонический вид уравнения параболического типа).

Уравнение эллиптического типа (∆

(канонический вид уравнения эллиптического типа). Здесь F и Ф — некоторые функции, зависящие от искомой функции и, ее первых производных и независимых переменных ξ, η. Вид функций F и Ф определяется исходным уравнением (1).

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

Как правило, приведениеуравнения(1) к каноническому виду путем замены независимых переменных имеет локальный характер, т. е. осуществимо лишь в некоторой достаточно малой окрестности рассматриваемой точки Mo(xo, уo).

Когда число п независимых переменных больше двух, также различают уравнения гиперболического, параболического и эллиптического типов. Например, при п = 4 простейшая каноническая форма таких уравнений имеет вид

Здесь и = и(х, у, z, t).

Замечание:

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

к каноническому виду возможно только в данной точке и невозможно в любой сколь угодно малой окрестности этой точки.

Мы ограничимся рассмотрением линейных дифференциальных уравнений 2-го порядка. К таким уравнениям приводит большое количество различных физических задач.

Так, колебательные процессы различной природы (колебания струн, мембран, акустические колебания газа в трубах, электромагнитные колебания и т. д.) описываются уравнениями гиперболического типа. Простейшим из таких уравнений является уравнение колебаний струны (одномерное волновое уравнение): (2)

Здесь х — пространственная координата, t — время, где Т — натяжение струны, р — ее линейная плотность.

Процессы теплопроводности и диффузии приводят к уравнениям параболического типа. В одномерном случае простейшее уравнение теплопроводности имеет вид
(3)

Здесь где р — плотность среды, с — удельная теплоемкость, k — коэффициент теплопроводности.

Наконец, установившиеся процессы, когда искомая функция не зависит от времени, определяются уравнениями эллиптического типа, типичным представителем которых является уравнение Лапласа
(4)

Непосредственной проверкой убеждаемся в том, что решением уравнения (2) является всякая функция и(х, t) вида

Можно показать, что решениями уравнения (3) являются функции вида

произвольные постоянные, А — числовой параметр). Интегрируя решение и(х, t; λ) = уравнения (3) по параметру λ в пределах от — ∞ до + ∞ , получим так называемое фундаментальное решение U(x, t) = уравнения теплопроводности.

Наконец, нетрудно убедиться, что действительнозначные функции Рn(х,у) и Qn(x, у), определяемые из соотношения

являются решениями уравнения Лапласа (4) для п = 0, 1, 2…..Этот последний результат есть частный, случай общего утверждения, что и действительная и мнимая части аналитической функции

f(z) = u(x, у) + iv(x, у)

комплексного переменного z = х + iy являются решениями уравнения Лапласа (4).

В силу линейности уравнения (4) ряды

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

Таким образом, для простейшей — канонической — формы уравнений гиперболического, параболического и эллиптического типов мы располагаем о решениях этих уравнений некоторой информацией.

Постановка основных задач для линейных дифференциальных уравнений второго порядка

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

Пример:

Общее решение уравнения

имеет вид и(х, у) = f(x) + g(y), где f(x) и g(y) — произвольные дифференцируемые функции. Поэтому чтобы выделить решение, описывающее данный физический процесс, необходимо задать дополнительные условия.

Различают три основных типа задач для дифференциальных уравнений с частными производными (число независимых переменных равно п):

а) задача Коши для уравнений гиперболического и параболического типов: задаются начальные условия, область Ω совпадает со всем пространством R n , граничные условия отсутствуют;

б) краевая задача для уравнений эллиптического типа: задаются граничные условия на границе S области Ω, начальные условия отсутствуют;

в) смешанная задача для уравнений гиперболического и параболического типов: задаются начальные и граничные условия, Ω ≠ R n

Решение заданий и задач по предметам:

Дополнительные лекции по высшей математике:

Образовательный сайт для студентов и школьников

Копирование материалов сайта возможно только с указанием активной ссылки «www.lfirmal.com» в качестве источника.

© Фирмаль Людмила Анатольевна — официальный сайт преподавателя математического факультета Дальневосточного государственного физико-технического института