- •Введение
- •Глава 1. Элементы векторного анализа
- •1. Векторы и действия над ними
- •2. Математическое понятие поля. Градиент
- •3. Дивергенция. Теорема Остроградского-Гаусса
- •4. Ротор. Теорема Стокса
- •5. Некоторые соотношения векторного анализа
- •6. Операции в криволинейных координатах
- •В цилиндрических координатах
- •В сферических координатах
- •7. О дифференциальных уравнениях с частными производными
- •Глава 2. Уравнения лапласа и пуассона
- •8. Дельта-функция Дирака
- •9. Интегрирование уравнения Пуассона
- •10. Граничные задачи для уравнения Лапласа
- •11. Метод разделения переменных
- •Г лава 3. Гармонические колебания и волны
- •12. Гармонические колебания и метод комплексных амплитуд
- •13. Волновые процессы и их математическое описание
- •14. Вращение декартовой системы координат
- •Глава 4. Решения волновых уравнений
- •15. Интегрирование неоднородного уравнения Гельмгольца и уравнения Даламбера
- •16. Уравнение Бесселя и цилиндрические функции
- •17. Решение однородного уравнения Гельмгольца методом разделения переменных
- •Глава 5. Краевые задачи электродинамики
- •18. Граничные задачи для уравнения Гельмгольца. Собственные функции и собственные значения
- •19. Ортогональные системы функций и ряды Фурье
- •20. Сведения из алгебры
- •21. Проекционные методы
- •Список литературы
- •Контрольные задания
19. Ортогональные системы функций и ряды Фурье
19.1. Некоторые свойства оператора Лапласа. Рассмотренные выше граничные задачи (18.1) и (18.11), а также задачи (18.41), (18,42) и ряд других можно кратко сформулировать в виде
(19.1)
где 2 - взятый с обратным знаком оператор Лапласа, параметр к - это k2 или χ 2 (в. зависимости от того, является ли задача трёхмерной или двумерной), а – векторные или скалярные ( = и) функции, подчинённые требуемым граничным условиям. Обычно граничные условия задаются вместе с оператором L; они определяют класс функций и, на которые может распространяться действие L, т. е. область определения оператора DL. Разумеется, функции , принадлежащие этому классу (кратко пишут: Dl ), должны допускать заданные операции дифференцирования. Говорят, что запись (19.1) выражает общую формулировку задачи на собственные значения оператора L. В частности если и = 0 на S для Dl, получаем задачу (18.1), если иτ = 0 и div = 0 на S для Dl - задачу (18.41), и т. д.
Введём интегральную величину, имеющую своё обозначение:
(19.2)
Здесь подразумевается область V трёхмерной задачи, а в случае задачи двумерной интегрирование производится по её области S; величина называется скалярным произведением функций , последние могут быть как векторами, так и скалярами. Заметим, что (v, и) = (и, v) *.
Оператор L, обладающий тем свойством, что для
(19.3)
является по определению симметрическим. Легко убедиться, что для всех поставленных в п.18 задач оператор 2 симметрический.
Действительно, в случае скалярных uиv на основании первой формулы Грина (5.13)
(19.4)
Как для первой граничной задачи (18.1) (и = 0, v = 0 на S), так и для второй задачи (18.11) поверхностный интеграл в (19.4) равен нулю. Меняя местами u и v*, мы оставляем правую часть (19.4) неизменной, а потому
(19.4а)
т. е. равенство (19.3) выполняется.
Взяв векторные , получаем:
(19.5)
(были использованы формулы векторного анализа (5.12), (5.10), (5.9) и (3.8)).
Если функции подчинены граничным условиям задачи (18.41) (uτ = 0, div = 0 и vτ = 0, div = 0 на S) или граничным условиям задачи (18.42) (uτ = 0, (rot )τ = 0 и vτ = 0, (rot )τ = 0 на S), то поверхностные интегралы в (19.5) равны нулю, и правая часть остается неизменной при замене и обратно. Поэтому
, (19.5a)
и выполняется равенство (19.3).
Установленное свойство (19.3) сохраняется и для двумерных областей, т. е. при замене V на S, а граничной поверхности S объема V на контур L поверхности S. Читателю рекомендуется проверить это в качестве упражнения.
Нетрудно убедиться, что собственные значения задачи (19.1) при симметрическом операторе 2 вещественны и неотрицательны. Образуя в (19.1) слева и справа скалярные произведения с , имеем:
(19.6)
Знаменатель положителен. В силу (19.3) , и в то же время *. Отсюда следует, что числитель - величина вещественная. Следовательно, вещественны и собственные значения
Рассматривая отдельно скалярные и векторные задачи, числитель в (19.6) преобразуем при помощи соотношений (19.4а) и (19.5а). Мы видим, что
(19.7)
и
(19 8)
Итак, собственные значения κ = κi задачи (19.1) неотрицательны и могут быть расположены в следующем порядке:
(19.9)
19.2. Ортогональные системы функций. Две функции и и v называются oртогональными, если их скалярное произведение равно нулю:
. (19.10)
Пусть (i = 1, 2, ...) - собственные функции задачи (19.1), которым соответствуют собственные значения κ; предположим, что все κi различны, или, как говорят, отсутствует вырождение. Взяв две любые собственные функции , перепишем формулировку (19.1):
Образуя скалярные произведения ( , ) и ( , ) и принимая во внимание вещественность собственных значений, получаем:
или, в силу симметричности L
= 0, (19.11)
а отсюда следует, что , т. е собственные функции иi ортогональны.
Задача на собственные значения (19.1) порождает, таким образом, ортогональные системы функций .
Ортогональная система всегда может быть нормирована, т. е. можно так подобрать постоянные коэффициенты в выражениях , что для всех i. Тогда получается ортонормированная система, для которой
(19.12)
где δik= 0 при i ≠ к и δik = 1 при i = k (символ Кронекера). Путём проверки нетрудно убедиться, что различные собственные функции, полученные в п.18, образуют ортогональные системы. Возьмём, например, систему функций итпр (18.14):
.
Скалярное произведение двух функций итпр и имеет вид:
(i и k надо понимать как совокупность чисел т, п, р и т', п', р' соответственно). Поскольку
и таковы же интегралы по у и z, то, как видно, скалярное произведение действительно обращается в нуль при i ≠ k (т. е. т ≠ т', п ≠ п' и р ≠ р'). Отсюда же следует, что система будет ортонормированной, если взять
при т ≠ 0, п ≠ 0, р ≠ 0. Когда среди чисел т, п и р имеются нули, и0 столько раз делится на , каково число нулей (один, два или три).
Предлагаю при помощи формул (16.24) и (16.26) проверить ортогональность и ввести нормировку функций ипт (18.29). При этом следует отдельно рассматривать функции с С ≠ 0, D = 0 и С = 0, D ≠ 0 (либо Q ≠ 0, Т = 0 и Q = 0,T ≠ 0).
19.3. Ортогональные ряды. Взяв ортонормированную систему функций и некоторую функцию , определённую в той же области, построим ряд
(19.13)
Он называется ортогональным рядом, или рядом Фурье функции , а ап - коэффициентами Фурье.
Отличительным свойством ряда Фурье является выполнение равенства
(19.14)
Действительно, составляя в (19.13) скалярное произведение с , справа получаем нуль во всех членах кроме k-го, который дает аk; равенство (19.13) превращается в выражение коэффициента Фурье аk.
Говорят, что ряд Фурье сходится в среднем к , а система {ип} полна (в этом смысле), если
0 (19.15)
Системы собственных функций оператора Лапласа обладают указанным свойством полноты при всевозможных функциях которые приходится рассматривать практически (например, с разрывами второго рода).
Легко убедиться, что обычные тригонометрические ряды Фурье дают частный пример ряда (19.13).
Возьмём, например, систему собственных функций одномерного оператора Лапласа из п.7.2. После нормировки в (7.9) , и система принимает вид:
По этой системе разложим некоторую функцию f(x), определённую на отрезке 0 ≤ х ≤ а. Согласно (19.13)
а это и есть обычный тригонометрический ряд Фурье по синусам.
В качестве второго примера рассмотрим ряд Фурье типа (12.22). В данном случае разлагается определённая на отрезке –T/2 ≤ t ≤ T/2 функция , по ортонормированной системе
Ряд (19.13) имеет вид
что совпадает с (12.22), (12.23). Заметим, что - это собственные функции одномерного оператора Лапласа при периодических граничных условиях и (-Т/2) = и(Т/2).