1
1
---
2
2
layout : single
3
- title : " (1) The Genesis of Fourier Analysis"
3
+ title : " 1 The Genesis of Fourier Analysis"
4
4
categories : mathematics
5
5
tags : [analysis, Fourier analysis, Fourier Series, differential equations, wave equation, heat equation, simple harmonic motion, Dirichlet problem, Laplacian, Fourier, Euler, d'Alembert]
6
6
use_math : true
@@ -74,7 +74,7 @@ BUTTON을 사용했다 (파랑색).
74
74
75
75
<div class =" notice--info " markdown =" 1 " >
76
76
1절에서 다루는 내용은 파동방정식(wave equation)이다.
77
- 먼저 가장 기본이 되는 simple harmonic motion이 언급된다.
77
+ 먼저 simple harmonic motion이 언급된다.
78
78
파동(wave), 조금 더 정확하게는, 양 끝이 고정된 줄의 모양은, 수평변위($x$)와 시간($t$)에 대한 수직변위($y$)의 함수 $y=u(x,t)$로 나타내어질 수 있다.
79
79
이 파동을 정상파(standing wave)와 traveling wave의 두 가지 타입으로 모델링하고 있는 것으로 보인다.
80
80
standing wave는 주어진 함수 $u(x,t)$의 변수 $x$, $t$를 분리하여 $u(x,t)=\phi(x)\psi(t)$로 두는 방식으로 파동이 주어진 파동모양을 기준으로 진폭만 시간에 따라 바뀐다는 생각에 기반한다.
@@ -134,9 +134,9 @@ This wave moves rightward if $c\gt0$ where $c$ can be thought of as the velocity
134
134
135
135
<div class =" notice--info " markdown =" 1 " >
136
136
1.1절에서는 파동방정식을 유도해본다.
137
- 줄을 공간적으로 infinitesimal하게 나누고, 시간적으로도 바로 다음 순간의 상태를 고려한다.
137
+ 줄을 공간적으로 infinitesimal하게 나누고, 시간적으로도 infinitesimal하게 바로 다음 순간의 상태를 고려한다.
138
138
줄의 한 부분(particle)이 받는 힘들을 계산하고 그 알짜힘에 의해 다음 순간에 바뀌는 위치를 계산하여 방정식을 만들 수 있다.
139
- 이때 공간적인 간격 $\Delta x$와 시간적인 간격 $\Delta t$를 모두 0으로 보냄으로써 파동방정식 $\frac{\partial^2u}{\partial t^2}=\frac{\partial^2u}{\partial x^2}$이 유도된다.
139
+ 이때 공간적인 간격 $\Delta x$와 시간적인 간격 $\Delta t$를 모두 0으로 보내면 파동방정식 $\frac{\partial^2u}{\partial t^2}=\frac{\partial^2u}{\partial x^2}$이 유도된다.
140
140
</div >
141
141
142
142
Consider a horizontal string of length $L$.
@@ -179,6 +179,8 @@ We analyze and solve the standard one and simply apply change of variables for t
179
179
180
180
Consider a physical problem of a string where $0\le x\le \pi$ :
181
181
182
+ <div class =" notice--success " markdown =" 1 " >
183
+
182
184
$$
183
185
\begin{aligned}
184
186
\frac{\partial^2 u}{\partial t^2}=\frac{\partial^2 u}{\partial x^2}
@@ -193,11 +195,13 @@ u(0,t)=u(\pi,t)=0
193
195
\tag a
194
196
$$
195
197
198
+ </div >
199
+
196
200
### 1.2 Solution to the wave equation
197
201
198
202
<div class =" notice--info " markdown =" 1 " >
199
203
1.2절에서는 파동방정식의 해를 구하는데 traveling wave와 정상파(standing wave)의 두 가지 방법을 사용해서 해본다.
200
- 조금 더 정확하게는 경계조건(boundary condition)과 초기조건(파동의 초기모양, 파동의 초기 속도)까지 포함된 문제에서의 해를 구하는데
204
+ 조금 더 정확하게는 경계조건(boundary condition)과 초기조건(initial condition, 파동의 초기모양, 파동의 초기 속도)까지 포함된 문제에서의 해를 구하는데
201
205
traveling wave와 정상파(standing wave)의 두 가지 방법을 사용해서 해본다.
202
206
<br ><br >
203
207
traveling wave의 방법은 $u(x,t)$가 적절한 함수 $F$, $G$에 대하여 $u(x,t)=F(x+t)+G(x-t)$로 표현될 수 있음에 기반한다.
@@ -207,7 +211,7 @@ traveling wave의 방법은 $u(x,t)$가 적절한 함수 $F$, $G$에 대하여 $
207
211
<br ><br >
208
212
정상파의 방법은 파동방정식의 양변을 각각 분리된 두 변수 $x$와 $t$로만 표현시키고 이를 simple harmonic motion에서 나온 미분방정식으로 만들어서 푼다.
209
213
이를 통해 보면 사인과 코사인 (혹은 지수함수쌍)으로 이루어진 함수들이 특정한 해인 것을 유도할 수 있고, 이 해들로 중첩한 것이 일반해인 것을 가정한다.
210
- 중첩에 사용되는 계수들(Fourier coefficient)은 초기조건으로부터 특정할 수 있는데, 이때 간단한 적분을 통해 얻을 수 있다.
214
+ 중첩에 사용되는 계수들(Fourier coefficient)은 초기조건으로부터 특정할 수 있는데, 간단한 적분을 통해 얻을 수 있다.
211
215
이러한 모든 이야기가 성립하기 위해서는, 어떤 함수 $f$가 주어졌을 때, 그 함수를 sine series와 cosine series 로 표현하는 것이 가능해야 한다.
212
216
이 문제는 이 책 전반에서 다루어질 중요한 질문이라는 점이 강조되어 언급된다.
213
217
이러한 정상파의 방법은 traveling wave의 방법보다 더 일반적이고 더 많은 활용을 가능하게 하며, 푸리에(Fourier)가 1807년에 발전시켰다.
@@ -222,16 +226,20 @@ We can get two kinds of solution to the wave equation ;
222
226
223
227
We claim that (a) has the general solution of the form
224
228
229
+ <div class =" notice--success " markdown =" 1 " >
230
+
225
231
$$
226
232
u(x,t) = \frac12\left[f(x+t)+f(x-t)\right]+\frac12\int_{x-t}^{x+t}g(y)\,dy.
227
233
\tag b
228
234
$$
229
235
236
+ </div >
237
+
230
238
Note that $f$ and $g$ were a function only on the interval $[ 0,\pi] $ in (a).
231
239
Now in the above solution (b), $f$ and $g$ are defined everywhere in $\mathbb R$.
232
240
That is, we conceive $f$ and $g$ as an extended versions in (b).
233
241
From (a), extend the domain $[ 0,\pi] $ to $[ -\pi, \pi] $ by defining $f(-x)=-f(x)$ and $g(-x)=-g(x)$ and then extend the domains to $\mathbb R$ by defining $f(x+2\pi)=f(x)$ and $g(x+2\pi)=g(x)$.
234
- That is, $f$ and $g$ are periodic functions of period $2\pi$.
242
+ That is, $f$ and $g$ are odd functions of period $2\pi$.
235
243
In order for $f$ and $g$ to be continuous, we need the conditions $f(0)=f(\pi)=0$ and $g(0)=g(\pi)=0$, which are guaranteed to be true from the original problem (a).
236
244
237
245
Before proving that (b) is a general solution of (a), i.e. (a) $\Rightarrow$ (b), let's check if (b) is a specific solution of (a), i.e. (b) $\Rightarrow$ (a).
@@ -259,6 +267,7 @@ For the boundary conditions, first note that
259
267
$f(\pi-s)=-f(-\pi+s)=-f(\pi+s)$
260
268
and
261
269
$g(\pi-s)=-g(-\pi+s)=-g(\pi+s)$ for real $s$.
270
+ That is, $y=f(x)$ is symmetric not only about the origin, but also about $(\pi,0)$
262
271
Then
263
272
264
273
$$
@@ -309,10 +318,11 @@ Equating $\frac{\partial^2 u}{\partial x^2}$ and $\frac{\partial^2 u} {\partial
309
318
310
319
$$ \frac{\partial^2 u}{\partial\xi\partial\eta}=0 $$
311
320
312
- Write $u(x,t)=v(\xi,\eta)$, this is
321
+ Write $u(x,t)=v(\xi,\eta)$, then
313
322
314
323
$$ \frac{\partial}{\partial\xi}\frac{\partial}{\partial\eta}v(\xi,\eta)=0 $$
315
324
325
+ for all $\xi$ and $\eta$
316
326
For any $\xi$, there exists $\zeta$ between $0$ and $\xi$ such that
317
327
318
328
$$
358
368
\begin{aligned}
359
369
F(x)
360
370
&=\int_0^x\frac12\left[f'(y)+g(y)\right]\,dy+C_1\\
361
- &=\frac12f(x)+\frac12\int_0^xg(y)\,dy+C_1' \\
371
+ &=\frac12f(x)+\frac12\int_0^xg(y)\,dy+C_1\\
362
372
G(x)
363
373
&=\int_0^x\frac12\left[f'(y)-g(y)\right]\,dy+C_2\\
364
- &=\frac12f(x)-\frac12\int_0^xg(y)\,dy+C_2' .
374
+ &=\frac12f(x)-\frac12\int_0^xg(y)\,dy+C_2.
365
375
\end{aligned}
366
376
$$
367
377
368
- Note that $C_1' +C_2' =0$ since $f(x)=F(x)+G(x)$.
378
+ Note that $C_1+C_2=0$ since $f(x)=F(x)+G(x)$.
369
379
Therefore,
370
380
371
381
$$
372
382
\begin{aligned}
373
383
u(x,t)
374
384
&=F(x+t)+G(x-t)\\
375
- &=\frac12f(x+t)+\frac12\int_0^{x+t}g(y)\,dy+C_1'
376
- +\frac12f(x-t)-\frac12\int_0^{x-t}g(y)\,dy+C_2' \\
385
+ &=\frac12f(x+t)+\frac12\int_0^{x+t}g(y)\,dy+C_1
386
+ +\frac12f(x-t)-\frac12\int_0^{x-t}g(y)\,dy+C_2\\
377
387
&=\frac12\left[f(x+t)+f(x-t)\right]+\frac12\int_{x-t}^{x+t}g(y)\,dy.
378
388
\end{aligned}
379
389
$$
382
392
383
393
We claim that (a) has the general solution of the form
384
394
395
+ <div class =" notice--success " markdown =" 1 " >
396
+
385
397
$$ u(x,t)=\sum_{m=1}^\infty\sin mx\left(A_m\cos mt+B_m\sin mt\right)\tag c $$
386
398
387
399
where
@@ -393,6 +405,7 @@ B_n&=\frac2{n\pi}\int_0^\pi g(x)\sin nx\,dx\\
393
405
\end{aligned}
394
406
$$
395
407
408
+ </div >
396
409
397
410
Separating variables $x$ and $t$, let
398
411
427
440
$$
428
441
429
442
by 3.6, where $m^2=-\lambda$.
430
- Apply the fourth condition of (a).
443
+ We exckyde the case when $\lambda\lt0$ since $e^{rx}$ and $e^{-rx}$ are not oscillatory (See Appendix).
444
+ Apply the boundary condition of (a).
431
445
Since the string has zero displacement at $x=0$,
432
446
433
447
$$ 0=\tilde A\left(A\cos mt + B\sin mt\right) $$
@@ -444,22 +458,27 @@ For every integer $m$,
444
458
445
459
$$ u_m(x,t)=\sin mx\left(A\cos mt + B\sin mt\right) $$
446
460
447
- is a solution for the wave equation.
461
+ is a solution to the wave equation with the boundary condition .
448
462
449
463
<!-- Let
450
464
451
465
$$u_m(x,t)=\sin mx\left(A_m\cos mt + B_m\sin mt\right)$$ -->
452
466
453
467
<div class =" notice--info " markdown =" 1 " >
454
468
<b >Remark </b > <br >
455
- Letting, for example, $A_m=1$, $B_m=0$,
469
+ Let, for example, $A_m=1$, $B_m=0$.
470
+
456
471
$$ u_1(x,t)=\cos t\sin x $$
472
+
457
473
is called the fundamental tone or the first harmonic.
474
+
458
475
$$ u_2(x,t)=\cos2t\sin2x $$
476
+
459
477
is called the first overtone or the second harmonic.
478
+
460
479
$$ u_3(x,t)=\cos3t\sin3x $$
461
- is called second overtone or the third harmonic, and so on.
462
480
481
+ is called second overtone or the third harmonic, and so on.
463
482
For fixed $t$, $u_m(x,t)=0$ if $x=\frac\pi mk$ for all $k\in\mathbb Z$.
464
483
For each $k$, $x=\frac\pi mk$ is called a node, while $x=\frac\pi m\left(k+\frac12\right)$ is called an antinode.
465
484
</div >
@@ -469,8 +488,8 @@ Thus,
469
488
470
489
$$ u(x,t)=\sum_{m=-\infty}^\infty\sin mx\left(A_m\cos mt+B_m\sin mt\right) $$
471
490
472
- also satisfy the wave equation.
473
- Suppose that this solution is a general one (we omit the proof here.)
491
+ also satisfy the wave equation with the boundary condition .
492
+ Suppose that this solution is the general one.
474
493
Since $u_m=0$ for $m=0$ and
475
494
476
495
$$
512
531
\end{aligned}
513
532
$$
514
533
515
- Thus, (Fourier sine coeffcient of $f$)
534
+ Thus,
516
535
517
536
$$
518
537
\begin{aligned}
@@ -521,9 +540,9 @@ B_n&=\frac2{n\pi}\int_0^\pi g(x)\sin nx\,dx.
521
540
\end{aligned}
522
541
$$
523
542
524
- And the claim is now proved.
543
+ These are called Fourier sine coeffcient of $f$ and the claim is now proved.
525
544
526
- Note that this argument assumes that a * reasonable* function $f$ on $[ 0,\pi] $ can be expressed as a linear combination of sine functions ;
545
+ Note that the above arguments assume that a * reasonable* function $f$ on $[ 0,\pi] $ can be expressed as a sine series ;
527
546
528
547
$$
529
548
f(x)=\sum_{m=1}^\infty A_m\sin mx\tag{$\ast$}
557
576
558
577
By Euler identity, (we can choose $a_m=\frac{A_m'-A_mi}2$ and $a_ {-m}=\frac{A_m'+A_mi}2$)
559
578
560
- <div class =" notice--success " markdown =" 1 " >
561
579
$$
562
580
F(x) = \sum_{m=-\infty}^\infty a_me^{imx}.
563
581
$$
564
582
565
- </div >
566
583
567
584
Multiplying $e^{-inx}$ and integrating over $[ -\pi, \pi] $,
568
585
@@ -972,13 +989,10 @@ $$u(r,\theta)=r^{|m|}\left(A_m e^{im\theta}+B_me^{-im\theta}\right)$$
972
989
is a solution to $\Delta u=0$ for each integer $m$.
973
990
Superposing this class of solutions we can write
974
991
975
- <div class =" notice--success " markdown =" 1 " >
976
992
$$
977
993
u(r,\theta)=\sum_{m=-\infty}^\infty r^ma_me^{im\theta}.
978
994
$$
979
995
980
- </div >
981
-
982
996
The boundary condition $f(\theta)=u(1,\theta)$ becomes
983
997
984
998
$$
0 commit comments