무한서고

점화식에서의 특성방정식(characteristic equation)

포스팅 개요

과거 피보나치 수열의 일반항 구하는 포스팅(>>피보나치 수열의 일반항과 비율의 극한(황금비)<<)을 작성하였다.

사실 작성하던 중에는 크게 못느꼈는데, 다른 사람에게 설명을 하던 중 특성방정식을 잠깐 빌려와서 근과 계수와의 관계로 풀어내는 과정에서 '왜 $ a_{n+2} = a_{n+1} + a_{n} \Leftrightarrow x^2 = x + 1 $ 인가?'에 대해서 너무나도 당연하게 받아들였다는 것을 깨닫고 추가로 더 공부해본 결과 이를 '특성방정식'이라고 한다는 것을 알게되고, 이에 포스팅을 작성한다.

점화식에서의 특성방정식(characteristic equation)

점화식을 풀 때 우리는 특성방정식(characteristic equation)을 이용해서 해결을 하게 된다.
(과거 고등학교 수학에서 나왔던 점화식의 해결법(계차의 등비수열로 해결)도 어떻게보면 특성방정식의 활용이다.)

그런데 '왜 특성방정식을 사용해서 점화식을 푸는가?'에 대해서 궁금하진 않은가? 그냥 된다니까 하기에는 조금 껄쩍지근하다.

한줄로 정의해보자면
'점화식 자체로는 뭔가를 찾기 힘드니까 본질이 같은 다른 것으로 바꾸어서 해를 찾자'이다.

여기서 좀 더 자세히 따져보자면
'본질이 같은'은 '같은 선형성을 가지는'이란 의미이고
'다른 것'은 '기저(basis)'를 뜻한다.
(갑자기 대수학에서 벡터공간에서 쓰는 '기저'라는 단어가? 싶기도 하겠지만, 사실 모든 '함수'는 벡터공간 안에서 구현이 가능하다는 점을 보면 이해할 수 있을 것이다.)

진짜 결국 그냥 '쉽게 구할 수 있는 걸로 변형하자!'이거다..

보통은 점화식 뿐만 아니라 미분방정식, 선형대수에 모두 사용가능하기 때문에 '특성방정식'이라고 검색하면 사실 점화식보다는 미분방정식이나 선형대수 관련한 벡터공간 관련 내용 더 나아가 고유값/고유벡터 등이 나오게 된다.(사실 고유값분해를 대충이라도 이해할 수 있다면 특성방정식이 뭐하는 놈인지는 쉽게 이해가 간다)

그러나 일단은 이렇게 복잡한 내용 이전에 여기서는 아주 간단하게 점화식에 대해서만 설명해보기로 한다.

이제부터 이해해야 하는 키워드는

1. 점화식의 구분
 *선형성
 *선형의 또다른 의미
2. 특성방정식이란
이다.

차례대로 알아보자

1. 선형성(linearity)이란 무엇인가?

선형성을 만족시키기 위해서는 두 가지 조건인 가산성(Additivity)와 동차(제차)성(Homogeneity)을 만족해야한다.

-가산성(Additivity)은 f(a)+f(b)=f(a+b)를 만족하는 함수를 말한다.
다른말로는 중첩의 원리(principle of superposition)라고도 한다.

-동차(제차)성(Homogeneity)이란 f(ax)=af(x)를 만족하는 함수를 동차성이 있다고 한다.(일부 서적의 번역으로는 제차성이라고한다)
쉽게말해 입력이 조정된 비율만큼 동일하게 결과도 조정된 값이 나온다는 것이다.
선형성과 마찬가지로 상수항이 없는 함수에 대하여 성립하며, 상수항(혹은 이에 준하는 상수를 출력하는 함수)이 있는 경우 동차성에 위배된다.

간단하게 선형성과 동차성에 대해 쉽게 알 수 있는 예제가 있다.

y=ax+b는 선형인가?

가산성을 먼저 살펴보자
f(x1)+f(x2)=f(x1+x2)이면 가산성이 있는 함수이다.
a(x1+x2)+2b $ \neq $ a(x1+x2)+b
가산성에 위배된다.

동차성은 어떤지 보자
f(k x1) = kf(x1)
akx1+b $ \neq $ k(ax1+b)
동차성에 위배된다.

결국 y=ax+b는 선형이 아니다.


근데 신기한건 엄밀한 선형이 아닌데도 '선형'이라는 단어를 붙이는 경우가 있다는 것이다.
아니 또 이것은 무엇인가?
선형이 아닌데 선형이라고?

아까 선형성의 정의에서 '가산성'을 보았다. 여기서 파생되어서 덧셈으로 연결되는 함수들을 보고 '함수들이 선형 결합을 한다'라고 말한다.

그리고 이게 줄어들어서 '선형'이 된 것이다.

그래서 다항함수는? 선형 함수이다. 만약 상수항이 없다면 엄밀한 '선형'함수이고, 상수항이 있다면 선형(결합을 한)함수이다.

따라서 특정 계수의 곱을 통한 덧셈으로 정의되는 일반적인 점화식의 경우 선형이다.

그러면 이렇게 선형이라는 말을 막쓰면, '엄밀한 선형'인지 '일반적 선형'인지는 또 어떻게 알 것인가?

그래서 선형이란 말 뒤에 '동차' 혹은 '비동차'라는 말을 써준다.

이렇게 되면 '선형 동차 점화식' 혹은 '선형 비동차 점화식'이라는 말이 생겨날 수 있는데, 결국 이 단어로 구분이 완벽하게 되는 것이다.

선형 동차 점화식: 모든 항이 선형 결합을 한 상태이며, 상수항(혹은 그에 준하는 상수를 출력하는 함수)이 없는 점화식
선형 비동차 점화식: 모든 항이 선형 결합을 한 상태이며, 상수항(혹은 그에 준하는 상수를 출력하는 함수)이 있는 점화식

자 이제 점화식의 종류에 대해서 알아보았다.
사실 '비'자가 붙으면 뭐든 어려워진다. 선형함수보다 비선형함수가 어렵고, 동차보다는 비동차가 어렵다.
현재 이 포스팅은 '쉽게!' 알아보려는게 목적이므로 앞으로는 '선형 동차'에 대해서만 써보려고한다.
이제 '특성방정식'이란게 뭔지 알아보러가자.

2. 특성방정식이란?

선형 동차 점화식은 결국 선형성을 띈다.
그리고 이것은 점화식의 어떠한 해 a, b에 대하여 이 둘에 특정 계수를 곱한 선형결합역시도 점화식을 만족시킨다는 것이다.
결국 어떠한 해를 찾고 이 해에 곱해지는 계수를 찾으면 점화식을 풀어낼 수 있다(일반항을 구할 수 있다)는 것이다.

그럼 일단 이 '어떠한 해'를 찾아야하는데, 이게 그냥 점화식만 뚫어져라 쳐다보면 툭 답이 나오는 것도 아니고.. 참 힘들다.
그래서 이 어떠한 해를 찾기위해 우리는 주어진 점화식을 변형할 것이다.
어떻게?
'같지만 다르게!'
같은 선형결합을 가지지만, 이 점화식을 다른 각도로 볼 수 있는 새로운 '틀'(=기저)을 찾아서 바꿔주면 될 것이다.
결국 이 '틀'은 진짜 오만가지 것이 다 되지만, 제일 간단하며 우리가 무언가를 찾아내기 쉬운 틀은 $ x^n $일 것이다.
x를 찾아내면 점화식의 어떠한 해를 찾아낸 것이며, 이 어떠한 해에 특정 계수를 곱한 선형결합이 점화식을 만족시킬 것이기 때문이다.
그리고 새로운 틀을 $ x^n $으로 정의했으니, 이 틀을 이용해 만든 변화시킨 점화식에서는(해가 두개라고 가정하면) $ k_1 * x_{1}^n + k_2 * x_{2}^n $이 일반적인 해를 나타낸다고 볼 수 있겠다.
그렇다면 결국 $ a_n $이라는 수열을 새로운 틀 $ x^n $으로 놓게 된 것이다.

좀 더 이해를 쉽게하기위해서 하나의 예시를 가지고 논지를 진행시켜보자

$ a_n = a_{n-1} + a_{n-2} $

라는 점화식이 있다고하자(상수항이 없는 이유는? 우리는 현재 선형 동차 점화식에 대해서 보고있다.)
이 점화식에서는 차수(order)가 2이다.
갑자기 뜬금없이 차수? 선형 동차 점화식이 몇개의 항으로 연결되었는가를 나타내는 것으로, 선형적인 다항함수의 차수와 그 의미가 같다.

이를 새로운 틀 $ x^n $으로 치환하면
$ x^n = x^{n-1} + x^{n-2} $
차수로 n을 제한하면(n의 최대값은 차수가 된다)
$ x^2 = x + 1 $ 이 된다.
(위에서도 논의했듯이 항이 두개로 연결된 점화식은 두개의 해를 가지며, 두개의 해를 가지기 위해서는 다항함수의 차수=수열의 항수 이다)

자, 여기까지의 논지가 바로 이전 포스팅(피보나치 수열의 일반항 구하기)에서 근과 계수와의 관계를 사용하기 위해 살짝 빌려왔던 개념되겠다. 모르고 그냥 받아들여도 크게 문제없지만 알면 더 신기한 그런거다.

각설하고, 더 논지를 진행시켜보자.

$ x^2 - x -1 = 0 $의 형태로 바꾸어 방정식을 만들어주고

여기서 해를 구하면
$ x = \frac{1+\sqrt{5}}{2} \mathrm{ or } \frac{1-\sqrt{5}}{2} $
이 나오며

이 두개의 해의 특정 상수배씩의 선형결합이 이 점화식의 일반항이 된다. (더 나아가서, 같은 말로 '특정한 벡터'를 고유값(eigen value, 아이젠 밸류)과 고유벡터(eigen vector, 아이젠 벡터)로 분해하였을 때, 두개의 고유벡터에 특정 상수배(고유값)씩의 선형결합이 '특정한 벡터'가 된다는 것과 동일하다.)

$ a_n = k_1 * x_{1}^n + k_2 * x_{2}^n $

이때 특정 상수배($ k_1,  k_2 $)를 구하는 방법은, 초기 조건을 가지고 구하면 된다.

결국 $ a_0 = 0 $, $ a_1 = 1 $ 을 이용하면

$ a_0 = 0 = k_1 + k_2 $
$ a_1 = 1 = k_1 * (\frac{1+\sqrt{5}}{2})^1 + k_2 * (\frac{1-\sqrt{5}}{2})^1 $

$ k_1 = \frac{1}{\sqrt{5}} $, $k_2 = -\frac{1}{\sqrt{5}} $

따라서 이 점화식의 일반항은

$ a_n = \frac{1}{\sqrt{5}} * (\frac{1+\sqrt{5}}{2})^n - \frac{1}{\sqrt{5}} * (\frac{1-\sqrt{5}}{2})^n $

어디서 많이 본 일반항 아닌가?

맞다. 피보나치 수열의 일반항이다.

전 포스팅에서는 특성방정식을 그대로 이용하지 않고, 특성방정식의 아주 일부분만 잠깐 빌려다가 쓰고(특성 방정식에서 '틀'을 바꿔 수열을 잠깐 다항함수의 방정식 형태로 바꾼 뒤 근과 계수와의 공식으로 수열의 계수를 바꾼 정도) 계차의 공비를 구하는 식으로 점화식을 풀어내었다.

솔직히 전 포스팅에서는 이정도로 자세하게 특성방정식의 개념과 그 활용을 사용하지 않았기에 '일단 받아들여보세요~'하고 진행하였지만(사실 그 부분만 참고 넘어가면 이후에 유도하는데는 전혀 문제가 없다) 이후 좀 더 자세한 설명이 필요해보이기에 추가 포스팅한다.

직선 두개로 원을 삼등분하기

원을 삼등분하는 방법은 정말 여러가지가 있다.

일본에서는 챌린지까지 열릴정도라고하니 열기가 대단함을 알 수 있다.

그 중에는 실용적인 것도 있지만, 오로지 수학적 아름다움만을 위해 있는 것들도 있고... 하지만, 이번에 소개하려는 방식은 그 챌린지에 없지만, 실용적으로 쓸 수 있는 직선 두개로 원을 삼등분 하는 방식이다.

물론 수치가 완전히 딱 떨어지게는 나오지 않기 때문에 실전에서 사용하려면 어느정도 비슷한 근사치로 봐야하지만, 그래도 빠르게 직선만으로 3등분을 내기에는 적합하지 않나 싶다.

그래서 우리의 목표는

위와같이 원을 세로로 한번, 가로로 한번 잘라 면적이 동일하게 나오게 하는 것이다. 면적이 동일하면 3등분!

그렇지만, 따져보면 왼쪽의 두 부분은 위아래 대칭으로 한쪽 넓이만 계산하면 자동으로 반대편 넓이가 나오므로, 여기서 가장 중요한 포인트는 첫번째 세로로 자를 때 주황부분과 초록부분이 넓이가 같아지는지가 제일중요하게 된다.

그러면 어느 x점에서 잘라야 하는지가 중요해지고, 이것만 계산할 수 있으면 원을 직선 두개로 삼등분 할 수 있게 된다.

자, 그러면 두 넓이가 같다는 식을 세워보자

음영처리된 부분의 넓이를 구하는게 키포인트인데, 일단 간단하게 음영처리된 부분의 넓이를 $ a $라고 하고 넓이구하는 식을 만들어 봅시다. 원은 단위원으로 반지름은 1입니다.

(노란부분) $ \frac{\pi}{4} + a = 2(\frac{\pi}{4} - a) $ (초록부분)

근데 식을 써보니 식이 a에 대해서 정리될 것 같습니다.

$ \frac{\pi}{4} + a = 2(\frac{\pi}{4} - a) $

$ 3a = \frac{\pi}{4} $

$ a = \frac{\pi}{12} $

즉, 전체 넓이를 구하는 식에서 음영처리된 넓이가 $ \frac{\pi}{12} $면 된다는 식으로 식이 더 간단해 졌군요.

그럼 이제 우리는 음영처리된 넓이만 구하면 되겠습니다.

1) 각도와 삼각함수로 구하기

이전 포스팅(링크)에서 구했던 방식인데, 여기서 x값으로 다시 나타내지 않고, 각도값 그대로 사용하여 푸는 방식입니다.

원과 각도가 주어졌고, 위와 같은 모양의 넓이면, 호의 넓이+삼각형 넓이 해서 구할수 있으며, 이 식은

$ \frac{1}{2} \theta $ (호의 넓이) + $ \frac{1}{2} sin\theta \cdot cos\theta $ (삼각형의 넓이) [r=1 생략]

그러면 우리가 구하고자 했던 $ \frac{\pi}{12} $로 방정식을 놓으면

$ \frac{1}{2} \theta + \frac{1}{2} sin\theta \cdot cos\theta = \frac{\pi}{12} $

정리하면

$ \theta + sin\theta \cdot cos\theta = \frac{\pi}{6} $

가 됩니다. 전에도 말씀드렸다시피 삼각함수의 일반각 계산은 사람이 할 수 없으므로 울프람 알파(https://www.wolframalpha.com/)를 돌려줍니다.

[위 식을 그대로 긁어서 붙이면 울프람에서 인식을 못합니다. θ + sin θ ⋅ cos θ = π/6 이거를 붙여 넣어주세요]

그러면 결과가

$ \theta = 0.268133... $

으로 나옵니다.

그러나 이거는 각도 값일 뿐, 우리가 원하는 것은 x좌표이므로, 위의 그림에서 놓은 $ \theta $ 방향을 보면, $ sin \theta $가 x좌표임을 알 수 있습니다.

다시 $ sin \theta $에 0.268133값을 넣어서 울프람 알파를 돌리면

$ sin 0.268133 = 0.264932... = x $

이와 같은 값을 얻고, 이는 바로 x좌표입니다.

2) 부정적분으로 구하기

이전 포스팅(링크)에서 주야장천 구했던 식을 가지고 푸는 방법입니다.

물론 계산은 울프람 알파가 해줄겁니다.

먼저 방정식을 세웁니다.

$ \frac{1}{2} \cdot x \cdot \sqrt{1-x^{2}} + \frac{1}{2} \cdot arcsin \, x = \frac{\pi}{12} $

정리하여

$ x \cdot \sqrt{1-x^{2}} + arcsin \, x = \frac{\pi}{6} $

으로 만들고 울프람을 돌리겠습니다. 마찬가지로 위 식을 그대로 붙이면 울프람에서 인식 못하니 x ⋅ √ (1 − x^2) + arcsin x = π/6로 붙여넣어 주세요.

그러면 이번에는 한번에 x 값을 알려줍니다. 결과는 1)번과 같습니다.

결론

위의 계산으로 단위원의 원점에서 약 0.265만큼 떨어진 위치에서 세로로 한번, 나머지 부분에 가로로 한번 잘라주면 3등분이 된다는 사실을 알았네요!

좀 더 실용적으로 말해보자면 어떤 크기의 원이든, 원점에서 약 1/4 지점보다 조금 더 나가서 T자로 자르면 3분할이 된다고 볼 수 있겠습니다!

원넓이의 부정적분 구하기 - 2) 부정적분으로 구해보기

intetral root(1-x^2) dx

굉장히 오랜만에 다시 써보는 포스팅이네요.

저번 시간에 원 내부의 사다리꼴과 같은 도형의 넓이를 구하는 방법을 가장 기본적인 공식(부채꼴 공식+삼각형 공식)을 가지고 구해보았습니다.

이번 시간에는 이것을 부정적분으로 x값을 가지고 바로 구하는 방법을 알아볼까 합니다.

저번에 각도 $\theta$를 부채꼴 부분으로 잡았는데, 이번에도 한번 이렇게 잡아서 수식을 전개해보려고 합니다.

일단 이렇게 특수한 상황에 가기 전에, 일반적으로 원의 넓이를 적분으로 어떻게 구하는지 다시 한번 살펴보도록하겠습니다.

자, 일단 원의 방정식은 $x^{2}+y^{2}=r^{2}$입니다. 저희는 $r$이 1인 단위원을 사용하기 때문에 식은 더욱 간단하게 $x^{2}+y^{2}=1$이 되겠네요.

이를 좀 더 보기 편하게 y에 대한 값으로 나타내면(x의 값에 따라 y값을 결정하는 방식) $y=\pm \sqrt{1-x^{2}}$으로 정리할 수 있습니다.

이 때 부호에 따라 양의 부호는 y축을 기준으로 0보다 위에 그려지는 반원을, 음의 부호는 아래쪽에 그려지는 반원을 의미합니다.

현재 저희는 위에 그려지는 반원 중에서도 1사분면 위의 사반원에 대해 구하려고 하고 있으므로, 이에 대한 적분 수식은 $\int \sqrt{1-x^2} \, dx$라고 볼 수 있습니다.

여기서 루트가 들어간 적분은 그냥 풀기에 너무 힘들기 때문에 x를 치환시켜 줄 것입니다.

예전에 고등학교 때 적분을 공부하면서 도대체 왜 치환하는지 의문을 가졌었는데, 실상은 치환해서 더 쉬운 형태로 만들어서 적분을 쉽게 만들기 위해서 하는 작업입니다 치환적분은!

루트를 없애줄 수 있으면서 적분 형태를 간단하게 해줄 수 있는 것이 무엇이 있나 한번 살펴보다보니, 언뜻 지나가는 공식이 있습니다.

$ sin^{2} \theta + cos^{2} \theta = 1 $이라는 공식이지요.(이 공식은 그냥 암기할 게 아니라, 너무 당연한 것을 표현한 것입니다. 아까 단위원의 방정식은 $x^{2}+y^{2}=1$이라고 했습니다. 이것은 원 위에서 무조건 성립하는 값입니다. 여기서 매개변수 표현법을 사용하면 $y=sin \theta, x=cos \theta$라고 했습니다. 즉, 단위원의 방정식에 매개변수 표현법을 사용하여 표현 방법만 x, y 변수가 아닌 $\theta$변수 로 바꿔준 것이 됩니다.)

이항해보면

$cos^{2} \theta =1-sin^{2} \theta$

제곱을 제거하면

$cos \theta = \sqrt{1-sin^{2} \theta}$

어디선가 많이 본 보양이지요?

즉, $x$를 $sin \theta$로 치환하면 자연스럽게 루트가 들어간 식이 정리되면서 적분이 가능한 형태로 바뀔 것 같습니다!

일단, $x$를 치환하면 $dx$도 같이 치환해 줘야 하므로 미분을 때려 봅시다.

$x=sin \theta$

$dx=cos \theta d\theta$

그럼 이렇게 준비된 x를 가지고 치환적분을 해보겠습니다.

$ \int \sqrt{1-x^{2}} dx $

$ \int \sqrt{1-sin^{2} \theta} cos \theta d\theta]_{x=sin \theta, dx=cos \theta d\theta} $

$ \int cos \theta \cdot cos \theta d\theta $

$ \int cos^{2} \theta \, d\theta $

여기서 다시 난관에 봉착합니다. $ cos^{2} \theta $를 적분하려면 많은 애로사항이 꽃핍니다.

일단 제곱을 떨어내야하는데, 어떻게 떨어내야할지 생각해봤더니... 배각공식을 역이용해서 떨어보겠습니다.

$ cos 2\theta \, = \, cos^{2} \theta - sin^{2} \theta $

참고로 배각공식은 삼각함수의 덧셈공식에서 온겁니다

$ cos (\alpha+\beta) = cos \alpha \cdot cos \beta - sin \alpha \cdot sin \beta $

자, 일단 $ cos^{2} \theta $를 $ cos 2\theta $로 바꿀 수 있는 실마리를 잡았는데, 뒤에 $ sin^{2} \theta $는 어떻게 없앨 수 있을까요?

여기서 삼각함수 무적의 공식 $ sin^{2} \theta + cos^{2} \theta = 1 $이 등장합니다.

$ sin^{2} \theta = 1 - cos^{2} \theta $로 만들고, 원 식에 대입하면

$ cos 2\theta = cos^{2} \theta - (1 - cos^{2} \theta) $

$ cos 2\theta = 2cos^{2} \theta - 1 $

우리는 $ cos^{2} \theta $를 바꿔야 하니 $ cos^{2} \theta $로 정리해보죠

$ cos^{2} \theta = \frac{cos 2\theta +1}{2} $

그럼 바로 대입하면

$ \int cos^{2} \theta \, d\theta $

$ \int \frac{cos 2\theta +1}{2} \, d\theta $

$ \frac{1}{2} (\int (cos 2\theta + 1) \, d\theta) $

$ \frac{1}{2} (\frac{1}{2}sin 2\theta + \theta) $ +C 생략

자, 드디어 적분을 완료해서 적분기호가 사라졌습니다.

그러나 $ sin 2\theta $는 뭔가 보기에 깔끔하지 않죠.. 똑같이 삼각함수 배각공식을 이용하여 단일 $ \theta $항으로 만들어줍시다.

$ sin 2\theta = 2 \cdot sin \theta \cdot cos\theta $

물론 이 배각공식도 덧셈공식에서 왔습니다.

$ sin (\alpha+\beta) = sin \alpha \cdot cos \beta + cos \alpha \cdot sin \beta $

따라서 원 식에 배각공식을 이용하여 풀어주면

$ \frac{1}{2} (\frac{1}{2}(2 \cdot sin \theta \cdot cos\theta) + \theta) $

$ \frac{1}{2} (sin \theta \cdot cos\theta + \theta) $

$ \frac{1}{2}sin \theta \cdot cos\theta + \frac{1}{2}\theta $

여기서

$ x = sin \theta $

$ \theta = arcsin x $

$ y = \sqrt{1-x^{2}} = cos \theta $

이므로, $ \theta $에 대한 식이 아닌, 원 x에 대한 식으로 바꿔주면

$ \frac{1}{2} \cdot x \cdot \sqrt{1-x^{2}} + \frac{1}{2} \cdot arcsin \, x $

이 나오고, 이는 이 전 포스팅의 결과 식과 완전히 같은 모양이 됩니다.

원넓이의 부정적분 구하기원넓이의 부정적분 구하기 - 1) 일반공식으로 구해보기

원넓이를 처음 배우는 것은 초등학교 때, $\pi$를 3.14 근사값으로 배우면서 공식 암기와 함께 시작한다.

이후 중학교 과정에서 수의 확장과 함께 무리수로 $\pi$를 배워 무리수가 들어간 공식으로 배우고, 고등학교에 이르러서는 적분을 통해 원의 넓이를 새삼스레 다시 구해본다.

결국 우리의 수학 교과과정은 원에 대해서 배우는 것이다 라고 말해도 과언이 아닐 정도이다.

여기서 고등학교에서 정적분으로 원의 넓이를 구할 때 기계적으로 치환, 공식대입, 정적분을 통해서 '아 그냥 그렇게 되는구나'라고 알고 넘어가는 사람들이 대다수 일 터.

이번에 뭔가 궁금증이 생겨서 다시 풀어보니, 예전엔 그냥 단순히 치환하고 공식을 대입해서 풀었던 여기에는 참 많은 이유들이 있다는 것을 알게되었다.

그리하여 부정적분을 통해서 왜 이렇게 치환하고 그것으로 어떻게 넓이를 구하는지 알아보고자 한다.

물론 원은 사분원 넓이의 네배이니까 적분으로 원넓이 공식을 유도할 때 처럼, 반지름이 1인 단위원을 기준으로 하여 사분원의 넓이를 구하는 식으로 진행한다.

1] 일반식으로 원 넓이 구하기

부정적분으로 넘어가기 전에, 우리가 아는 일반 공식으로 $x$축에 대한 원의 넓이를 구할 수 있다.

딱 위의 그림에서와 같이 한번에 구하려면 왠지 적분을 써야 할 것 같지만,

위 그림과 같이 부채꼴과 삼각형으로 나눠서 구한다면 쉽게 원의 부분 넓이를 구할 수 있다.

부채꼴의 넓이는 호도법으로 전체 각도($2\pi$(360도))에 대한 원 넓이 $\pi r^2$을 전체 각도에 대한 부분각도의 비 만큼 곱해주면( $\frac{\theta}{2\pi}$ ) 부채꼴의 넓이($\frac{\theta}{2\pi}*\pi r^2 = \frac{1}{2} r^2\theta$)가 나온다.

물론 여기서는 r(반지름)을 1로 놨으니 r 변수는 사라질 것이다.

부채꼴의 각도를 기준으로 놨으니, 이제 삼각형도 계산할 수 있다. 삼각형의 넓이는 가로*세로/2이다.

xy좌표축의 x와 y의 값을 $\theta$로 표현하면 단위원의 매개변수 표현법에 의거 $y=sin\theta, x=cos\theta$로 표현할 수 있다. 다만 여기서는 부채꼴의 각도를 기준으로 표현을 했으니 우리가 쓰는 좌표축의 $y$값이 $cos\theta$가 될것이며 $x$값은 $sin\theta$가 될 것이다. 뭔가 두 길이가 달라진 것 같지만, 사실상 직사각형에서의 삼각형이니 두개는 대칭이다.(엄밀히 매개변수 표현법으로 $y=sin\theta, x=cos\theta$이니, 위의 예에서 $y=sin(90-\theta), x=cos(90-\theta)$가 되고 각각, $y=sin(90-\theta)=cos\theta$, $x=cos(90-\theta)=sin\theta$의 관계가 성립한다.)

수식으로 표현하면 $\frac{1}{2}sin\theta cos\theta$가 삼각형의 넓이가 될 것이다.

그리하여 두 식을 더한 $\frac{1}{2} \theta + \frac{1}{2}sin\theta cos\theta$ 값이 저 사다리꼴과 같은 도형의 부분 넓이가 된다는 것을 알았다.

여기까지 우리는 각도를 알면 그 각도에 해당하는 원의 사다리꼴과 같은 도형의 부분넓이를 구할 수 있게 되었으나, 반드시 각도를 알아야 한다는 단점이 있다.

여기서 x축 값 만으로 이 넓이를 구하려면, arcsin값만 알면 된다. 위의 그래프에서 x값과 같은 값을 나타내는 것은 $\theta$ 각도를 가지고 구한 $sin$값과 같다는 것을 알 수 있다. 그렇다면, 반대로 x값을 sin함수의 역함수인 arcsin에 넣어주면, 그 값에 해당하는 각도가 구해질 것이고, 그 각도값으로 부채꼴의 넓이 공식에 적용하면 부채꼴의 넓이를 알 수 있으므로 arcsin만 써주면 해결이다.

그렇다면 삼각형 부분은 어떻게 해결할 것인가?

사실 이 삼각형 부분은 x축 값이 밑변, 원의 방정식에서 x값을 대입한 값이 y값이다. 즉, y값은 $y = \sqrt{1-x^2}$이다.

이렇게 되면, 우리는 일반식으로 원에서의 사다리꼴과 같은 형태의 넓이를 x값에 따라 얻을 수 있는 일반식을 만들 수 있다.

$\frac{1}{2} arcsin\, x + \frac{1}{2}\cdot x \cdot \sqrt{1-x^2}$

다음 포스팅에서는 이를 부정적분으로 구해보는 시간을 가져볼 예정이다.

피보나치 수열의 일반항과 비율의 극한(황금비)

피보나치 수열하면 모르는 사람이 없을 정도로 아주 간단한 규칙을 가진 수열이다.

바로 앞의 두 숫자를 더하면 다음 숫자가 나오는 수열이다.

여기서 앞의 두 숫자는 1, 1 이다.

그러면 바로 아래와 같은 수열이 나오게 된다.

1 1 2 3 5 8 13 ...

물론 이 수열의 극한은 무한대로 발산할 것이 분명하지만, 이 수열의 두 항의 '비율'의 극한은 수렴할까? 수렴한다면 어디로 수렴할까? 한번 확인해보자.

여기서 수열의 극한을 확인하려면 항상 일반항이 있어야 한다. 그러나 피보나치 수열은 '앞의 두 수를 더하면 다음 숫자가 된다'는 점화식만 있는 형태이다. 그러면 이 점화식을 통해서 일단 피보나치 수열의 일반항을 구해보도록 하자.

피보나치 수열의 일반항 구하기

1. 피보나치 수열의 점화식을 써보자.

  피보나치 수열은 이 전의 두 항을 더하면 다음 항이 되는 수열이다.

  $ a_{n+2} = a_{n+1} + a_{n} $

  이러한 형태 점화식만 있는 상태로 등차, 등비, 멱급수 등등등 그 어떤 수열의 형태도 아니다.

2. 일반식으로 확장

  이 수열의 상태만으로는 우리가 뭔가 찝쩍거릴 건덕지가 없으니까, 일반적인 일반식으로 확장한 뒤 근과 계수와의 관계(Vieta's formulas, 두 근을 $ \alpha, \; \beta $로 놓으면 $ px^2+qx+r=0 $의 방정식에서 $ \alpha + \beta = - \frac{q}{p}, \alpha \beta = \frac{r}{p} $의 관계가 생긴다는 공식)를 활용하여 근을 활용한 일반식으로 변화시켜 볼 것이다. 참고로 수열에서 항수는 차수가 다른 방정식과 동일하게 볼 수 있다.(더 자세한 내용은 >>점화식에서의 특성방정식(characteristic equation)<<에서 확인할 수 있다.)

  $ a_{n+2} = a_{n+1} + a_{n} \Rightarrow x^2 = x + 1 \Leftrightarrow x^2 -x -1 = 0 $와 같이 쓴뒤, $ px^2+qx+r=0 $의 일반식으로 변환시켜주면, $ p = 1, q = -1, r = -1 $이 되고, 근과 계수와의 관계에서 $ \alpha+\beta=-\frac{q}{p}=1, \; \alpha \beta = \frac{r}{p} = -1 $이다.

  이는 다시 쓰면, $ p $가 기본적으로 1이기 때문에 $ \alpha+\beta = -q, \; \alpha \beta = r $이라고 놓을 수 있다.

  그래서 일반식을 다시 근과 계수와의 관계를 이용하여 계수가 아닌 근의 형태로 표현해주면

  $ x^2-(\alpha+\beta)x+\alpha \beta = 0 $

  이를 다시 수열의 항을 통해서 표현해주면

  $ a_{n+2} = (\alpha+\beta)a_{n+1} - \alpha \beta a_n $과 같은 근을 활용한 일반식으로 확장이 되었다.

  이때, $ a_1 = 1, \; a_2 = 1, \; \alpha + \beta = 1, \; \alpha \beta = -1 $이다.

3. 반복되는 형태를 만들어서 계산가능하게 만들자

  과거 >>https://omnil.tistory.com/172<<포스팅에서 감마함수를 팩토리얼로 변환하는 과정과 같이 등식의 좌변과 우변이 반복되는 형태를 만들어주게 되면 계산이 되지 않을 것 같은 등식도 계산이 된다. 특히 최종단계를 우리가 직접 계산해서 값을 알 수 있다면 더더욱이 말이다. 참고로 감마함수는 n=1일때 값이 1이며, 우리는 뭔가 이런단계를 거치면 1항이 1, 2항이 1이라는 것을 통해서 값을 구할 수 있을 것이다.

  $ a_{n+2} = (\alpha+\beta)a_{n+1} - \alpha \beta a_n $

  $ a_{n+2} = \alpha a_{n+1} + \beta a_{n+1} - \alpha \beta a_n $

  $ a_{n+2}-\alpha a_{n+1} = \beta a_{n+1} - \alpha \beta a_n $

  $ a_{n+2}-\alpha a_{n+1} = \beta (a_{n+1} - \alpha a_n) $

  이렇게 변환하면 등식의 좌변과 우변의 공동되는 부분의 한 항 차이가 $ \beta $배라는 것을 알 수 있다. 바로 이것으로 우리가 아는 $ a_2 $와 $ a_1 $를 가지고 계산할 수 있는 형태로 반복계산이 가능하다.

  $ a_{n+2}-\alpha a_{n+1} = \beta (a_{n+1} - \alpha a_n) $

  $ a_{n+1}-\alpha a_{n} = \beta (a_{n} - \alpha a_{n-1}) $

  $ \Rightarrow a_{n+2}-\alpha a_{n+1} = \beta^2 (a_{n} - \alpha a_{n-1}) $

  이런식으로 $ \beta $배씩 곱해주면 우항을 a2와 a1항으로 계산할 수 있는 형태로 만들어줄 수 있다.

  이 때, $ \beta $가 몇개 생기는지는 항 수를 보고 생각하면 된다.

  우변의 맨 오른쪽항이 a2항에서 a1항으로 떨어지게 되면, $ \beta $는 한개가 생길 것이다. 즉, an항에서 a1항으로 떨어지면 (n-1)개의 $ \beta $가 생성될 것이다.

  $ a_{n+2}-\alpha a_{n+1} = \beta \cdot \beta^{n-1} \cdot (a_{2} - \alpha a_1) $

  $ a_{n+2}-\alpha a_{n+1} = \beta \cdot \beta^{n-1} \cdot (1 - \alpha \cdot 1) \leftarrow \because a_2=1,\; a_1=1 $

  $ a_{n+2}-\alpha a_{n+1} = \beta \cdot \beta^{n-1} \cdot \beta \leftarrow \because \alpha + \beta = 1 $

  $ a_{n+2}-\alpha a_{n+1} = \beta^{n+1} $

  즉  $ a_{n+2}-\alpha a_{n+1} $는 $ \beta $를 $ n+1 $번 곱한 것이니 항수 만큼 $ \beta $를 곱해주는 횟수가 된다는 것을 알 수 있다. 그렇다면 우리가 알고싶은 $ a_n $을 기준으로 하는 식으로 바꿔주면

  $ a_{n}-\alpha a_{n-1} = \beta^{n-1} \cdots (1)$

  이 되고, 이는 $ \alpha $ 변수와 $ \beta $ 변수를 바꾸어도 변수위치만 바뀐 동일한 식이 나온다.

  $ a_{n}-\beta a_{n-1} = \alpha^{n-1} \cdots (2)$

4. 연립하여 $ a_n $에 대한 일반항으로 풀어준다.

  변수 두개에 식이 두개가 나왔으니 연립방정식으로 풀 수 있다.

  (2)식에 $ \frac{\alpha}{\beta} $배를 해준 뒤 (1)-(2)식을 해줘서 $ a_{n-1} $항을 소거하여 $ a_n $의 일반항을 얻을 수 있다.

  $ a_{n}-\alpha a_{n-1} = \beta^{n-1} \cdots (1)$

  $ \frac{\alpha}{\beta}a_{n}-\alpha a_{n-1} = \frac{\alpha^{n}}{\beta} \cdots (2)$

  $ (1)-(2) $

  $ a_{n}-\frac{\alpha}{\beta}a_n = \beta^{n-1}-\frac{\alpha^n}{\beta} $

  $ \beta a_{n}-\alpha a_n = \beta^{n}-\alpha^n $

  $ (\beta -\alpha) a_n = \beta^{n}-\alpha^n $

  $ \therefore a_n = \frac{\beta^{n}-\alpha^n}{\beta -\alpha} $

  일반항 겟!!

  이제 일반항에 값만 대입해주면 진짜 n에 몇번째 항인지만 대입해주면 거기에 해당하는 값이 나오는 일반항이 된다.

5. $ \alpha $와 $ \beta $의 값 구하여 일반항에 대입하기

  여기서 $ \alpha $와 $ \beta $는 사실 $ x^2 -x -1 = 0 $의 두 근과 같기 때문에 근의 공식을 통하여 바로 값을 구할 수 있다.

  $ ax^2+bx+c = 0 $에서 두 근은 $ \frac{-b\pm \sqrt{b^2-4ac}}{2a} $식으로 구할 수 있다.
  $ \frac{1\pm \sqrt{5}}{2}, \; a=1, \: b=-1, \: c=-1 $

  $ \beta = \frac{1 + \sqrt{5}}{2}, \; \alpha = \frac{1 - \sqrt{5}}{2} $

  $ a_n = \frac{\left(\frac{1+\sqrt{5}}{2}\right)^n-\left(\frac{1-\sqrt{5}}{2}\right)^n}{\frac{1+\sqrt{5}}{2}-\frac{1-\sqrt{5}}{2}} $

  $ \therefore a_n = \frac{1}{\sqrt{5}}\left(\left(\frac{1+\sqrt{5}}{2}\right)^n-\left(\frac{1-\sqrt{5}}{2}\right)^n\right) $

  이렇게 피보나치 수열의 일반항을 구했다!!

  근데, 유리수의 합으로 나타나는 피보나치 수열에서 일반항에 무리수가 들어가는 것이 신기하지 않은가!

피보나치 수열의 비율의 극한

이렇게 일반항을 구했으면 비율의 극한도 쉽게 구할 수 있다.

여기서는 더 큰수를 더 작은수로, 즉 $ \frac{a_{n+1}}{a_n} $의 비를 구할 것이다.

이번엔 비율을 구할 것이기 때문에, 숫자까지 들어간 일반항 보다는 문자로 표현된 더 한눈에 보기 간편한 일반항을 사용하여 극한을 구해볼 것이다.

1. 비율 식 구하기

  $ a_{n+1} = \frac{\beta^{n+1}-\alpha^{n+1}}{\beta -\alpha} $

  $ a_{n} = \frac{\beta^{n}-\alpha^{n}}{\beta -\alpha} $

  $ \frac{a_{n+1}}{a_n} = \frac{\frac{\beta^{n+1}-\alpha^{n+1}}{\beta -\alpha}}{\frac{\beta^{n}-\alpha^{n}}{\beta -\alpha}} $

  $ \frac{a_{n+1}}{a_n} = \frac{\beta^{n+1}-\alpha^{n+1}}{\beta^{n}-\alpha^{n}} $

  $ \frac{a_{n+1}}{a_n} = \frac{\beta-\alpha \left(\frac{\alpha}{\beta}\right)^n}{1-\left(\frac{\alpha}{\beta}\right)^{n}} $

2. 극한 씌워주기

  $ \lim_{n \to \infty} \frac{a_{n+1}}{a_n} = \lim_{n \to \infty} \frac{\beta-\alpha \left(\frac{\alpha}{\beta}\right)^n}{1-\left(\frac{\alpha}{\beta}\right)^{n}} $

  여기서 $ \beta = \frac{1 + \sqrt{5}}{2}, \; \alpha = \frac{1 - \sqrt{5}}{2} $이고, $ \beta $가 $ \alpha $보다 크기 때문에 $ \left(\frac{\alpha}{\beta}\right)^n $항은 $ n $이 무한대로 갈 때 값이 0으로 수렴한다.

  참고로 실제 값을 대입해서 계산해본 $ \left(\frac{\alpha}{\beta}\right) $ 값은 $ \frac{\sqrt{5}-3}{2} $이며, 그 값은 약 -0.382이다. 즉, 이 값을 무한대로 제곱할 경우 양과 음을 반복 진동하며 수렴한다.

  즉, 극한을 취한 뒤의 값은

  $ \lim_{n \to \infty} \frac{a_{n+1}}{a_n} = \frac{\beta-\alpha 0}{1-0} $

  $ \lim_{n \to \infty} \frac{a_{n+1}}{a_n} = \beta $

  $ \therefore \beta $

  이며, 이 $ \beta $값은  $ \frac{1 + \sqrt{5}}{2} $이므로, 피보나치 수열의 비율의 극한 값은 $ \frac{1 + \sqrt{5}}{2} $이 된다.

  그러면 이 값은 과연 무엇일까

황금비

  인생을 살면서 '황금비'라는 단어를 한번은 들어본다.

  황금비는 1: 1.618로써 근사하면 5:8정도의 비율을 나타내는 것을 황금비라고 한다.

  이것은 우리가 어떤 비율을 봤을 때 가장 아름답다고 생각하는 비율이라고 하는데, 이 1.618이라는 값은

  $ \frac{1 + \sqrt{5}}{2} $을 계산하면 나오는 값이다.

  즉, 피보나치 수열의 비율을 극한으로 가져가면 황금비를 가진다는 사실!