본문 바로가기
Physics/수리물리

푸리에 해석

by hydrogendeuteride 2021. 1. 29.

테일러 급수를 아는 사람들은 알겠지만, 어떤 임의의 함수는 멱급수 형태로 전개할 수 있다.

그럼 반대로 생각해보자 어떤 임의의 주기함수는 사인함수의 합으로 이루어질 수 있는가?

 

어떤 함수 f(x)f(x)

12a0+a1cosx+a2cos2x+a3cos3x...+b1sinx+b2sin2x+b3sin3x...12a0+a1cosx+a2cos2x+a3cos3x...+b1sinx+b2sin2x+b3sin3x...

이렇게 주기함수의 합으로 나타낼 수 있다고 가정해보자

 

이런 식을 구하기 위해서는 한 주기에 대한 삼각함수의 평균값에 대해 알아야 한다.


sinmxcosnxsinmxcosnx의 평균값은

12πππsinmxcosnxdx=012πππsinmxcosnxdx=0

 

sinmxsinnxsinmxsinnx의 평균값은

12πππsinmxsinnxdx=0(mn),12(m=n0),0(m=n=0)12πππsinmxsinnxdx=0(mn),12(m=n0),0(m=n=0)

 

cosmxcosnxcosmxcosnx의 평균값은

12πππcosmxcosnxdx=0(mn),12(m=n0),1(m=n=0)12πππcosmxcosnxdx=0(mn),12(m=n0),1(m=n=0)

 

여기서 두 주기함수의 곱이 1212인 것을 직교한다고 한다.(위의 것을 함수의 내적아라고도 한다 나중에 르장드르 다항식 설명할 기회가 되면 그때 설명하겠다)


 

그럼 여기서 a0a0의 값을 구할 수 있다.

12a0+a1cosx+a2cos2x+a3cos3x...+b1sinx+b2sin2x+b3sin3x...12a0+a1cosx+a2cos2x+a3cos3x...+b1sinx+b2sin2x+b3sin3x...

여기서 사인이나 코사인같은 삼각함수의 주기에 대한 평균값은 0이기 때문에 12a012a0을 제외한 보든 항은 0이다.

 

에서 양변의 (π,π)에서의 평균값을 구하기 위해 양변을 적분하면 된다.

f(x)=12a0+a1cosx+a2cos2x+a3cos3x...+b1sinx+b2sin2x+b3sin3x...

 

12πππf(x)dx=

12a012πππdx+a112πππcosxdx+a212πππcos2xdx+a312πππcos3xdx...

+b112πππsinxdx+b212πππsin2xdx+b312πππsin3xdx...

 

그럼 a0/212πππf(x)dx 이다.

 

이걸 다른 항들에 적용해보자.

 

f(x)cos(x)를 곱하고 평균을 구해보자

12πππf(x)cosxdx=

12a012πππcosxdx+a112πππcos2xdx+a212πππcos2xcosxdx+a312πππcos3xcosxdx...

+b112πππsinxcosxdx+b212πππsin2xcosxdx+b312πππsin3xcosxdx...

이다.

 

위에서 한 주기에 대한 주기함수의 곱의 평균값을 이용하여 위 식의 값을 구해보면

12πππf(x)cosxdx=a112πππcos2xdx

 

a1에 대해 풀면

1πππf(x)cosxdx=a1

이다.

 

다른 항들에 대해 똑같은 과정은 거치면 다음을 알 수 있다.

12πππf(x)cosnxdx=an,12πππf(x)sinnxdx=bn

 

이것이 푸리에 급수이다.

 

복소수 형태로는 위와 비슷하게 전개할 수 있다.

f(x)=cneinx

cn=12πππf(x)einxdx


Dirichlet 조건

어떤 주기함수를 항상 푸리에 급수로 전개할 수 있는 것은 아니다.

1.한 주기 내의 모든 함숫값이 (1,1)이고

2.최댓값과 최솟값이 유한개이고

3.불연속점이 유한 개여야 한다.


Parseval의 정리

f(x)의 평균값은 (12a0)2+121a2n+121b2n이다.

따라서 위의 식에서 한 개의 항이라도 빠진다면 평균값이 달라지게 되고, 푸리에 급수로 표현 불가능해진다.


푸리에 변환

푸리에 급수를 복소수 형태로 쓰면

f(x)=cneiπnx l

인데 여기서 항들의 합을 적분으로 바꿀 수 있다

 

그러면

f(x)=cneiπnx l

*여기서 cn=12lllf(x)einπx ldx은 푸리에 급수의 계수로 위해서 구하기 위해 유도했었고, 복소수 형태로 다시 썼다.

 

그러면

f(x)=g(α)eiaxdα

g(α)=12πf(x)eiαxdx

(g(a)는 푸리에 계수인cn이다)

 

그러면 이제 적분으로 바뀌었기 때문에 특정 구간에서만이 아니라 전체 구간에서의 푸리에 계수를 구할 수 있다.


아주아주 간단히 설명하고 넘어간 부분이 있지만, 대충 이정도면 된 것 같다.

예시를 들어보자

 

f(x)=1(π<x<0),1(0<x<π)인 함수를 전개하면

f(x)=4π(sinx1+sin3x3+sin5x5...)이다. (위의 푸리에 급수에 넣어서 한번 해 보자)

 

이걸  그려보면

대충 이렇게 되겠다.

항의 수가 늘어나면 더 직선에 가까워지지만, 여기서는 잘 보이게 하기 위해 11개의 항만 사용하였다.

 

동영상을 봐 보자 코드는 파이썬으로 작성하였다. 곧 공개하겠다.

'Physics > 수리물리' 카테고리의 다른 글

변분법  (0) 2021.01.15