Методы математической физики.-1
.pdfМинистерство образования и науки Российской Федерации Федеральное государственное бюджетное образовательное учреждение
высшего профессионального образования «Томский государственный университет систем управления и
радиоэлектроники»
Кафедра электронных приборов
МЕТОДЫ МАТЕМАТИЧЕСКОЙ ФИЗИКИ
Методические указания к практическим занятиям для студентов направления «Фотоника и оптоинформатика» и
«Электроника и микроэлектроника»
2012
Гейко Павел Пантелеевич
Методы математической физики: методические указания к практическим занятиям для студентов направления «Фотоника и оптоинформатика» и «Электроника и микроэлектроника» /П.П. Гейко; Министерство образования и науки Российской Федерации, Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования Томский государственный университет систем управления и радиоэлектроники, Кафедра электронных приборов. - Томск: ТУСУР, 2012. -
31с.
Впособии рассматриваются основные понятия и определения, связанные с уравнениями с частными производными и вопросы приведения к каноническому виду линейных уравнений второго порядка. Излагаются вопросы, относящиеся к аналитическим методам решения основных уравнений математической физики (гиперболических, параболических и эллиптических) и численным методам их решений.
Пособие предназначено для студентов очной и заочной форм, обучающихся по направлению «Фотоника и оптоинформатика» и «Электроника и микроэлектроника» по дисциплине «Методы математической физики».
© Гейко Павел Пантелеевич , 2012
Министерство образования и науки Российской Федерации Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования
«Томский государственный университет систем управления и радиоэлектроники»
Кафедра электронных приборов
УТВЕРЖДАЮ Зав.кафедрой ЭП
_____________С.М. Шандаров «___» _____________ 2012 г.
МЕТОДЫ МАТЕМАТИЧЕСКОЙ ФИЗИКИ
Методические указания к практическим занятиям для студентов направления «Фотоника и оптоинформатика» и
«Электроника и микроэлектроника»
Разработчик
_________ П.П. Гейко «____»____________2012 г
2012
Содержание |
|
Введение |
4 |
1. Классификация линейных уравнений второго порядка |
4 |
2.Приведение линейных уравнений второго порядка к 6 канонической форме
3.Канонические формы линейных уравнений с постоянными 7 коэффициентами
4.Нахождение общего решения линейного однородного 10 уравнения 1-го порядка.
5.Краевая задача для однородного уравнения теплопроводности. 11
6. |
Краевая задача для однородного волнового уравнения. |
12 |
7. |
Краевая задача для неоднородного волнового уравнения. |
14 |
8.Формула Даламбера решения задачи Коши для волнового 15 уравнения
9. Краевые задачи для уравнения Лапласа |
17 |
10.Уравнение Лапласа в цилиндрических координатах. Решение 18 задачи Дирихле для кольца
11. Решение задачи Дирихле для уравнения Лапласа в круге |
19 |
12.Численные методы решения задач по уравнениям 21 математической физики
13. |
Варианты индивидуальных заданий |
22 |
14. |
Задачи к экзамену |
28 |
Литература |
29 |
ВВЕДЕНИЕ Математическая физика – это математический аппарат изучения физических полей –
одного из центральных объектов современной физики и инженерии. Только привлекая рассмотрение физических полей и соответствующий математический аппарат, удается наиболее полно описать физические явления, а в целом ряде случаев без такого привлечения даже не удается сформулировать первоначальные понятия и простейшие утверждения. Поэтому знание тех или иных разделов математической физики оказывается необходимым каждому современному инженеру.
Термин "математическая физика" имеет и более узкий, "классический" смысл. Он относится к уравнениям в частных производных, являющимся теоретическим аппаратом гидромеханики, теории теплопроводности и диффузии, теории упругости, "классической" части теории электромагнитного поля. Поля, рассматриваемые в этих классических разделах, оказывается возможным трактовать как механические системы с бесконечным числом степеней свободы, что и обусловило общность соответствующего математического аппарата.
1. Классификация линейных уравнений второго порядка
Будем рассматривать уравнения с частными производными второго порядка, линейные относительно старших производных, т.е. имеющие вид
a11 uxx + 2a12 uxy + a22 uyy |
+ F ( x, y, u, ux , uy ) = 0, |
(1.1) |
где a11 ,a12 ,a22 являются функциями x и y . |
|
|
|
|
|
С помощью преобразования переменных |
|
|
ξ = ϕ( x, y), η =ψ( x, y), |
|
допускающего обратное преобразование (для этого достаточно потребовать, чтобы функциональный определитель
D = |
|
ϕx |
ψ x |
||
|
|||||
|
|
ϕ |
y |
ψ |
y |
|
|
|
|
был отличен от нуля), можно получить уравнение, эквивалентное исходному. Нас будет интересовать вопрос: как выбрать новые переменные ξ и η, чтобы относительно них уравне-
ние имело наиболее простой (канонический) вид. Перейдя к новым переменным, будем иметь
ux = uξξx + uηηx , uy = uξξy + uηηy ,
u |
xx |
= (u ξ |
x |
+ u η |
x |
) |
ξ |
ξ |
x |
+ (u ξ |
x |
+ u η |
x |
) |
η |
η |
x |
= u ξ 2 |
+ u (ξ |
x |
) |
ξ |
ξ |
x |
+ |
|
|||||||||||||||||
|
|
|
ξ |
|
η |
|
|
|
|
|
ξ |
|
|
η |
|
|
|
|
|
ξξ x |
|
ξ |
|
|
|
|
|
|
|||||||||||||||
+ u η |
x |
ξ |
x |
+ u (η |
x |
) |
ξ |
ξ |
x |
+ u ξ |
x |
η |
x |
+ u (ξ |
x |
) η |
x |
+ u η2 |
+ u (η |
x |
) |
η |
η |
x |
= |
||||||||||||||||||
|
|
ηξ |
|
|
η |
|
|
|
|
ξη |
|
ξ |
|
|
|
|
η |
|
ηη |
x |
η |
|
|
|
|
|
|
||||||||||||||||
= uξξξx2 |
+ 2uξηξxηx + uηηηx2 + uξ |
[(ξx )ξ ξx + (ξx )ηηx ]+ uη [(ηx )ξ ξx + (ηx )ηηx ]= |
= uξξξx2 + 2uξηξxηx + uηηηx2 + uξξxx + uηηxx .
Аналогично получаем
uxy = uξξξxξy + uξη (ξxηy +ξ yηx ) + uηηηxηy + uξξxy + uηηxy ,
uyy = uξξξy2 + 2uξηξyηy + uηηη2y + uξξyy + uηηyy .
(1.2)
(1.3)
После подстановки полученных производных в (1.1) получим уравнение a11 uξξ + 2a12 uξη + a22 uηη + F1 (ξ,η, uξ , uη , u) = 0, (1.4)
где
4
a |
11 |
= a ξ 2 |
+ 2a ξ |
ξ |
y |
+a ξ 2 , |
|
|
11 x |
12 x |
|
22 y |
|
||
a12 |
= a11ξxηx + a12 (ξxηy +ηxξy ) +a22ξyηy , |
(1.5) |
|||||
a |
22 |
= a η2 |
+ 2a η η |
y |
+ a η2 . |
|
|
|
11 x |
12 x |
22 y |
|
Заметим, что если исходное уравнение линейно, т.е.
F ( x, y, u, ux , uy ) = b1 ux + b2 uy + cu + f ,
то F1 имеет вид
F1 (ξ,η, uξ , uη , u) = β1 uξ + β2 uη +γu +δ .
Таким образом, уравнение в этом случае снова получается линейным.
Попытаемся выбрать переменную ξ =ϕ( x, y) так, чтобы коэффициент a11 в уравнении (6.4) был равен нулю. Для этого необходимо, чтобы ξ =ϕ( x, y) было решением уравнения
a ξ2 |
+2a ξ |
ξ |
y |
+a ξ2 |
= 0. |
(1.6) |
11 x |
12 x |
|
22 y |
|
Уравнение (6.6) можно записать в виде произведения
(a11ξx − (− a12 + a122 − a11a22 )ξy )(a11ξx − (− a12 − a122 − a11a 22 )ξy ).
Таким образом, решение уравнения (1.6) свелось к решению двух линейных однородных уравнений первого порядка
a11ξx − (− a12 ± |
|
a122 |
− a11a22 )ξy = 0. |
(1.7) |
|||
для решения уравнений (1.7) надо найти общий интеграл каждого из |
|
||||||
уравнений |
|
|
|
|
|
|
|
dy |
= |
a12 ± |
a122 −a11a22 |
. |
(1.8) |
||
dx |
|
a11 |
|
||||
На вид решений уравнений (1.8) существенно влияет знак подкоренного выражения |
|||||||
a122 −a11a22 . По знаку этого выражения определяется тип уравнения (1.1). |
|
||||||
Будем называть уравнение (1.1) в точке |
M |
|
|
||||
гиперболического типа, если a122 − a11a22 |
> 0 , |
|
|
||||
эллиптического типа, если a122 |
− a11a22 < 0 , |
|
|
|
|||
параболического типа, если a122 − a11a22 |
= 0 . |
|
|
|
|||
Можно убедиться в справедливости равенства |
|
|
a122 −a11a22 = (a122 −a11a22 )D2 ,
из которого следует, что тип уравнения не меняется при преобразовании переменных. Следует отметить также, что тип уравнения зависит от точки M и в разных точках мо-
жет быть разным.
Пример 1.1. Рассмотрим уравнение |
|
uxx + xuyy = 0, |
(1.9) |
здесь a11 = 1, a12 = 0 и a22 = x , следовательно, |
|
a122 − a11a22 = −x. |
|
Тем самым при x < 0 уравнение (1.9) гиперболического типа, при x = 0 |
– параболиче- |
ского типа, а при x > 0 – эллиптического типа. |
|
5
2. Приведение линейных уравнений второго порядка к канонической форме
Уравнение
a11dy2 −2a12dxdy +a22dx2 = 0 |
(2.1) |
будем называть характеристическим для уравнения (6.1), а его интегралы – характеристиками. Уравнение (2.1) распадается на два уравнения (1.8) и играет основную роль в задаче приведения к каноническому виду уравнения (1.1). Затем, что для уравнения гиперболического типа характеристики действительные и различные, для уравнений эллиптического типа – комплексные и различные, а для уравнений параболического типа обе характеристики действительные и совпадают.
Разберем каждый из этих случаев в отдельности.
1. Для уравнений гиперболического типа a122 − a11a22 > 0 и правые части уравнений (1.8) действительные и различные. Общие интегралы их ϕ( x, y) = C и ψ( x, y) = C опреде-
ляют семейства характеристик. Положим
|
ξ = ϕ( x, y), η =ψ( x, y), |
тогда коэффициенты |
a11 и a22 (1.5) обратятся в нуль и уравнение (1.4) после деления на |
коэффициент при uξη |
приведется к виду |
|
uξη = Ф(ξ,η, u, uξ , uη ), |
где Ф = − |
F1 |
. Это – так называемая каноническая форма уравнений гиперболического ти- |
||||||||
|
||||||||||
|
2a12 |
|
|
|
|
|
|
|||
па. Часто пользуется другой канонической формой. Положим |
||||||||||
|
|
|
|
α = |
ξ +η |
, |
β = |
ξ −η |
, |
|
где α и β - новые переменные. Тогда |
2 |
|
2 |
|
|
|||||
|
|
|
|
|
|
|||||
|
|
uξ = |
1 |
(uα + uβ ), |
uη = (uα − uβ ), |
|
uξη = (uαα − uββ ) |
|||
|
|
|
|
|||||||
2 |
|
|
|
|
|
|
|
|||
и уравнение (2.2) примет вид |
uαα − uββ |
= Ф1 , |
|
|
||||||
|
|
|
|
|
|
|
||||
где Ф1 = 4Ф. |
|
|
|
|
|
|
||||
2. Для уравнений параболического типа a122 |
− a11a22 |
= 0, и уравнение (1.8) дает один |
||||||||
общий интеграл ϕ( x, y) = C. Положим в этом случае |
|
|
||||||||
|
|
|
|
ξ =ϕ( x, y), |
η =ψ( x, y), |
где ψ( x, y) - любая функция, допускающая вместе с ϕ( x, y) обратное преобразование переменных. Тогда
a11 |
= a11ξx |
+ 2a12ξxξy + a12ξy |
|
ξx + |
a12 |
|
2 |
= a11 |
ξy |
= 0 , |
|||||
|
2 |
2 |
|
|
|
|
|
|
|
|
|
a11 |
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
a12 |
|
a12 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
a12 |
= a11 |
ξx + |
|
ξ ηx + |
|
ηy |
= 0, |
|
|
|
|
|
|
a11 |
a11 |
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
||
т.к. a 2 |
= a |
11 |
a |
22 |
. После деления уравнения 6.4) на коэффициент при u |
получим канониче- |
|||||||
12 |
|
|
|
|
|
|
|
|
|
ξη |
скую форму для уравнения параболического типа
uηη = Ф(ξ,η, u, uξ , uη ),
где Ф = − F11 .
a22
6
3. Для уравнения эллиптического типа a122 −a11a12 < 0 и правые части уравнения (1.8)
комплексно сопряженные, поэтому общие интегралы этих уравнений будут также комплексно сопряженными
ϕ( x, y) = C, ϕ ( x, y) = C .
Чтобы не иметь дела с комплексными переменными и функциями, введем новые, уже вещественные, переменные α и β
|
|
α = |
ϕ + |
ϕ |
, |
β = |
ϕ − |
ϕ |
|
, |
2 |
|
2i |
|
|||||||
|
|
|
|
|
|
|
|
|||
т.е. ϕ =α + iβ, |
ϕ |
=α − iβ. Таким образом, |
|
|
|
|
|
|
|
0 = a ϕ2 |
+ 2a ϕ |
ϕ |
y |
+a ϕ2 = a |
11 |
(α |
x |
+ iβ |
x |
)2 |
+ 2a |
12 |
(α |
x |
+ iβ |
x |
)(α |
y |
+ iβ |
y |
) + |
||||||||||||||||||
|
|
|
|
11 |
x |
|
|
12 |
x |
|
|
22 y |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
+a |
22 |
(α |
y |
+ iβ |
y |
)2 |
= (a α |
2 + 2a α |
α |
y |
+a α2 ) −(a |
11 |
β 2 |
+ 2a |
12 |
β |
x |
β |
y |
+a |
22 |
β 2 ) + |
|||||||||||||||||
|
|
|
|
|
|
|
|
|
11 |
x |
12 |
x |
|
|
22 |
|
y |
|
|
x |
|
|
|
|
|
|
|
y |
||||||||||||
|
+ 2i(a11αx βx +a12 (αx βy +αy βx ) +a22αy βy ). |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
Отсюда следует (см. (6.5) ), что |
|
a11 = a22 |
и a12 |
= 0. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
Уравнение (1.4) после деления на коэффициент при uαα |
принимает вид |
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
uαα + uββ = Ф(α, β, u, uξ , uβ ), |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
где Ф = − |
F1 |
. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
a11 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Итак, в зависимости от знака выражения |
|
a122 − a11a22 (т.е. от типа уравнения) получаем |
||||||||||||||||||||||||||||||||||||||
следующие канонические формы уравнения (1.1): |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
1. |
|
|
Гиперболический тип: |
uξη |
|
= Ф или uαα − uββ |
= Ф. |
|
|
|
|
|
|
|
2.Параболический тип: uηη = Ф.
3.Эллиптический тип: uαα + uββ = Ф.
3. Канонические формы линейных уравнений с постоянными коэффициентами
Линейное уравнение второго порядка с постоянными коэффициентами имеет вид
a11 uxx + 2a12 uxy + a22 uyy + b1 ux + b1 uy + cu + f ( x, y) = 0. |
(3.1) |
Решая уравнение (1.8) получаем характеристики, которые будут прямыми линиями y = λ1 x +C1 , y = λ2 +C2 ,
где λ1 и λ2 корни уравнения (его удобно в данном случае тоже называть характеристическим)
|
|
|
a11λ2 −2a12λ +a22 = 0. |
|
|
(3.2) |
||
Теперь с помощью соответствующего преобразования переменных, а именно: |
||||||||
1. |
Если λ1 и λ2 вещественные и различные (гиперболический тип) |
|
|
|||||
|
|
ξ = y − λ1 x, η = y − λ2 x или |
ξ = y − |
λ1 + λ2 |
x, η = |
λ2 − λ1 |
x; |
|
|
|
|
2 |
|||||
2. |
Если |
λ1 = λ2 = λ (параболический тип) |
2 |
|
|
|||
|
|
|
|
|
||||
|
|
|
ξ = y − λx, η = x; |
|
|
|
||
3. |
Если |
λ1,2 = a ± ib |
(b ≠ 0) (эллиптический тип) |
|
|
|
||
|
|
|
ξ = y − ax, |
η = bx; |
|
|
|
|
уравнение (1.1) приводится к одному из следующих видов: |
|
|
|
|||||
1. |
uξη +Ф = 0 или uξξ |
− uηη +Ф = 0; |
|
|
|
|
|
7
2.uηη +Ф = 0;
3.uξξ − uηη +Ф = 0;
здесь Ф = β1 uξ + β2 uη + cu + f .
Пример 3.1. Привести к каноническому виду уравнение |
(3.3) |
|
uxx + 3uxy + 2uyy = 0. |
||
|
||
Напишем характеристическое уравнение (3.2) |
|
|
λ2 − 3λ + 2 = 0. |
|
Следовательно, λ1 = 1, λ2 = 2 . Уравнение - гиперболического типа, поэтому делаем замену
ξ = y − x, |
η = y − 2 x или ξ = y − |
3 |
x, |
η = x. |
|
||||
|
2 |
|
|
|
После замены переменных уравнение имеет вид uξη = 0 или uξξ |
− uηη = 0. |
|||
Заметим, что решение уравнения uξη |
= 0 рассматривалось в примере 1.3. Тем самым мы мо- |
жем выписать общее решение уравнения (3.3)
u = ϕ(ξ) +ψ(η) = ϕ( y − x) +ψ( y − 2 x).
Пример 3.2. Привести к каноническому виду уравнение
uxx + 2uxy + uyy + ux + uy = 0. |
(3.4) |
Решая характеристическое уравнение
λ2 −2λ +1 = 0,
получаем λ1 = λ2 = 1. Следовательно, уравнение (3.4) параболического типа. Делаем замену ξ = y − x, η = x . Имеем
ux = uξξx + uηηx = uη − uξ , uy = uξξy + uηηy = uξ ,
uxx = (uη − uξ )ξ ξx + (uη − uξ )ηηx = −(uηξ − uξξ ) + (uηη − uξη ) = uξξ − 2uξη + uηη , uyy = (uξ )ξ ξy + (uξ )ηηy = uξξ ,
uxy = (uη − uξ )ξ ξy + (uη − uξ )ηηy = uξη − uξξ .
Подставляя полученные выражения в уравнение (3.4) и приводя подобные члены, будем иметь
uηη + uη = 0.
Заметим, что мы получили уравнение, которое можно рассматривать и как обыкновенное дифференциальное уравнение, зависящее от параметра ξ. Решая его, получаем
u = C1 (ξ) + C2 (ξ)e−η = C1 ( y − x) + C2 ( y − x)e−x .
Пример 3.3. Привести к каноническому виду уравнение
uxx − 2uxy + 2uyy + ux + uy = 0. |
(3.5) |
Для корней характеристического уравнения
λ2 + 2λ + 2 = 0
имеем λ1,2 = −1 ± i . Уравнение – эллиптического типа, поэтому делаем замену
ξ = x + y, η = x.
Подставляя выражения
8
ux |
= uξξx + uηηx |
= uξ + uη , |
uy |
= uξξy + uηηy |
= uξ , |
uxx = uξξ + 2uξη + uηη , |
uxy = uξξ + uξη , |
uyy |
= uξξ , |
в уравнение (3.5), получаем |
|
(3.6) |
uξξ + uηη + 2uξ + uη |
= 0. |
Для линейных уравнений второго порядка с постоянными коэффициентами возможны дальнейшие упрощения канонической формы уравнений. Введем для этого вместо u новую функцию ϑ :
u = eλξ +μηϑ,
где λ и μ - некоторые постоянные. Тогда
uξ = eλξ +μη (ϑξ + λϑ), |
|
uη = eλξ +μη (ϑη + μϑ), |
|
uξξ = eλξ +μη (ϑξξ + 2λϑξ + λ2ϑ), |
(3.7) |
uξη = eλξ +μη (ϑξη + λϑη + μϑξ + λμϑ), |
|
uηη = eλξ +μη (ϑηη + 2μϑη + μ2ϑ). |
|
Подставляем эти выражения, например, в уравнение
uξη + β1 uξ + β1 uη + cu + f = 0
и сокращая затем на eλξ +μη , получим
ϑξη +(μ + β1 )ϑξ +(λ + β2 )ϑη +(λμ + β1λ + β2 μ +C )ϑ + f1 = 0.
Выберем параметры λ и μ так, чтобы коэффициенты при первых производных обратились
в нуль (λ = −β2 , μ = −β1 ) . В результате получим
ϑξη +γϑ + f1 = 0, где γ = λμ + β1λ + β2 μ + c = c − β1 β2 , f1 = fe−(λξ +μη) .
Аналогично упрощения проводятся и для других канонических форм. Окончательно приходим к следующим каноническим формам для уравнений с постоянными коэффициентами:
1. Гиперболический тип
ϑξη +γϑ + f1 = 0 или ϑξξ −ϑηη −γϑ + f1 = 0.
2. Параболический тип
ϑηη + β1ϑξ + f1 = 0.
3. Эллиптический тип
ϑξξ +ϑηη +γϑ + f1 = 0.
Пример 3.4. Упростим уравнение (3.6), полученное в примере 3.3. Подставляя выражения для производных (8.7) и сокращая на e λξ +μη будем иметь
ϑξξ +ϑηη + 2(λ +1)ϑξ + (2μ +1)ϑη + (λ2 + μ2 − 2λ + μ)ϑ = 0.
Выберем λ = −1, μ = − 12 . Тогда уравнение будет иметь вид
9