Fourier series for beginners. Fourier series: the history and influence of the mathematical mechanism on the development of science
transcript
1 MINISTRY OF EDUCATION AND SCIENCE OF THE RUSSIAN FEDERATION NOVOSIBIRSK STATE UNIVERSITY FACULTY OF PHYSICS R. K. Belkheeva FOURIER SERIES IN EXAMPLES AND TASKS Tutorial Novosibirsk 211
2 UDC BBK V161 B44 B44 Belkheeva R. K. Fourier series in examples and problems: Textbook / Novosib. state un-t. Novosibirsk, s. ISBN B study guide the basic information about Fourier series is presented, examples are given for each topic studied. An example of applying the Fourier method to solving the problem of transverse vibrations of a string is analyzed in detail. Illustrative material is given. There are tasks for independent solution. Designed for students and teachers Faculty of Physics NSU. Published according to the decision of the Methodological Commission of the Faculty of Physics of NSU. Reviewer Dr. phys.-math. Sciences. V. A. Aleksandrov ISBN c Novosibirsk State University, 211 c Belkheeva R. K., 211
3 1. Decomposition 2π- periodic function in a Fourier series Definition. The Fourier series of the function f(x) is the functional series a 2 + (a n cosnx + b n sin nx), (1) where the coefficients a n, b n are calculated by the formulas: a n = 1 π b n = 1 π f(x) cosnxdx, n = , 1,..., (2) f(x) sin nxdx, n = 1, 2,.... (3) Formulas (2) (3) are called the Euler Fourier formulas. The fact that the function f(x) corresponds to the Fourier series (1) is written as a formula f(x) a 2 + (a n cosnx + b n sin nx) (4) and they say that the right side of formula (4) is a formal series Fourier functions f(x). In other words, formula (4) means only that the coefficients a n, b n are found by formulas (2), (3). 3
4 Definition. A 2π-periodic function f(x) is called piecewise smooth if the interval [, π] contains a finite number of points = x< x 1 <... < x n = π таких, что в каждом открытом промежутке (x j, x j+1) функция f(x) непрерывно дифференцируема, а в каждой точке x j существуют конечные пределы слева и справа: f(x j) = lim h + f(x j h), f(x j +) = lim h + f(x j + h), (5) f(x j h) f(x j) f(x j + h) f(x j +) lim, lim. h + h h + h (6) Отметим, что последние два предела превратятся в односторонние производные после замены предельных значений f(x j) и f(x j +) значениями f(x j). Теорема о представимости кусочно-гладкой функции в точке своим рядом Фурье (теорема о поточечной сходимости). Ряд Фурье кусочно-гладкой 2π-периодической функции f(x) сходится в каждой точке x R, а его сумма равна числу f(x), если x точка непрерывности функции f(x), f(x +) + f(x) и равна числу, если x точка разрыва 2 функции f(x). ПРИМЕР 1. Нарисуем график, найдем ряд Фурье функции, заданной на промежутке [, π] формулой, f(x) = x, предполагая, что она имеет период 2π, и вычислим суммы 1 1 числовых рядов (2n + 1) 2, n 2. n= Решение. Построим график функции f(x). Получим кусочно-линейную непрерывную кривую с изломами в точках x = πk, k целое число (рис. 1). 4
5 Fig. 1. Graph of the function f(x) nx + π n n 2 = 2 π (1) n 1 n 2 = b n = 1 π π = 2 π f(x) cosnxdx = cos nx cos n 2 = 4 πn2, for odd n, for even n, f(x ) sin nxdx = because the function f(x) is even. We write the formal Fourier series for the function f(x): f(x) π 2 4 π k= 5 cos (2k + 1)x (2k + 1) 2.
6 Find out whether the function f(x) is piecewise smooth. Since it is continuous, we calculate only the limits (6) at the end points of the interval x = ±π and at the break point x = : and f(π h) f(π) π h π lim = lim h + h h + h = 1, f(+ h) f(+) + h () lim = lim h + h h + h f(+ h) f(+) + h lim = lim = 1, h + h h + h = 1, f(h) f () h () lim = lim = 1. h + h h + h The limits exist and are finite, hence the function is piecewise smooth. By the pointwise convergence theorem, its Fourier series converges to the number f(x) at each point, i.e., f(x) = π 2 4 π k= cos (2k + 1) + x (2k + 1) 2 = = π 2 4 (cosx + 19 π cos 3x) cos 5x (7) Figures 2 and 3 show the character of the approximation of the partial sums of the Fourier series S n (x), where S n (x) = a n 2 + (a k coskx + b k sin kx), k=1, to the function f(x) in the interval [, π] . 6
7 Fig. Fig. 2. Graph of the function f(x) with superimposed graphs of partial sums S (x) = a 2 and S 1(x) = a 2 + a 1 cos x 3. Graph of the function f (x) with a partial sum graph superimposed on it S 99 (x) \u003d a 2 + a 1 cos x + + a 99 cos 99x 7
8 Substituting in (7) x = we get: = π 2 4 π k= 1 (2k + 1) 2, from where we find the sum of the number series: = π2 8. Knowing the sum of this series, it is easy to find the following sum We have: S = ( ) S = ()= π S, hence S = π2 6, that is, 1 n = π The sum of this famous series was first found by Leonhard Euler. It is often found in mathematical analysis and its applications. EXAMPLE 2. Draw a graph, find the Fourier series of the function given by the formula f(x) = x for x< π, предполагая, что она имеет период 2π, и вычислим суммы числовых (1) n) рядов + n= ((2n + 1,) (k k + 1) Решение. График функции f(x) приведен на рис. 4. 8
9 Fig. 4. Graph of the function f(x) The function f(x) is continuously differentiable on the interval (, π). At the points x = ±π, it has finite limits (5): f() =, f(π) = π. In addition, there are finite limits (6): f(+ h) f(+) lim = 1 and h + h f(π h) f(π +) lim = 1. h + h Hence, f(x) is piecewise smooth function. Since the function f(x) is odd, then a n =. The coefficients b n are found by integration by parts: b n = 1 π f(x) sin πnxdx= 1 [ x cosnx π πn + 1 n = 1 πn [(1)n π + (1) n π] = 2(1)n+ one. n Let us compose the formal Fourier series of the function 2(1) n+1 f(x) sin nx. n 9 cosnxdx ] =
10 According to the pointwise convergence theorem for a piecewise smooth 2π-periodic function, the Fourier series of the function f(x) converges to the sum: 2(1) n+1 sin nx = n f(x) = x if π< x < π, = f(π) + f(π +) 2 =, если x = π, (8) f() + f(+) =, если x =. 2 На рис. 5 8 показан характер приближения частичных сумм S n (x) ряда Фурье к функции f(x). Рис. 5. График функции f(x) с наложенным на него графиком частичной суммы S 1 (x) = a 2 + a 1 cos x 1
11 Fig. Fig. 6. Graph of the function f(x) with the graph of the partial sum S 2 (x) superimposed on it. 7. Graph of the function f(x) with the graph of the partial sum S 3 (x) 11 superimposed on it
12 Fig. 8. Graph of the function f(x) with the graph of the partial sum S 99 (x) superimposed on it. We use the obtained Fourier series to find the sums of two numerical series. We put in (8) x = π/2. Then 2 () +... = π 2, or = n= (1) n 2n + 1 = π 4. We easily found the sum of the well-known Leibniz series. Putting x = π/3 in (8), we find () +... = π 2 3, or (1+ 1) () (k) 3π +...= 3k
13 EXAMPLE 3. Draw a graph, find the Fourier series of the function f(x) = sin x, assuming that it has a period of 2π, and 1 calculate the sum of the number series 4n 2 1. Solution. The graph of the function f(x) is shown in fig. 9. Obviously, f(x) = sin x is a continuous even function with period π. But 2π is also the period of the function f(x). Rice. 9. Graph of the function f(x) Let's calculate the Fourier coefficients. All b n = because the function is even. Using trigonometric formulas, we calculate a n for n 1: a n = 1 π = 1 π sin x cosnxdx = 2 π sin x cosnxdx = (sin(1 + n)x sin(1 n)x) dx = = 1 () π cos( 1 + n)x cos(1 n)x + = 2 () 1 + (1) n = π 1 + n 1 n π 1 n 2 ( 4 1 if n = 2k, = π n 2 1 if n = 2k
14 This calculation does not allow us to find the coefficient a 1 because at n = 1 the denominator goes to zero. Therefore, we calculate the coefficient a 1 directly: a 1 = 1 π sin x cosxdx =. Since f(x) is continuously differentiable on (,) and (, π) and at the points kπ, (k is an integer), there are finite limits (5) and (6), the Fourier series of the function converges to it at every point: = 2 π 4 π sinx = 2 π 4 π cos 2nx 4n 2 1 = (1 1 cos 2x cos 4x + 1) cos 6x 1. Graph of the function f(x) with the graph of the partial sum S(x) superimposed on it 14
15 Fig. Fig. 11. Graph of the function f(x) with the graph of the partial sum S 1 (x) superimposed on it. Fig. 12. Graph of the function f(x) with the graph of the partial sum S 2 (x) superimposed on it. 13. Graph of the function f(x) with the graph of the partial sum S 99 (x) 15 superimposed on it
16 1 Calculate the sum of the number series. To do this, we put 4n 2 1 in (9) x =. Then cosnx = 1 for all n = 1, 2,... and Therefore, 2 π 4 π 1 4n 2 1 =. 1 4n 2 1 = = 1 2. EXAMPLE 4. Let us prove that if a piecewise smooth continuous function f(x) satisfies the condition f(x π) = f(x) for all x (i.e., is π-periodic) , then a 2n 1 = b 2n 1 = for all n 1, and vice versa, if a 2n 1 = b 2n 1 = for all n 1, then f(x) is π-periodic. Solution. Let the function f(x) be π-periodic. Let us calculate its Fourier coefficients a 2n 1 and b 2n 1: = 1 π (a 2n 1 = 1 π f(x) cos(2n 1)xdx + f(x) cos(2n 1)xdx =) f(x) cos (2n 1)xdx. In the first integral we make the change of variable x = t π : f(x) cos(2n 1)xdx = f(t π) cos(2n 1)(t + π) dt. 16
17 Using the fact that cos(2n 1)(t + π) = cos(2n 1)t and f(t π) = f(t), we get: a 2n 1 = 1 π (f(x) cos(2n 1)x dx+) f(x) cos(2n 1)x dx =. It is proved similarly that b 2n 1 =. Conversely, let a 2n 1 = b 2n 1 =. Since the function f(x) is continuous, then, by the theorem on the representability of a function at a point by its Fourier series, we have Then f(x π) = f(x) = (a 2n cos 2nx + b 2n sin 2nx). (a2n cos 2n(x π) + b 2n sin 2n(x π)) = (a2n cos 2nx + b 2n sin 2nx) = f(x), which means that f(x) is a π-periodic function. EXAMPLE 5. Let us prove that if a piecewise smooth function f(x) satisfies the condition f(x) = f(x) for all x, then a = and a 2n = b 2n = for all n 1, and vice versa, if a = a 2n = b 2n =, then f(x π) = f(x) for all x. Solution. Let the function f(x) satisfy the condition f(x π) = f(x). Let us calculate its Fourier coefficients: 17
18 = 1 π (a n = 1 π f(x) cos nxdx + f(x) cosnxdx =) f(x) cosnxdx. In the first integral we make the change of variable x = t π. Then f(x) cosnxdx = f(t π) cosn(t π) dt. Using the fact that cos n(t π) = (1) n cosnt and f(t π) = f(t), we obtain: a n = 1 π ((1) n) f(t) cosnt dt = if n even, = 2 π f(t) cos nt dt, if n is odd. π It is proved similarly that b 2n =. Conversely, let a = a 2n = b 2n =, for all n 1. Since the function f(x) is continuous, then, by the theorem on the representability of a function at a point, its Fourier series satisfies the equality f(x) = (a 2n 1 cos ( 2n 1)x + b 2n 1 sin (2n 1)x). eighteen
19 Then = f(x π) = = = f(x). EXAMPLE 6. Let us study how to extend the function f(x) integrable on the interval [, π/2] to the interval [, π], so that its Fourier series has the form: a 2n 1 cos(2n 1)x. (1) Solution. Let the graph of the function have the form shown in Fig. 14. Since in series (1) a = a 2n = b 2n = for all n, it follows from Example 5 that the function f(x) must satisfy the equality f(x π) = f(x) for all x. This observation gives a way to extend the function f(x) to the interval [, /2] : f(x) = f(x+π), fig. 15. From the fact that series (1) contains only cosines, we conclude that the continued function f (x) must be even (i.e., its graph must be symmetrical about the Oy axis), Fig.
20 Fig. 14. Graph of the function f(x) 15. Graph of the continuation of the function f(x) on the interval [, /2] 2
21 So, the desired function has the form shown in fig. 16. Fig. 16. Graph of the continuation of the function f(x) on the interval [, π] Summing up, we conclude that the function should be continued as follows: f(x) = f(x), f(π x) = f(x), that is interval [π/2, π], the graph of the function f(x) is centrally symmetric about the point (π/2,), and on the interval [, π], its graph is symmetric about the Oy axis. 21
22 GENERALIZATION OF EXAMPLES 3 6 Let l >. Consider two conditions: a) f(l x) = f(x); b) f(l + x) = f(x), x [, l/2]. FROM geometric point From the point of view, condition (a) means that the graph of the function f(x) is symmetric about the vertical line x = l/2, and condition (b) that the graph f(x) is centrally symmetric about the point (l/2;) on the x-axis. Then the following statements are true: 1) if the function f(x) is even and condition (a) is satisfied, then b 1 = b 2 = b 3 =... =, a 1 = a 3 = a 5 =... = ; 2) if the function f(x) is even and condition (b) is satisfied, then b 1 = b 2 = b 3 =... =, a = a 2 = a 4 =... = ; 3) if the function f(x) is odd and condition (a) is satisfied, then a = a 1 = a 2 =... =, b 2 = b 4 = b 6 =... = ; 4) if the function f(x) is odd and condition (b) is satisfied, then a = a 1 = a 2 =... =, b 1 = b 3 = b 5 =... =. PROBLEMS In problems 1 7 draw graphs and find the Fourier series for the functions, (assuming they have a period of 2π: if< x <, 1. f(x) = 1, если < x < π. 1, если < x < /2, 2. f(x) =, если /2 < x < π/2, 1, если π/2 < x < π. 3. f(x) = x 2 (< x < π). 4. f(x) = x 3 (< x < π). { π/2 + x, если < x <, 5. f(x) = π/2 x, если < x < π. 22
23 ( 1 if /2< x < π/2, 6. f(x) = 1, если π/2 < x < 3π/2. {, если < x <, 7. f(x) = sin x, если < x < π. 8. Как следует продолжить интегрируемую на промежутке [, π/2] функцию f(x) на промежуток [, π], чтобы ее ряд Фурье имел вид: b 2n 1 sin (2n 1)x? Ответы sin(2n 1)x sin(2n + 1)x. π 2n 1 π 2n + 1 n= 3. 1 (1) n () 12 3 π2 + 4 cosnx. 4. (1) n n 2 n 2π2 sin nx. 3 n 5. 4 cos(2n + 1)x π (2n + 1) (1) n cos(2n + 1)x. π 2n + 1 n= n= 7. 1 π sin x 2 cos 2nx. 8. Функцию следует продолжить следующим образом: f(x) = f(x), f(π x) = f(x), π 4n 2 1 то есть на промежутке [, π], график функции f(x) будет симметричен относительно вертикальной прямой x = π/2, на промежутке [, π] ее график центрально симметричен относительно точки (,). 23
24 2. Expansion of a function given in the interval [, π] only in terms of sines or only in terms of cosines Let a function f be given in the interval [, π]. In order to expand it in this interval into a Fourier series, we first extend f into the interval [, π] in an arbitrary way, and then we use the Euler Fourier formulas. The arbitrariness in the continuation of a function leads to the fact that for the same function f: [, π] R we can obtain different Fourier series. But it is possible to use this arbitrariness in such a way as to obtain an expansion only in sines or only in cosines: in the first case, it suffices to continue f in an odd way, and in the second, in an even way. Solution algorithm 1. Continue the function in an odd (even) way on (,), and then periodically with a period of 2π continue the function to the entire axis. 2. Calculate the Fourier coefficients. 3. Compose the Fourier series of the function f(x). 4. Check the conditions for the convergence of the series. 5. Specify the function to which this series will converge. EXAMPLE 7. Expand the function f(x) = cosx,< x < π, в ряд Фурье только по синусам. Решение. Продолжим функцию нечетным образом на (,) (т. е. так, чтобы равенство f(x) = f(x) выполнялось для всех x (, π)), а затем периодически с периодом 2π на всю ось. Получим функцию f (x), график которой приведен на рис
25 Fig. 17. Graph of the continued function Obviously, the function f (x) is piecewise smooth. Let's calculate the Fourier coefficients: a n = for all n because the function f (x) is odd. If n 1, then b n = 2 π f(x) sin πnxdx = 2 π cosx sin nxdx = = 2 π dx = = 2 π cos (n + 1) x cos (n 1) x + = π n + 1 n 1 = 1 (1) n (1)n 1 1 = π n + 1 n 1 = 1 if n = 2 k + 1, (1)n+1 (n 1) + (n + 1) = π ( n + 1)(n 1) 2 2n if n = 2k. π n 2 1 For n = 1 in the previous calculations, the denominator vanishes, so the coefficient b 1 can be calculated directly.
26 Essentially: b 1 = 2 π cosx sin xdx =. Compose the Fourier series of the function f (x) : f (x) 8 π k=1 k 4k 2 1 sin 2kx. Since the function f (x) is piecewise smooth, then, by the pointwise convergence theorem, the Fourier series of the function f (x) converges to the sum cosx if π< x <, S(x) =, если x =, x = ±π, cosx, если < x < π. В результате функция f(x) = cosx, заданная на промежутке (, π), выражена через синусы: cosx = 8 π k=1 k 4k 2 1 sin 2kx, x (, π). Рис демонстрируют постепенное приближение частичных сумм S 1 (x), S 2 (x), S 3 (x) к разрывной функции f (x). 26
27 Fig. Fig. 18. Graph of the function f (x) with the graph of the partial sum S 1 (x) superimposed on it. 19. Graph of the function f(x) with the graph of the partial sum S 2 (x) superimposed on it 27
28 Fig. Fig. 2. Graph of the function f (x) with the graph of the partial sum S 3 (x) superimposed on it. 21 shows graphs of the function f (x) and its partial sum S 99 (x). Rice. 21. Graph of the function f (x) with a graph of the partial sum S 99 (x) 28 superimposed on it
29 EXAMPLE 8. Let us expand the function f(x) = e ax, a >, x [, π], in a Fourier series only in cosines. Solution. We continue the function in an even way to (,) (i.e., so that the equality f(x) = f(x) holds for all x (, π)), and then periodically with a period of 2π to the entire real axis. We obtain the function f (x), the graph of which is shown in Fig. 22. Function f (x) at points 22. The graph of the continued function f (x) x = kπ, k is an integer, has kinks. Let us calculate the Fourier coefficients: b n =, since f (x) is even. Integrating by parts, we get 29
30 a n = 2 π a = 2 π = 2 cosnxd(e ax) = 2 πa e ax dx = 2 π a (eaπ 1), f(x) cos πnxdx = 2 π πa eax cosnx = 2 πa (eaπ cosnπ 1 ) + 2n πa 2 π e ax cos nxdx = + 2n e ax sin nxdx = πa sin nxde ax = = 2 π a (eaπ cos n π 1) + 2n π sin nx π a 2eax 2n2 e ax cos nxdx = 2 π a 2 π a (eaπ cos n π 1) n2 a a n. 2 Therefore, a n = 2a e aπ cos n π 1. π a 2 + n 2 Since f (x) is continuous, according to the pointwise convergence theorem, its Fourier series converges to f (x). Hence, for all x [, π] we have f(x) = 1 π a (eaπ 1)+ 2a π k=1 e aπ (1) k 1 a 2 + k 2 coskx (x π). Figures demonstrate the gradual approximation of the partial sums of the Fourier series to a given discontinuous function. 3
31 Fig. 23. Graphs of functions f (x) and S (x) 24. Graphs of functions f (x) and S 1 (x) 25. Graphs of functions f (x) and S 2 (x) 26. Graphs of functions f (x) and S 3 (x) 31
32 Fig. 27. Graphs of functions f (x) and S 4 (x) 28. Graphs of the functions f (x) and S 99 (x) PROBLEM 9. Expand the function f (x) = cos x, x π, in a Fourier series only in cosines. 1. Expand the function f (x) \u003d e ax, a >, x π, in a Fourier series only in terms of sines. 11. Expand the function f (x) \u003d x 2, x π, in a Fourier series only in sines. 12. Expand the function f (x) \u003d sin ax, x π, in a Fourier series in terms of cosines only. 13. Expand the function f (x) \u003d x sin x, x π, in a Fourier series only in sines. Answers 9. cosx = cosx. 1. e ax = 2 [ 1 (1) k e aπ] k sin kx. π a 2 + k2 k=1 11. x 2 2 [ π 2 (1) n 1 π n + 2 ] n 3 ((1)n 1) sin nx. 32
33 12. If a is not an integer, then sin ax = 1 cosaπ (1 + +2a cos 2nx ) + π a 2 (2n) 2 +2a 1 + cosaπ cos(2n 1)x π a 2 (2n 1) 2; if a = 2m is an even number, then sin 2mx = 8m cos(2n 1)x π (2m) 2 (2n 1) 2; if a = 2m 1 is a positive odd number, then sin(2m 1)x = 2 ( cos 2nx ) 1 + 2(2m 1). π (2m 1) 2 (2n) π 16 n sin x sin 2nx. 2 π (4n 2 1) 2 3. Fourier series of a function with an arbitrary period Assume that the function f(x) is defined in the interval [ l, l], l >. By substituting x = ly, y π, we obtain the function g(y) = f(ly/π) defined in the interval π [, π]. This function g(y) corresponds to the (formal) Fourier series () ly f = g(y) a π 2 + (a n cosny + b n sin ny), whose coefficients are found by the Euler Fourier formulas: a n = 1 π g(y) cosny dy = 1 π f (ly π) cos ny dy, n =, 1, 2,..., 33
34 b n = 1 π g(y) sinny dy = 1 π f () ly sin ny dy, n = 1, 2,.... π l, we obtain a slightly modified trigonometric series for the function f(x): where f(x) a 2 + a n = 1 l b n = 1 l l l l l (a n cos πnx l f(x) cos πnx l f(x) sin πnx l + b n sin πnx), (11) l dx, n =, 1, 2,..., (12) dx, n = 1, 2,.... (13) Formulas (11) (13) are said to define expansion in a Fourier series of a function with an arbitrary period. EXAMPLE 9. Find the Fourier series of the function given in the interval (l, l) by the expression ( A if l< x, f(x) = B, если < x < l, считая, что она периодична с периодом 2l. Решение. Продолжим функцию периодически, с периодом 2l, на всю ось. Получим функцию f (x), кусочно-постоянную в промежутках (l + 2kl, l + 2kl), и претерпевающую разрывы первого рода в точках x = lk, k целое число. Ее коэффициенты Фурье вычисляются по формулам (12) и (13): 34
35 a = 1 l l f(x) dx = 1 l A dx + 1 l l B dx = A + B, l l a n = 1 l l l f(x) cos πnx l dx = = 1 l = 1 l l A cos πnx l = A + B π n l b n = 1 l dx + 1 l l B cos πnx l sin πn = if n, l l A sin πnx l f(x) sin πnx l dx + 1 l l dx = B sin πnx l = B A (1 cosπn). πn Compose the Fourier series of the function f (x) : f(x) A + B π (B A Since cosπn = (1) n, then n dx = dx = (1 cosπn) sin πnx). l for n = 2k we get b n = b 2k =, for n = 2k 1 b n = b 2k 1 = 35 2(B A) π(2k 1).
36 Hence f(x) A + B (B A) π (sin πx + 1 3πx sin + 1 5πx sin +... l 3 l 5 l According to the pointwise convergence theorem, the Fourier series of the function f(x) converges to the sum A, if l< x, S(x) = A + B, если x =, x = ±l, 2 B, если < x < l. Придавая параметрам l, A, B конкретные значения получим разложения в ряд Фурье различных функций. Пусть l = π, A =, B = 3π. На рис. 29 приведены графики первых пяти членов ряда, функции f (x) и частичной суммы S 7 (x) = a 2 + b 1 sin x b 7 sin 7x. Величина a является средним значением функции на промежутке. Обратим внимание на то, что с возрастанием ча- 2 стоты гармоники ее амплитуда уменьшается. Для наглядности графики трех высших гармоник сдвинуты по вертикали. На рис. 3 приведен график функции f(x) и частичной суммы S 99 (x) = a 2 + b 1 sin x b 99 sin 99x. Для наглядности на рис. 31 приведен тот же график в другом масштабе. Последние два графика иллюстрируют явление Гиббса. 36).
37 Fig. 29. Graph of the function f (x) with superimposed graphs of the harmonics S (x) = a 2 and S 1 (x) = b 1 sinx. For clarity, the graphs of the three higher harmonics S 3 (x) \u003d b 3 sin 3πx, S l 5 (x) \u003d b 5 sin 5πx l and S 7 (x) \u003d b 7 sin 7πx are shifted vertically up l 37
38 Fig. Fig. 3. Graph of the function f(x) with the graph of the partial sum S 99 (x) superimposed on it. 31. Fragment of fig. 3 in another scale 38
39 PROBLEMS In problems, expand the specified functions in Fourier series in given intervals. 14. f(x) = x 1, (1, 1). 15. f(x) = ch2x, (2, 2] f(x) = x (1 x), (1, 1]. 17. f(x) = cos π x, [ 1, 1] f(x ) = sin π x, (1, 1).( 2 1 if 1< x < 1, 19. f(x) = 2l = 4., если 1 < x < 3; x, если x 1, 2. f(x) = 1, если 1 < x < 2, 2l = 3. { 3 x, если 2 x < 3;, если ωx, 21. f(x) = 2l = 2π/ω. sin ωx, если ωx π; Разложить в ряды Фурье: а) только по косинусам; б) только по синусам указанные функции в заданных промежутках (, l) { 22. f(x) = { 23. f(x) = ax, если < x < l/2, a(l x), если l/2 < x < l. 1, если < x 1, 2 x, если 1 x 2. Ответы 14. f(x) = 4 cos(2n 1)πx. π 2 (2n 1) f(x) = sh sh4 (1) n nπx cos 16 + π 2 n f(x) = cos 2nπx. π 2 n f(x) = 2 π + 8 π (1) n n 1 4n 2 cosnπx. 39
40 18. f(x) = 8 (1) n n sin nπx. π 1 4n (1) n 2n + 1 cos πx. π 2n πn 2πnx π 2 sin2 cos n π sin ωx 2 cos 2nωx π 4n 2 1. (l 22. a) f(x) = al 4 2) 1 (4n 2)πx cos, π 2 (2n 1) 2 l b) f(x) = 4al (1) n 1 (2n 1) πx sin. π 2 (2n 1) 2 l 23. a) f(x) = (cos π π 2 2 x 2 2 cos 2π 2 2 x cos 3π 2 2 x cos 5π), 2 2 x... b) f( x) = 4 (sin π π 2 2 x 1 3 sin 3π)+ 2 2 x (sin π π 2 x cos 2π) 2 x Complex form of the Fourier series Decomposition f(x) = c n e inx, where c n = 1 2π f (x)e inx dx, n = ±1, ±2,..., is called the complex form of the Fourier series. The function expands into a complex Fourier series under the same conditions under which it expands into a real Fourier series. four
41 EXAMPLE 1. Find the Fourier series in the complex form of the function given by the formula f(x) = e ax in the interval [, π), where a is a real number. Solution. Let us calculate the coefficients: = c n = 1 2π f(x)e inx dx = 1 2π e (a in)x dx = 1 ((1) n e aπ (1) n e aπ) = (1)n sh aπ. 2π(a in) π(a in) The complex Fourier series of the function f has the form f(x) sh aπ π n= (1) n a in einx. Let us verify that the function f(x) is piecewise smooth: in the interval (, π) it is continuously differentiable, and at the points x = ±π there are finite limits (5), (6) lim h + ea(+h) = e aπ, lim h + ea(π h) = e aπ, e a(+h) e a(+) lim h + h = ae aπ e a(π h) e a(π), lim h + h = ae aπ. Therefore, the function f(x) can be represented by a Fourier series sh aπ π n= (1) n a in einx, which converges to the sum: ( e S(x) = ax if π< x < π, ch a, если x = ±π. 41
42 EXAMPLE 11. Find the Fourier series in the complex and real form of the function given by the formula f(x) = 1 a 2 1 2a cosx + a2, where a< 1, a R. Решение. Функция f(x) является четной, поэтому для всех n b n =, а a n = 2 π f(x) cosnxdx = 2 (1 a2) π cos nxdx 1 2a cosx + a 2. Не будем вычислять такой сложный интеграл, а применим следующий прием: 1. используя формулы Эйлера sin x = eix e ix 2i = z z 1, cosx = eix + e ix 2i 2 = z + z 1, 2 где z = e ix, преобразуем f(x) к рациональной функции комплексной переменной z; 2. полученную рациональную функцию разложим на простейшие дроби; 3. разложим простейшую дробь по формуле геометрической прогрессии; 4. упростим полученную формулу. Итак, по формулам Эйлера получаем = f(x) = 1 a 2 1 a(z + z 1) + a 2 = (a 2 1)z (z a)(z a 1) = a z a az. (14) 42
43 Recall that the sum of an infinite geometric progression with denominator q (q< 1) вычисляется по формуле: + n= q n = 1 1 q. Эта формула верна как для вещественных, так и для комплексных чисел. Поскольку az = a < 1 и a/z = a < 1, то az = + a n z n = a n e inx, a z a = a z 1 1 a/z = a z n= + n= a n z = + n n= n= a n+1 z = + a n+1 e i(n+1)x. n+1 После замены переменной (n + 1) = k, < k < 1, получим: 1 a z a = a k e ikx. Следовательно, f(x) + n= k= c n e inx, где c n = n= { a n, если n, a n, если n <, то есть c n = a n. Поскольку функция f(x) непрерывна, то в силу теоремы о поточечной сходимости имеет место равенство: f(x) = + n= a n e inx. Тем самым мы разложили функцию f(x) в ряд Фурье в комплексной форме. 43
44 Now let's find the Fourier series in real form. To do this, we group the terms with numbers n and n for n: a n e inx + a n e inx = 2a neinx + e inx Since c = 1, then 2 = 2a n cos nx. f(x) = 1 a 2 1 2a cosx + a = a n cosnx. 2 This is a Fourier series in the real form of the function f(x). Thus, without calculating a single integral, we found the Fourier series of the function. In doing so, we calculated a hard integral depending on the parameter cos nxdx 1 2a cosx + a = 2 π an 2 1 a2, a< 1. (15) ПРИМЕР 12. Найдем ряд Фурье в комплексной и вещественной форме функции, заданной формулой a sin x f(x) = 1 2a cosx + a2, a < 1, a R. Решение. Функция f(x) является нечетной, поэтому для всех n a n = и b n = 2 π f(x) sin nxdx = 2a π sin x sin nxdx 1 2a cosx + a 2. Чтобы записать ряд Фурье нужно вычислить сложные интегралы или воспользоваться приемом, описанным выше. Поступим вторым способом: 44
45 a(z z 1) f(x) = 2i (1 a(z z 1) + a 2) = i 2 + i (a + a 1)z 2 2 (z a)(z a 1) = = i 2 + i () a 2 z a + a 1. z a 1 We expand each of the simple fractions according to the geometric progression formula: + a z a = a 1 z 1 a = a a n z z n, n= z a 1 z a = az = a n z n. n= This is possible because az = a/z = a< 1. Значит + ia n /2, если n <, f(x) c n e inx, где c n =, если n =, n= ia n /2, если n >, or, more briefly, c n = 1 2i a n sgnn. Thus, the Fourier series in complex form is found. Grouping terms with numbers n and n, we obtain the Fourier series of the function in real form: = f(x) = + a sin x 1 2a cosx + a + 2 (1 2i an e inx 1 2i an e inx n= +) = c n e inx = a n sin nx. Again, we managed to calculate the following complex integral: sin x sin nxdx 1 2a cosx + a 2 = π an 1. (16) 45
46 PROBLEM 24. Using (15), calculate the integral cos nxdx 1 2a cosx + a 2 for real a, a > Using (16), calculate the integral sin x sin nxdx for real a, a > a cosx + a2 In problems, find the series Fourier in complex form for functions. 26. f(x) = sgn x, π< x < π. 27. f(x) = ln(1 2a cosx + a 2), a < 1. 1 a cosx 28. f(x) = 1 2a cosx + a2, a < Докажите, что функция f, определенная в промежутке [, π], вещественнозначна, если и только если коэффициенты c n ее комплексного ряда Фурье связаны соотношениями c n = c n, n =, ±1, ±2, Докажите, что функция f, определенная в промежутке [, π], является четной (т. е. удовлетворяет соотношению f(x) = f(x)), если и только если коэффициенты c n ее комплексного ряда Фурье связаны соотношениями c n = c n, n = ±1, ±2, Докажите, что функция f, определенная в промежутке [, π], является нечетной (т. е. удовлетворяет соотношению f(x) = f(x)), если и только если коэффициенты c n ее комплексного ряда Фурье связаны соотношениями c n = c n, n =, ±1, ±2,.... Ответы 1 2π 24. a n a π a n i + e 2inx, где подразумевается, что слагаемое, соответствующее n =, пропущено. π n n= a n n cosnx. 28. a n cosnx. n= 46
47 5. Lyapunov's equality Theorem (Lyapunov's equality). Let a function f: [, π] R be such that f 2 (x) dx< +, и пусть a n, b n ее коэффициенты Фурье. Тогда справедливо равенство, a (a 2 n + b2 n) = 1 π называемое равенством Ляпунова. f 2 (x) dx, ПРИМЕР 13. Напишем равенство Ляпунова для функции { 1, если x < a, f(x) =, если a < x < π и найдем с его помощью суммы числовых рядов + sin 2 na n 2 и + Решение. Очевидно, 1 (2n 1) 2. 1 π f 2 (x) dx = 1 π a a dx = 2a π. Так как f(x) четная функция, то для всех n имеем b n =, a = 2 π f(x) dx = 2 π a dx = 2a π, 47
48 a n = 2 π f(x) cosnxdx = 2 π a cos nxdx = 2 sin na πn. Therefore, the Lyapunov equality for the function f(x) takes the form: 2 a 2 π + 4 sin 2 na = 2a 2 π 2 n 2 π. From the last equality for a π we find sin 2 na n 2 = a(π a) 2 Assuming a = π 2, we obtain sin2 na = 1 for n = 2k 1 and sin 2 na = for n = 2k. Therefore, k=1 1 (2k 1) 2 = π2 8. EXAMPLE 14. Let's write the Lyapunov equality for the function f(x) = x cosx, x [, π], and use it to find the sum of the number series (4n 2 + 1) 2 (4n 2 1) 4. 1 π Solution. Direct calculations give = π π f 2 (x) dx = 1 π x 2 cos 2 xdx = 1 π x sin 2xdx = π π x cos x = π x 21 + cos 2x dx = 2 π 1 4π cos 2xdx =
49 Since f(x) is an even function, then for all n we have b n =, a n = 2 π = 1 π 1 = π(n + 1) = f(x) cosnxdx = 2 π 1 cos(n + 1)x π (n + 1) 2 x cosxcosnxdx = x (cos(n + 1)x + cos(n 1)x) dx = 1 π sin(n + 1)xdx sin(n 1)xdx = π(n 1) π π 1 + cos(n 1)x = π(n 1) 2 1 (= (1) (n+1) 1) 1 (+ (1) (n+1) 1) = π(n + 1) 2 π(n 1) 2 () = (1)(n+1) 1 1 π (n + 1) + 1 = 2 (n 1) 2 = 2 (1)(n+1) 1 n k π (n 2 1) = π (4k 2 1) 2 if n = 2k, 2 if n = 2k + 1. The coefficient a 1 must be calculated separately, since in the general formula for n = 1 the denominator of the fraction vanishes. = 1 π a 1 = 2 π f(x) cosxdx = 2 π x(1 + cos 2x)dx = π 2 1 2π 49 x cos 2 xdx = sin 2xdx = π 2.
50 Thus, the Lyapunov equality for the function f(x) has the form: 8 π + π (4n 2 + 1) 2 π 2 (4n 2 1) = π 2 1) = π π PROBLEM 32. Write the Lyapunov equality for the function ( x f(x) = 2 πx if x< π, x 2 πx, если π < x. 33. Напишите равенства Ляпунова для функций f(x) = cos ax и g(x) = sin ax, x [, π]. 34. Используя результат предыдущей задачи и предполагая, что a не является целым числом, выведите следующие классические разложения функций πctgaπ и (π/ sin aπ) 2 по рациональным функциям: πctgaπ = 1 a + + 2a a 2 n 2, (π) = sin aπ (a n) 2. n= 35. Выведите комплексную форму обобщенного равенства Ляпунова. 36. Покажите, что complex form The Lyapunov equation is valid not only for real-valued functions, but also for complex-valued functions. 5
51 π (2n + 1) = π sin 2απ 2απ = 2sin2 απ α 2 π 2 Answers + 4 sin2 απ π 2 α 2 (α 2 n 2) 2; sin 2απ 1 2απ = απ n 2 4sin2 π 2 (α 2 n 2) 2. 1 π 35. f(x)g(x) dx= c n d n, where c n is the Fourier coefficient 2π of f(x), and d n is the Fourier coefficient functions g(x). 6. Differentiation of Fourier series Let f: R R be a continuously differentiable 2π-periodic function. Its Fourier series has the form: f(x) = a 2 + (a n cos nx + b n sin nx). The derivative f (x) of this function will be a continuous and 2π-periodic function, for which a formal Fourier series can be written: f (x) a 2 + (a n cos nx + b n sin nx), where a, a n, b n, n = 1 , 2,... Fourier coefficients of the function f (x). 51
52 Theorem (on term-by-term differentiation of Fourier series). Under the above assumptions, the equalities a =, a n = nb n, b n = na n, n 1 are true. EXAMPLE 15. Let a piecewise-smooth function f(x) be continuous in the interval [, π]. Let us prove that when the condition f(x)dx = is satisfied, the inequality 2 dx 2 dx, called Steklov's inequality, holds, and we verify that equality in it is realized only for functions of the form f(x) = A cosx. In other words, Steklov's inequality gives conditions under which the smallness of the derivative (in rms) implies the smallness of the function (in rms). Solution. Let us extend the function f(x) to the interval [, ] evenly. Denote the extended function by the same symbol f(x). Then the continued function will be continuous and piecewise smooth on the interval [, π]. Since the function f(x) is continuous, then f 2 (x) is continuous on the interval and 2 dx< +, следовательно, можно применить теорему Ляпунова, согласно которой имеет место равенство 1 π 2 dx = a () a 2 n + b 2 n. 52
53 Since the continued function is even, then b n =, a = by condition. Consequently, the Lyapunov equality takes the form 1 π 2 dx = a 2 π n. (17) Let us make sure that f (x) satisfies the conclusion of the theorem on term-by-term differentiation of the Fourier series, that is, that a =, a n = nb n, b n = na n, n 1. Let the derivative f (x) undergo breaks at the points x 1, x 2,..., x N in the interval [, π]. Denote x =, x N+1 = π. Let us divide the integration interval [, π] into N +1 intervals (x, x 1),..., (x N, x N+1), on each of which f(x) is continuously differentiable. Then, using the additivity property of the integral and then integrating by parts, we get: b n = 1 π = 1 π = 1 π f (x) sin nxdx = 1 π N f(x) sin nx j= N f(x) sin nx j= x j+1 x j x j+1 x j n n π N j= x j+1 x j x j+1 x j f (x) sin nxdx = f(x) cosnxdx = f(x) cosnxdx = = 1 π [(f(x 1) sin nx 1 f(x) sin nx) + + (f(x 2) sinnx 2 f(x 1) sin nx 1)
54 + (f(x N+1) sin nx N+1 f(x N) sin nx N)] na n = = 1 π na n = = 1 π na n = na n. x j+1 a = 1 f (x)dx = 1 N f (x)dx = π π j= x j = 1 N x j+1 f(x) π = 1 (f(π) f()) = . x j π j= Similarly, we get a n = nb n. We have shown that the theorem on term-by-term differentiation of Fourier series for a continuous piecewise-smooth 2π-periodic function whose derivative in the interval [, π] undergoes discontinuities of the first kind is true. So f (x) a 2 + (a n cosnx + b n sin nx) = (na n)sin nx, since a =, a n = nb n =, b n = na n, n = 1, 2,.... Because 2dx< +, то по равенству Ляпунова 1 π 2 dx = 54 n 2 a 2 n. (18)
55 Since each term of the series in (18) is greater than or equal to the corresponding term of the series in (17), then 2 dx 2 dx. Recalling that f(x) is an even continuation of the original function, we have 2 dx 2 dx. Which proves the Steklov equality. Now let us examine for which functions equality holds in Steklov's inequality. If for at least one n 2, the coefficient a n is nonzero, then a 2 n< na 2 n. Следовательно, равенство a 2 n = n 2 a 2 n возможно только если a n = для n 2. При этом a 1 = A может быть произвольным. Значит в неравенстве Стеклова равенство достигается только на функциях вида f(x) = A cosx. Отметим, что условие πa = f(x)dx = (19) существенно для выполнения неравенства Стеклова, ведь если условие (19) нарушено, то неравенство примет вид: a a 2 n n 2 a 2 n, а это не может быть верно при произвольном a. 55
56 PROBLEMS 37. Let a piecewise-smooth function f(x) be continuous on the interval [, π]. Prove that under the condition f() = f(π) = the inequality 2 dx 2 dx, also called Steklov's inequality, holds, and make sure that equality in it holds only for functions of the form f(x) = B sin x. 38. Let a function f be continuous in the interval [, π] and have in it (with the possible exception of only a finite number of points) a square-integrable derivative f(x). Prove that if the conditions f() = f(π) and f(x) dx = are satisfied, then the inequality 2 dx 2 dx, called the Wirtinger inequality, holds, and the equality in it takes place only for functions of the form f(x ) = A cosx + B sinx. 56
57 7. Application of Fourier series for solving partial differential equations When studying a real object (natural phenomena, production process, control system, etc.), two factors turn out to be significant: the level of accumulated knowledge about the object under study and the degree of development of the mathematical apparatus. On the present stage scientific research the following chain was developed: phenomenon physical model mathematical model. The physical formulation (model) of the problem is as follows: the conditions for the development of the process and the main factors influencing it are identified. The mathematical formulation (model) consists in describing the factors and conditions chosen in the physical formulation in the form of a system of equations (algebraic, differential, integral, etc.). A problem is said to be well-posed if, in a certain functional space, the solution of the problem exists, uniquely and continuously depends on the initial and boundary conditions. The mathematical model is not identical to the object under consideration, but is its approximate description. Derivation of the equation of free small transverse vibrations of the string. We will follow the textbook. Let the ends of the string be fixed, and the string itself be taut. If the string is taken out of equilibrium (for example, by pulling or striking it), then the string will start 57
58 hesitate. We will assume that all points of the string move perpendicular to its equilibrium position (transverse vibrations), and at each moment of time the string lies in the same plane. Let us take a system of rectangular coordinates xou in this plane. Then, if at the initial time t = the string was located along the axis Ox, then u will mean the deviation of the string from the equilibrium position, that is, the position of the string point with the abscissa x at an arbitrary time t corresponds to the value of the function u(x, t). For each fixed value of t, the graph of the function u(x, t) represents the shape of the vibrating string at time t (Fig. 32). At a constant value of x, the function u(x, t) gives the law of motion of a point with the abscissa x along a straight line parallel to the Ou axis, the derivative u t is the speed of this motion, and the second derivative 2 u t 2 is the acceleration. Rice. 32. Forces applied to an infinitely small section of a string Let's write an equation that the function u(x, t) must satisfy. To do this, we make some more simplifying assumptions. We will assume that the string is absolutely flexible.
59 coy, that is, we will assume that the string does not resist bending; this means that the stresses arising in the string are always directed tangentially to its instantaneous profile. The string is assumed to be elastic and subject to Hooke's law; this means that the change in the magnitude of the tension force is proportional to the change in the length of the string. Let us assume that the string is homogeneous; this means that its linear density ρ is constant. We neglect external forces. This means that we are considering free vibrations. We will study only small vibrations of a string. If we denote by ϕ(x, t) the angle between the abscissa axis and the tangent to the string at the point with the abscissa x at the time t, then the condition for small oscillations is that the value of ϕ 2 (x, t) can be neglected in comparison with ϕ (x, t), i.e., ϕ 2. Since the angle ϕ is small, then cos ϕ 1, ϕ sin ϕ tg ϕ u, therefore, the value (u x x,) 2 can also be neglected. It immediately follows from this that in the process of oscillation we can neglect the change in the length of any section of the string. Indeed, the length of a piece of string M 1 M 2 projected into the interval of the x-axis, where x 2 = x 1 + x, is equal to l = x 2 x () 2 u dx x. x Let us show that, under our assumptions, the value of the tension force T will be constant along the entire string. To do this, we take some part of the string M 1 M 2 (Fig. 32) at time t and replace the action of the discarded parts
60 kov by the tension forces T 1 and T 2. Since, according to the condition, all points of the string move parallel to the Ou axis and there are no external forces, the sum of the projections of the tension forces on the Ox axis must be equal to zero: T 1 cosϕ(x 1, t) + T 2 cosϕ(x 2, t) =. Hence, due to the smallness of the angles ϕ 1 = ϕ(x 1, t) and ϕ 2 = ϕ(x 2, t), we conclude that T 1 = T 2. Denote general meaning T 1 \u003d T 2 through T. Now we calculate the sum of the projections F u of the same forces onto the Ou axis: F u \u003d T sin ϕ (x 2, t) T sin ϕ (x 1, t). (2) Since for small angles sin ϕ(x, t) tg ϕ(x, t), and tg ϕ(x, t) u(x, t)/ x, equation (2) can be rewritten as F u T (tan ϕ(x 2, t) tan ϕ(x 1, t)) (u T x (x 2, t) u) x (x 1, t) x x T 2 u x 2(x 1, t) x . Since the point x 1 is chosen arbitrarily, then F u T 2 u x2(x, t) x. After all the forces acting on the section M 1 M 2 are found, we apply Newton's second law to it, according to which the product of mass and acceleration is equal to the sum of all acting forces. The mass of a piece of string M 1 M 2 is equal to m = ρ l ρ x, and the acceleration is equal to 2 u(x, t). Newton's t 2 equation takes the form: 2 u t (x, t) x = u 2 α2 2 x2(x, t) x, where α 2 = T ρ is a constant positive number. 6
61 Reducing by x, we get 2 u t (x, t) = u 2 α2 2 x2(x, t). (21) As a result, we have obtained a linear homogeneous partial differential equation of the second order with constant coefficients. It is called the string vibration equation or the one-dimensional wave equation. Equation (21) is essentially a reformulation of Newton's law and describes the motion of a string. But in the physical formulation of the problem, there were requirements that the ends of the string are fixed and the position of the string at some point in time is known. We will write these conditions in equations as follows: a) we will assume that the ends of the string are fixed at the points x = and x = l, i.e., we will assume that for all t the relations u(, t) =, u(l, t ) = ; (22) b) we will assume that at the time t = the position of the string coincides with the graph of the function f(x), i.e., we will assume that for all x [, l] the equality u(x,) = f( x); (23) c) we will assume that at the time t = the point of the string with the abscissa x is given speed g(x), i.e., we will assume that u (x,) = g(x). (24) t Relations (22) are called boundary conditions, and relations (23) and (24) are called initial conditions. Mathematical model of free small transverse 61
62 string vibrations is that it is necessary to solve equation (21) with boundary conditions (22) and initial conditions (23) and (24) Solution of the equation of free small transverse vibrations of the string by the Fourier method< t <, удовлетворяющие граничным условиям (22) и начальным условиям (23) и (24), будем искать методом Фурье (называемым также методом разделения переменных). Метод Фурье состоит в том, что частные решения ищутся в виде произведения двух функций, одна из которых зависит только от x, а другая только от t. То есть мы ищем решения уравнения (21), которые имеют специальный вид: u(x, t) = X(x)T(t), (25) где X дважды непрерывно дифференцируемая функция от x на [, l], а T дважды непрерывно дифференцируемая функция от t, t >. Substituting (25) into (21), we get: X T = α 2 X T, (26) or T (t) α 2 T(t) = X (x) X(x). (27) It is said that there has been a separation of variables. Since x and t do not depend on each other, the left side in (27) does not depend on x, but the right side does not depend on t, and the total value of these ratios is 62
63 must be constant, which we denote by λ: T (t) α 2 T(t) = X (x) X(x) = λ. From this we get two ordinary differential equations: X (x) λx(x) =, (28) T (t) α 2 λt(t) =. (29) In this case, the boundary conditions (22) take the form X()T(t) = and X(l)T(t) =. Since they must be fulfilled for all t, t >, then X() = X(l) =. (3) Let us find solutions to equation (28) satisfying boundary conditions (3). Let's consider three cases. Case 1: λ >. Denote λ = β 2. Equation (28) takes the form X (x) β 2 X(x) =. Its characteristic equation k 2 β 2 = has roots k = ±β. Consequently, common decision equation (28) has the form X(x) = C e βx + De βx. We must choose the constants C and D so that the boundary conditions (3) are met, i.e. X() = C + D =, X(l) = C e βl + De βl =. Since β, then this system of equations has a unique solution C = D =. Hence X(x) and 63
64 u(x, t). Thus, in case 1 we have obtained a trivial solution, which we will not consider further. Case 2: λ =. Then equation (28) takes the form X (x) = and its solution is obviously given by the formula: X(x) = C x+d. Substituting this solution into the boundary conditions (3), we obtain X() = D = and X(l) = Cl =, hence C = D =. Hence X(x) and u(x, t), and we again have a trivial solution. Case 3: λ<. Обозначим λ = β 2. Уравнение (28) принимает вид: X (x)+β 2 X(x) =. Его характеристическое уравнение имеет вид k 2 + β 2 =, а k = ±βi являются его корнями. Следовательно, общее решение уравнения (28) в этом случае имеет вид X(x) = C sin βx + D cosβx. В силу граничных условий (3) имеем X() = D =, X(l) = C sin βl =. Поскольку мы ищем нетривиальные решения (т. е. такие, когда C и D не равны нулю одновременно), то из последнего равенства находим sin βl =, т. е. βl = nπ, n = ±1, ±2,..., n не равно нулю, так как сейчас мы рассматриваем случай 3, в котором β. Итак, если β = nπ (nπ) 2, l, т. е. λ = то существуют l решения X n (x) = C n sin πnx, (31) l C n произвольные постоянные, уравнения (28), не равные тождественно нулю. 64
65 In what follows, we will assign to n only positive values n = 1, 2,..., since for negative n, solutions of the same form (nπ) will be obtained. The values λ n = are called eigenvalues, and the functions X n (x) = C n sin πnx eigenfunctions of differential equation (28) with boundary conditions (3). Now let's solve equation (29). For him, the characteristic equation has the form k 2 α 2 λ =. (32) l 2 Since we found out above that nontrivial solutions X(x) of Eq. (28) exist only for negative λ equal to λ = n2 π 2, it is these λ that we will consider below. The roots of equation (32) are k = ±iα λ, and the solutions of equation (29) have the form: T n (t) = A n sin πnαt + B n cos πnαt, (33) l l where A n and B n are arbitrary constants. Substituting formulas (31) and (33) into (25), we find particular solutions of equation (21) that satisfy boundary conditions (22): (u n (x, t) = B n cos πnαt + A n sin πnαt) C n sin pnx. l l l Entering the factor C n in brackets and introducing the notation C n A n = b n and B n C n = a n, we write u n (X, T) as (u n (x, t) = a n cos πnαt + b n sin πnαt) sin pnx. (34) l l l 65
66 The vibrations of the string corresponding to the solutions u n (x, t) are called natural vibrations of the string. Since equation (21) and boundary conditions (22) are linear and homogeneous, then a linear combination of solutions (34) (u(x, t) = a n cos πnαt + b n sin πnαt) sin πnx (35) l l l will be a solution to equation (21 ) satisfying the boundary conditions (22) with a special choice of the coefficients a n and b n, which ensures the uniform convergence of the series. Now we choose the coefficients a n and b n of solution (35) so that it satisfies not only the boundary conditions, but also the initial conditions (23) and (24), where f(x), g(x) are given functions (moreover, f() = f (l) = g() = g(l) =). We assume that the functions f(x) and g(x) satisfy the Fourier expansion conditions. Substituting the value t = into (35), we obtain u(x,) = a n sin πnx l = f(x). Differentiating series (35) with respect to t and substituting t =, we obtain u t (x,) = πnα b n sin πnx l l = g(x), and this is the expansion of the functions f(x) and g(x) into Fourier series. Therefore, a n = 2 l l f(x) sin πnx l dx, b n = 2 l g(x) sin πnx dx. πnα l (36) 66
67 Substituting the expressions for the coefficients a n and b n into series (35), we obtain a solution to equation (21) that satisfies boundary conditions (22) and initial conditions (23) and (24). Thus, we have solved the problem of free small transverse vibrations of a string. Let us clarify the physical meaning of the eigenfunctions u n (x, t) of the problem of free vibrations of a string, defined by formula (34). Let us rewrite it as where u n (x, t) = α n cos πnα l α n = a 2 n + b2 n, (t + δ n) sin πnx, (37) l πnα δ n = arctg b n. l a n Formula (37) shows that all points of the string perform harmonic oscillations with the same frequency ω n = πnα and phase πnα δ n. The oscillation amplitude depends on l l the abscissa x of the string point and is equal to α n sin πnx. With such an oscillation, all points of the string simultaneously reach their l maximum deviation in one direction or another and simultaneously pass the equilibrium position. Such oscillations are called standing waves. A standing wave will have n + 1 fixed points given by the roots of the equation sin πnx = in the interval [, l]. The fixed points are called the nodes of the standing wave. In the middle between the nodes - l mi are the points at which the deviations reach a maximum; such points are called antinodes. Each string can have its own oscillations of strictly defined frequencies ω n = πnα, n = 1, 2,.... These frequencies are called natural frequencies of the string. The lowest l tone that a string can produce is determined by itself 67
68 low natural frequency ω 1 = π T and is called the fundamental tone of the string. The remaining tones corresponding to l ρ frequencies ω n, n = 2, 3,..., are called overtones or harmonics. For clarity, we will depict the typical profiles of a string emitting the fundamental tone (Fig. 33), the first overtone (Fig. 34) and the second overtone (Fig. 35). Rice. Fig. 33. Profile of the string that emits the fundamental tone. Fig. 34. Profile of a string emitting the first overtone. Fig. 35. Profile of a string emitting a second overtone If the string performs free vibrations determined by the initial conditions, then the function u(x, t) is represented, as can be seen from formula (35), as a sum of individual harmonics. Thus arbitrary oscillation 68
The 69th string is a superposition of standing waves. In this case, the nature of the sound of the string (tone, sound strength, timbre) will depend on the ratio between the amplitudes of individual harmonics. Strength, pitch and timbre of the sound A vibrating string excites air vibrations perceived by the human ear as a sound emitted by a string. The strength of sound is characterized by the energy or amplitude of vibrations: the greater the energy, the greater the strength of the sound. The pitch of a sound is determined by its frequency or period of oscillation: the higher the frequency, the higher the sound. The timbre of sound is determined by the presence of overtones, the distribution of energy over harmonics, i.e., the method of excitation of vibrations. The amplitudes of the overtones are, generally speaking, less than the amplitude of the fundamental, and the phases of the overtones can be arbitrary. Our ear is not sensitive to the phase of oscillations. Compare, for example, the two curves in Fig. 36, borrowed from . This is a recording of sound with the same fundamental tone, extracted from the clarinet (a) and the piano (b). Both sounds are not simple sinusoidal oscillations. The fundamental frequency of the sound in both cases is the same and this creates the same tone. But the curve patterns are different because different overtones are superimposed on the fundamental tone. In a sense, these drawings show what timbre is. 69
Equations of hyperbolic type. Vibrations of an infinite and semi-infinite string. Fourier method Fourier method Standing waves 4 Lecture 4.1 Hyperbolic type equations. Fluctuations of infinite and semi-infinite
MOSCOW STATE TECHNICAL UNIVERSITY OF CIVIL AVIATION V.M. Lyubimov, E.A. Zhukova, V.A. Ukhova, Yu.A. Shurinov
MINISTRY OF EDUCATION AND SCIENCE OF RUSSIA Federal State Budgetary Educational Institution of Higher Professional Education MATI Russian State Technological University named after K. E. Tsiolkovsky
Ministry of Education of the Republic of Belarus Vitebsk State Technological University Topic. "Rows" Department of Theoretical and Applied Mathematics. developed by Assoc. E.B. Dunina. Main
Federal Agency for Education Federal State Educational Institution of Higher Professional Education SOUTHERN FEDERAL UNIVERSITY R. M. Gavrilova, G. S. Kostetskaya Methodical
Topic Fourier series Practical lesson Fourier series in orthogonal systems of functions Space of piecewise continuous functions Generalized Fourier series 3 Bessel inequality and convergence of Fourier series Space
THEORY OF SERIES The theory of series is the most important component of mathematical analysis and finds both theoretical and numerous practical applications. Distinguish between numerical and functional series.
CONTENTS Fourier series 4 The concept of a periodic function 4 Trigonometric polynomial 6 3 Orthogonal systems of functions 4 Trigonometric Fourier series 3 5 Fourier series for even and odd functions 6 6 Decomposition
Federal Agency for Education Moscow State University of Geodesy and Cartography (MIIGAiK) METHODOLOGICAL INSTRUCTIONS AND TASKS FOR INDEPENDENT WORK on the course HIGHER MATHEMATICS
Lecture 4. Harmonic analysis. Fourier series Periodic functions. Harmonic analysis In science and technology, one often has to deal with periodic phenomena, i.e., those that repeat through
TOPIC V FOURIER SERIES LECTURE 6 Expansion of a periodic function in a Fourier series Many processes occurring in nature and technology have the properties to repeat at certain intervals Such processes
METHODOLOGICAL INSTRUCTIONS FOR CALCULATION TASKS ON THE COURSE OF HIGHER MATHEMATICS "ORDINARY DIFFERENTIAL EQUATIONS SERIES DOUBLE INTEGRALS" PART III THEME SERIES Contents Series Numerical series Convergence and divergence
6 Fourier series 6 Orthogonal systems of functions Fourier series in terms of an orthogonal system of functions The functions ϕ () and ψ (), defined and integrable on the segment [, ], are called orthogonal on this segment if
DEFINITE INTEGRAL. Integral Sums and Definite Integral Let a function y = f () defined on the segment [, b ], where< b. Разобьём отрезок [, b ] с помощью точек деления на n элементарных
5 Power series 5 Power series: definition, domain of convergence Function series of the form (a + a) + a () + K + a () + K a) (, (5) numbers are called power series Numbers
BELARUSIAN STATE UNIVERSITY FACULTY OF APPLIED MATHEMATICS AND INFORMATION SCIENCE Department of Higher Mathematics Teaching aid for students of the Faculty of Applied Mathematics and Informatics
Let's look at some examples. Example. Let's find the sum of an infinite geometric progression The formula for the common term of this series is a+aq+...+aq n +... (a). a n = aq n. Let us calculate its partial sums. If q =, then
Task 1.1. Find solutions y = y(x) of the differential equation that are non-identically zero in the indicated area and satisfy the given boundary conditions (the Sturm-Liouville problem) Solution: Consider
Mathematical analysis Topic: Definite integral Improper integrals Lecturer Pakhomova E.G. 2017 CHAPTER II. Definite integral and its applications 1. Definite integral and its properties 1. Tasks,
Lecture 8 4 Sturm-Liouville problem
Explanations to the text: the sign is read as "equivalent" and means that the equations to the right of the sign and to the left of the sign have the same set of solutions, the sign IR denotes the set of real numbers, the sign IN
82 4. Section 4. Functional and power series 4.2. Lesson 3 4.2. Lesson 3 4.2.. Taylor expansion of a function DEFINITION 4.2.. Let the function y = f(x) be infinitely differentiable in some neighborhood
MINISTRY OF EDUCATION AND SCIENCE OF RUSSIA FEDERAL STATE BUDGETARY EDUCATIONAL INSTITUTION OF HIGHER PROFESSIONAL EDUCATION "SAMARA STATE TECHNICAL UNIVERSITY" Department of Applied Mathematics
Federal Agency for Railway Transport Ural State University of Railway Transport Department "Higher and Applied Mathematics" N. P. Chuev Elements of Harmonic Analysis Methodical
Lecture 3 Taylor and Maclaurin series Application of power series Expansion of functions into power series Taylor and Maclaurin series For applications, it is important to be able to expand a given function into a power series, those functions
S A Lavrenchenko wwwwrckoru Lecture Fourier transform Concept of integral transformation The method of integral transformations is one of the powerful methods of mathematical physics and is a powerful solution
Integrability of a function (according to Riemann) and a definite integral Examples of problem solving 1. The constant function f(x) = C is integrable on , since for any partitions and any choice of points ξ i the integral
I course, task. Prove that the Riemann function, if 0, m m R(), if, m, m 0, and the fraction is irreducible, 0, if irrational, is discontinuous at every rational point and continuous at every irrational one. Solution.
1 2 Table of contents 1 Fourier series 5 1.1 Trigonometric Fourier series .................. 5 1.2 Only sin & cos ......................... 7 1.3 Fourier series in complex form............. 11 1.4 f(x) = c k?...............
EQUATIONS OF MATHEMATICAL PHYSICS 1. Partial differential equations
Lecture 4. Wave equations 1. Derivation of the equation of string vibrations 2. Equation of longitudinal vibrations of a rod 3. Initial conditions, boundary conditions 4. Problem statement 1. Derivation of the equation of string vibrations
1. Electrostatics 1 1. Electrostatics Lesson 6 Separation of variables in Cartesian coordinates 1.1. (Problem 1.49) The z = plane is charged with density σ (x, y) = σ sin (αx) sin (βy), where σ, α, β are constants.
Module Topic Function sequences and series Properties of uniform convergence of sequences and series Power series Lecture Definitions of function sequences and series Uniformly
Equations of parabolic type. Method of separation of variables Homogeneous boundary value problem Source function Inhomogeneous heat equation 7 Lecture 7.1 Equations of parabolic type. Separation method
Lecture Numerical series Signs of convergence Number series Signs of convergence An infinite expression of a numerical sequence + + + +, composed of members of an infinite one, is called a numerical series
35 7 Trigonometric Fourier series Fourier series for periodic functions with period T. Let f(x) be a piecewise continuous periodic function with period T. Consider the basic trigonometric system
Faculty of Metallurgy Department of Higher Mathematics
Department of Mathematics and Informatics Elements of Higher Mathematics Educational and methodological complex for students of secondary vocational education studying using distance technologies Module Differential calculus Compiled by:
9. Antiderivative and indefinite integral 9.. Let the function f() be given on the interval I R. The function F () is called the antiderivative function f() on the interval I, if F () = f() for any I, and the antiderivative
DIFFERENTIATION OF FUNCTIONS OF A ONE VARIABLE The concept of a derivative, its geometric and physical meaning Problems leading to the concept of a derivative Definition of Tangent S to the line y f (x) at point A x ; f(
Equations of hyperbolic type. Vibrations of an infinite and semi-infinite string. d'Alembert's method Infinite string. d'Alembert formula Semi-infinite string 3 Lecture 3.1 Hyperbolic type equations.
Title Introduction. Basic concepts.... 4 1. Volterra integral equations... 5 Homework options.... 8 2. Resolvent of the Volterra integral equation. 10 Homework options.... 11
ROWS. Number lines. Basic definitions Let an infinite sequence of numbers be given The expression (infinite sum) a, a 2,..., a n,... a i = a + a 2 + + a n +... () i= is called a number series. Numbers
8. Power series 8.. A functional series of the form c n (z) n, (8.) n= where c n is a numerical sequence, R is a fixed number, and z R is called a power series with coefficients c n. By changing the variables
~ ~ Indefinite and definite integrals The concept of antiderivative and indefinite integral. Definition: A function F is called an antiderivative with respect to a function f if these functions are related as follows
3724 SERIES OF MULTIPLE AND CURVILINEAR INTEGRALS 1 WORKING PROGRAM OF SECTIONS "SERIES OF MULTIPLE AND CURVILINEAR INTEGRALS" 11 Number series The concept of a number series Properties of number series A necessary criterion for convergence
EAT. ORE MATHEMATICAL ANALYSIS. NUMERICAL AND FUNCTIONAL SERIES NOVOSIBIRSK 200 2 RUSSIAN MINISTRY OF EDUCATION AND SCIENCE SEI HPE "NOVOSIBIRSK STATE PEDAGOGICAL UNIVERSITY" E.M. Rudoy MATHEMATICAL ANALYSIS.
LECTURE N 7 .Power
QUADRATIC EQUATIONS
SECTION OF TASKS WITH PARAMETERS Comment Tasks with parameters are traditionally complex tasks in the USE structure, requiring the applicant not only to master all the methods and techniques for solving various
Differential calculus Introduction to mathematical analysis Sequence and function limit. Disclosure of uncertainties within. Function derivative. Differentiation rules. Application of the derivative
Fourier series Orthogonal systems of functions From the point of view of algebra, the equality where are functions of a given class and are coefficients from R or C simply means that the vector is a linear combination of vectors B
1. Definite integral 1.1. Let f be a bounded function defined on the segment [, b] R. A partition of the segment [, b] is a set of points τ = (x, x 1,..., x n 1, x n ) [, b] such that = x< x 1 < < x n 1
Ch Power series a a a A series of the form a a a a a () is called a power series, where, a, are constants, called the coefficients of the series. Sometimes a power series of a more general form is considered: a a (a) a (a) a (a) (), where
Fourier series is a representation of an arbitrarily taken function with a specific period as a series. In general terms, this solution is called the decomposition of an element in an orthogonal basis. The expansion of functions in a Fourier series is a fairly powerful tool for solving various problems due to the properties of this transformation when integrating, differentiating, as well as shifting an expression in an argument and convolution.
A person who is not familiar with higher mathematics, as well as with the works of the French scientist Fourier, most likely will not understand what these “series” are and what they are for. Meanwhile, this transformation has become quite dense in our lives. It is used not only by mathematicians, but also by physicists, chemists, physicians, astronomers, seismologists, oceanographers and many others. Let us also take a closer look at the works of the great French scientist, who made a discovery ahead of his time.
Man and the Fourier Transform
Fourier series is one of the methods (along with analysis and others) This process occurs every time a person hears any sound. Our ear automatically transforms elementary particles in an elastic medium, they are decomposed into rows (along the spectrum) of successive values of the volume level for tones of different heights. Next, the brain turns this data into sounds familiar to us. All this happens in addition to our desire or consciousness, by itself, but in order to understand these processes, it will take several years to study higher mathematics.
More on the Fourier Transform
The Fourier transform can be carried out by analytical, numerical and other methods. Fourier series refer to the numeral way of decomposing any oscillatory processes - from ocean tides and light waves to cycles of solar (and other astronomical objects) activity. Using these mathematical techniques, it is possible to analyze functions, representing any oscillatory processes as a series of sinusoidal components that go from minimum to maximum and vice versa. The Fourier transform is a function that describes the phase and amplitude of sinusoids corresponding to a specific frequency. This process can be used to solve very complex equations that describe dynamic processes that occur under the influence of thermal, light or electrical energy. Also, Fourier series make it possible to isolate the constant components in complex oscillatory signals, which made it possible to correctly interpret the obtained experimental observations in medicine, chemistry and astronomy.
History reference
The founding father of this theory is the French mathematician Jean Baptiste Joseph Fourier. This transformation was subsequently named after him. Initially, the scientist applied his method to study and explain the mechanisms of heat conduction - the spread of heat in solids. Fourier suggested that the original irregular distribution can be decomposed into the simplest sinusoids, each of which will have its own temperature minimum and maximum, as well as its own phase. In this case, each such component will be measured from minimum to maximum and vice versa. The mathematical function that describes the upper and lower peaks of the curve, as well as the phase of each of the harmonics, is called the Fourier transform of the temperature distribution expression. The author of the theory reduced the general distribution function, which is difficult to describe mathematically, to a very convenient series of cosine and sine, which sum up to give the original distribution.
The principle of transformation and the views of contemporaries
The scientist's contemporaries - the leading mathematicians of the early nineteenth century - did not accept this theory. The main objection was Fourier's assertion that a discontinuous function describing a straight line or a discontinuous curve can be represented as a sum of sinusoidal expressions that are continuous. As an example, consider Heaviside's "step": its value is zero to the left of the gap and one to the right. This function describes the dependence of the electric current on the time variable when the circuit is closed. The contemporaries of the theory at that time had never encountered such a situation, when a discontinuous expression would be described by a combination of continuous, ordinary functions, such as an exponential, sinusoid, linear or quadratic.
What confused French mathematicians in Fourier theory?
After all, if the mathematician was right in his statements, then by summing the infinite trigonometric Fourier series, one can obtain an exact representation of the stepwise expression even if it has many similar steps. At the beginning of the nineteenth century, such a statement seemed absurd. But despite all the doubts, many mathematicians have expanded the scope of the study of this phenomenon, taking it beyond the scope of studies of thermal conductivity. However, most scientists continued to be tormented by the question: "Can the sum of the sinusoidal series converge to the exact value of the discontinuous function?"
Fourier Series Convergence: An Example
The question of convergence is raised whenever it is necessary to sum infinite series of numbers. To understand this phenomenon, consider a classic example. Can you ever reach the wall if each successive step is half the size of the previous one? Suppose you are two meters from the goal, the first step brings you closer to the halfway point, the next one to the three-quarters mark, and after the fifth step you will cover almost 97 percent of the way. However, no matter how many steps you take, you will not achieve the intended goal in a strict mathematical sense. Using numerical calculations, it can be shown that in the end it is possible to approach an arbitrarily small given distance. This proof is equivalent to demonstrating that the total value of one-half, one-fourth, etc. will tend to one.
A Question of Convergence: The Second Coming, or Lord Kelvin's Appliance
This question was raised again at the end of the nineteenth century, when Fourier series were tried to be used to predict the intensity of ebb and flow. At this time, Lord Kelvin invented a device, which is an analog computing device that allowed sailors of the military and merchant fleet to track this natural phenomenon. This mechanism determined the sets of phases and amplitudes from a table of tide heights and their corresponding time moments, carefully measured in a given harbor during the year. Each parameter was a sinusoidal component of the tide height expression and was one of the regular components. The results of the measurements were entered into Lord Kelvin's calculator, which synthesized a curve that predicted the height of the water as a function of time for the next year. Very soon similar curves were drawn up for all the harbors of the world.
And if the process is broken by a discontinuous function?
At that time, it seemed obvious that a tidal wave predictor with a large number of counting elements could calculate a large number of phases and amplitudes and thus provide more accurate predictions. Nevertheless, it turned out that this regularity is not observed in those cases when the tidal expression to be synthesized contained a sharp jump, that is, it was discontinuous. In the event that data is entered into the device from the table of time moments, then it calculates several Fourier coefficients. The original function is restored thanks to the sinusoidal components (according to the found coefficients). The discrepancy between the original and restored expression can be measured at any point. When carrying out repeated calculations and comparisons, it can be seen that the value of the largest error does not decrease. However, they are localized in the region corresponding to the discontinuity point, and tend to zero at any other point. In 1899, this result was theoretically confirmed by Joshua Willard Gibbs of Yale University.
Convergence of Fourier series and the development of mathematics in general
Fourier analysis is not applicable to expressions containing an infinite number of bursts in a certain interval. In general, Fourier series, if the original function is the result of a real physical measurement, always converge. Questions of the convergence of this process for specific classes of functions have led to the emergence of new sections in mathematics, for example, the theory of generalized functions. It is associated with such names as L. Schwartz, J. Mikusinsky and J. Temple. Within the framework of this theory, a clear and precise theoretical basis was created for such expressions as the Dirac delta function (it describes an area of a single area concentrated in an infinitely small neighborhood of a point) and the Heaviside “step”. Thanks to this work, Fourier series became applicable to solving equations and problems in which intuitive concepts appear: a point charge, a point mass, magnetic dipoles, and also a concentrated load on a beam.
Fourier method
Fourier series, in accordance with the principles of interference, begin with the decomposition of complex forms into simpler ones. For example, a change in heat flow is explained by its passage through various obstacles made of irregularly shaped heat-insulating material or a change in the surface of the earth - an earthquake, a change in the orbit of a celestial body - the influence of planets. As a rule, similar equations describing simple classical systems are elementarily solved for each individual wave. Fourier showed that simple solutions can also be summed to give solutions to more complex problems. Expressed in the language of mathematics, Fourier series is a technique for representing an expression as the sum of harmonics - cosine and sinusoids. Therefore, this analysis is also known as "harmonic analysis".
Fourier series - the ideal technique before the "computer age"
Before the creation of computer technology, the Fourier technique was the best weapon in the arsenal of scientists when working with the wave nature of our world. The Fourier series in a complex form allows solving not only simple problems that can be directly applied to the laws of Newton's mechanics, but also fundamental equations. Most of the discoveries of Newtonian science in the nineteenth century were made possible only by Fourier's technique.
Fourier series today
With the development of computers, Fourier transforms have risen to a qualitatively new level. This technique is firmly entrenched in almost all areas of science and technology. An example is a digital audio and video signal. Its realization became possible only thanks to the theory developed by a French mathematician at the beginning of the nineteenth century. Thus, the Fourier series in a complex form made it possible to make a breakthrough in the study of outer space. In addition, this influenced the study of the physics of semiconductor materials and plasma, microwave acoustics, oceanography, radar, and seismology.
Trigonometric Fourier series
In mathematics, a Fourier series is a way of representing arbitrary complex functions as a sum of simpler ones. In general cases, the number of such expressions can be infinite. Moreover, the more their number is taken into account in the calculation, the more accurate the final result is. Most often used as the simplest trigonometric functions cosine or sine. In this case, the Fourier series are called trigonometric, and the solution of such expressions is called the expansion of the harmonic. This method plays an important role in mathematics. First of all, the trigonometric series provides a means for the image, as well as the study of functions, it is the main apparatus of the theory. In addition, it allows solving a number of problems of mathematical physics. Finally, this theory contributed to the development and brought to life a number of very important sections of mathematical science (the theory of integrals, the theory of periodic functions). In addition, it served as a starting point for the development of the following functions of a real variable, and also marked the beginning of harmonic analysis.
Fourier series expansion of even and odd functions expansion of a function given on a segment into a series in terms of sines or cosines Fourier series for a function with an arbitrary period Complex representation of the Fourier series Fourier series in general orthogonal systems of functions Fourier series in an orthogonal system Minimal property of Fourier coefficients Bessel's inequality Equality Parseval Closed systems Completeness and closedness of systems
Fourier series expansion of even and odd functions The function f(x), defined on the segment \-1, where I > 0, is called even if the Graph of the even function is symmetrical about the y-axis. The function f(x) defined on the segment J, where I > 0, is called odd if the Graph of the odd function is symmetrical with respect to the origin. Example. a) The function is even on the segment |-jt, jt), since for all x e b) The function is odd, since the Fourier series expansion of even and odd functions is the expansion of a function given on the segment in a series of sines or cosines Fourier series for a function with an arbitrary period Complex notation of the Fourier series Fourier series in general orthogonal systems of functions Fourier series in an orthogonal system Minimum property of Fourier coefficients Bessel inequality Parseval equality Closed systems Completeness and closedness of systems c) Function f(x)=x2-x, where does not belong neither to even nor to odd functions, since Let the function f(x) satisfying the conditions of Theorem 1 be even on the segment x|. Then for all i.e. /(g) cos nx is even function, and f(x)sinnx is odd. Therefore, the Fourier coefficients of an even function /(x) will be equal Therefore, the Fourier series of an even function has the form 00 If f(x) - odd function on the segment [-mr, ir|, then the product f(x) cosnx will be an odd function, and the product f(x) sin nx will be an even function. Therefore, we will have Thus, the Fourier series of an odd function has the form We have Applying integration by parts twice, we get that Hence, the Fourier series of this function looks like this: or, in expanded form, This equality is valid for any x €, since at the points x = ±ir the sum of the series coincides with the values of the function f(x ) = x2, since the graphs of the function f(x) = x and the sums of the resulting series are given in fig. Comment. This Fourier series allows you to find the sum of one of the convergent numerical series, namely, for x \u003d 0, we get that The function /(x) satisfies the conditions of Theorem 1, therefore it can be expanded into a Fourier series, which, due to the oddness of this function, will have the form Integrating by parts, we find the Fourier coefficients Therefore, the Fourier series of this function has the form This equality holds for all x В points x - ±tg the sum of the Fourier series does not coincide with the values of the function / (x) = x, since it is equal to Outside the segment [- *, n-] the sum of the series is a periodic continuation of the function / (x) \u003d x; its graph is shown in Fig. 6. § 6. Expansion of a function given on an interval into a series in terms of sines or cosines. Let a bounded piecewise monotonic function / be given on an interval . The values of this function on the interval 0| can be defined in various ways. For example, it is possible to define the function / on the segment mc] in such a way that /. In this case it is said that) "is extended to the segment 0] in an even way"; its Fourier series will contain only cosines. If, however, the function /(x) is defined on the segment [-x, mc] so that /(, then an odd function is obtained, and then we say that / "is extended to the segment [-*, 0] in an odd way"; in this In this case, the Fourier series will contain only sines.So, each bounded piecewise-monotone function /(x), defined on the segment , can be expanded into a Fourier series both in terms of sines and cosines.Example 1. Expand the function in a Fourier series: a) by cosines; b) along the sines. M This function, with its even and odd extensions to the segment |-x, 0) will be bounded and piecewise monotonic. a) We continue / (z) into the segment 0) a) We continue j \ x) into the segment (-m, 0 | in an even way (Fig. 7), then its Fourier series i will have the form P \u003d 1 where the Fourier coefficients are equal, respectively for Therefore, b) Let's continue /(z) in the segment [-x,0] in an odd way (Fig. eight). Then its Fourier series §7. Fourier Series for a Function with an Arbitrary Period Let the function fix) be periodic with a period of 21.1 ^ 0. To expand it into a Fourier series on the interval where I > 0, we make a change of variable by setting x = jt. Then the function F(t) = / ^tj will be a periodic function of the argument t with a period and it can be expanded on a segment in a Fourier series Returning to the variable x, i.e., setting, we obtain , remain valid for periodic functions with an arbitrary period 21. In particular, sufficient sign expandability of a function in a Fourier series. Example 1. Expand in a Fourier series a periodic function with a period of 21, given on the segment [-/,/] by the formula (Fig. 9). Since this function is even, its Fourier series has the form Substituting the found values of the Fourier coefficients into the Fourier series, we obtain Note one thing important property periodic functions. Theorem 5. If a function has a period T and is integrable, then for any number a the equality m holds. i.e. the integral on a segment whose length is equal to the period T has the same value regardless of the position of this segment on the real axis. Indeed, We make a change of variable in the second integral, assuming This gives and therefore, Geometrically, this property means that in the case of the area shaded in Fig. 10 areas are equal to each other. In particular, for a function f(x) with a period, we obtain at the Fourier series expansion of even and odd functions the expansion of a function given on a segment into a series in terms of sines or cosines Fourier series for a function with an arbitrary period Complex representation of the Fourier series Fourier series in general orthogonal systems functions Fourier series in an orthogonal system Minimal property of Fourier coefficients Bessel inequality Parseval equality Closed systems Completeness and closedness of systems that the Fourier coefficients of a periodic function f(x) with a period of 21 can be calculated using the formulas where a is an arbitrary real number(note that the cos - and sin functions have a period of 2/). Example 3. Expand in a Fourier series a function given on an interval with a period of 2x (Fig. 11). 4 Find the Fourier coefficients of this function. Putting in the formulas we find that for Therefore, the Fourier series will look like this: At the point x = jt (discontinuity point of the first kind) we have §8. Complex notation of the Fourier series In this section, some elements of complex analysis are used (see chapter XXX, where all the operations performed here with complex expressions , are strictly justified). Let the function f(x) satisfy sufficient conditions for expansion into a Fourier series. Then on the segment x] it can be represented by a series of the form Using the Euler formulas Substituting these expressions into the series (1) instead of cos nx and sin xy we will have We introduce the following notation Then the series (2) takes the form Thus, the Fourier series (1) is presented in the complex form (3). Let us find expressions for the coefficients in terms of integrals. We have Similarly, we find Finally, the formulas for с„, с_п and с can be written as follows: . . The coefficients cn are called the complex Fourier coefficients of the function For a periodic function with a period), the complex form of the Fourier series takes the form values w if limits exist Example. Expand the period function into a complex Fourier series This function satisfies sufficient conditions for expansion into a Fourier series. Let Find the complex Fourier coefficients of this function. We have for odd for even n, or, in short. Substituting the values), we finally obtain Note that this series can also be written as follows: Fourier series in general orthogonal systems of functions 9.1. Orthogonal Systems of Functions Denote by the set of all (real) functions that are square-defined and integrable on the interval [a, 6], i.e., those for which there exists an integral. In particular, all functions f(x) that are continuous on the interval [a , 6], belong to 6], and the values of their Lebesgue integrals coincide with the values of the Riemann integrals. Definition. The system of functions, where, is called orthogonal on the interval [a, b\, if Condition (1) assumes, in particular, that none of the functions is identically equal to zero. The integral is understood in the sense of Lebesgue. and we call the quantity the norm of a function. If in an orthogonal system for any n we have, then the system of functions is called orthonormal. If the system (y>n(x)) is orthogonal, then the system Example 1. A trigonometric system is orthogonal on a segment. The system of functions is an orthonormal system of functions on, Example 2. The cosine system and the sine system is orthonormal. Let us introduce the notation that they are orthogonal on the segment (0, f|, but not orthonormal (for I ↦ 2). Since their norms are COS that the functions form an orthonormal system of functions on a segment Let us show, for example, the orthogonality of the Legendre polynomials. Let m > n. In this case, integrating n times by parts, we find since for the function t/m = (z2 - I)m all derivatives up to order m - I inclusive vanish at the ends of the interval [-1,1). Definition. The system of functions (pn(x)) is called orthogonal on the interval (a, b) by overhang p(x) if: 1) there are integrals for all n = 1,2,... Here it is assumed that the weight function p(x) is defined and positive everywhere on the interval (a, b), with the possible exception of a finite number of points where p(x) can vanish. After performing differentiation in formula (3), we find. It can be shown that the Chebyshev-Hermite polynomials are orthogonal on the interval Example 4. The system of Bessel functions (jL(pix)^ is orthogonal on the interval of zeros of the Bessel function system Let an orthogonal system of functions in the interval (a, 6) and let the series (cj = const) converge on this interval to the function f(x): By virtue of the orthogonality of the system, we obtain that This operation has, generally speaking, a purely formal character. However, in some cases, for example, when the series (4) converges uniformly, all functions are continuous and the interval (a, 6) is finite, this operation is legal. But it is the formal interpretation that is important for us now. So let's say a function is given. We form the numbers c * according to the formula (5) and write The series on the right side is called the Fourier series of the function f (x) with respect to the system (^n (n)) - The numbers Cn are called the Fourier coefficients of the function f (x) in this system. The sign ~ in formula (6) only means that the numbers Cn are related to the function f(x) by formula (5) (in this case, it is not assumed that the series on the right converges at all, much less converges to the function f(x)). Therefore, the question naturally arises: what are the properties of this series? In what sense does it "represent" the function f(x)? 9.3. Average Convergence Definition. A sequence converges to an element ] on average if the norm is in space Theorem 6. If a sequence ) converges uniformly, then it also converges on average. M Let the sequence ()) converge uniformly on the segment [a, b] to the function f(x). This means that for any, for all sufficiently large n, we have Hence, from which our assertion follows. The converse is not true: the sequence () can converge on average to /(x), but not be uniformly convergent. Example. Let us consider the sequence nx It is easy to see that But this convergence is not uniform: there exists e, for example, such that no matter how large n is, on the segment Fourier series for a function with an arbitrary period Complex representation of the Fourier series Fourier series in general orthogonal systems of functions Fourier series in an orthogonal system Minimal property of Fourier coefficients Bessel inequality Parseval equality Closed systems Completeness and closedness of systems and let ) in the orthonormal system b Consider a linear combination where n ^ 1 is a fixed integer, and find the values of the constants for which the integral takes its minimum value. Let us write it in more detail Integrating term by term, due to the orthonormality of the system, we obtain The first two terms on the right side of equality (7) are independent, and the third term is nonnegative. Therefore, the integral (*) takes on a minimum value at ak = sk. The integral is called the root-mean-square approximation of the function f(x) as a linear combination of Tn(x). Thus, the root-mean-square approximation of the function /\ takes on a minimum value when. when Tn(x) is the 71st partial sum of the Fourier series of the function /(x) in the system (. Setting ak = ck, from (7) we obtain Equality (9) is called the Bessel identity. Since its left side is non-negative, then from it Bessel's inequality follows Since i is arbitrary here, Bessel's inequality can be represented in a strengthened form, i.e., for any function /, the series of squared Fourier coefficients of this function in an orthonormal system ) converges. Since the system is orthonormal on the segment [-x, r], then inequality (10) translated into the usual notation of the trigonometric Fourier series gives the relation do valid for any function f(x) with an integrable square. If f2(x) is integrable, then due to necessary condition convergence of the series on the left side of inequality (11), we obtain that. Parseval's equality For some systems (^n(x)) the inequality sign in formula (10) can be replaced (for all functions f(x) 6 x) by an equals sign. The resulting equality is called the Parseval-Steklov equality (completeness condition). The Bessel identity (9) allows us to write condition (12) in an equivalent form by the space norm 6]. Definition. An orthonormal system ( is called complete in b2[ay b] if any function can be approximated with any accuracy on the average by a linear combination of the form with a sufficiently large number of terms, i.e. if for any function f(x) ∈ b2[a, b\ and for any e > 0 there is natural number nq and numbers a\, a2y..., such that No The above reasoning implies Theorem 7. If, by orthonormalization, the system ) is complete in space, the Fourier series of any function / in this system converges to f(x) on the average, i.e., by the norm It can be shown that the trigonometric system is complete in space. This implies the assertion. Theorem 8. If a function /0 its trigonometric Fourier series converges to it on the average. 9.5. closed systems. Completeness and closedness of systems Definition. An orthonormal system of functions \, is called closed if in the space Li\a, b) there is no non-zero function orthogonal to all functions. In the space L2\a, b\ the concepts of completeness and closedness of orthonormal systems coincide. Exercises 1. Expand the function in the Fourier series in the interval (-i-, x) 2. Expand the function in the Fourier series in the interval (-r, r) 3. Expand the function in the Fourier series in the interval (-r, r) 4. Expand in a Fourier series in the interval (-jt, r) function 5. Expand in a Fourier series in the interval (-r, r) the function f (x) \u003d x + x. 6. Expand in a Fourier series in the interval (-jt, r) the function n 7. Expand in a Fourier series in the interval (-r, x) the function / (x) \u003d sin2 x. 8. Expand in a Fourier series in the interval (-m, jt) the function f(x) = y 9. Expand in a Fourier series in the interval (-mm, -k) the function f(x) = | sinx|. 10. Expand in a Fourier series in the interval (-x-, r) the function f(x) = g. 11. Expand in a Fourier series in the interval (-r, r) the function f (x) \u003d sin §. 12. Expand in a Fourier series the function f (x) = n -2x, given in the interval (0, x), continuing it in the interval (-x, 0): a) in an even way; b) in an odd way. 13. Expand in a Fourier series in terms of sines the function / (x) \u003d x2, given in the interval (0, x). 14. Expand in a Fourier series the function / (x) \u003d 3-x, given in the interval (-2,2). 15. Expand in a Fourier series the function f (x) \u003d |x |, given in the interval (-1,1). 16. Expand in a Fourier series in terms of sines the function f (x) \u003d 2x, specified in the interval (0,1).
One of the types of functional series is the trigonometric series
The task is to choose the coefficients of the series so that it converges to a function given in the interval [-π, π]; in other words, it is required to expand the given function into a trigonometric series. A sufficient condition for the solvability of this problem is that the function be piecewise continuous and piecewise differentiable in the interval [-π, π], i.e., that the interval [-π, π] can be divided into a finite number of partial intervals, in each of which the given function is continuous and has a derivative (at the ends of the partial intervals, the function must have finite one-sided limits and one-sided derivatives, in the calculation of which its one-sided limit is taken as the value of the function at the end of the partial interval). The condition of piecewise differentiability can be replaced by the condition of piecewise monotonicity of the function, i.e., the requirement that the function be monotonic in each of the partial intervals. A sufficient condition for the expansion of a function in the interval [-π, π] into a trigonometric series is also the requirement that the function has a limited change in this interval. By the definition of the function f(x) has a bounded change in an interval if, for any division of this interval into a finite number of intervals
magnitude
bounded above by the same number.
It is with such functions that one has to deal with in solving practical problems.
When performing any of the three sufficient conditions the function f(x) is represented in the interval [-π, π] by a trigonometric series whose coefficients are determined by the formulas
With such coefficients, the trigonometric series is called near Fourier. This series converges to f(x) at every point of its continuity; at breakpoints, it converges to the arithmetic mean of the left and right limit values, i.e. k, if x is a breakpoint (Fig. 1); at the boundaries of the segment, the series converges to .
Picture 1.
The function expressed by the Fourier series is a periodic function, and therefore the series compiled for the function given on the segment [-π, π] converges outside this segment to a periodic continuation of this function (Fig. 2).
Figure 2.
If the Fourier series represents the function f(x), given in an arbitrary interval [α, α+2π] of length 2π, then the coefficients of the series a 0 , a k , b k (Fourier coefficients) can be determined by the indicated formulas, in which the limits of integration are replaced by α and α+2π. In general, since the formulas for a 0 , a k , b k contain functions with a period of 2π, integration can be carried out over any interval with a length of 2π.
The Fourier series can be used for an approximate representation of the function, namely: the function f(x) is replaced by the sum s n (x) of the first few terms of the Fourier series, which is approximately equal to it:
The expression s n (x), where a 0 , a k , b k are the Fourier coefficients of the function f(x), in comparison with other expressions of the same form with the same value of n, but with different coefficients, leads to the minimum standard deviation s n (x ) of f(x), which is defined as
Some simplifications are possible depending on the kind of symmetry of the function. If the function is even, i.e. f(-x)=f(x), then
and the function expands into a series in cosines. If the function is odd, i.e. f(-x)=-f(x), then
and the function expands into a series in terms of sines. If the function satisfies the condition f(x+π)=-f(x), i.e., the curve referring to half of the segment of length 2π is a mirror image of the other half of the curve, then
The function can be defined not only on a segment of length 2π, but also on a segment of any length 2l. If it satisfies the above conditions on this segment, then it can be expanded into a Fourier series of the following form:
where the coefficients of the series are calculated by the formulas
In table. 1 expansions of some functions are given.
Table 1.
The trigonometric series can also be written in the following form:
The Fourier series of the function f(x) converges the faster, the smoother the function is. If the function f (x) and its derivatives f "(x), f" (x), ..., f k -1 (x) are everywhere continuous, and f (k) (x) allows only discontinuity points of the 1st kind in a finite number, then the Fourier coefficients a n , b n of the function f (x) will be
The symbol denotes a value such that
The expansion into a trigonometric series is called harmonic analysis, and the trigonometric functions included in this series are called harmonics. The calculation of the component harmonics is called harmonic synthesis.
When calculating structures, it is often necessary to expand in a Fourier series various functions given by graphs, and above all, representing the load. In table. 2 and 3, expansions are given for some functions characteristic of loads, including series corresponding to concentrated forces.
Table 2.
Function Graph |
Fourier series |
n |