Математика. Двумерное преобразование Фурье и его приложения. Учебное пособие
.pdfПереходим к полярным координатам. Для точек окружности имеем
R cos , |
R sin ; |
|
|
|
|||||||||
d R sin d , |
|
|
d R cos d ; |
||||||||||
g1 , |
|
|
1 |
|
sin |
|
|
|
|||||
|
|
|
|
|
|
|
2 |
|
; |
||||
|
4 R sin |
|
|||||||||||
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
g2 , |
|
1 |
|
|
|
cos |
|
|
|||||
|
|
|
|
|
2 |
|
. |
||||||
|
4 R sin |
|
|
||||||||||
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
Внося эти значения в уравнения (2.48), (2.49), после надлежащих алгебраических преобразований, приходим к такой паре преобразований Гилберта на окружности:
|
|
|
|
|
1 2 |
|
|
|
|
|
|||||||
u R, A |
|
|
v R, ctg |
|
|
2 |
d ; |
(2.67) |
|||||||||
2 |
|
|
|||||||||||||||
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
||
|
|
|
|
|
1 2 |
|
|
|
|
|
|||||||
v R, B |
|
u R, ctg |
|
|
|
d , |
(2.68) |
||||||||||
2 |
|
|
2 |
||||||||||||||
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
||
где |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
A |
1 |
2 u |
|
R, |
|
d , |
B |
1 |
2 v |
|
R, d . |
(2.69) |
|||||
|
|
||||||||||||||||
|
2 |
|
|
|
|
|
|
2 |
|
|
|
|
|||||
|
0 |
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|||
41
ГЛАВА 3
3.1. Е- и H-волны в слоисто-неоднородной среде
Рассмотрим двумерные преобразования Фурье, сводящие систему уравнений Максвелла в произвольной слоистонеоднородной среде к двум независимым системам обыкновенных дифференциальных уравнений, описывающих распространение E- и H-волн.
Пусть электромагнитные параметры среды произвольным образом изменяются вдоль некоторого фиксированного направления, задаваемого осью z прямоугольной системы координат. В плоскостях, перпендикулярных к этому направлению, среда неограниченна и однородна.
Рассмотрим комплексную форму системы уравнений Максвелла для монохроматических колебаний:
rot |
H |
i a (z) |
E |
, |
rot |
E |
i a (z) |
H |
. |
(3.1) |
По поперечным координатам х, y векторы напряженности электрического и магнитного полей представим в виде двумерных интегралов Фурье:
|
|
|
|
|
|
|
|
|||
|
|
(x,y,z) |
E |
( 1, 2,z)ei 1x i 2yd 1d 2; |
|
|||||
E |
(3.2) |
|||||||||
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
i 1x i 2y |
|
|
|||
|
|
|
|
|
|
,z)e |
d 1d 2. |
|
||
|
|
|
|
|
||||||
H(x,y,z) H( 1, 2 |
|
|
||||||||
Подставляя данные разложения в систему уравнений Максвелла (3.1) и приравнивая векторы спектральных плотностей в левых и правых частях равенств, получаем для них следующую векторную форму системы обыкновенных дифференциальных уравнений:
42
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
k |
H |
|
|
|
|
|
|
|
|||||||||||||
|
z |
i H i a |
(z)E; |
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(3.3) |
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
E |
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
k |
|
|
|
|
|||||||||||||||||
|
|
|
|
i E i a |
(z)H, |
||||||||||||||||
|
|
|
|
|
z |
|
|
|
|
|
|
|
|
|
|
|
|
||||
где i , j , k – единичные базисные векторы прямоугольной системы координат, а 1 i 2 j .
Введем новую правую тройку базисных векторов с помощью равенств
e1 , e2 k , e3 k
и разложим векторы E и H по этому базису:
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
; |
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
E E1 e1 |
E2 e2 |
E3 e3 |
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
H |
H1 e1 |
H2 e2 |
H3 e3. |
|
|
|
|
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
, е |
, е |
|||||
Переход от компонент векторов E и H в базисе е |
||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
2 |
3 |
||||
к их компонентам в исходном базисе i , j , k осуществляется по формулам:
|
1 |
|
|
2 |
|
|
Ex |
|
E1 |
|
|
E2 |
; |
|
2 |
|
|
1 |
|
|
Ey |
|
E1 |
|
|
E2 |
; |
|
|
|
|
|
|
|
|
Ez E3. |
|
|
|||
Проектируем уравнения системы (3.3) на векторы е1 , е2 , е3 . При этом наша система распадается на две независимые системы из трех уравнений:
|
|
|
|
|
|
H2 |
|
|
|
||
|
z |
i a E1; |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
E1 |
|
; |
(3.4) |
||
z |
i E3 |
i a H2 |
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
i H2 i a E3. |
|
|
||
43
|
|
|
|
|
|
E2 |
|
|
|
||
|
z |
i a H1; |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(3.5) |
H1 |
|
; |
|||
z |
i H3 |
i a E2 |
|
||
|
|
|
|
|
|
|
|
|
|
|
|
i E2 |
i a H3. |
|
|
||
Электромагнитное поле, представляющее собой результат решения системы (3.4), называется E-волной, что подчеркивает наличие в нем продольной компоненты электрического поля ( E3 = Ez 0) и отсутствие таковой компоненты в магнитном поле ( H3 = Hz = 0). Соответственно, поле, описываемое системой уравнений (3.5), называется H-волной, так как в ней Ez = 0, Hz 0.
Поскольку продольные составляющие E3 и H3 не дифференцируются, их удобно исключить из рассматриваемых уравнений. В результате получаются системы из двух уравнений, по форме полностью совпадающие с соответствующими системами для плоских волн. А именно, для E-волны
|
|
|
|
|
|
|
|
|
|
H2 |
|
|
|
|
|
|
|
||
z |
|
i a E1; |
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
E1 |
|
|
|
|
|
|
; |
(3.6) |
|
z |
i a H2 |
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
E3 |
|
|
H2. |
|
|
|
|||
a |
|
|
|
||||||
Для H-волны |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
E2 |
|
|
|
|
|
|
|
||
z |
|
i a H1; |
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
H1 |
|
|
|
|
; |
(3.7) |
|||
z |
|
i a E2 |
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
H3 |
|
E2. |
|
|
|
||||
a |
|
|
|
||||||
44
Здесь были использованы такие обозначения:
1 |
2 |
, |
k2 2 a a . |
|
k2 |
||||
|
|
|
Внося спектральные плотности, отвечающие уравнениям (3.6), (3.7), в интегралы Фурье (3.2), получим следующие представления векторов электромагнитных полей, физически существующих независимо друг от друга E- и H-волн.
Поле E-волны
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
i 1x |
i 2y |
d 1d 2; |
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
Eэ E1 e1 |
E3 e3 e |
|
|
|
|
|
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(3.8) |
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
i 1x i 2y |
d 1d 2. |
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
Hэ H2 e2e |
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
Поле H-волны |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
e |
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
i 1x i 2y |
d 1d 2 |
; |
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
Hм H1 e1 |
H3 e3 |
|
|
|
|
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(3.9) |
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
i 1x i 2y |
d 1d 2. |
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
Eм E2 e2e |
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
Спектральные плотности в |
этих интегралах зависят |
||||||||||||||||||||||||||||
от волнового параметра |
|
|
2 2 . По отдельности волно- |
||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
2 |
|
|
|
|
|
|||
вые числа 1, 2 входят только в определение базисных век-
торов e1 и e2 . Поэтому в полярной системе координат интегрирование по углу может быть произведено в явном виде. Введем на плоскости волновых чисел 1, 2 и на плоскости физических переменных х, y полярные системы координат по формулам:
1 cos , |
2 sin ; |
x cos , |
y sin . |
45
Тогда
e1 cos i sin j; e2 sin i cos j ,
и, внося полярные координаты в интегралы (3.8), (3.9), приводим их к виду
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
d , |
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
Eэ E1T E3 N |
Hэ k H2 T d ; |
|||||||||||||||||||||||||||||
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
Eм |
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
Hм H1T H3 N d , |
|
|
k E2 T d , |
||||||||||||||||||||||||||||
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|||||||||||||
где |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|||||||
|
|
(cos i sin |
|
)ei cos( )d , |
|
|
|
|
ei cos( )d . |
|||||||||||||||||||||||||
T |
j |
N |
k |
|||||||||||||||||||||||||||||||
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
||||||||||
Функции Бесселя нулевого и первого порядка имеют та- |
||||||||||||||||||||||||||||||||||
кие представления: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
J (z) |
1 |
eiz costdt, |
|
J (z) |
1 |
|
cos t eiz costdt, |
|||||||||||||||||||||||||
|
|
|
|
i |
|
|||||||||||||||||||||||||||||
0 |
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
||
посредством которых, как нетрудно показать, введенные векторы T и N выражаются так:
|
|
|
T |
2 iJ1 ( ) |
e |
, |
N |
|
2 J0 ( ) |
k |
, |
|
|
|||||||
где |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
cos i sin |
|
, |
|
|
|
|
|
|
sin i cos |
|
. |
|||||||
e |
j |
e |
k |
e |
j |
|||||||||||||||
Таким образом, поля независимых E- и H-волн в слоистонеоднородной среде представляются однократными интегралами типа Фурье–Бесселя.
46
Поле E-волны:
Eэ 2 i E1 ( ,z)J1 ( ) d e 2 E3 ( ,z)J0 ( ) d k;
0 |
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
( ,z)J1 ( ) d e . |
|
|
|
|||||||||
|
|
|
Hэ 2 i H2 |
|
|
|
|||||||||||
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
||
Поле H-волны: |
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
( ,z)J1 ( ) d e |
( ,z)J0 |
( ) d k; |
|||||||||||||
Hм 2 i H1 |
2 H3 |
||||||||||||||||
0 |
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
( ,z)J1 |
( ) d e . |
|
|
|
||||||||
|
|
|
Eм 2 i E2 |
|
|
|
|||||||||||
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
||
47
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1.Анго А. Математика для электро- и радиоинженеров. М.: Наука,1965.
2.Кошляков Н.С., Глинер Э.Б., Смирнов М.М. Уравнения в частных производных математической физики. М.: Высшая школа. 1970.
3.Романовский П.И. Ряды Фурье. Теория поля. Аналитические и специальные функции. Преобразование Лапласа. М.: Наука, 1973.
4.Треногин В.А., Недосекина И.С. Методы математической физики. М.: Изд. Дом НИТУ «МИСиС», 2012.
5.Цапенко Н.Е. Аналитические функции и интегральные преобразования. М.: Физматлит, 2012.
6.Цапенко Н.Е. Интеграл Фурье и его приложения. М.: Изд. Дом НИТУ «МИСиС», 2021.
48
Учебное издание
Цапенко Николай Евгеньевич
МАТЕМАТИКА
Двумерное преобразование Фурье и его приложения
Учебное пособие
Корректор В.В. Демидова
Технический редактор Т.В. Суханова Верстальщик Ю.Б. Пашкова
Подписано в печать 30.05.23 |
Уч.-изд. л. 3,06 |
Формат 60 90 1/16
Издательский Дом НИТУ МИСИС, 119049, Москва, Ленинский пр-кт, д. 4, стр. 1
Издательский Дом НИТУ МИСИС, 119049, Москва, Ленинский пр-кт, д. 2А Тел. 8 (495) 638-44-06
Отпечатано в типографии Издательского Дома НИТУ МИСИС,
119049, Москва, Ленинский пр-кт, д. 4А Тел. 8 (495) 638-44-16, 8 (495) 638-44-43
