Text
                    АКАДЕМИЯ ФЕДЕРАЛЬНОЙ СЛУЖБЫ ОХРАНЫ
РОССИЙСКОЙ ФЕДЕРАЦИИ
В. М. Радыгин, С. И. Федорова
Тригонометрические ряды
и интеграл Фурье
Пособие
Орел 2010
. - м	омам
-: й 6 Я А Я
*отв«л Г
•3W-A4L-
00117013

УДК 517.443 ББК 22.161.1 Р159 Радыгин, В. М. Р159 Тригонометрические ряды и интеграл Фурье : пособие / В. М. Рады- гин, С. И. Федорова. — Орел : Академия ФСО России, 2010. - 70 с. Пособие написано в соответствии с программами по дисциплинам "Математи- ка" и "Математический анализ" для подготовки курсантов по специальностям 210404-210406, 230102, 090106. Помимо теоретического материала оно содержит большое число примеров, которые призваны оказать помощь курсантам в овладении навыками в разложении функций в тригонометрические ряды. Значительное место в пособии отводится гармоническому анализу и спектральному представлению наи- более типичных сигналов. Обсуждено и рекомендовано к изданию на заседании кафедры № 4, протокол № 7 от 10 июня 2009 г. УДК 517.443 ББК 22.161.1 © Академия ФСО России, 2010
СОДЕРЖАНИЕ § 1. Периодические функции и некоторые их свойства..........4 § 2. Гармонические колебания................................6 § 3. Сложение гармонических колебаний. Тригонометрический ряд.9 § 4. Ортогональность функций тригонометрической системы....10 § 5. Коэффициенты Фурье. Ряд Фурье..........................13 § 6. Достаточные условия разложения функции в ряд Фурье. Теорема Дирихле.............................................17 § 7. Ряд Фурье для функции, заданной на промежутке длины Т..22 § 8. Ряды Фурье четных и нечетных функций .................26 Г Т § 9. Ряды Фурье функции, заданной на промежутке 0; — ......31 2 § 10. Ряд Фурье в комплексной форме........................36 §11. Гармонический анализ периодических функций............40 § 12. Формула Парсеваля-Ляпунова............................47 § 13. Обобщенный ряд Фурье............................... 49 § 14. Интеграл Фурье. Теорема Фурье........................52 § 15. Интеграл Фурье четных и нечетных функций.............57 § 16. Комплексная форма интеграла Фурье....................59 § 17. Преобразования Фурье.................................65 БИБЛИОГРАФИЧЕСКИЙ СПИСОК...................................69
§ 1. ПЕРИОДИЧЕСКИЕ ФУНКЦИИ И НЕКОТОРЫЕ ИХ СВОЙСТВА Многие важные процессы, с которыми приходится сталкивать- ся при изучении различных явлений в физике и технике совершают- ся циклически, т. е. воспроизводятся в прежнем виде через опреде- ленный промежуток времени Т, называемый периодом рассматри- ваемого процесса. Периодическими, в частности, являются многочисленные ко- лебательные процессы, происходящие в электрических и радиотех- нических цепях. Различные величины, характеризующие такие про- цессы, представляют собой периодические функции времени t. Так, сигналы, с которыми приходится иметь дело в современной теории связи, как правило, являются изменяющимися во времени токами или напряжениями и весьма часто представляют собой периодиче- ские функции времени. Напомним, что функция /(х) называется периодической, если существует действительное число Т тН), такое, что при любом х из области определения данной функции числа х-Т и х + Т принад- лежат этой области и выполняется неравенство /(ЛТ) = /(Л) = ./'(Л- + Т). (I) Число Т называют периодом функции. По определению, для того чтобы получить представление о функции /(х), достаточно изучить ее на любом интервале длины Т. За такой интервал чаще всего принимают один из двух интервалов (0; Т) или а - 2Т а — Т 0 а Рис.1 4
Примерами периодических функций могут служить тригоно- метрические функции j = sinx, j = cosx с периодом Т = 2л, j = tgx, у = ctgx с периодом л. Периодической функцией будет и функция у = с (с - константа), ее периодом является любое дей- ствительное число Т. В дальнейшем будут широко использованы следующие свой- ства периодических функций: 1. Если Т - период функции /(х), то любое число кТ, где к - целое число, также является периодом функции. 2. Сумма, разность, произведение и частное двух периодиче- ских функций периода Т есть периодическая функция периода Т. 3. Сумма сходящегося ряда функций периода Т также есть функция периода Т. 4. Если функция /(х) имеет период Т, то функция Т ф(х) = f (сох) имеет период —. со Действительно, = /(сох + Т) - f (сох) = ф(х). 5. Интеграл от периодической функции /(х) по любому про- межутку [а; ос + 77], длина которого равна периоду Т, имеет одно и то же постоянное значение: а+Т Т \f{x)dx = J/(x)c/x. а О (2) Геометрически равенство (2) выражает равенство площадей заштрихованных криволинейных трапеций с основаниями [0; Г] и ос; ос + г] (рис. 2). 5
Рис. 2 § 2. ГАРМОНИЧЕСКИЕ КОЛЕБАНИЯ Простейший периодический процесс - гармоническое колеба- ние. Опишем этот процесс математически. Предположим, что точка z движется по окружности радиуса А с постоянной угловой скоростью со (рис. 3). Рис. 3 Ее проекции zx и zv на ось Ох и Оу соответственно будут совершать при этом ко- лебательные движения. Най- дем законы движения точек Z И Z у . Пусть в начальный мо- мент времени t = 0 точка z занимает положение zq, ее ра- диус-вектор составляет с осью Ох угол (pg. По прошествии времени t радиус-вектор точки z составит с осью Ох угол со/1 + (pg • Абсцисса точки z изменяется по закону х = zx = A cos(co^ + фо), а ордината - по закону у = z у = A sin(co/ + Фо) • (3) (4) 6
Колебания, описываемые функциями (3) и (4), называются гармони ч еским г i. Итак, при равномерном движении точки z по окружности ра- диуса А ее проекции zx и zy на оси координат совершают гармо- нические колебания. При этом А называют амплитудой, (pg - на- чальной фазой, а со - угловой частотой гармонического колебания. Период гармонического колебания равен промежутку времени со за который точка z, двигаясь по окружности, возвратится вновь в свое исходное положение zq. В соответствии с этим угловая частота гармонического коле- бания 2п (0 =— Т равна тому углу, на который поворачивается радиус-вектор Oz дви- жущейся точки z за одну секунду. Любую из гармоник (3) и (4) всегда можно записать в виде a cos coz + b sin со/, где а и b - некоторые постоянные числа. Например, A sin (со/ + (pg) — A sin cpg cos со/ + A cos фд sin со/ = - a cos со/ + b sin со/, где а = A sin (pg, b = Л cos (pg. Справедливо и обратное утверждение: всякая функция вида a cos со/ + b sin со/ 7
/2 2 представляет собой гармонику с амплитудой А = уа +Ь , угловой частотой (pg, определяемой из условий а = Язтфф, b = Лсозфф: a cos coz + /?sin со/ = Л5ш(со/ + фо). При описании закона движения точки z по окружности вместо двух действительных функций (3) и (4) удобно рассматривать одну комплексную функцию действительного переменного: ^(cos(co/ + 9o) + /sin(cD/ + 9o))= (5) Эту функцию называют комплексной гармоникой. Действительные гармоники (3) и (4), выражающие законы из- менения координат вектора (9z, можно представить в виде действи- тельной и мнимой части функции (5): Л cos (со/ + фо) = Re(y4ez(wz+<*>°)), Л8т(сп/ + фо)= Комплексную гармонику (5) можно представить в виде Z = = Ае1™, где А = Ае1(?° - так называемая комплексная амплитуда гармоники. Комплексная амплитуда А геометрически изображается ра- диусом-вектором Ozq начального положения zq движущейся точки z по окружности. Этот вектор вращается с постоянной угловой ско- ростью (0. 8
§ 3. СЛОЖЕНИЕ ГАРМОНИЧЕСКИХ КОЛЕБАНИИ. ТРИГОНОМЕТРИЧЕСКИЙ РЯД Задача исследования периодических процессов, описывающих работу современной аппаратуры связи, значительно облегчается, ес- ли эти процессы представимы в виде наложения гармонических ко- лебаний, а величины, их характеризующие, в виде суммы гармоник. Рассмотрим следующую последовательность гармоник: = А\ sin(cox + cp| )= <7} coscox + Z>i sin сох, у2 — ^2 sin(2cox + <р2)= cos 2сох + ^2 s^n 2cox, yn = An sin(/?cox + <pw) - an cos ncox + bn sin /7 cox. (6) Независимая переменная здесь обозначена через х. Нулевая гармоника - постоянная А. Она же обозначена и через —— для еди- 2 нообразия приводимых в дальнейшем формул. Последующие гар- моники имеют частоты, кратные со: со, 2со, ..., исо. Периоды этих гармоник: 5 2л n(i) 2л Очевидно, что число Т = — будет общим периодом для всех со гармоник. Сумма этих гармоник Aq + ^Afc sin (Acox + <р£) = —— + а к cos Acox + sin k(dx (7) к=] 2 А=1 9
есть, согласно свойству периодических функций (3), периодическая . 2л функция с периодом Т = —. со Замечено, что суммы гармоник приводят к самым разнообраз- ным периодическим функциям, отнюдь не похожим на гармоники. Естественно возникает вопрос: нельзя ли произвольную перио- дическую функцию представить в виде суммы гармоник или, как говорят, разложить на гармоники? Оказалось, что этого, как прави- ло, сделать нельзя, если ограничиться конечной суммой гармоник. Если же привлечь к рассмотрению ряды, членами которых являются гармоники, то практически любую периодическую функцию можно разложить на гармоники. Итак, рассмотрим функциональный ряд ОО OQ Aq + ^Ап sin(«COX + Фн) = — + ^r//?COS/?01Y + ^Sin/7(D.r, (8) гг=1 2 п=1 составленный из гармоник (6). Этот ряд называется тригонометрическим рядом, а числа 6/q,с/] ,6/2носят название коэффициентов триго- нометрического ряда. Если этот ряд сходится, то его сумма S(x) бу- дет периодической функцией периода Т = 2л со § 4. ОРТОГОНАЛЬНОСТЬ ФУНКЦИЙ ТРИГОНОМЕТРИЧЕСКОЙ СИСТЕМЫ В основании любого тригонометрического ряда — + 2L ап cos п^х + sin «сох 2 /7 = 1 лежит тригонометрическая система функций 10
—, cos cor, sin cox, cos 2a)x, sin 2cox,..., cos ий)х, sin и nix... (9) Эта система обладает одним замечательным свойством, кото- рое служит стержнем всей теории тригонометрических рядов. Прежде чем сформулировать это свойство, введем: Определение. Система функций Ф1(х),ф2(х),...,фл(х)... назы- вается ортогональной на отрезке [«;/?], если Ь /фл (х)фш(х)с/х = а О прии /и; п,т = 1, 2, 3... с > 0 при/? = /и; п = 1, 2, 3... (10) Убедиться в том, что тригонометрическая система функций ор- тогональна, нетрудно. Для этого надо показать, что для функций системы выполняются условия (10). В самом деле 2 j — cosncnxt/x = -----sin /7 (lit 2тссо — sin исо— = — sin ил = 0, ий) 2 ий) f — sinncoxcZx — 0, так как подынтегральная функция нечетна. Применяя известные формулы тригонометрии: cos а • cos р = - (cos(oc - р) + cos(a + р)), cos а • sin Р = - (sin(a - р) + sin(a + Р)),
sin ос • sin p = — (cos(ot - P) - cos(ot + P)), получим при натуральных к и I f cos kaxx cos Iwxdx — jcos(£ — /)сох<7х + fcos(A: + l)(j)xdx — О 2 f cos Acox sin laxxdx T k*l 2 js in(Zc - l)oxxdx + — 2 Jsin(£ + l)(jixdx — 0, T 2 2 f sin kcax sin ktoxdx = k*l — jcos(£ - l)(tixdx — jcos(£ + /)сох<Ух = 0. Кроме того, г 2 I cos и(охах = L 2 j. 2 2 j(l + cos 2n($x)dx = 2 2 1 2r sin nmxdx = — J(1 - cos 2hcdx)<7x = 12
Заметим, что в силу периодичности подынтегральных функций полученные равенства справедливы для любого другого отрезка длины Т. Следовательно, тригонометрическая система функций ортого- нальна на любом отрезке длины Т. § 5. КОЭФФИЦИЕНТЫ ФУРЬЕ. РЯД ФУРЬЕ Пусть тригонометрический ряд правильно сходится на всей числовой оси к периодической функции /(х) с периодом Г, т. е. оо -----1- 2, ап cos исох + bn sin яшг = /(х). 2 /7-1 Выразим коэффициенты ряда а$,ап,Ьп через сумму ряда /(х). Для нахождения коэффициента проинтегрируем в пределах от----до — равенство (11): OQ 2 2 2 Jcos nttixdx + bn Jsin ясохс/х = [f (x)<7x. T T T В силу ортогональности тригонометрической системы все сла- гаемые, стоящие под знаком суммы, равны нулю, поэтому 13
Учитывая, что Та^ ---, получаем ^0=- (12) Определим коэффициент ап (п - целое положительное чис- ло). Для этого умножим обе части равенства (И) на costar и про- интегрируем полученное равенство по отрезку [ f(x) cos k(f>xdx = f 00 ^0 7 ,7 V4 — cos kttixdx + [cos «arc cos koxxdx + 2 + bn [sin лсох cos kwxdx. > £ I у? В силу ортогональности тригонометрической системы все ин- тегралы в правой части равны нулю, кроме одного - при коэффици- енте который равен —. Будем иметь f (x) cos kdixdx = Ofc 2 или, заменив к на и, 14
2 J f (x) cos nayxdx. T (13) 2 Для отыскания коэффициента bn обе части равенства (11) ум- ножим на sin/xox и проинтегрируем в пределах от-до —: 2 2 J/(х)sinkcaxdx - J — sinktiixdx + ^an (cosrztDxsinkayxdx + T _T 2 n=X _T 2 2 2 T 2 + bn (sin «coxsin AclhyA. T Все интегралы в правой части равны нулю, кроме интеграла, Т являющегося коэффициентом при Ьп, равного —: J f (х) sin kitixdx — b^ ~ или 2 Т 2 (/(х) sin natxdx. Т (14) Определение. Коэффициентами Фурье функции f (х) называ- ются числа ао,апфп, вычисленные по формулам 15
а0 2 \f(x)dx; ап 2 Т 2 J/(x) cos n(dxdx; bn T 2 T 2 {/(x)sin nfjyxdx, T где T = 2л co Тригонометрический ряд, коэффициенты которого являются коэффициентами Фурье для функции /’(х), называется рядом Фурье для функции /(х). Замечание. В силу известного свойства интеграла от периоди- ческой функции /(х) при вычислении коэффициентов Фурье вме- сто промежутка можно брать любой другой промежуток длины Т: 2 а+Т «0 = — |/(х)с/х J ос+Г f / (х) cos noxxdx; а 2 Т а+Т f/(x)sin ntnxdx. а Пусть функция /(х) такова, что для нее можно найти коэффи- циенты Фурье. Тригонометрический ряд оо cos п сох + bn sin п сох, (15) 16
где a§,an,bn - коэффициенты Фурье функции /(х), будет рядом Фурье для этой функции независимо от того, сходится он или рас- ходится, и независимо от того, совпадает ли его сумма 5’(х) с функ- цией / (х), для которой он составлен, или нет. § 6. ДОСТАТОЧНЫЕ УСЛОВИЯ РАЗЛОЖЕНИЯ ФУНКЦИИ В РЯД ФУРЬЕ. ТЕОРЕМА ДИРИХЛЕ Какие же требования надо предъявлять к функции /(х), чтобы составленный для нее ряд Фурье сходился и его сумма 5(х) совпа- дала с/(х)? Для этого, очевидно, достаточно предположить, чтобы функ- ция /(х) была интегрируема на отрезке . Ведь для любой интегрируемой, в частности, непрерывной или кусочно- непрерывной функции /(х) можно вычислить коэффициенты Фурье ац,ап,Ьп и формально записать соответствующий ей ряд Фурье: f (х) ~ —— + 22 ап cos п сох + bn sin п сох. 2 И = 1 (16) Однако нельзя утверждать, что формально записанный ряд Фурье для функции /(х) сходится. Он может расходиться для всех х. Более того, даже в случае сходимости ряда его сумма £(х) может не совпадать с / (х). В этой связи возникает вопрос о том, какие условия следует наложить на функцию /(х), чтобы составленный для нее ряд Фурье, во-первых, сходился и, во-вторых, его сумма S(x) совпадала с /(х). Только в этом случае можно считать решенной задачу о разложении функции /(х) в ряд Фурье и в соотношении (16) символ соответ- ствия заменить знаком равенства: ./«=? оо + 22 ап cos сох + sin и шх. Л? = 1 (17) 17
На поставленный вопрос отвечает теорема Дирихле, которую ниже приводим без доказательства. Но прежде дадим определение. Определение. Функция /(х) кусочно-монотонна на отрезке [«; />], если этот отрезок можно разбить конечным числом точек XI, Х2,..., хн-| на интервалы (а; х0, (х^ х2; b) так, что на каждом из интервалов функция монотонна, т. е. либо невозрастаю- щая, либо неубывающая. Из определения следует, что если функция /’(х) монотонна и ограничена на отрезке [о; b], то она может иметь только точки раз- рыва первого рода. Теорема Дирихле. Пусть периодическая функция /(х) с пе- риодом Т на каком-либо отрезке длины Т монотонная и ограничен- ная. Тогда ряд Фурье, составленный для функции / (х), сходится на всей числовой прямой, причем его сумма S(x) равна функции /(х) во всех точках непрерывности этой функции 5(л) = Дх) и равна среднему арифметическому пределов слева и справа во вся- кой точке х разрыва функции S(.v) = Дх - 0) + Дх + 0) Замечание. Условия монотонности и ограниченности называ- ются условиями Дирихле. Пример 1. Осциллограмма напряжения, зарегистрирован- ная на экране осциллографа, представляет периодическую последо- вательность импульсов г/(х) = < 0 при - 2 < х < 0 1 при 0 < х < 2 с периодом Т = 4. Разложить в ряд Фурье функцию z/(x). 18
Решение. Условия теоремы Дирихле для функции w(x) вы- полнены. Период функции г/(х) Т = 4, поэтому со = — = —- = —. Вычислим коэффициенты Фурье: «О = \и(х)с1х = — ]м(х)б/х = пл -2 пл , cos — хах = 1 . пл —sin — пл 2 ---(sin 71/7-0) = 0 прии = 1, 2, 3... пл ш(х) sin nuxxdx — плх , 1 ---dx = — 2 2 пл Г1 . пл . I IJ • sin — xdx =------cos плх О 2 пл =-----(cos пл — I) = о при п четном, пл при п нечетном. Следовательно, в точках непрерывности 71 . Зя .571 -x + sin—x + sin — 2 2 2 и т J 1 гр 2 п 5 О Т 2 п 2 О о О 19
В точках разрыва х = 2п при п = 0, ± 1, ± 2... сумма ряда равна полусумме левого и правого пределов: S(2n) = 2 На рисунке 4 изображены графики частичных сумм получен- ного ряда: 2 2 . Зя —sin — С увеличением числа входящих в них гармоник суммы S|(x), (х), *?з(х) начинают все точнее и точнее воспроизводить осцил- лограмму импульсного напряжения. Пример 2. Разложить в ряд Фурье функцию 1,если 1 с периодом Т - 2. 20
Решение. Построим график функции (рис. 5). /(*) Рис. 5 Заметим, что данная функция на промежутке [0; 2] удовлетво- ряет условиям теоремы Дирихле. Период функции Т = 2, поэтому 2л 0) = — = л. Т Найдем коэффициенты Фурье: 2 2 1 1 3 а0 = ~ J~ ff (x)dx = [xdx + jl-dx = - + l = -, 1 0 20 0 1 2 2 2T 22 1 an =— J/(x) cos ntnxdx = — J/(x) cos mixdx = Jx cos ллхт/х + T 0 20 0 ИЛ COS/7HX JO — sin mix rm 0 (cos Л2Л — 1) = при n четном при n нечетном, т. e. U2n 0, 1 — 21
2Т 22 1 bn = — J/(х) sin ntaxdx = — J j (x) sin nnxdx = x • sin mtxdx + T 0 20 0 ---COSHflX + mt Sin /7ЛХ — cos лях mt ПТ1 Следовательно, во всех точках непрерывности функции /‘(х) имеем разложение 3 00 4 Л?=1 cos(2h - 1)лх + — sin лях. Ttn В точках разрыва х = 2£, к = 0, ± 1, ± 2... сумма ряда равна по- лусумме левого и правого пределов S(2k) = f(2k - 0) + f(2k + 0) _ 0 + 1 § 7. РЯД ФУРЬЕ ДЛЯ ФУНКЦИИ, ЗАДАННОЙ НА ПРОМЕЖУТКЕ ДЛИНЫ Т В приложениях часто возникает задача о разложении в триго- нометрический ряд функции /(х), заданной лишь на некотором промежутке длины Т, например [oc;a + 7j (рис. 6), и удовлетво- ряющей на этом промежутке условиям Дирихле. О периодичности /(х) здесь, следовательно, не может быть и речи. 22
--• а-Т О ос + Т Рис. 6 Поступим следующим образом. Будем рассматривать функцию j(x) как "участок" некоторой периодической функции ср(х) с пе- риодом Т, которая в промежутке [а; а + Т) совпадает с функцией <р(х + Г) = <р(х), если хе — со* ф(х) = /(х), если хе [а; а + Г]. (18) Такую функцию ф(х) естественно называть периодическим продолжением функции /(х) на всю ось Ох (рис. 7). Рис. 7 Вычислим коэффициенты Фурье ф(х): tf0=- |ф(х)<7х, а 23
а+Т j ср(х) cos ntsyxdx, 01 2 T a+T jcp(x) sin nwxdx. a Запишем ряд Фурье для функции ср(х): со СЮ 1 —— + £_,ап cos п сох + bn sin и сох. 2 п=\ Этот ряд во всех точках непрерывности промежутка (ос; ос + Т) имеет суммой значение функции /(х). Посмотрим теперь, чему равна сумма ряда на концах промежутка [ос; ос + 7]. Здесь возможны два случая: 1) /(ос) = /(ос + Т) - в этом случае точки а + кТ (к = 0, ±1, ±2...) будут точками непрерывности функции ср(/), а по- этому сумма ряда в этих точках будет /(ос) или / (a + Т); 2) /(ос) Ф /(а + Т) - в этом случае точки а + кТ окажутся точ- ками разрыва 1-го рода, а поэтому сумма ряда в этих точках будет l(/(a) + /(a + T)). Таким образом, ряд Фурье функции /(х), заданной на проме- жутке длины Т, есть ряд Фурье для функции, полученной из /(х) периодическим продолжением на всю числовую ось. Пример 3. Разложить в ряд Фурье функцию /« = 3 1-- если хе [0; 2]. Решение. График функции /(х) вместе с ее периодическим продолжением изображен на рисунке 8. 24
На промежутке [0; 2] функция /(х) ограничена и монотонна. Длину промежутка, на котором задана функция /(х), примем за пе- 2л риод Т, т. е. Т = 2, тогда частота со = — = л. Найдем коэффициенты Фурье ац,ап,Ьп: cos rmxdx = 0, — sin rmxdx =---- 2) rm COSMTEX Sin /77LXJq = 3 rm Отсюда в промежутке (0; 2) имеем разложение 25
На концах промежутка [О; 2] ряд Фурье сходится к среднему арифметическому пределов функции слева и справа: 5(0) = 5(2) = /(0 + 0) + /(2-0) _ 3 + 0 _ 3 2 2 ~ 2 § 8. РЯДЫ ФУРЬЕ ЧЕТНЫХ И НЕЧЕТНЫХ ФУНКЦИЙ Функция /(х), определенная на симметричном относительно начала координат промежутке [-а; а|, называется четной, если /(~х) = /(х), Vx G [— а; а , и нечетной, если f (—х) = -/’(х), Vx е [— сг, а . (19) График четной функции симметричен относительно оси Оу (рис. 9), а график нечетной функции симметричен относительно на- чала координат (рис. 10). Рис. 9 Рис. 10 Нетрудно убедиться в том, что произведение двух четных и двух нечетных функций есть функция четная, а произведение двух 26
функций, из которых одна четная, а другая нечетная, есть функция нечетная. Кроме того, если /(х) интегрируема в промежутке [-о; а], то а 2 |/(х)б/х если f(x) четная, О О, если /(х) нечетная. Пусть функция /(х) является четной на промежутке 2’ 2 и удовлетворяет в нем условиям Дирихле. Так как cos мох - функция четная, a sin max - нечетная, то произведение /(х) cos «сох будет функцией четной, а /(x)sinHCOx - нечетной. На основании предыдущего замечания для коэффициентов Фурье четной функции /(х) будем иметь формулы 2 4 2 J/(x)dx = - J/(x)t/x; T 1 о ап - — j/ (х) cos исохб/х = j f (х) cos noxxdx; О п 2 — J/(x)sinHcav6/x = 0. Соответственно этому разложение в ряд Фурье четной функ- ции /(х) в точках ее непрерывности имеет вид 27
оо ДД = ~^ + Е ап COS MDX . п=\ В случае нечетной функции /(х) получим: «О = — f f(x)dx - 0; п f (x) cos nttixdx = 0; T Т п j f (х) sin ntisxdx = — J/(х) sin ntoxdx и поэтому для нечетной функции /(х) в точках непрерывности бу- дем иметь оо f W = X bn sin «сох. /7 — 1 Итак, как говорят, четная функция раскладывается в ряд Фурье "по косинусам", а нечетная - "по синусам". Пример 4. Разложить в ряд Фурье функцию ДД = 1 1, если хе (-2; 0), -1, если хе [0; 2). Решение. График функции изображен на рисунке 11. 28
Рис. 1 1 Функция f(x) кусочно-монотонна, ограничена на промежутке (-2; 2) и является нечетной. 2л 71 Примем Т = 4, поэтому со = — = —. Найдем коэффициенты Фурье функции /(х). В силу четности функции f(x) получим: ^0 Q @ п 0 ’ Л71 при п четном, при п нечетном . Таким образом, при хе (-2; 0)и(0; 2) будет .71 1.3л sin—x + -sin— 2 3 2 1 . 5л -sin — х 5 2 В точке х = 0 сумма ряда Фурье равна 0 и не совпадает со зна- чением функции /(х) в этой точке. 2 Пример 5. Разложить в ряд Фурье функцию у~х , хе [-1; 1]. Решение. График функции изображен на рисунке 12. 29
Принимая Т = 2 (длину промежутка, на которой задана функ- ция), найдем коэффициенты Фурье «О 2 1 2 j f (%) cos nisxxdx = 2 Jxz T О cos mixdx = sin mix + mi 2x cos mix - sin mix bn = 0, так как f(x) - четная функция. Следовательно, внутри промежутка (—1; 1) имеем разложение cos их +—cos 2лх + — cos 3 лх +... 4 9 Ввиду того что /(-!) = /(1), полученное разложение справед- ливо и на концах промежутка. 30
§ 9. РЯДЫ ФУРЬЕ ФУНКЦИИ, ЗАДАННОЙ НА ПРОМЕЖУТКЕ 0;- Пусть функция /(х) задана на промежутке 0; — и удовле- творяет на нем условиям Дирихле. При помощи какой-либо функ- ции <?(х), не нарушающей условий Дирихле, "продолжим" функцию /(х) на промежуток---; 0 (рис. 13). В результате получим функцию q(x), если хе--------;0 (р(х) = f (х), если хе 0; — 2 определенную на промежутке длины Т (рис. 13) и удовлетворяю- щую там условиям Дирихле. Ряд Фурье функции <р(х) на отрезке 0; во всех внутренних точках непрерывности этой функции будет сходиться к заданной функции /(х). 31
Поскольку существует бесконечное множество продолжений д(х) для функции /(х) на промежуток то и рядов Фурье для функции /(х), заданной на промежутке 0; можно полу- чить бесконечно много. На практике чаще всего используют четное (рис. 14) или не- четное (рис. 15) продолжение. Рис. 14 Рис. 15 Ряд Фурье для четного продолжения функции /(х) будет со- держать только косинусы. Но так как в промежутке четное продолжение совпадает с функцией /(х), то во всех точках непре- рывности промежутка будем иметь разложение ,/(х) в ряд Фурье только по косинусам «0 /(*) = оо + COS И сох. /1=1 32
Если же функцию f(x), определенную в промежутке продолжить нечетным образом на промежуток то рассуж- дая аналогично, можно получить ее разложение в ряд Фурье по си- нусам сю /(х) = bn sin исох, И = 1 справедливое во всех точках непрерывности внутри промежутка Пример 6. Функцию /(х) =---------, заданную в промежутке 0; л], разложить в ряд Фурье по косинусам. Решение. График функции f (х) вместе с ее периодическим продолжением изображен на рисунке 16. Жирной линией обозначен график функции /(х), а тонкой - график ее четного периодического продолжения с периодом 2л. Функция /(х) монотонна и непрерывна в промежутке [0; я], следо- вательно, удовлетворяет условиям Дирихле. Учитывая, что 2л 2л , О) = — = — = 1 Т 2л , вычислим коэффициенты Фурье: 33
т = “ \f\x)dx - J О J f (x) cos nwxdx = — J 0 Ko cos nxdx ~ cos nx при n четном, при n нечетном, bn —0. 0 Таким образом, разложение функции f(x) на отрезке [0; п имеет вид Л -X л 4 Ъх + 25 cos5x + ... 2 Здесь учтено, что периодическое продолжение функции /(х) является непрерывной на всей числовой оси функцией. Поэтому ряд Фурье /(х) на концах промежутка [0; л] в точках х-0 и х = л, сходится к значению функции /(х) в соответствующих точках. л — х пр и м е р 7. Функцию f (х) =----, заданную в промежутке [0; л], разложить в ряд Фурье по синусам. Решение. График функции /(х) вместе с ее нечетным пе- риодическим продолжением изображен на рисунке 17. Жирной ли- нией обозначен график самой функции /(х). 34
Рис. 17 Как и в предыдущем примере, /(х) на промежутке [0; л] удов- летворяет условиям Дирихле, а периодическое продолжение имеет гр >-р 2л 2л период 7 = 2л. Тогда со = — = — = 1. Найдем коэффициенты Фурье: п = — j/(x) sin nuxxdx — ~f Т О sin nxdx = О — cos пх----— sin пх " л2 О ” Запишем разложение функции f (х) в ряд Фурье по синусам sin 2х sin3x которое имеет место во всех внутренних точках промежутка [0; л]. Выясним, будет ли справедливо это разложение на концах промежутка [О; л]. Точка х = л является точкой непрерывности периодического продолжения, поэтому 5(л) = /(л). В точке х = 0 периодическое продолжение терпит разрыв, поэтому сумма ряда находится как по- 35
лусумма левого и правого пределов функций в этой точке. Отсюда 5(0) = 0, что не совпадает со значением f (0). Итак, найденное разложение функции /(х) в ряд Фурье по си- нусам верно для всех точек промежутка (0; л . Одну и ту же функцию f(x) =-----, заданную на промежутке 0; л], мы представили и рядом Фурье по косинусам, и рядом Фурье по синусам. § 10. РЯД ФУРЬЕ В КОМПЛЕКСНОЙ ФОРМЕ Пусть функция f(x) в промежутке удовлетворяет условиям Дирихле. Тогда в точках непрерывности этого промежутка справедливо разложение / (х) = — + 22 ап cos + bn s*n n<jXx ’ 2 где ап 2 J f (х) cos naxxdx 5 п j/(x)sin ntitxdx По формулам Эйлера cos исох —in(Ax in(f)x 36
1 i moir -zTzgix | 1 I in(Ax . —in(AK smn(Ax =—e -e = — -e + e 2Г 1 2v Подставив эти выражения для coshCDx и sin/7 сох в (17), получим «О оо ШОХ —inox in (АХ —inox _"0 со in ОХ -1П(АХ (20) п п п п Введем обозначения (21) Тогда разложение (17) можно записать в виде оо / / ,х . . V ЭПОХ /(х) — cq + Л — 1 оо /7=—оо (22) Правая часть этого равенства называется комплексной формой ряда Фурье для функции f (х), а коэффициенты сп (п = 0, ± 1, ±2...), вычисляемые по формулам (21), называются комплексными коэф- фициентами Фурье для функции / (х). Очевидно, что коэффициенты являются комплексно- сопряженными коэффициентам сп: Для коэффициентов сц, сп, с_п можно получить интегральные представления. 37
Действительно, с0 = J/(x)t7x n ~ j/(x) cos nttixdx — i j/(x)sin/?0)хб/х = Г f (x)(cos n($x — i sin n($x)dx = ~z7?cox “n —n ~ n J f (x) cos nwxdx 4- i \f (x) sin z?G)xc/x — 2 1 2 \f (x)(cos «cox + i sin n(f)x)dx = — \f (x)ein(i)xdx Эти формулы можно объединить в одну формулу L 1 2 сп=- \f(x)e~m(,ixdx, (И = о,+1,±2,±3...), (23) 1 JL 7 38
в которой п принимает все целочисленные значения (положитель- ные, отрицательные и нулевое). П р и м е р 8. Разложить функцию /(V) = -1, если - л < х < О, 1, если0 <х <п с периодом Т - 2л в ряд Фурье в комплексной форме. Решение. Функция удовлетворяет условиям Дирихле (рис. 18). Рис. 18 Учитывая, что Т = 2тг, со = 1 согласно формулы (23), найдем комплексные коэффициенты Фурье 39
inn - inn inn —inn 2nni nni —(1 — cos ил) = nni nm Следовательно oo inx 711 Это разложение справедливо для всех х, кроме х = пк, где § 11. ГАРМОНИЧЕСКИМ АНАЛИЗ ПЕРИОДИЧЕСКИХ ФУНКЦИЙ Большинство периодических функций, встречающихся при исследовании периодических процессов, особенно в теории связи, таковы, что их можно представить в виде суммы тригонометриче- ского ряда „ 00 + У (24) Членами тригонометрического ряда являются гармоники cos «сох + bn sin ждх = Ап sin(//cox + ) fl fl- fl- с г с амплитудами А 2 2л п , угловыми частотами со - — 4л 2/77Г Пв) =------ начальными фазами (рл, определяемыми из условии п -Sin<pH = П COS (pw 40
В соответствии с этим, разложение (24) можно записать в виде оо /'(х) = Aq + Ап sin (//сох 4- <р/7). п=1 (25) Физический смысл равенства (25) состоит в том, что периоди- ческая функция /’(х), быть может довольно сложная по своей при- роде, оказалась представленной в виде суммы бесконечного числа более простых гармонических колебаний. В этой связи представление периодической функции /(х) в виде суммы тригонометрического ряда (24) или ряда (25) называется гармоническим анализом периодической функции Постоянная составляющая, или "нулевая гармоника", Т = -у- = 7 \fO)dx 2 1 т 2 является средним значением функции ,/(х) в промежутке Первую гармонику Л] sin(wx + (pi) называют основной гармо- никой, а ее угловую частоту со = 2л Т - основной частотой. Угловые частоты /до последующих гармоник Ап sin(ncox + срл) кратны угло- вой частоте со основной гармоники. Очевидно, что разложение (25) функции /(х) на составляю- щие гармоники полностью определяется набором чисел Aq,A[,^2,...,Ап... - амплитуд гармоник и чисел ср|,ср2,---,(Р/7 ••• _ начальных фаз гармоник. Совокупность амплитуд Aq, Ai , Л2,..., Ап ... гармоник, входя- щих в разложение (25), называется амплитудным спектром, а сово- купность их начальных фаз ср|,(р2,---,(Рл _ фазовым спектром пе- риодической функции /(х). 41
Спектры функции обычно изображают в виде графиков: на го- ризонтальной оси откладывают частоты соп=/7со, на вертикаль- ной - соответствующие амплитуды Ап или начальные фазы ф/7 гар- моник. На рисунке 19 изображен амплитудный спектр периодиче- ской функции /(х). Амплитуды отдельных гармоник представлены величинами соответствующих вертикальных отрезков, которые на- зывают спектральными линиями. Рис. 19 Ниже будет показано, что значения А„ с ростом п стремятся к нулю. Аналогичный характер имеет и фазовый спектр (рис. 20). Рис. 20 В связи с этим говорят, что периодическая функция /(х) имеет дискретные линейчатые амплитудный и фазовый спектры. 42
В тех же самых условиях функция f (х) может быть представ- лена рядом Фурье в комплексной форме: 4-00 /(х)= П = — оо (26) с комплексными коэффициентами Если наряду с частотами со, 2со,...,/?со... комплексных гармоник 1П(ИХ П п ~ спе рассматривать отрицательные частоты - со, - 2со,...,-/7С0..., то шкалу частот нужно дополнить отрицательной полуосью. Совокупность модулей сп при п = 0, ± 1, ± 2, ±3... комплекс- ных амплитуд сл, отложенных против соответствующих частот сон - nd) при п = 0, ± 1, ± 2, ± 3..., составит линейчатый амплитуд- ный спектр функции /(%). Ввиду того что и точки соп = исо находятся на оди- наковом расстоянии со друг от друга, амплитудный спектр будет располагаться симметрично относительно оси ординат (рис. 21). Рис. 21 43
Коэффициенты сп гармоник могут, в частности, оказаться дей- ствительными числами. В этом случае, при построении амплитудно- го спектра функции f(x), можно не переходить к модулям коэффи- циентов откладывая соответствующие спектральные линии с г t учетом знаков чисел сп. Для построения соответствующего фазового спектра нужно найти аргументы ф„ комплексных коэффициентов сп. Ввиду со- пряженности сп и с_п их аргументы отличаются лишь знаками. По- этому линейчатый фазовый спектр функции ,/(х) будет распола- гаться симметрично относительно начала координат (рис. 22). Рис. 22 Таким образом, в зависимости от представления периодиче- ской функции f{x) рядом Фурье в тригонометрической или ком- плексной форме получаются спектры: либо на положительной полу- оси, либо на обеих полуосях частот. С увеличением периода Т функции /(х) уменьшается рас- 2л стояние со = —- между соседними спектральными линиями. При Т —> сю расстояние между спектральными линиями будет безгранич- но сокращаться. Естественно ожидать, что амплитудный и фазовый спектры, все более и более сгущаясь, перейдут в пределе в непре- рывные или, как говорят, выродятся в непрерывные амплитудный и фазовый спектры. 44
При исследовании процессов, характеризующих работу совре- менных средств связи часто приходится иметь дело с сигналами, яв- ляющимися периодическими функциями времени. Для того чтобы изучить структуру и свойства сигнала, прибегают к его гармониче- скому анализу. Пример 9. Найти амплитудный спектр периодического сиг- нала f (/), график которого изображен на рисунке 23. Рис. 23 Решение. Аналитически сигнал можно представить в виде ./(О = с при — Предположим, что период Т сигнала значительно больше дли- тельности т прямоугольного импульса. Величина — = N называется скважностью импульса. Условия теоремы Дирихле для сигнала /(/) выполнены. Тогда в точках непрерывности имеем разложение 45
сю Найдем п ~~ '7 -inoM Tino ТпО inco— —ino— 2 2П . ПОП ---sin--- Тпо 2 . пот sin —— non Спектральные линии обращаются нуль при пот , ----= ЛА', 2л —intfrf 2 2 /7(0 = на шкале частот. Эти точки не зависят от величины периода Т, так как определяются лишь длительностью т импульса. ( 2л В промежутке 0; между началом координат и первым ну- 2л 2л лем — располагается п — — т тсо = N спектральных линий , 2л (включая нулевую линию, приходящуюся на точку — При Т —> оо или, что одно и то же, при увеличении скважности f.. 2л импульса число спектральных линий в промежутке 0; как и в других подобных промежутках, безгранично растет, а расстояние 46
между ними со = — безгранично сокращается. Дискретный спектр сп сигнала сгущается, вырождаясь в пределе в непрерывный спектр. § 12. ФОРМУЛА ПАРСЕВАЛЯ-ЛЯПУНОВА Пусть периодическая функция /’(х) периода Т непрерывна на всей числовой оси и на любом отрезке длины Т кусочно-монотонна. Тогда по теореме Дирихле при всех х имеет место равенство со f (х) = — + cos + bn sin исйх) . 2 /7 = 1 Умножая обе части этого равенства на — /(х) и интегрируя полученное равенство по промежутку [0; Г], получим I J./'2 = "0'7 J/(х№х + £ "/7 7 J7 W cos ntoxdx + 7 0 7 0 /7 = 1у 1 о 2Г . . + Ьп —- J/(х) sin /7(Дх<Тх . 7 0 ) Учитывая формулы (12, 13, 14), для коэффициентов а$,ап,Ьп будем иметь т \f\x)dx 0 И 2 Т (27) Полученное равенство называется равенством Парсеваля- Ляпунова. Если учесть, что 47
где Ап (п = О, 1, 2...) - амплитуда л-й гармоники, то равенство Пар- севаля-Ляпунова можно записать в другом виде 1 | оо J/2(.v)dx = ^02+-XA,2 о 2и=1 (28) Левая часть этого равенства есть конечная величина, следова- тельно, ряд в правой его части сходится. Тогда, согласно необходи- мому признаку сходимости числовых рядов при д—его общий член А„ —> 0. Это означает, что амплитуда Ап с ростом п (что то же самое при со/7 —» «=) стремится к нулю (заметим, что при п —> о© коэффициенты ап и Ьп также стремятся к нулю). Рассмотрим пример использования равенства Парсеваля- Ляпунова. Пусть сигнал описывается напряжением г/(0, которое является 2 периодической функцией с периодом Т. Тогда p(t) = u (/) есть мгновенная мощность, выделяемая на сопротивление в 1 Ом. Энер- гия колебания Е на отрезке [0; Т] определяется формулой Е = \p2(t)dt = \u2dt. О О Средняя мощность колебания на промежутке [0; Т] есть Из равенства Парсеваля-Ляпунова следует, что 48
Полученное равенство дает представление о распределении средней мощности сигнала по диапазону частот и позволяет опреде- лить среднюю мощность по известному амплитудному спектру. Если, например, функция описывает сигнал-помеху, то можно выделить полосы частот, в основном образующие данный сигнал. Построив тот или иной фильтр, не пропускающий в прием- ник эти полосы частот, можно свести к минимуму действие помехи. § 13. ОБОБЩЕННЫЙ РЯД ФУРЬЕ Рассмотренный ранее ряд Фурье функции / (х) со cos п(х)х + b„ sin п(йх f I- в основе своей имеет основную тригонометрическую систему функ- ций —, cosoix, sin ок, cos2oir, sin 2сох..., ортогональную на отрезке длины Т = со Поэтому точнее такой ряд следует называть триго- нометрическим рядом Фурье, ибо ряд Фурье можно строить и на ос- нове других систем функций. Пусть на отрезке [а; 6] задана система непрерывных функций Ф1(*)> Ф2(Д-> ФмМ- (29) Нормой какой-либо функции этой системы называется число Г Фи(Х) = {фл? (х)с/х (и = 1,2...). V и 49
Определение. Система функций {ф/7(х)} называется ортонор- мированной, если нормы всех функций этой системы равны 1 и лю- бые две ее функции ортогональны, т. е. при к Ф т tyk(x)'tym(x)dx = 0. Теорема. Если функция f(x) представлена в виде суммы функционального ряда /(*) = Хап<?п(х) /7 = 1 (30) где фл(х) (и = 1, 2, 3...) образуют ортогональную на отрезке [я; b систему, то такое разложение единственно и коэффициенты ап мо- гут быть определены по формуле Ь |/(х)ф„(х)б7х (31) [фя(х)<7х а Коэффициенты, вычисленные по формуле (31) называются ко- эффициентами Фурье, а ряд с такими коэффициентами - обобщен- ным рядом Фурье данной функции /(х), соответствующим системе функций {фи(х)}. Доказательство. Умножим обе части равенства (30) на ФА- (х): /(х)Фа (х) = Х"пФ«(х)Фа(х). /7 = 1 50
Проинтегрируем обе части этого равенства на промежутке [б/; b], при этом предположим, что ряд можно почленно интегриро- вать: h оо b Е 1<рп (х) 9 к (x)dx а п—\ а В силу ортогональности системы функций ср] (х),(р2 (*)•• все слагаемые в правой части при пФ к обращаются в ноль, поэтому Ь Ь (xXv = ак J<P£2 (x)dx> а а откуда и вытекает формула (31). Отметим, что знаменателем этой формулы является квадрат нормы функции (х) (тем самым эта величина отлична от нуля и коэффициент определяется однозначно). Обобщенный ряд Фурье широко используется в теории пере- дачи сигналов, при этом используются разнообразные системы функций: полиномы Чебышева, Эрмита, Лаггера, Лежандра, функ- ции Уолша и др. В 1933 г. выдающийся отечественный ученый В. А. Котельни- ков, получивший мировое признание за свои работы, доказал сле- дующую теорему. Теорема Котельникова. Если колебание w(/) задается на ин- тервале (- оо; + ©о) и имеет спектр, ограниченный сверху частотой F, то его можно представить обобщенным рядом Фурье: w(0 = ~ , д ч sin 2nF(t - кМ) У г/(ЛА/)------------------, X 2tcF(z - W) (32) где AZ = — 2F Разложение в данном случае выполнено по ортогональной (но ненормированной) системе функций: 51
<?к (0 = sin 2nF(t - кМ) 2nF(t-kM) ’ называемых функциями отсчета. Равенство (32) позволяет восстановить непрерывное колебание u(t) по его отсчетам, взятым через интервалы А/. Таким образом, можно передавать не весь сигнал, а лишь его мгновенные значения, взятые через промежутки времени А/. Такой процесс в теории пере- дачи сообщений называют дискретизацией сигнала по времени t. § 14. ИНТЕГРАЛ ФУРЬЕ. ТЕОРЕМА ФУРЬЕ Пусть функция f (х) удовлетворяет условиям Дирихле на вся- ком конечном промежутке и абсолютно интегрируема на всей числовой оси, т. е. существует несобственный интеграл -рОО j|/(x)| dx = Q. — оо По теореме Дирихле в каждой точке х ряд Фурье этой функ- ции сходится и его сумма ST (х) ~ + X (ап cos п^х + s^n ясох). 2 W = 1 Подставим сюда выражения для коэффициентов Фурье: 2 2 //('Ж 1 _т 2 52
т и получим ап 2 J/(Z) cos noMdt, Т 2 [/(/) sin noMdt Т сю 2 Т 2 J f (/) cos nortdt cos a? cox + Г 2 J ff (7) sin nortdt sin aicox . (33) Рассмотрим теперь переменную частоту (0Z2 = лсо гармоник, из которых составлен ряд (17). В промежутке (0; оо) она принимает равноотстоящие значения (О] = со, 0)2 = 2со,сон ~ 77(0..., получая каждый раз приращение Д(О„ = со/7 — (0л_| = со = 53
Отсюда = Дсо„ • -. Г t л В соответствии с этим равенство (33) можно записать в виде оо И п • (34) Посмотрим, что произойдет с этой формулой при Т —> «>. Пер- вое слагаемое правой части этого равенства стремится к нулю, так как /(О dt = — оо Ввиду того, что точка х считается фиксированной, а интегри- рование ведется по переменной /, интегралы, стоящие под знаком суммы в равенстве (34), будут функциями со„. Обозначим: Т j 2 — ]/(/) COS — £7у (оЭ/7 ), 71 _т_ 2 54
2 f/(Z) sin (f)ntdt = b? ((f) n), т тогда сумму в правой части равенства (34) можно записать в виде OQ cos солх + h'p ((f) п )sin cd^x^co^. /7 = 1 (35) Эта сумма напоминает интегральную сумму для функции fl((o)cos сох + b((f))sin сох на промежутке [0; оо). Поэтому можно допустить, что сумма (34) при Т -> оо будет интегралом" СО /(х) - J(c/((o)coscox + Z?(co)sincox)c/co, О (36) — со Z?(co) = — —со f/'(Z) sin (f)tdt. — со Интеграл /(х) называется интегралом Фурье. Заметим, что формула (36) получена из разложения функции /’(х) в ряд Фурье в результате предельного перехода при Т —> со, поэтому в точках непрерывности функции /(х) значение /(х) равно Строгое доказательство предельного соотношения (36) содержится в подробных курсах анализа. 55
/(х), а в точках разрыва 7(х) = | (/(х - 0) +/(х + 0)), т. е. имеет место формула оо I(x) = J(<i(co)cos (fix + b(co)sin сох)с/со = 0 /(х), если х - точка непрерывности, 1 (37) + /(х + 0)), еслих — точка непрерывности. Формула (37) называется интегральной формулой Фурье. Итак, вышеописанное позволяет сформулировать теорему Фурье, которую приводим без доказательства. Теорема Фурье. Если функция /(х) удовлетворяет условиям Дирихле на всяком промежутке и абсолютно интегрируема на всей числовой оси, то оо j(<7(co)cOS сох + /фсо)sin (ox)Jco = 0 /(x), если x - точка непрерывности, 1 (/(X - 0) +/(x + ОД если х - точка непрерывности. Интеграл Фурье можно представить и в другом виде. Учиты- вая, что я (со) cos (fix + Z?(co)sin сох = /i((o)sin(cfix — (pg), где Л((о) - амплитуда и фо - начальная фаза, соответствующие час- 7 2 2 а (со) + b (со), а фд определяется из ус- ловий <7(со) = Л(со)с08фо, 6(С0) = Л(со)5тфо- Интеграл Фурье дает разложение, вообще говоря, непериоди- ческой функции на непрерывную "сумму" гармоник 56
.1 (co)sin(шл' - фф), частоты которых также изменяются непрерывно. Зависимость амплитуды Л (со) от частоты (амплитудный спектр) изображается непрерывной линией. Анализируя эту зависимость можно выделить те полосы частот, которые играют наибольшую роль в формировании данной функции. § 15. ИНТЕГРАЛ ФУРЬЕ ЧЕТНЫХ И НЕЧЕТНЫХ ФУНКЦИЙ Случаи четной и нечетной функций рассматриваются анало- гично ряду Фурье. Пусть /(х) - четная функция, тогда cz(co) = — J / (х) cos сохб/х = — [/ (х) cos сохс/х — сю О = ~ J/(x)sincoxc/x = 0. —оо (38) (39) Следовательно, интеграл Фурье функции /(х) имеет вид оо [<7(CO)COS СОХб/CD. О (40) Аналогично, если /(х) - нечетная функция, то интеграл Фурье примет вид оо fZ?(co) sin сохс/со, о (41) где Z>(cd) = оо [ / (х) sin ош/х. 0 (42) ), то доопределяя ее Если f (х) определена на полуоси [0; + <*> произвольным образом (с соблюдением требуемых ограничений) на 57
промежутке (- 0) получим представление ее интегралом Фурье (естественно, оно не будет однозначным). Практически удобно брать либо четное, либо нечетное продолжение. Пример 10. Представить функцию, график которой изо- бражен на рисунке 24 интегралом Фурье. Построить амплитудный спектр данной функции. Решение. Данная функция имеет следующее аналитическое выражение: /« = I 1, если - 2 < х < 2, 0, если х <-2, х > 2. Эта функция отвечает условиям теоремы Фурье, поэтому пред- ставима интегралом Фурье, а так как она четная, то 2sinuix 2sin 2со тгсо JLCD 0 Рис. 24 Рис. 25 Интеграл Фурье имеет вид /(х) = - к COS Он сходится при всех х и при этом 58
1, если /(л) = 2 О, если х Амплитуда Л (со) в данном случае равна я(со) так как (/?(со) = 0). График этой функции - амплитудный спектр (рис. 25). § 16. КОМПЛЕКСНАЯ ФОРМА ИНТЕГРАЛА ФУРЬЕ В электротехнике и теории цепей широко используется ком- плексная форма интеграла Фурье. Аналогично ряду Фурье исполь- зуем формулы Эйлера cos сох = - (е/Сйг 4- е-/03х ) 2\ /’ sin сох = (e/01Y - e-/Oir). 2V ' Подставим их в интеграл Фурье: ОО Z J / 7(х)= j (a(co))--(e/OXY 0< 2 1 ОО / = - Д(с/(со) - /7>(со)>/ок 2 о + (я(со) +/£>(со))е ZOir)(7co. Обозначим <?(со) -/Ь(со) = S(co), тогда j т*00 । Н-оо 5(со) = — J/(x)coscoxJx-z — f/(x)sincox<7x = — оо —оо | Ч“°° J ~|-оо ~ — [/(x)(cosсох - i sin cox)cZx = — [ f (x)e-/cox dx. л л J — co —oo 59
Соответственно S (со) = «(со) + /Л (со), тогда 1 Ч- °° 5(<о) = - [ f(x)eiVKdx = S(- со). 71 — оо Теперь интеграл Фурье примет вид /(х) = - °f (s(co)e'“¥ + S(-W)e_'“v )<Ло = 2 О = - [$Хю>'&к<Ло + -pI-coK itaxdw. 2 о 2 о Во втором интеграле сделаем замену: -co=z, d($=-dz и по- лучим 1 00 1 о 1 00 /(х) = - p(C0)e'“vc7w— fS(z)e'“zdz = - JS(w)e'“'Jco. 2 2 2 Итак, интеграл Фурье представлен в виде 1 + оо /(х) = - f S(<o)eiwcdu, — оо (43) где 1 т °° 5(со) = - ( л — оо /(х)е"'“хг/х. (44) Форма интеграла Фурье (43) называется комплексной формой, а функция 5(со) - спектральной функцией (или спектральной плот- ностью). Отметим, что S(co) = (со) + b (со) = Л (со) 60
т. е. модуль спектральной функции равен амплитуде гармоники, со- ответствующей частоте со. Поэтому при нахождении амплитудного спектра удобно пользоваться именно комплексной формой интегра- ла Фурье. Пример 11. Найти амплитуду A((i)) для функции /(X) = при х < О, при х > 0. Р е ш е н и е. По формуле (44) найдем : Тогда Я(со) = 5(со) = В теории цепей достаточно часто приходится иметь дело с им- пульсными непериодическими сигналами. Амплитудные и фазовые спектры сигналов позволяют всесторонне изучить их свойства и на этой основе конструировать более совершенную аппаратуру связи для приема и передачи информации. Пример 12. Найти спектральную плотность прямоугольно- го импульса высоты h и длительности т (рис. 26) /(О = 0, при t 61
Рис. 26 Решение. По формуле (44) найдем т к J л J_ л/со 1 — оо I 2 2 /СОТ /СОТ 2/7 Л (О 2/? . (от ----sin — лсо 2 . сот sm — — fa----2_ Л(ОТ Амплитудный спектр сигнала изображен на рисунке 27. £((0) Рис. 27 Полученный непрерывный спектр одиночного прямоугольного импульса полезно сравнить с дискретным спектром периодического импульса той же формы, что и у рассмотренного в примере 9. 62
Пример 13. Импульс напряжения на выходе дифференци- рующей цепочки изменяется по закону (рис. 28) г/вых ае~ш О при при t > 0 (а > 0), КО. 11 в,х Рис. 28 Найти спектральную характеристику, амплитудный и фазовый спектры импульса. Решение. Находим спектральную плотность (характеристи- ку) импульса: 5(ю) = 4ZCD (a+ild)t 0 63
Амплитудный спектр 5(а>) а 7l(<7 + /СО) а Фазовый спектр ср(со) = argS(co) = arctg со со = -arctg— а Графики амплитудного и фазового спектров изображены на рисунках 29 и 30. 64
§ 17. ПРЕОБРАЗОВАНИЯ ФУРЬЕ Представим интеграл Фурье в комплексной форме в следую- щем виде: • е/Сйгс7со. —сю Обозначим 1 F(co) = -—= f f(x)e~i<Mdx. а/2тс j (45) Функция ^(co) называется образом /(x) по Фурье (от спек- тральной функции 5(со) она отличается только постоянным множи- телем). Переход от / (х) к F(co) называют преобразованием Фурье. Здесь видна аналогия с преобразованием Лапласа: /(х) - оригинал, F(w) - изображение. Предполагая, что функция /(х) непрерывна, можно записать J 4~ОЭ ,/(л) = ^= J F(w)t>'Ao. (46) Эта формула представляет обратное преобразование Фурье. Она позволяет по образу F(yri) найти оригинал /(х). Аналогично вводятся косинус и синус преобразования Фурье. Пусть /(х) - четная функция (либо задана на полуоси [0; ©о)), тогда 65
f(x)= j(7(to) COS G)X<7CO = j — j/(x) COS (llWx COS GJt/(D - 0 > о f (x) COS (DXc/x V 0 COSGX/O). Выражение Fc (CD) = + oo J/(x) coscmz/x 0 (47) называется косинус-образом Фурье функции /(х), а переход от /(х) к Fc(cd) - косинус-преобразованием функции /(х). Формула 4“ОО fF(o))coscoxJco о (48) определяет обратное косинус-преобразование Фурье. Аналогично, если функция /(х) - нечетна либо задана на по- луоси [0; с» ), то формулы Fc (со) = — [ f (х) sin сшх/х VTt 0 (49) г^"+со f (х) - Л — fF(co) sin oixz/co Мтг 0 (50) определяют прямое и обратное синус-преобразование Фурье. Эти преобразования непосредственно используются при реше- нии дифференциальных и интегральных уравнений, уравнений ма- 66
тематической физики. Кроме того, они нашли широкое применение в различных вопросах радиотехники, автоматического регулиро- вания. Идея заключается в следующем. На вход некоторого преобра- зователя В (например, усилителя) подается известная функция /(х) (например, напряжение переменного тока). Преобразователь В пре- вращает функцию /(х) в другую функцию ср(х) = £?(/(х)). При этом известна спектральная характеристика преобразователя - функция В(со), которая описывает воздействие преобразователя на функцию входа, точнее на ее спектральную функцию 5(со) или амплитуду И((о). Требуется найти функцию ср(х) по заданным функциям /(х) и /?(со). Эта задача решается по такому плану: 1. По функции /(х) находим ее образ по Фурье F(co) (либо, если это удобнее, Fc(co) или Fs,(co)). 2. Зная воздействие преобразователя на 5(со) (функция F(co), как известно, отличается от £(со) лишь постоянным множителем), * * находим 5 (со) (и соответственно F (со)), являющуюся спектраль- ной функцией функции ср(х) на выходе. * 3. С помощью обратного преобразования Фурье по F (со) на- ходим <р(х). П р и м е р 14. На вход преобразователя поступает функция /(*) = Преобразователь действует таким образом, что все колебания с частотами со> 1 гасятся, а при 0 < со< 1 амплитуда колебания 5(со) умножается на множитель ———. Найти функцию <р(х) на выходе, sin лсо Решение. Данная функция четная, поэтому будем использо- вать косинус-преобразование 67
z ч |2V , |2 sin сох K.(co)= — 2-coscoxax= — 2------- VTCq h CO n 2 sinrcco = 2J---------- V я co Эта функция отличается от амплитуды Л(со) (в полосе О < со < 1) только постоянным множителем, поэтому функция * Fc (со) имеет вид («•) = со sinrcco при 0 < со< 1, при (О> 1, т. е. Fc* (со) = 2 — при сое [О; 1]. V я Для нахождения функции ф(х) на выходе сделаем обратное ко- синус-преобразование ср(х) = 1~2+°° [2 (со) cos сохб/со = J— I 2 J— cos coxJco =
библиографический список 1. Бугров, Я. С. Высшая математика (в 3-х томах), т. 3 / Я. С. Бугров, С. Н. Никольский. - М. : Дрофа, 2004. - 512 с. 2. Гаврилина, Н. П. Сборник задач по математике, ч. IV / Н. П. Гаври- лина, С. И. Федорова. - Орел : ВИПС, 1994. 3. Ильин, В. А. Основы математического анализа / В. А. Ильин, Э. Г. Позняк. - М. : Наука, 2000.-447 с. 4. Кудрявцев, Л. Д. Математический анализ / Л. Д. Кудрявцев. - М. : Высшая школа, 1973.-481 с. 5. Колмогоров, А. И. Элементы теории функций и функционального анализа / А. Н. Колмогоров, С. В. Фомин. - М. : Наука, 1989. - 623 с. 6. Письменный, Д. Т. Конспект лекций по высшей математике, ч. 2 / Д. Т. Письменный. - М. : Рольф, 2000. - 256 с. 7. Титчмарш, Э. Н. Введение в теорию интегралов Фурье / Э. Н. Тит- чмарш. - М. : Ком Книга, 2007. - 480 с. 8. Федорова, С. И. Ряды Фурье. Интеграл Фурье / С. И. Федорова. - Орел : ВИПС, 1994.- 50 с. 9. Харди, Г. Г. Ряды Фурье / Г. Г. Харди, В. В. Рогозинский. - М. : Ком Книга, 2006.- 158 с.
Учебное издание Радыгин Владимир Михайлович Федорова Светлана Ивановна Тригонометрические ряды и интеграл Фурье Пособие Технический редактор М. Г. Лобанова Корректор О. Н. Горбачева Выпускающий Е. Н. Оброткина Компьютерная верстка М. М. Мухиной Подписано в печать 08.12.2009 г. Формат 60x84/16. Печать офсетная. Бумага офсетная. Гарнитура «Times New Roman». Усл. печ. л. 4,06. Уч.-изд. л. 3,20. Тираж 102 экз. Заказ № 419. Отпечатано в типографии Академии ФСО России 302034, г. Орел, ул. Приборостроительная, 35.