Ряды Фурье: Теория и методы разложения

Ортогональность тригонометрической системы функций на интервале [a; a + 2π]

Добро пожаловать в курс «Ряды Фурье». Это первая и фундаментальная статья, на которой строится вся теория гармонического анализа. Прежде чем мы научимся раскладывать сложные сигналы и функции на простые составляющие, нам нужно понять, как эти «простые составляющие» взаимодействуют друг с другом.

Ключевым понятием здесь является ортогональность. Без этого свойства ряды Фурье были бы просто красивой, но бесполезной теорией. В этой статье мы разберем математическую природу тригонометрической системы и строго докажем её свойства.

Аналогия с векторами: что такое ортогональность?

Чтобы понять ортогональность функций, вспомним школьную геометрию и векторы. Два вектора и называются перпендикулярными (или ортогональными), если их скалярное произведение равно нулю:

где — скалярное произведение, а и — длины векторов.

В функциональном анализе функции рассматриваются как «векторы» в бесконечномерном пространстве. Аналогом скалярного произведения для двух функций и на интервале является определенный интеграл от их произведения.

Определение: Две функции и называются ортогональными на интервале , если интеграл от их произведения равен нулю:

где и — границы интервала интегрирования, и — исследуемые функции.

Согласно scask.ru, это свойство позволяет «выделять» отдельные гармоники из сложного сигнала, подобно тому как мы находим проекции вектора на оси координат.

Тригонометрическая система функций

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

Запишем эту систему:

Здесь: * — постоянная составляющая (можно представить как ). * — натуральное число (), указывающее на частоту колебания. * — независимая переменная.

Наша задача — доказать, что любые две различные функции из этого набора ортогональны друг другу на любом интервале длины . Обычно рассматривают интервал или , но в общем виде мы возьмем интервал .

Доказательство ортогональности

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

Нам нужно проверить три типа комбинаций функций:

Косинус на косинус (, где ).

Синус на синус (, где ).

Синус на косинус (, для любых и ).

Также отдельно стоит проверить ортогональность любой функции с единицей (первым элементом системы).

1. Ортогональность с единицей

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

где — целое число (), — начало интервала.

Подставим пределы интегрирования:

Так как синус имеет период , то добавление к аргументу не меняет значения функции: . Следовательно, числитель обращается в ноль:

Аналогично доказывается для синуса: . Это означает, что все синусы и косинусы ортогональны константе .

2. Произведение косинусов (при )

Рассмотрим интеграл от произведения двух разных косинусов. Используем формулу:

где и — различные натуральные числа.

Интеграл примет вид:

Так как , то и — это целые ненулевые числа. Как мы доказали в пункте 1, интеграл от косинуса (или синуса) целой кратности по полному периоду всегда равен нулю. Значит, оба слагаемых равны нулю.

Вывод: при .

3. Произведение синусов (при )

Используем формулу:

Логика абсолютно та же, что и для косинусов. Интеграл распадается на сумму двух интегралов от косинусов с целыми коэффициентами и . Оба интеграла равны нулю на интервале длиной .

Вывод: при .

4. Произведение синуса и косинуса

Здесь нам нужно рассмотреть произведение . Важно отметить, что здесь результат будет нулем даже если . Используем формулу:

Интегрируем:

* Если , то мы интегрируем синусы с целыми частотами по полному периоду. Результат — 0. * Если , то первое слагаемое превращается в , а второе — в , интеграл от которого по периоду тоже равен 0.

Вывод: Синусы и косинусы всегда ортогональны друг другу, независимо от их индексов и .

Норма функций (Квадрат нормы)

Мы выяснили, что интеграл от произведения разных функций системы равен нулю. А что будет, если умножить функцию саму на себя? Это называется квадратом нормы функции. Это значение понадобится нам в следующих статьях для вычисления коэффициентов ряда.

Для синусов и косинусов ()

Найдем значение интеграла для квадрата косинуса:

где — обозначение квадрата нормы.

Используем формулу понижения степени:

Подставим в интеграл:

Разберем слагаемые:

Первый интеграл: .

Второй интеграл: интеграл от косинуса по полному периоду равен 0.

Таким образом:

Аналогичный расчет для синуса через формулу дает тот же результат:

Для единицы

Для первой функции системы (константы 1) расчет самый простой:

Обратите внимание: квадрат нормы для единицы () отличается от квадрата нормы для синусов и косинусов (). Именно поэтому в формуле ряда Фурье первое слагаемое часто записывают как , чтобы компенсировать эту разницу («двойку») и использовать общие формулы.

Итоги

Мы строго доказали свойства, которые лежат в основе теории рядов Фурье. Краткое резюме:

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

Универсальность интервала: Результат верен для любого интервала , так как подынтегральные функции периодичны.

Нормировка:

* Интеграл от квадрата любого или равен . * Интеграл от квадрата единицы равен .

Смысл: Это свойство позволяет вычислять коэффициенты ряда Фурье по отдельности, не решая бесконечных систем уравнений.

Разложение функции на основном интервале в ряд Фурье и формулы коэффициентов

В предыдущей статье мы разобрали фундамент гармонического анализа — ортогональность тригонометрической системы функций. Мы выяснили, что синусы и косинусы с разными частотами «не видят» друг друга при интегрировании: их скалярное произведение равно нулю.

Теперь мы переходим к главной магии рядов Фурье. Мы научимся брать произвольную функцию (сигнал) и раскладывать её на сумму простейших гармонических колебаний. Это похоже на то, как призма раскладывает белый свет на спектр цветов, или как музыкант слышит в аккорде отдельные ноты.

Наша цель — получить формулы, которые позволят вычислить «вес» (амплитуду) каждой ноты в общем звучании функции. Эти веса называются коэффициентами Фурье.

Тригонометрический ряд Фурье

Предположим, что у нас есть периодическая функция с периодом , или же функция, заданная на интервале . Мы хотим представить её в виде бесконечной суммы:

Где: * — исходная функция. * — коэффициенты ряда (числа), которые нам предстоит найти. * — номер гармоники (), определяющий частоту колебания. * и — базисные функции тригонометрической системы.

Обратите внимание на первое слагаемое . Часто возникает вопрос: почему делим на 2? Это сделано для удобства, чтобы формула для выглядела так же, как и общая формула для , и нам не пришлось запоминать лишнее исключение. Математическая причина кроется в норме функции, которую мы вычисляли в прошлой статье (норма единицы была , а косинусов — ).

Вывод формул Эйлера-Фурье

Как найти неизвестные коэффициенты , если нам дана только ? Нам поможет свойство ортогональности. Идея метода состоит в том, чтобы «вырезать» из бесконечной суммы нужный нам коэффициент, обнулив все остальные слагаемые с помощью интегрирования.

1. Поиск коэффициента (Постоянная составляющая)

Проинтегрируем обе части уравнения ряда по интервалу :

Где: * — определенный интеграл по периоду.

Вспомним свойства из прошлой статьи:

Интеграл от или по полному периоду равен 0. Значит, вся огромная сумма превращается в ноль.

Остается только первое слагаемое.

Отсюда выражаем :

Это формула для вычисления среднего значения функции на интервале (с поправкой на коэффициент 2 в ряде).

2. Поиск коэффициентов (при косинусах)

Чтобы найти конкретный коэффициент (где — фиксированный номер), умножим обе части равенства ряда на и снова проинтегрируем по .

Разберем, что происходит слагаемыми справа:

Первое слагаемое: (интеграл от косинуса по периоду).

Слагаемые с синусами: для любых и (ортогональность синусов и косинусов).

Слагаемые с косинусами: . Согласно свойству ортогональности, этот интеграл равен 0 во всех случаях, кроме одного: когда . В этом случае интеграл равен .

Таким образом, из бесконечной суммы выживает только одно слагаемое, где :

Выражаем (заменив индекс обратно на для общности):

3. Поиск коэффициентов (при синусах)

Аналогично умножаем уравнение ряда на и интегрируем. Все слагаемые с косинусами и свободный член обнуляются. Остается только интеграл от квадрата синуса при , который равен .

Результат:

Итоговые формулы коэффициентов

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

Согласно Wikipedia, для функции , интегрируемой на промежутке , коэффициенты определяются так:

Где: * — искомые коэффициенты. * — разлагаемая функция. * — индекс суммирования (номер гармоники).

Пример разложения: «Прямоугольная волна»

Лучший способ понять формулы — применить их на практике. Рассмотрим простую, но важную для электроники функцию — прямоугольный импульс (или функцию знака), заданную на интервале .

Пусть функция задана следующим образом:

Это «ступенька»: слева от нуля она равна -1, справа равна 1.

Шаг 1: Вычисляем

Вычисляем интегралы: * Первая часть: . * Вторая часть: . * Сумма: .

Это логично: среднее значение функции, которая половину времени равна -1, а половину +1, равно нулю.

Шаг 2: Вычисляем

Косинус — функция четная (симметричная), а наша — нечетная (антисимметричная). Произведение нечетной и четной функции дает нечетную функцию. Интеграл от нечетной функции по симметричному интервалу всегда равен нулю.

Мы могли бы это посчитать честно через первообразные, но использование свойств симметрии экономит время. В ряде Фурье для нечетной функции будут только синусы.

Шаг 3: Вычисляем

Здесь нечетная и нечетная. Их произведение — функция четная. Значит, интеграл по левой половине равен интегралу по правой. Можно просто удвоить интеграл по правой части:

Интегрируем синус:

Подставляем пределы:

Упрощаем выражение в скобках: * . * (это чередование: -1 при , 1 при и т.д.).

Проанализируем результат: * Если — четное (), то , и скобка . Все четные коэффициенты равны нулю. * Если — нечетное (), то , и скобка . Тогда .

Итоговый ряд

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

Или в развернутом виде:

Согласно mathprofi.ru, этот результат показывает, как прямоугольный сигнал можно собрать из суммы синусоид. Первая гармоника () задает основной тон, а добавление третьей, пятой и последующих гармоник делает края волны всё более крутыми, приближая форму к «ступеньке».

Итоги

Суть метода: Любую периодическую (или заданную на конечном отрезке) функцию можно представить как сумму синусов и косинусов с разными амплитудами.

Роль интеграла: Интеграл выступает как «фильтр», который благодаря ортогональности выделяет из функции амплитуду конкретной частоты.

Формулы: Коэффициенты и находятся путем интегрирования произведения исходной функции на и соответственно, деленного на .

Симметрия: Если функция нечетная (симметрична относительно начала координат), её ряд состоит только из синусов (). Если четная (симметрична относительно оси Y) — только из косинусов ().