Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Методы математической физики.-1

.pdf
Скачиваний:
7
Добавлен:
05.02.2023
Размер:
589.02 Кб
Скачать

Министерство образования и науки Российской Федерации Федеральное государственное бюджетное образовательное учреждение

высшего профессионального образования «Томский государственный университет систем управления и

радиоэлектроники»

Кафедра электронных приборов

МЕТОДЫ МАТЕМАТИЧЕСКОЙ ФИЗИКИ

Методические указания к практическим занятиям для студентов направления «Фотоника и оптоинформатика» и

«Электроника и микроэлектроника»

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)ex .

Пример 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

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]