Processing math: 100%
 

오일러의 공식 (1) 멱급수, 매클로린 급수

반응형

오일러에 대해서 찾아보다가 멱급수를 연구하였다고 하기에 정리해보고자 합니다.

 

etc-image-0
위키백과

저 빨간 박스 안의 식이 어떻게 유도되는지 보겠습니다.

 


멱급수

 

ax2+bx+c

이 식은 많이 보던 다항식입니다.

이러한 다항식을

 

a0+a1x+a2x2+a3x3+...

이와 같이 다항식을 일반화 한 급수를 멱급수(power series)라고 합니다.

* 급수: 수열의 합

 

예를 들어,

1+2x+3x2+4x3+...

이 식은 멱급수입니다.

 


매클로린 급수

 

다음 함수를 생각해봅니다.

f(x)=a0+a1x+a2x2+a3x3+...

 

위 함수의 x에 0을 대입하면

f(0)=a0

 

f(x)를 미분하면

f(x)=a1+2a2x+3a3x2+4a4x3+...

 

위 함수의 x에 0을 대입하면

f(0)=a1

 

f(x)를 미분하면

f(x)=2a2+32a3x+43a4x2+54a5x3+...

 

위 함수의 x에 0을 대입하면

f(0)=2a2

a2=f(0)2

 

f(x)를 미분하면

f(x)=32a3+432a4x+543a5x2+654a6x3+...

 

위 함수의 x에 0을 대입하면

f(0)=32a3

a3=f(0)32

 

a0=f(0),a1=f(0),a2=f(0)2,a3=f(0)32 이므로

a4=f(4)(0)432 일 것이라 알 수 있습니다.

 

이것을 일반화시키면

0!=1,1!=1이므로

an=f(n)(0)n!입니다.

 

따라서 f(x)는 다음과 같은 전개식을 가질 수 있습니다.

 

f(x)=f(0)0!+f(0)1!x+f(0)2!x2+f(0)3!x3+...

 

이 식을 f에 대한 매클로린 급수라고 합니다.

 


sinx을 매클로린 급수로 나타내기

 

sinx를 매클로린 급수로 나타내봅시다.

 

f(x)=f(0)0!+f(0)1!x+f(0)2!x2+f(0)3!x3+...

 

f(x)=sinxf(0)=0

f(x)=cosxf(0)=1

f(x)=sinxf(0)=0

f(x)=cosxf(0)=1

f(4)=sinxf(4)(0)=0

이므로

 

sinx=00!+11!x+02!x2+13!x3+04!x4+15!x5+06!x6+17!x7+...

 

sinx=x13!x3+15!x517!x7+...

위 식으로 나타낼 수 있습니다.

 

 

그래프로 확인을 해보겠습니다.

(그래프를 그리는 프로그램은 옥타브(Octave)를 사용하였습니다.)

 

f(x)=x

 

x = [-2*pi:0.1:2*pi]
y = sin(x)
z1 = x
plot(x,y,x,z1,'r')
legend('sin(x)','x','Location','northwest')
axis([-6, 6, -20, 20])

 

z1.png

 

 

f(x)=xx33!

 

x = [-2*pi:0.1:2*pi]
y = sin(x)
z2 = x-x.^3./factorial(3)
plot(x,y,x,z2,'r')
legend('sin(x)','x-(x^3/3!)','Location','northwest')
axis([-6, 6, -20, 20])

 

z2.png

 

f(x)=xx33!+x55!

 

x = [-2*pi:0.1:2*pi]
y = sin(x)
z3 = x-x.^3./factorial(3)+x.^5./factorial(5)
plot(x,y,x,z3,'r')
legend('sin(x)','x-(x^3/3!)+(x^5/5!)','Location','northwest')
axis([-6, 6, -20, 20])

 

z3.png

 

 

f(x)=xx33!+x55!x77!

 

x = [-2*pi:0.1:2*pi]
y = sin(x)
z4 = x-x.^3./factorial(3)+x.^5./factorial(5)-x.^7./factorial(7)
plot(x,y,x,z2,'r')
legend('sin(x)','x-(x^3/3!)+(x^5/5!)-(x^7/7!)','Location','northwest')
axis([-6, 6, -20, 20])

 

z4.png

 

점점 sinx 그래프에 가까워짐을 알 수가 있습니다.

 

x = [-2*pi:0.1:2*pi]
y = sin(x)
z1 = x
z2 = x-x.^3./factorial(3)
z3 = x-x.^3./factorial(3)+x.^5./factorial(5)
z4 = x-x.^3./factorial(3)+x.^5./factorial(5)-x.^7./factorial(7)
plot(x,y,x,z1,'r-',x,z2,'g--',x,z3,'m:',x,z4,'k-.')
legend('sin(x)','x','x-(x^3/3!)','x-(x^3/3!)+(x^5/5!)','x-(x^3/3!)+(x^5/5!)-(x^7/7!)','Location','northwest','NumColumns',2)
axis([-6, 6, -20, 20])

 

화면 캡처 2022-12-14 220712.png


같은 방법으로 매클로린 급수로 

cosxex를 나타내면

 

cosx=1x22!+x44!x66!+...

ex=1+x+x22!+x33!+x44!+...

 

로 나타낼 수 있습니다.

 

 

다음 글에서는 위에서 설명한 수식들을 이용하여 '세상에서 가장 아름다운 수식'이라 불리는

 

eiπ+1=0

 

오일러 항등식을 유도해보겠습니다.

반응형