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

Математика. Двумерное преобразование Фурье и его приложения. Учебное пособие

.pdf
Скачиваний:
0
Добавлен:
09.08.2026
Размер:
474 Кб
Скачать

№ 4751 МИНИСТЕРСТВО НАУКИ И ВЫСШЕГО ОБРАЗОВАНИЯ РФ

УНИВЕРСИТЕТ НАУКИ И ТЕХНОЛОГИЙ МИСИС

ИНСТИТУТ БАЗОВОГО ОБРАЗОВАНИЯ

Кафедра математики

Н.Е. Цапенко

МАТЕМАТИКА

Двумерное преобразование Фурье и его приложения

Учебное пособие

Рекомендовано редакционно-издательским советом университета

Москва 2023

УДК 51 Ц17

Р е ц е н з е н т ы :

канд. физ.-мат. наук, доц. кафедры математического и компьютерного моделирования В.П. Григорьев (Институт информационных

и вычислительных технологий НИУ «МЭИ»)

Цапенко, Николай Евгеньевич.

Ц17 Математика. Двумерное преобразование Фурье и его приложения : учеб. пособие / Н.Е. Цапенко. – Москва : Издательский Дом НИТУ МИСИС, 2023. – 49 с.

ISBN 978-5-907560-83-3

Исходя из формул Грина и их следствий, записаны интегралы Фурье от функций, заданных в различных плоских областях, как ограниченных, так и неограниченных. Выведены Фурьепредставления частных производных первого и высших порядков, в которые входят граничные значения как самой функции, так и её производных. С помощью этих представлений записаны решения граничных задач для двумерного уравнения Лапласа в случае нескольких плоских областей. Также рассмотрена система уравнений Максвелла, описывающая распространение электромагнитных волн в слоисто-неоднородной среде.

Предназначено для студентов старших курсов и магистров, специализирующихся по направлению «Прикладная математика».

 

УДК

51

 

Н.Е. Цапенко, 2023

 

ISBN 978-5-907560-83-3

НИТУ МИСИС, 2023

СОДЕРЖАНИЕ

ПРЕДИСЛОВИЕ. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .4

ГЛАВА 1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .5 1.1. Формулы Грина . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5 1.2. Двумерные преобразования Фурье . . . . . . . . . . . . . . . . 7 1.3. Функция, заданная в бесконечной полосе . . . . . . . . . 11 1.4. Функция, заданная в бесконечном квадрате. . . . . . . . . . . . . . . . . . . . . . . . . . . . 13 1.5. Круговая область.

Преобразование Фурье–Бесселя. . . . . . . . . . . . . . . . . . . . . 15

ГЛАВА 2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .24 2.1. Гармонические функции . . . . . . . . . . . . . . . . . . . . . . . 24 2.2. Интегральные уравнения для граничных значений гармонических функций . . . . . . . . . . . . . . . . . . 30 2.3. Решение уравнений (2.48), (2.49)

в случае бесконечной полосы . . . . . . . . . . . . . . . . . . . . . . . 33 2.4. Решение уравнений (2.48), (2.49)

в круговой области. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 40

ГЛАВА 3 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .42 3.1. Е- и H-волны в слоисто-неоднородной среде . . . . . . . 42

Библиографический список . . . . . . . . . . . . . . . . . . . . . . . . .48

3

ПРЕДИСЛОВИЕ

В нашем университете существует такой курс, как «Уравнения в частных производных математической физики». Он читается для студентов, обучающихся по специальности «Прикладная математика» в течение 4-го учебного семестра.

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

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

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

4

ГЛАВА 1

1.1. Формулы Грина

Пусть в некоторой плоской области D с границей L заданы две непрерывные и дифференцируемые в каждой точке этой области функции P(x, y) и Q(x, y). Имеют место следующие интегральные соотношения:

 

 

 

 

P

 

 

 

 

 

 

 

dxdy Pdx ,

(1.1)

 

 

y

 

 

 

D

 

 

 

 

L

 

 

 

 

 

 

Q

 

 

 

 

 

 

 

dxdy Qdy

(1.2)

 

 

 

x

 

 

 

D

L

 

и, соответственно,

 

 

 

 

 

 

 

 

 

 

Q

 

P

 

 

 

x

 

 

y dxdy Pdx Qdy .

(1.3)

D

 

 

 

 

 

 

L

 

Полагая в формулах (1.1), (1.2) P Q u v , получим двумерные аналоги обычной формулы интегрирования по частям:

v

u

 

 

v

y dxdy uvdx u y dxdy ,

D

 

L

D

 

 

 

u

 

 

v

v x dxdy uvdy u

 

dxdy .

x

D

 

 

D

 

 

(1.4)

(1.5)

Если теперь в формуле (1.4) функцию u(x, y) заменить на ее производную по y, а в формуле (1.5) эту же функцию u(x, y) поменять на ее производную по x и получившиеся таким образом формулы сложить, то придем к такой формуле:

 

 

 

 

 

u v

 

u v

 

 

u

u

 

v udxdy

 

 

 

 

dxdy v

x dy

y dx , (1.6)

x x

y y

D

 

 

D

 

 

 

 

 

 

L

 

 

 

где

2

 

2

–двумерный оператор Лапласа.

 

 

x2

y2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5

В последнем равенстве (1.6) мы можем поменять роли букв u и v . Тогда

 

 

 

u v

 

u v

 

 

 

v

v

 

u vdxdy

 

 

 

 

dxdy u

 

dy

y dx . (1.7)

x x

y y

x

D

D

 

 

 

 

 

 

L

 

 

 

 

 

Наконец, вычитая из (1.6) (1.7), получим такую важную формулу, тоже известную, как формула Грина:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

v u u v dxdy

 

 

 

 

 

D

 

 

 

 

 

 

 

 

(1.8)

 

 

u

u

 

 

 

v

v

 

 

 

 

v

x dy

y dx

u

 

dy

y dx .

 

x

 

L

 

 

 

 

L

 

 

 

 

 

 

Эту формулу можно записать без привязки к прямоугольной системе координат, если использовать понятие производной по внешней нормали.

Введем для номинальных векторов плоской кривой (рис. 1.1) следующие обозначения:

– единичный вектор касательной к данной кривой;

n – единичный вектор внешней нормали к данной кривой; dl dxi dyj – вектор элемента дуги плоской кривой;

 

dl

dl

 

dx2 dy2

– величина элемента дуги плоской

кривой;

 

 

 

 

 

 

 

 

 

 

 

 

k

n

 

– единичный вектор, перпендикулярный к пло-

скости данной кривой, то есть к плоскости рисунка.

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

Рис. 1.1

 

 

 

 

 

 

 

 

6

 

 

 

Вэтих обозначениях имеем

dl dx i dy j ;

dl dl dl

n k dy i dx j. dl dl

Производная функции u = u(x, y) по направлению внешней нормали записывается так:

u u dy u dx .n x dl y dl

Поэтому

 

 

 

 

u dy

u dx

 

u

dl ,

 

 

x

y

n

и, соответственно, формула Грина (1.8) может быть представлена в такой инвариантной форме

 

 

v u u v ds

 

v

 

u

 

u

 

v

 

dl.

(1.9)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

n

 

D

 

 

L

 

 

 

 

 

 

 

 

 

 

1.2. Двумерные преобразования Фурье

Пусть в некоторой односвязной плоской области D с границей L задана функция двух переменных f(x, y). Прямым преобразованием Фурье этой функции называется двойной интеграл:

F 1, 2

 

1

e i 1x i 2yf x,y dxdy .

(1.10)

 

2

 

2 D

 

Функция F 1, 2

от двух параметров 1, 2 называет-

ся спектральной характеристикой функции f x,y

(или, ее

спектральной плотностью), а сами параметры 1, 2

называ-

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

7

 

 

 

 

 

 

 

 

 

 

 

 

 

f

x,y

 

 

ei 1x i 2y F

,

d d

2

.

(1.11)

 

 

 

 

1

2

1

 

 

 

 

 

 

 

 

 

 

 

 

 

С помощью преобразования Фурье можно решать многие

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

Выразим спектральные характеристики

производных

функции f(x, y) через спектральную характеристику самой этой функции и ее граничные значения. Имеем

Fx1 1, 2

1

 

e i 1x i 2y

 

f

dxdy ;

2

 

 

 

 

2

 

D

 

x

Fy1 1, 2

 

1

 

e i 1x i 2y

f

dxdy.

2

 

 

 

 

2

 

D

 

y

Применяя к этим интегралам формулы (1.4), (1.5), получим

Fx1 1, 2 i 1F 1, 2

1

e i 1x i 2yfdy;

(1.12)

2

 

 

2

L

 

Fy1 1, 2 i 2F 1, 2

 

1

e i 1x i 2yfdx.

(1.13)

 

2

 

 

2

L

 

Для спектральных характеристик вторых производных

2

1, 2

 

1

 

 

e

i x i y 2f

Fxx

 

 

 

 

1

2

 

 

 

 

dxdy;

 

2

 

 

x

2

 

 

 

 

 

 

2

 

D

 

 

 

 

 

 

 

 

2

1, 2

 

1

 

 

e

i x i y 2f

Fyy

 

 

 

 

1

2

 

 

 

dxdy;

 

2

 

 

y

2

 

 

 

 

 

2

 

D

 

 

 

 

 

 

 

 

2

1, 2

 

1

 

e

i x i y

 

2f

Fxy

 

 

 

 

1

2

 

 

 

 

 

dxdy,

 

2

 

 

 

 

 

 

 

 

 

 

 

2

 

D

 

 

 

x y

аналогично, по формулам (1.4), (1.5) (или по уже написанным формулам (1.12), (1.13)), получаются такие выражения:

8

Fxx2 1, 2 12F 1, 2

 

 

 

 

 

 

1

 

 

e

i x i y

 

 

 

f

 

 

 

 

 

 

 

 

1

2

i 1f

 

 

 

dy;

 

 

 

2

 

x

 

 

 

 

 

 

 

2

L

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Fyy2 1, 2 22F 1, 2

 

 

 

 

 

 

 

1

 

 

e

i x i y

 

 

f

 

 

 

 

 

 

 

 

 

 

 

 

 

1

2

i 2f

 

 

 

 

dx;

 

 

 

 

 

 

2

 

 

 

y

 

 

 

 

 

 

2

 

L

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Fxy2 1, 2 1 2F 1 2

 

 

 

 

 

1

e

i x i y

 

 

 

 

f

 

 

 

 

 

 

 

 

 

 

 

 

 

1

2

i 2fdy

 

 

dx

 

 

 

 

 

 

 

2

 

 

x

 

2

L

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 2F 1, 2

 

 

 

 

 

 

 

 

 

 

 

1

 

e

i x i

y

 

 

 

 

 

f

 

 

 

 

 

 

 

 

 

 

 

1

2

i 1fdx

 

 

dy .

 

 

 

 

2

 

 

y

 

 

 

2

L

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1.14)

(1.15)

(1.16)

Сумма преобразований (1.14), (1.15) даст преобразование от двумерного оператора Лапласа:

 

1

 

e i 1x i 2y fdxdy 12

22 F 1 2

 

 

 

 

2

 

 

2

D

 

 

 

 

 

 

 

 

 

(1.17)

 

 

1

 

 

e

i x i y

 

f

 

f

 

 

 

 

 

 

 

1 2

i 1f

 

dy i 2f

 

dx .

 

 

2

 

x

y

 

2

L

 

 

 

 

 

 

 

 

 

Если функция f(x, y) удовлетворяет уравнению Лапласа, т.е. если

f 2f 2f 0,x2 y2

то ее спектральная плотность определяется через ее граничные значения так:

9

F 1, 2

1

 

 

1

e

i x i y

 

 

 

 

 

 

1

2

2

 

2

2

 

 

2

1

2

L

 

 

 

 

(1.18)

 

f

 

 

 

 

f

 

 

 

 

 

i 2f

 

 

dx i 1f

 

 

 

dy .

 

 

 

 

 

 

 

y

 

 

 

 

x

 

 

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

f x,y ei 1x i 2y F 1, 2 d 1d 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

ei 1x i 2y

 

 

 

i

x i

y

 

 

 

 

 

 

 

 

 

e

 

1

2

 

 

 

2

 

2

2

 

 

 

2

 

 

1

2

 

L

 

 

 

 

 

 

 

f

 

 

 

 

 

f

 

 

 

 

i 2f

 

 

 

dx i 2f

 

 

 

dy d 1d 2.

y

 

 

 

 

 

 

 

 

 

 

x'

 

 

 

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

 

 

f

f

 

 

f x,y

G x x ,y y

 

dy

 

dx

 

x

y

L

 

 

 

 

 

 

f x ,y Gx x x ,y y dy Gy x x ,y y dx .

L

Здесь использованы обозначения:

G ,

 

 

1

 

 

 

 

 

e

i 1 i 2

 

 

 

1

 

ln 2 2 ;

 

 

 

 

 

d 1d 2

 

 

 

2

 

2

2

 

 

 

 

 

 

 

 

2

 

 

 

 

1

2

 

 

 

4

 

 

 

 

 

 

G

 

 

 

1

 

 

 

i ei 1 i 2

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

d 1d 2

 

 

 

 

 

 

 

 

;

 

 

 

2

 

2

2

 

 

 

 

2

2

 

 

2

 

 

 

 

1

2

 

 

 

 

2

 

 

G

 

 

 

1

 

 

i ei 1 i 2

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

d 1d 2

 

 

 

 

 

 

 

.

 

 

 

2

 

2

2

 

 

 

 

 

2

2

 

 

2

 

 

 

 

1

2

 

 

 

 

2

 

 

10

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