bilim.techО проекте

Ряд Фурье

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

01 Что это такое

Ряд Фурье раскладывает любой периодический сигнал на синусоиды кратных частот.

Вибрация вала, ток в сети, звук мотора повторяются с некоторым периодом, но по форме бывают далеко не синусоидой. Оказывается, любую такую форму можно сложить из простых синусоид: одна с основной частотой, другие в 2, 3, 4 и больше раз чаще. У каждой своя высота. Набор этих высот называют спектром: он показывает, сколько каждой гармоники сидит в сигнале.

сигнал во временивремя tмеандрсумма трех синусоидспектр: высота каждой гармоники1234567номер гармоники k
Три синусоиды с частотами 1, 3 и 5 уже дают почти прямоугольный сигнал. Справа те же три слагаемые в виде спектра: чем выше столбик, тем сильнее гармоника.
k
Номер гармоники. Во сколько раз гармоника чаще основной частоты. Первая гармоника идет с частотой самого сигнала, третья в 3 раза чаще.
bₖ
Амплитуда гармоники. Высота k-й синусоиды в сумме. Все bₖ вместе и есть спектр сигнала.
K
Сколько гармоник взяли. Сумма первых K слагаемых называется частичной суммой. Чем больше K, тем точнее она повторяет сигнал.
02 Опыт
Добавляй гармоники по одной. Смотри, как из гладких волн получается угол. РАСЧЕТ В РЕАЛЬНОМ ВРЕМЕНИ
РЯДЫ ФУРЬЕ

Из каких синусоид сложен сигнал?

Рис. 6.3
СИГНАЛ
Ненулевых слагаемых
Первая гармоника, b₁
Ошибка суммы

ⓘ Коэффициенты bₖ для меандра, пилы и треугольника найдены по формуле с интегралом заранее и записаны готовыми выражениями. Сигнал вибрации здесь задан самой суммой гармоник, поэтому ряд для него точный уже при K = 11.

03 Формулы по буквам

Сумма синусоид, меандр и коэффициенты

Первая формула говорит, как сигнал собирается из гармоник. Вторая показывает ряд для меандра. Третья объясняет, откуда берется высота каждой гармоники. Под каждой формулой числа из опыта в выбранной точке t.

Из чего состоит сигнал?

f(t) = Σk=1∞ bₖ · sin(k·ω·t)

Берем синусоиду с частотой сигнала, потом в 2 раза чаще, в 3 раза чаще и так далее. Каждую умножаем на свою амплитуду bₖ и складываем. На практике берут конечное число слагаемых K.

k
номер гармоники: 1, 2, 3 и далее
bₖ
амплитуда k-й гармоники, высота столбика в спектре
ω
основная угловая частота: ω = 2π/T, где T период сигнала

Как сложить прямоугольник из волн?

4π(sin ωt + sin 3ωt3 + sin 5ωt5 + …)

Меандр прыгает между +1 и −1. В его ряду только нечетные гармоники, а их высота падает как 1/k. Поэтому ряд сходится медленно: острый угол требует очень высоких частот.

4/π
амплитуда первой гармоники, около 1,27: она выше самого меандра
1/k
каждая следующая гармоника слабее во столько раз, какой у нее номер
…
дальше sin 7ωt / 7, sin 9ωt / 9 и так без конца

Как найти высоту гармоники?

bₖ = 2T ∫0T f(t) · sin(kωt) dt

Умножаем сигнал на пробную синусоиду k-й частоты и усредняем за период. Если сигнал похож на эту синусоиду, произведение почти всегда положительно и среднее большое. Если не похож, плюсы и минусы гасят друг друга.

T
период сигнала
∫0T … dt
интеграл за период: площадь под произведением
2/T
делим на период, получаем среднее, и умножаем на 2, потому что среднее sin² равно 1/2

Почему гармоники не мешают друг другу

Перемножим две синусоиды разных частот, например sin t и sin 3t, и посмотрим на среднее за период. Произведение то положительное, то отрицательное, и площади плюсов и минусов в точности равны. Среднее равно нулю. А произведение синусоиды самой на себя, sin² t, всегда не меньше нуля, и его среднее равно 1/2.

Поэтому, когда мы умножаем сигнал на sin(kωt) и усредняем, все гармоники с другими номерами дают ноль, и остается только bₖ · 1/2. Отсюда множитель 2 в формуле. Каждую гармонику можно найти отдельно, не зная остальных.

У меандра только нечетные гармоники из-за симметрии: вторая половина периода повторяет первую с обратным знаком. Синусоида четного номера за полпериода делает целое число колебаний и в обеих половинах идет одинаково. Сигнал в половинах противоположный, вклады половин взаимно уничтожаются, и bₖ для четных k равен нулю.

среднее sin kt · sin mt
равно 0, если k ≠ m: разные гармоники не похожи друг на друга
среднее sin² kt
равно 1/2: гармоника полностью похожа на себя
f(t + T/2) = −f(t)
симметрия меандра и треугольника: остаются только нечетные k
bₖ ~ 1/k
у сигнала со скачком, как меандр и пила
bₖ ~ 1/k²
у сигнала с изломом без скачка, как треугольник: сходится быстрее
sin xsin 3x / 3sin 5x / 5сложим по вертикали:сумма почти прямоугольная
Каждая следующая волна чаще и ниже. Вместе они выравнивают верх и делают крутые фронты.
04 Где и зачем это применяют

Везде, где сигнал повторяется

Допущения модели+

В опыте все сигналы нечетные: f(−t) = −f(t), поэтому в ряду только синусы. В общем случае к ним добавляются косинусы с коэффициентами aₖ и постоянная составляющая a₀/2, а спектр показывает полную амплитуду √(aₖ² + bₖ²). Сигнал считается строго периодическим и бесконечным. Реальный прибор берет конечный кусок записи и считает спектр быстрым преобразованием Фурье по отсчетам, поэтому частоты в спектре идут с шагом 1/(длительность записи), а частота отсчетов должна быть хотя бы вдвое выше самой высокой гармоники. Выброс около 9% у скачка (эффект Гиббса) не исчезает при любом конечном K, он лишь становится уже.

Еще в направлении «Математика»

Все темы направления «Математика» →