к оглавлению

Непрерывные преобразования Фурье и Лапласа

  1. Интеграл Фурье
  2. Тригонометрическая форма интеграла Фурье
  3. Преобразование Лапласа
  4. Обобщенный ряд Фурье
  5. Быстрое преобразование Фурье
  6. Преобразование Фурье в среде MathCAD
  7. Преобразование Лапласа в среде MathCAD
  8. Преобразование Фурье в среде Maple
  9. Литература по преобразованию Фурье

Интеграл Фурье

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

Зададим периодическую последовательность импульсов и разложим импульс на одном периоде Т в ряд Фурье (формула 4.1.2). Не меняя положения импульса на интервале Т, увеличим значение Т в два раза (продлеваем интервал нулями), при этом выражение (4.1.2) для вычисления спектра остается без изменения, но по ней рассчитывается в 2 раза большее количество гармоник с уменьшением в 2 раза частоты первой гармоники и шага Dw=2p/T. Увеличение интервала Т не влияет на результаты вычисления интеграла функции (4.1.2), т.к. интервал продления заполнен нулевыми значениями сигнала.

По существу, при увеличении периода Т без изменения финитного сигнала форма спектра по оси частот остается без изменения, изменяется только шаг дискретизации спектра и, за счет множителя 1/Т, в 2 раза уменьшаются значения спектра. Новые гармоники располагаются в интервалах между гармониками первого ряда. Пример изменения спектра при увеличении периода Т в 2 раза приведен на рис. 4.2.1.

Рис. 4.2.1.

Процесс можно продолжить дальнейшим последовательным увеличением периода, при этом спектр будет приближаться к непрерывной функции. В пределе, при T ® Ґ , периодическая последовательность импульсов заменяется одиночным финитным сигналом, дискретные частоты nDw обращаются в непрерывные текущие значения w (Dw = 2p/T ® 0), суммирование амплитудных значений заменятся интегрированием, а сами значения спектра становятся бесконечно малыми (1/Т® 0). Для исключения последнего множитель 1/Т из расчетной формулы S(w) исключается, и интегральное преобразование Фурье приобретает следующий вид:

s(t) = (1/2p)S(w)exp(jwt) dw, (4.2.1)

S(w) =s(t)exp(-jwt) dt. (4.2.2)

Формулу (4.2.2) обычно называют формулой прямого преобразования Фурье, а формулу (4.2.1) – обратного преобразования Фурье. Этими выражениями устанавливается взаимно однозначная связь сигнала и его спектра, а точнее – плотности спектра сигнала в последовательной полосе малых (стремящихся к нулю) полосах частот. Эту величину называют спектральной плотностью сигнала. Спектральные функции содержат ровно столько информации, сколько и исходный сигнал.

При преобразовании сигнала в пространство гармонических частот и обратно формулы прямого и обратного преобразований Фурье тождественны за исключением знака аргументов экспоненты:

s(t) =S(f)exp(j2pft) df, (4.2.1')

S(f) =s(t)exp(-j2pft) dt. (4.2.2')

На рис. 4.2.2 сплошной кривой приведен пример непрерывного сигнала s(t), энергия которого сосредоточена на конечном интервале T = (0,25). Если нас не интересует форма данного сигнала за пределами интервала Т, то спектр сигнала в виде ряда Фурье можно определить по формуле (4.1.2). При обратном преобразовании Фурье по формуле (4.1.1), т.е. при восстановлении сигнала по его спектру, в интервале Т будет восстановлен исходный сигнал s(t). Но если интервал для восстановления будет задан больше интервала Т, например равным 0-2Т, то за пределами этого интервала начнется периодическое повторение исходного сигнала, как это показано пунктиром на рис. 4.2.2. Если такой процесс нежелателен и за пределами интервала Т должны быть сохранены нулевые значения сигнала, то необходимо использовать интегральное преобразование Фурье (4.2.1, 4.2.2). При этом следует учитывать особенности интегрального преобразования.

Спектральная функция S(w) представляет собой комплексную спектральную плотность сигнала, непрерывную на частотном интервале от - Ґ до Ґ . Если s(t) – вещественная функция, то спектр этой функции является сопряжено симметричным относительно нулевой частоты

S(-w) = S*(w)

и содержит четную действительную и нечетную мнимую части:

S(w) = A(w) - jB(w), (4.2.3)

A(w) =s(t)cos(wt) dt, (4.2.4)

B(w) =s(t)sin(wt) dt. (4.2.5)

Как и в случае рядов Фурье, вещественные четные функции имеют вещественный четный спектр, представленный только спектральной функцией A(w), а вещественные нечетные – нечетный и только мнимый спектр, представленный спектральной функцией B(w).

Пример спектральной функции S(f) для сигнала s(t) на рис. 4.2.2 приведен на рис. 4.2.3. Как правило, графическое отображение спектральных функций выполняется в виде модуля и аргумента спектральной функции (амплитудного и фазового спектра), приведенных на рис. 4.2.4.

Рис. 4.2.3. Рис. 4.2.4.

Такое представление аналогично (4.1.3'):

R(w) = , (4.2.6)

j(w) = arctg(-B(w)/A(w)), (4.2.7)

но в отношении функции модуля также имеет смысл спектральной плотности модуля.

Заметим также, что сопряженная симметричность спектральной функции позволяет в формулах (4.2.1)-(4.2.2) менять местами знаки аргументов в экспонентах, при этом изменяется только знак мнимой части и аргумента спектра.

Еще раз подчеркнем различие между спектрами и спектральными функциями сигналов. При практическом использовании формулы (4.2.2) для вычисления спектральных функций конечных сигналов, заданных на определенном интервале Т, пределы интегрирования обычно устанавливаются по границам интервала Т, так как нет необходимости выполнять интегрирование в бесконечных пределах, если за пределами интервала Т мы имеем нулевые (или незначимые) значения сигнала. Однако при сравнении формулы (4.2.2) с выражением (4.1.2) можно наглядно видеть, что значения интеграла (4.2.2) не нормируются на величину интервала Т. Отсюда следует, что числовые отсчеты значений модуля функции S(w) для определенных значений wi не являются амплитудными значениями соответствующих гармонических колебаний с частотой wi. Значения S(w) по сравнению со значениями функции S(nDw) по (4.1.2) при nDw = wi завышены на множитель Т. Это можно объяснить тем, что обратное преобразование Фурье по (4.1.1) представляет собой прямое суммирование гармоник с соответствующими амплитудами колебаний, в то время как интегрирование по (4.2.1) представляет собой предельное суммирование значений S(wi)Ч dwi, где dw = 2p/T (или, в обычном частотном представлении, df = 1/T) при Т Ю Ґ .

Что касается спектра фазовых углов, то значения по (4.2.7) и по (4.1.3') при nDw = wi полностью совпадают, так как их вычисление производится по отношению мнимой и действительной части спектра, наличие (или отсутствие) постоянного множителя в которых не меняет значение отношения.

Тригонометрическая форма интеграла Фурье

Тригонометрическая форма интеграла Фурье (при объединении комплексно сопряженных частей спектральных функций):

s(t) = (1/2p) [A(w)cos(wt)+B(w)sin(wt)] dw. (4.2.8)

s(t) = (1/2p)R(w)cos(wt - j(w)) dw. (4.2.8')

Прямое и обратное преобразование Фурье подобны. Любая теорема, доказанная для прямого преобразования Фурье, справедлива и для обратного преобразования, и наоборот. Это непосредственно следует из выражений прямого и обратного преобразования Фурье, которые различаются только знаком в экспоненте. Особенно наглядно (см. рис. 4.2.5) это видно для четных сигналов (заданных функциями, симметричными относительно t = 0), для которых В(w) = 0 и, соответственно, фазовый спектр равен нулю:

s(t) = 2S(f)cos(2pft)df, S(f) = 2s(t)cos(2pft)dt.

В математическом анализе для упрощения записей используют символическую форму обозначения преобразования Фурье:

s(t) Ы S(f), s(t) Ы S(w),

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

Для физических сигналов и их достаточно корректных математических моделей преобразование Фурье, как правило, всегда существует. С чисто математических позиций сигналу s(t) можно сопоставить спектральную плотность S(w), если существует интеграл:

|s(t)| dt < Ґ . (4.2.9)

Преобразование Лапласа

Если условие (4.2.9) не выполняется, то определенные приближения спектральных плотностей вычисляются с использованием специальных методов, одним из которых является преобразование Лапласа.

Допустим, что функция s(t) задана на интервале (0, Ґ ), равна нулю при t<0, а интеграл спектральной функции (4.2.2) расходится. Умножим s(t) на экспоненциальную функцию exp(-ct), где с - положительная константа, и выберем значение 'с' таким, чтобы произведение u(t) = s(t)Ч exp(-ct) удовлетворяло условию абсолютной интегрируемости. Сущность данной операции хорошо видна на рис. 4.2.6. Интегрируемость функции u(t) может быть установлена для любой функции s(t) соответствующим выбором коэффициента 'c'. При этом спектральная плотность функции u(t) может быть вычислена по формуле (4.2.2):

U(w,c) =[s(t) exp(-ct)] exp(-jwt) dt.

После объединения экспоненциальных функций это выражение можно переписать следующим образом:

U(c+jw) =s(t) exp[-(c+jw)t] dt. (4.2.10)

Соответствующее обратное преобразование Фурье функции U(c+jw):

(1/2p)S(c+jw) exp(jwt) dw = s(t) exp(-ct).

Для восстановления функции s(t) достаточно умножить обе части данного выражения на exp(ct), объединить экспоненциальные множители под интегралом и заменить переменную интегрирования w на c+jw:

s(t) = (1/2p)S(c+jw) exp[(c+jw)t] d(c+jw). (4.2.11)

Обозначим комплексную переменную c+jw в выражениях (4.2.10,4.2.11) через р и получим общепринятую форму прямого и обратного преобразования Лапласа:

S(p) =s(t) exp[-pt] dt. (4.2.10')

s(t) = (1/2p)S(p) exp(pt) dp. (4.2.11')

Сигнальную функцию s(t) в преобразованиях Лапласа обычно называют оригиналом, а ее спектральную функцию S(p) - изображением оригинала. Пример спектральной функции Лапласа для оригинала - сложного и неограниченного во времени сигнала, состоящего из каузальной суммы трех гармоник, приведен на рис. 4.2.7. По спектральной функции Лапласа можно выделить эти три основных частоты сигнала и оценить соотношение их амплитуд. Ширина пиков спектральной функции при выделении "чистых" гармоник зависит от значения коэффициента 'c' и уменьшается при его уменьшении.

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

Обобщенный ряд Фурье

Тригонометрические функции не является единственно возможными функциями разложения сигналов. В общем случае разложение сигнала s(t) на интервале (a, b) в ряд вида ckjk(t) может быть выполнено по произвольным функциям jk(t). При задании минимальной погрешности приближения

Ds =[s(t) - ckjk(t)]2 dt

коэффициенты ck могут быть найдены из системы линейных уравнений:

=[s(t) - ckjk(t)]2 dt = 0, k = 0,1,2,…N.

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

jm(t) jn(t) dt =,

то процесс нахождения коэффициентов ck оказывается наиболее простым:

ck =s(t) jk(t) dt,

и для принятого значения N погрешность приближения Ds является минимальной. Если при N ® Ґ имеет место Ds ® 0, система функций jk(t) называется базисной системой координат пространства сигналов L2[a, b] . При этом имеет место равенство:

s(t) =ckjk(t).

Разложение по ортонормированной системе базисных функций называется обобщенным рядом Фурье, а набор коэффициентов ck представляет собой спектр функции s(t) в соответствующем базисе.  В зависимости от специфики решаемых задач применяются различные системы базисных функций. В частности, используются разложения по полиномам Лежандра, Чебышева, Лагерра, Эрмита, функциям Хаара и Уолша и т.п.

Литература по преобразованию Фурье

  1. Бендат Дж., Пирсол А. Прикладной анализ случайных данных. – М.: Мир, 1989. – 540 с.
  2. Гольденберг Л.М. и др. Цифровая обработка сигналов: Учебное пособие для вузов. - М.: Радио и связь, 1990.- 256 с.
  3. Канасевич Э.Р. Анализ временных последовательностей в геофизике. - М.: Недра, 1985.- 300 с.
  4. Никитин А.А. Теоретические основы обработки геофизической информации: Учебник для вузов. - М.: Недра, 1986.- 342 с.
  5. Оппенгейм А.В., Шафер Р.В. Цифровая обработка сигналов. – М.: Связь, 1979. – 416 с.
  6. Рапопорт М.Б. Вычислительная техника в полевой геофизике: Учебник для вузов. - М.: Недра, 1993.- 350 с.
к оглавлению

Copyright ©2006 Davydov А.V.

Александр Пасечник: У кемеровской трагедии есть имя и фамилия

ВСЕ ВИДЕО
Знаете ли Вы, что cогласно релятивистской мифологии "гравитационное линзирование - это физическое явление, связанное с отклонением лучей света в поле тяжести. Гравитационные линзы обясняют образование кратных изображений одного и того же астрономического объекта (квазаров, галактик), когда на луч зрения от источника к наблюдателю попадает другая галактика или скопление галактик (собственно линза). В некоторых изображениях происходит усиление яркости оригинального источника." (Релятивисты приводят примеры искажения изображений галактик в качестве подтверждения ОТО - воздействия гравитации на свет)
При этом они забывают, что поле действия эффекта ОТО - это малые углы вблизи поверхности звезд, где на самом деле этот эффект не наблюдается (затменные двойные). Разница в шкалах явлений реального искажения изображений галактик и мифического отклонения вблизи звезд - 1011 раз. Приведу аналогию. Можно говорить о воздействии поверхностного натяжения на форму капель, но нельзя серьезно говорить о силе поверхностного натяжения, как о причине океанских приливов.
Эфирная физика находит ответ на наблюдаемое явление искажения изображений галактик. Это результат нагрева эфира вблизи галактик, изменения его плотности и, следовательно, изменения скорости света на галактических расстояниях вследствие преломления света в эфире различной плотности. Подтверждением термической природы искажения изображений галактик является прямая связь этого искажения с радиоизлучением пространства, то есть эфира в этом месте, смещение спектра CMB (космическое микроволновое излучение) в данном направлении в высокочастотную область. Подробнее читайте в FAQ по эфирной физике.

НОВОСТИ ФОРУМАФорум Рыцари теории эфира
Рыцари теории эфира
 20.04.2018 - 20:37: СОВЕСТЬ - Conscience -> РУССКИЙ МИР - Карим_Хайдаров.
28.03.2018 - 18:15: СОВЕСТЬ - Conscience -> Проблема государственного терроризма - Карим_Хайдаров.
22.03.2018 - 09:33: ВОСПИТАНИЕ, ПРОСВЕЩЕНИЕ, ОБРАЗОВАНИЕ - Upbringing, Inlightening, Education -> Просвещение от Ю.Ю. Болдырева - Карим_Хайдаров.
04.10.2017 - 15:26: ЭКОНОМИКА И ФИНАНСЫ - Economy and Finances -> ПРОБЛЕМА КРИМИНАЛИЗАЦИИ ЭКОНОМИКИ - Карим_Хайдаров.
04.10.2017 - 05:02: Беседка - Chatter -> "Зенит"ы с "Протон"ами будут падать - Карим_Хайдаров.
03.10.2017 - 18:16: ВОСПИТАНИЕ, ПРОСВЕЩЕНИЕ, ОБРАЗОВАНИЕ - Upbringing, Inlightening, Education -> Просвещение от О.Н. Четвериковой - Карим_Хайдаров.
03.10.2017 - 07:24: ЦИТАТЫ ЧУЖИХ ФОРУМОВ - Outside Quotings -> ЗА НАМИ БЛЮДЯТ - Карим_Хайдаров.
03.10.2017 - 05:48: Беседка - Chatter -> WHO IS WHO - КТО ЕСТЬ КТО - Карим_Хайдаров.
02.10.2017 - 19:04: АСТРОФИЗИКА - Astrophysics -> Апериодическая комета C/2014 Q2 Lovejoy - Карим_Хайдаров.
02.10.2017 - 14:57: СОВЕСТЬ - Conscience -> РАСЧЕЛОВЕЧИВАНИЕ ЧЕЛОВЕКА. КОМУ ЭТО НАДО? - Карим_Хайдаров.
01.10.2017 - 13:58: Беседка - Chatter -> ЭПИСТОЛЯРНАЯ ФИЗИКА - Карим_Хайдаров.
01.10.2017 - 07:23: СОВЕСТЬ - Conscience -> НАСАтые астропиндосы - Карим_Хайдаров.
Bourabai Research Institution home page

Боровское исследовательское учреждение - Bourabai Research Bourabai Research Institution