Ряд фурье и волновое уравнение (6 видео)

Содержание

Лекция 3. Метод Фурье
Общая схема метода Фурье.
Основная лемма метода Фурье.
Решение первой начально-краевой задачи для волнового уравнения.
Метод Фурье
Метод Фурье
Первое из граничных условий
Практическое применение преобразования Фурье для анализа сигналов. Введение для начинающих
1. Преобразование Фурье и спектр сигнала
2. Непрерывная функция и представление её рядом Фурье
3. Дискретные сигналы и дискретное преобразование Фурье
Некоторые итоги
🔥 Видео

Видео:Метод Фурье для волнового уравненияСкачать

Лекция 3. Метод Фурье

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

Общая схема метода Фурье.

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

При применении метода Фурье удобно использовать следующую лемму.

Основная лемма метода Фурье.

Если в прямоугольнике R плоскости XOY:

для некоторых функций выполняется тождество

то в этом случае

Доказательство. Предположим противное, т.е. что

Тогда существуют значения такие, что

Рассмотрим точки (x₁,y) и (x₂,y), принадлежащие прямоугольнику R. На R справедливо тождество (8), а поэтому

Сравнивая эти равенства, приходим к противоречию с нашим предположением. Следовательно X(x) = const, а тогда Y(y)=const.

Решение первой начально-краевой задачи для волнового уравнения.

Рассмотрим волновое уравнение

Граничные условия первого рода

И начальные условия

Решим эту задачу методом Фурье.

Шаг 1. Представим функцию U(x,t) в виде

Найдем частные производные U_xx и U_tt и подставим в уравнение (9):

В полученном уравнении левая часть зависит только от x, а правая- только от t. Используя основную лемму, заключаем:

Из граничных условий (10) получим

Шаг 2. Решим задачу Штурма-Лиувилля

Она имеет собственные значения и собственные функции

Шаг 3. Подставим найденные значения λ_n в уравнение а) и решим его:

Шаг 4. Выпишем частные решения уравнения (9):

Для волнового уравнения эти решения называются собственными колебаниями. В лекции 6 мы изучим их подробнее. В силу линейности и однородности уравнения (9) линейная комбинация этих решений

Замечание 1. Здесь мы предполагаем, что полученный функциональный ряд равномерно сходится и его можно дважды почленно дифференцировать по x и по t в области 0 0. Об условиях, при которых это можно сделать, будет рассказано в лекции 5.

Шаг 5. Определим коэффициенты A_nи B_n в формуле (12), используя начальные условия (11). Из первого начального условия получим

Равенство (13) означает, что начальная функция φ(x) разлагается в ряд Фурье по синусам, которые в данном случае являются собственными функциями X_n(x) задачи Штурма-Лиувилля.

Коэффициенты Фурье вычисляются по формулам

Из второго начального условия находятся коэффициенты B_n.

Вычислив коэффициенты A_n и B_n для конкретных начальных функций и подставив их значения в (12), мы получим решение первой начально-краевой задачи.

Замечание 2. Используя формулу (12), можно получить решение первой начально-краевой задачи для уравнения колебания струны: Для этого проведем замену переменной τ=at и получим

При этом начальное условие не изменится, а условие преобразуется к виду Тогда решение задачи в переменных (x,τ) будет иметь вид

Возвращаясь к переменным (x,t), получим

Видео:5. Решение волнового уравнения на отрезке методом ФурьеСкачать

Метод Фурье

Содержание:

По этой ссылке вы найдёте полный курс лекций по математике:

Метод Фурье, или метод разделения переменных, является одним из наиболее распространенных методов решения уравнений с частными производными. Рассмотрим этот метод, обратившись к простейшей задаче о свободных колебаниях однородной струны длины i, закрепленной на концах. §4. Свободные колебания однородной струны, закрепленной на концах Задача о свободных колебаниях однородной струны с закрепленными концами сводится к решению уравнения при граничных условиях и начальных условиях.

Метод Фурье Задачу (1 )-(3) называют смешанной: она содержит и начальные и граничные условия. Решение задачи начнем с поиска частных решений уравнения (1) вида При этом будем предполагать, что каждое из них удовлетворяет граничным условиям (2), но не равно нулю тождественно. Подставляя функцию и<х, t) в форме (4) в уравнение (1), получаем ИЛИ Последнее равенство (его левая часть зависит только от а правая — только от х) возможнолишь втом случае, если обе его части не зависят ни от ty ни от х,т.е. равны одной и той же постоянной.

Обозначим эту постоянную (разделения) через (-А), Из равенства (5) получаем два обыкновенных дифференциальных уравнения Граничные условия (2) дают откуда (T(t) £ 0) следует, что функция Х(х) должна удовлетворять граничным условиям Чтобы получить нетривиальные решения tt(x, t) вида (4), удовлетворяющие граничным условиям (2), необходимо найти нетривиальные решения уравнения удовлетворяющие граничным условиям.

Таким образом, мы приходим к следующей задаче: найти значения параметра А, при которых существуют нетривиальные решения задачи (7)-(8), а также сами эти решения. Такие значения параметра А называются собственными значениями, а соответствующие им нетривиальные решения — собственными функциями задачи (7)-(8). Сформулированную таким образом задачу называют задачей Штурма—Лиувилля. Найдем собственные значения и собственные функции задачи (7)-(8).

Рассмотрим отдельно три случая, когда 1.

При общее решение уравнения (7) имеет вид Потребовав выполнения граничных условий (8), получим (6) (7) Так как определитель системы (9) отличен от нуля, то . Следовательно, Х(х) = 0, т. е. при нетривиальных решений задачи не существует. (9) 2. При А = 0 общее решение уравнения (7) имеет вид Граничные условия (8) дают откуда С, = С2 = 0, и следовательно, при А = 0 нетривиальных решений задачи (7)-(8) также не существует. 3.

При Л > 0 общее решение уравнения (7) имеет вид Потребовав выполнение граничных условий (8), получим Система (10) имеет нетривиальные решениятогда и толькотогда, когда определитель системы равен нулю, Метод Фурье будут собственными функциями задачи. Собственные функции определены с точностью до постоянного множителя, который мы выбрали равным единице. При А = А* общее решение у равнения (6) имеетвид ктга кчга где Аки Bk — произвольные постоянные. Таким образом, функции удовлетворяют уравнению (1) и граничным условиям (2) при любых Ак и Вку В силу линейности и однородности уравнения (1) всякая коневая сумма решений будет также решением уравнения (1).

То же справедливо и для ряда если он сходится равномерно и его можно дважды почленно дифференцировать по х и по t. Поскольку каждое слагаемое в ряде (11) удовлетворяет граничным условиям (2), то этим условиям будет удовлетворять и сумма u(s, t) этого ряда. Остается определить в формуле (11) постоянные .4* и Вк так, чтобы выполнялись и начальные условия (3). Продифференцируем формально ряд (11) по t.

Имеем Полагая в соотношениях (l 1) и (12) t = 0, в силу начальных условий (3) получим Формулы (13) представляют собой разложения заданных функций вряд Фурье по синусам в интервале Коэффициенты разложений (13) вычисляются по известным формулам / I Теорема 2. Если и удоъчетворяет условиям и удовлетворяет условию то сумма tx(x, £) ряда (11), где -А* и В* опредыяются формулами (14), имеет в области непрерывные частные производные до второго порядка включительно по каждому из аргументов, удовлетворяет уравнению (1), граничным условиям (2) и начальным условиям (3), т. е. является решением задачи (1 )-(3).

Пример. Найти закон свободных колебаний однородной струны длины I, закрепленной на концах, если в начальный момент t = 0 струна имеет форму параболы — const), а начальная скорость отсутствует. 4 Задача сводится к решению уравнения при граничных условиях и начальных условиях.

Возможно вам будут полезны данные страницы:

Метод Фурье

Применяя метод Фурье, ищем нетривиальные решения уравнения (1), удовлетворяющие граничным условиям (2), в виде Подставляя «(*,*) в форме (4) в уравнение (1) и разделяя переменные, получим откуда причем в силу (2) Как было установлю но выше, собственные значения задачи (7)-(8) а соответствующие собственные функции Для А = Ащ общее решение уравнения (6) имеет вид пяа ижа Будем иска тъ решение исходной задачи в виде ряда Для определен ия коэффициентов -4Я и Z?„ воспользуемся начальными условия ми (3).

Имеем Из формулы (II) срезу

получаем, что 2?„ = 0 для любог о п, а из (10) Метод Фурье откуда, интегрируя по частям дважды . находи м . Подставляя наеденные значения А, и в ряд (9), получим решение поставленной задачи , Замечание. Если начальные фукхдда не удовлетворяют условиям теоремы 2, то дважды непрерывно дифференцируемого решения смешанной задачи (1)-(3) может и не существовать.

Однако если , то ряд (II) сходетс* равномерно при и любом t и определяет непрерывную функюао u(xtt). В этом случае можно говорить лишь об обобщенная решении задачи. Каждая из функций определяет так называемые собств енные колебания струны, закрепленной на концах. При собственных колебаниях, отвечающих к = 1, струна издает основной, самый низкий тон.

При колебаниях, соответствующих ббльшим Л.она издает более высокие тоны, обертоны. Записав *) в виде заключаем, что собственные колебания струны — стоячие волны, при которых точки струны совершают гармонические колебания с амплитудой Нк sin частотой Метод Фурье Мы рассмотрели случай свободных колебаний однородной струны, закрепленной на концах. Рассмотрим теперьслуч ай других граничных условий.

Пусть, например, левый конец струны закреплен, u(0, t) = 0, а правый конец х — 1 упругосвязан со своим положением равновесия, что соответствует условию . Нетривиальное решение u(x, t) уравнения (1), удовлетворяющее поставленным граничным условиям, будем опять искать в виде В результате подстановки в уравнение (1) приходим к следующей задаче о собственных значениям: найти такие значения параметра Л, для которых дифференциальное уравнение при граничных условиях имеет нетривиальные решения Х(х). Общее решение уравнения (15) имеет вид (А > 0)

Первое из граничных условий

Первое из граничных условий (16) дает С = 0, так что функциями Х(х) с точностью до постоянного множителя являются sin у/Хх. Из второго граничного условия Положим А = ir. Тогда Для отыскания и получаем трансцендентное уравнение. Корни этого уравнения можно найти графически, взяв в плоскости (f, z) сечения последовательных ветвей кривой z = tg(i//) прямой линией z = (рис. 7).

Обе части уравнения (18) — нечетные функции относительно р, поэтому каждому положительному корню i/fc соответствует равный ему по абсолютной величине отрицательный корень. Поскольку изменение знака Uk не влечет за собой появления новых собственных функций (они только изменят знак, что несущественно), достаточно ограничиться положительными корнями уравнения (18).

В результате опять получается последовательность собственных значений и отвечающие им последовательности собственных функций и собственных колебаний Кстати, для n-ой собственной частоты ип получается асимптотическое соотношение в частности, для I = т имеем Если правый конец струны х = I свободен, получаем cos vl = 0. Отсюда ul = § + тиг, так что в случае свободного конца собственные значения и собственные функции соответственно равны

Присылайте задания в любое время дня и ночи в ➔

Официальный сайт Брильёновой Натальи Валерьевны преподавателя кафедры информатики и электроники Екатеринбургского государственного института.

Все авторские права на размещённые материалы сохранены за правообладателями этих материалов. Любое коммерческое и/или иное использование кроме предварительного ознакомления материалов сайта natalibrilenova.ru запрещено. Публикация и распространение размещённых материалов не преследует за собой коммерческой и/или любой другой выгоды.

Сайт предназначен для облегчения образовательного путешествия студентам очникам и заочникам по вопросам обучения . Наталья Брильёнова не предлагает и не оказывает товары и услуги.

Видео:Уравнение колебаний струны. Метод разделения переменных. Метод ФурьеСкачать

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

1. Преобразование Фурье и спектр сигнала

Во многих случаях задача получения (вычисления) спектра сигнала выглядит следующим образом. Имеется АЦП, который с частотой дискретизации Fd преобразует непрерывный сигнал, поступающий на его вход в течение времени Т, в цифровые отсчеты — N штук. Далее массив отсчетов подается в некую программку, которая выдает N/2 каких-то числовых значений (программист, который ~~утянул из инета~~ написал программку, уверяет, что она делает преобразование Фурье).

Чтобы проверить, правильно ли работает программа, сформируем массив отсчетов как сумму двух синусоид sin(10*2*pi*x)+0,5*sin(5*2*pi*x) и подсунем программке. Программа нарисовала следующее:

рис.1 График временной функции сигнала

рис.2 График спектра сигнала

На графике спектра имеется две палки (гармоники) 5 Гц с амплитудой 0.5 В и 10 Гц — с амплитудой 1 В, все как в формуле исходного сигнала. Все отлично, программист молодец! Программа работает правильно.

Это значит, что если мы подадим на вход АЦП реальный сигнал из смеси двух синусоид, то мы получим аналогичный спектр, состоящий из двух гармоник.

Итого, наш реальный измеренный сигнал, длительностью 5 сек, оцифрованный АЦП, то есть представленный дискретными отсчетами, имеет дискретный непериодический спектр.

С математической точки зрения — сколько ошибок в этой фразе?

Теперь ~~начальство решило~~ мы решили, что 5 секунд — это слишком долго, давай измерять сигнал за 0.5 сек.

рис.3 График функции sin(10*2*pi*x)+0,5*sin(5*2*pi*x) на периоде измерения 0.5 сек

рис.4 Спектр функции

Что-то как бы не то! Гармоника 10 Гц рисуется нормально, а вместо палки на 5 Гц появилось несколько каких-то непонятных гармоник. Смотрим в интернетах, что да как…

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

рис.5 Добили нулей до 5 сек

рис.6 Получили спектр

Все равно не то, что было на 5 секундах. Придется разбираться с теорией. Идем в Википедию — источник знаний.

2. Непрерывная функция и представление её рядом Фурье

Математически наш сигнал длительностью T секунд является некоторой функцией f(x), заданной на отрезке (X в данном случае — время). Такую функцию всегда можно представить в виде суммы гармонических функций (синусоид или косинусоид) вида:

(1), где:

k — номер тригонометрической функции ( номер гармонической составляющей, номер гармоники)
T — отрезок, где функция определена (длительность сигнала)
Ak — амплитуда k-ой гармонической составляющей,
θk- начальная фаза k-ой гармонической составляющей

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

(Более строго, среднеквадратичное отклонение ряда от функции f(x) будет стремиться к нулю, но несмотря на среднеквадратичную сходимость, ряд Фурье функции, вообще говоря, не обязан сходиться к ней поточечно. См. https://ru.wikipedia.org/wiki/Ряд_Фурье.)

Этот ряд может быть также записан в виде:

(2),
где , k-я комплексная амплитуда.

(3)

Связь между коэффициентами (1) и (3) выражается следующими формулами:

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

Итого:
Математической основой спектрального анализа сигналов является преобразование Фурье.

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

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

Например, на графике рис.7 исходная функция определена на отрезке , а ряд Фурье представляет периодическую функцию, определенную на всей оси х.

Это происходит потому, что синусоиды сами являются периодическими функциями, соответственно и их сумма будет периодической функцией.

рис.7 Представление непериодической исходной функции рядом Фурье

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

Периоды гармонических составляющих кратны величине отрезка , на котором определена исходная функция f(x). Другими словами, периоды гармоник кратны длительности измерения сигнала. Например, период первой гармоники ряда Фурье равен интервалу Т, на котором определена функция f(x). Период второй гармоники ряда Фурье равен интервалу Т/2. И так далее (см. рис. 8).

рис.8 Периоды (частоты) гармонических составляющих ряда Фурье (здесь Т=2π)

Соответственно, частоты гармонических составляющих кратны величине 1/Т. То есть частоты гармонических составляющих Fk равны Fk= кТ, где к пробегает значения от 0 до ∞, например к=0 F0=0; к=1 F1=1T; к=2 F2=2T; к=3 F3=3T;… Fk= кТ (при нулевой частоте — постоянная составляющая).

Пусть наша исходная функция, представляет собой сигнал, записанный в течение Т=1 сек. Тогда период первой гармоники будет равен длительности нашего сигнала Т1=Т=1 сек и частота гармоники равна 1 Гц. Период второй гармоники будет равен длительности сигнала, деленной на 2 (Т2=Т/2=0,5 сек) и частота равна 2 Гц. Для третьей гармоники Т3=Т/3 сек и частота равна 3 Гц. И так далее.

Шаг между гармониками в этом случае равен 1 Гц.

Таким образом сигнал длительностью 1 сек можно разложить на гармонические составляющие (получить спектр) с разрешением по частоте 1 Гц.
Чтобы увеличить разрешение в 2 раза до 0,5 Гц — надо увеличить длительность измерения в 2 раза — до 2 сек. Сигнал длительностью 10 сек можно разложить на гармонические составляющие (получить спектр) с разрешением по частоте 0,1 Гц. Других способов увеличить разрешение по частоте нет.

Существует способ искусственного увеличения длительности сигнала путем добавления нулей к массиву отсчетов. Но реальную разрешающую способность по частоте он не увеличивает.

3. Дискретные сигналы и дискретное преобразование Фурье

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

Обычная схема измерения и оцифровки сигнала выглядит следующим образом.

рис.9 Схема измерительного канала

Сигнал с измерительного преобразователя поступает на АЦП в течение периода времени Т. Полученные за время Т отсчеты сигнала (выборка) передаются в компьютер и сохраняются в памяти.

рис.10 Оцифрованный сигнал — N отсчетов полученных за время Т

Какие требования выдвигаются к параметрам оцифровки сигнала? Устройство, преобразующее входной аналоговый сигнал в дискретный код (цифровой сигнал) называется аналого-цифровой преобразователь (АЦП, англ. Analog-to-digital converter, ADC) ( Wiki).

Одним из основных параметров АЦП является максимальная частота дискретизации (или частота семплирования, англ. sample rate) — частота взятия отсчетов непрерывного во времени сигнала при его дискретизации. Измеряется в герцах. (( Wiki))

Согласно теореме Котельникова, если непрерывный сигнал имеет спектр, ограниченный частотой Fмакс, то он может быть полностью и однозначно восстановлен по его дискретным отсчетам, взятым через интервалы времени , т.е. с частотой Fd ≥ 2*Fмакс, где Fd — частота дискретизации; Fмакс — максимальная частота спектра сигнала. Другими слова частота оцифровки сигнала (частота дискретизации АЦП) должна как минимум в 2 раза превышать максимальную частоту сигнала, который мы хотим измерить.

А что будет, если мы будем брать отсчеты с меньшей частотой, чем требуется по теореме Котельникова?

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

Рис. 11. Появление ложного сигнала низкой частоты при недостаточно высокой частоте дискретизации

Чтобы избежать эффекта алиасинга перед АЦП ставят специальный антиалиасинговый фильтр — ФНЧ (фильтр нижних частот), который пропускает частоты ниже половины частоты дискретизации АЦП, а более высокие частоты зарезает.

Для того, чтобы вычислить спектр сигнала по его дискретным отсчетам используется дискретное преобразование Фурье (ДПФ). Отметим еще раз, что спектр дискретного сигнала «по определению» ограничен частотой Fмакс, меньшей половине частоты дискретизации Fd. Поэтому спектр дискретного сигнала может быть представлен суммой конечного числа гармоник, в отличие от бесконечной суммы для ряда Фурье непрерывного сигнала, спектр которого может быть неограничен. Согласно теореме Котельникова максимальная частота гармоники должна быть такой, чтобы на нее приходилось как минимум два отсчета, поэтому число гармоник равно половине числа отсчетов дискретного сигнала. То есть если в выборке имеется N отсчетов, то число гармоник в спектре будет равно N/2.

Рассмотрим теперь дискретное преобразование Фурье (ДПФ).

Сравнивая с рядом Фурье

видим, что они совпадают, за исключением того, что время в ДПФ имеет дискретный характер и число гармоник ограничено величиной N/2 — половиной числа отсчетов.

Формулы ДПФ записываются в безразмерных целых переменных k, s, где k – номера отсчетов сигнала, s – номера спектральных составляющих.
Величина s показывает количество полных колебаний гармоники на периоде Т (длительности измерения сигнала). Дискретное преобразование Фурье используется для нахождения амплитуд и фаз гармоник численным методом, т.е. «на компьютере»

Возвращаясь к результатам, полученным в начале. Как уже было сказано выше, при разложении в ряд Фурье непериодической функции (нашего сигнала), полученный ряд Фурье фактически соответствует периодической функции с периодом Т. (рис.12).

рис.12 Периодическая функция f(x) с периодом Т0, с периодом измерения Т>T0

Как видно на рис.12 функция f(x) периодическая с периодом Т0. Однако из-за того, что длительность измерительной выборки Т не совпадает с периодом функции Т0, функция, получаемая как ряд Фурье, имеет разрыв в точке Т. В результате спектр данной функции будет содержать большое количество высокочастотных гармоник. Если бы длительность измерительной выборки Т совпадала с периодом функции Т0, то в полученном после преобразования Фурье спектре присутствовала бы только первая гармоника (синусоида с периодом равным длительности выборки), поскольку функция f(x) представляет собой синусоиду.

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

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

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

Рис.13 Пример функции и спектра сигнала кинематической погрешности редуктора

При меньшей длительности картина будет выглядеть «хуже»:

Рис.14 Пример функции и спектра сигнала вибрации ротора

На практике бывает сложно понять, где «реальные составляющие», а где «артефакты», вызванные некратностью периодов составляющих и длительности выборки сигнала или «скачками и разрывами» формы сигнала. Конечно слова «реальные составляющие» и «артефакты» не зря взяты в кавычки. Наличие на графике спектра множества гармоник не означает, что наш сигнал в реальности из них «состоит». Это все равно что считать, будто число 7 «состоит» из чисел 3 и 4. Число 7 можно представить в виде суммы чисел 3 и 4 — это правильно.

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

Некоторые итоги

1. Реальный измеренный сигнал, длительностью T сек, оцифрованный АЦП, то есть представленный набором дискретных отсчетов (N штук), имеет дискретный непериодический спектр, представленный набором гармоник (N/2 штук).

2. Сигнал представлен набором действительных значений и его спектр представлен набором действительных значений. Частоты гармоник положительны. То, что математикам бывает удобнее представить спектр в комплексной форме с использованием отрицательных частот не значит, что «так правильно» и «так всегда надо делать».

3. Сигнал, измеренный на отрезке времени Т определен только на отрезке времени Т. Что было до того, как мы начали измерять сигнал, и что будет после того — науке это неизвестно. И в нашем случае — неинтересно. ДПФ ограниченного во времени сигнала дает его «настоящий» спектр, в том смысле, что при определенных условиях позволяет вычислить амплитуду и частоту его составляющих.

Использованные материалы и другие полезные материалы.