김 용묵의 절대공간 - 블로그

이제 우리는 오일러 방정식 말고 다른 방식으로도 사이클로이드의 최단 강하 성질을 유도해 보고자 한다.
이것은 뉴턴 말고 동시대의 다른 천재 과학자이던 요한 베르누이가 최초로 발견했다.

물리 과목이 고전역학과 유체역학, 열역학, 전자기학을 넘어 파동과 입자 파트로 넘어갈 즈음에, 우리는 스넬의 법칙이란 걸 접하게 된다.
빛이 한 매질 속에서 v1이라는 속도로 나아가다가 속도가 다른(v2) 매질로 진입했을 때 잘 알다시피 굴절이 발생하는데.. 그 굴절 각도 사이에

sin θ1 : v1 = sin θ2 : v2

비율이 동일하게 유지된다는 법칙이다.
물리학이라는 게 딱딱한 물체의 가시적인 움직임만 기술하는 줄 알았는데.. 물 속에 비치는 물체가 실제보다 얕게 보이는 이유와 정확하게 얼마나 얕게 보이는지까지 수식으로 알려준다는 것이 경이롭다. 이건 게임에다 비유하자면 말 그대로 오브젝트들의 운동 역학뿐만 아니라 그래픽 렌더링에까지 관여한다고 볼 수 있다.

뭐, 0도~90도 사이에서 sin x는 x 자체와 별 차이 없을 정도로 단조증가만 하니, 위의 법칙이 아주 특이한 결과를 만들지는 않는다. 빛의 속도가 빠른 곳에서는 각도가 작아지고, 느린 곳에서는 각도가 대체로 커진다. 단지 고각이 될수록 그 정도에는 선형적이지 않은 차이가 발생할 뿐이다.
(낮은 각도에서 x와 sin(x)는 그리 크게 차이가 나지 않는다는 점은 단진자 주기의 근사값을 구할 때도 쓰인다. 이에 대해서는 곧 다시 다루어질 것이므로 참고하시라.)

그런데 빛은 애초에 질량이 없는 물건이고 구슬처럼 데구르르 구르지도 않을 텐데 저 법칙이 사이클로이드 내지 최단 강하 곡선과 무슨 관계가 있는 걸까..??
고지대는 매질이 진해서(...;) 빛의 속도가 제일 느리고, 아래로 갈수록 점진적으로, 선형적으로 매질의 농도가 옅어져서 빛이 더 빠르게 진행할 수 있다고 가정해 보자. 이게 빛에게 중력 역할을 하는 거나 마찬가지이다.

다시 말해, 공기 아니면 물 같은 이산적인 흑백이 아니라, 그러데이션을 생각하면 된다.
이런 상황을 만들어 놓고 위에서 아래로 빛을 조금이라도 비스듬하게 쏘면.. 빛 역시 단순무식하게 직선 형태로 나아가지 않게 된다. 최고 속도를 낼 수 있는 방향으로 쭈루룩 휘어지고, 그 궤적이 놀랍게도 사이클로이드가 된다.

즉, 궤적의 어느 점에서나, 어느 높이에서나 sinθ / v의 값은 일정해야 한다.
출처에 따라서 sinθ / sqrt(h) = C(상수) 라고 표기하기도 한다. 앞서 값을 구한 바와 같이, 속도 v는 sqrt(2*g*y)로, 높이의 제곱근에 비례하기 때문이다. y축이나 높이 h나 동일한 개념인 것은 다들 아실 테고..

그리고 여기서 각도 θ라는 건 보다시피 y축(세로선) 기준의 값이다. 이를 x축 기준으로 환산하면 90도-θ나 마찬가지이다. 즉, y축 기준 입사각 90도는 x축 기준 0도와 같다는 뜻.. 그러므로 y축 기준으로 sin θ는 우리 입장에서는 cos θ와 같은데..;;

여기서 놀라운 발상의 전환이 가능하다.
cosθ는 1/sqrt(tan(θ)^2 + 1) 로 형태를 바꿔치기 할 수 있다! (sinθ^2 + cosθ^2 = 1이니까)
여기서 탄젠트라는 건? 고맙게도 x축 대비 y의 변화량.. dy/dx, 다시 말해 y' 그 자체이다.

그러므로 sinθ/v에서 분자는 1/sqrt(1+y'^2)로 깔끔하게 나왔고, 분모 v에다가는 sqrt(2gy)를 집어넣으면..
이게 곧 미분 방정식

1/( sqrt(2gy)*sqrt(1+ y'^2)) = C

이 나오는데, 이건 앞서 오일러-라그랑주 방정식에서 d/dy 안에다 대입해서 집어넣었던 ∂F/∂x'와 같은 형태이다. 그때는 x를 먼저 구했지만 지금은 y를 먼저 구할 뿐... 그리고 그때도 지금처럼 y를 먼저 구하는 식으로 풀 수 있다.

잡다한 상수와 근호를 치우고 식을 정리하면 y(1+y'^2) = C 를 얻는다.
제곱이 좀 압박스럽게 보이긴 하지만, 얘는 dx = dy sqrt( y/(C-y) ) 로 변형 가능하며 이건.. 역시 이전 방식에서 얻었던 x'의 형태와 같다.
그러니 y를 삼각함수의 제곱으로 치환해서 계속 풀어 나가면 된다. 이상.. 더 이상의 자세한 설명은 다른 참고 사이트에서 열람하시기 바란다. (☞ 참고 1, 참고 2)

사이클로이드는 공대 1학년의 기초 필수 코스인 벡터 미적분학 교육과정에서 정말 최고로 적합한 교보재인 것 같다.
스넬의 법칙 덕분에 풀이의 일부가 좀 간소화된 듯한 느낌인데.. 사실 스넬의 법칙도 직접 유도해 보면 물리에 앞서 오일러-라그랑주 방정식과 동급의 수학 원리가 이미 깔끔하게 담긴 것이나 마찬가지이다.

빛이 속도를 최대화하기 위해 사이클로이드 궤적 형태로 굴절된다니, 이건 우주 로켓 발사와도 동일하지는 않지만 비슷한 구석이 있어 보인다.
공기가 가장 짙고 공기 저항도 극심한 지상에서는 이 구간을 빨리 빠져나가는 게 중요하므로 로켓이 닥치고 수직 상승한다. 하지만 공기가 충분히 옅어지면 이제 수평으로도 속도를 내서 지구 궤도에 진입할 채비를 하니까 말이다.

그렇게 각 고도별로 수직과 수평 동력을 잘 배분하여 연료를 가장 적게 쓰고 궤도에 진입할 수 있다. 그리고 그 경로를 짜는 건 굉장한 노하우가 필요하다. 여기에도 사이클로이드 궤적이 쓰이는지 궁금하다.

2. 등시 곡선

아이고 힘들다.. 최단 강하 곡선 얘기가 이렇게 길어질 줄 몰랐네.. 그런데 이게 사이클로이드의 역학적 특성의 끝이 아니다.
얘는 최단 강하 곡선인 동시에 놀랍게도 ‘등시 곡선’이기도 하다.
무슨 말이냐 하면.. 공을 곡선의 맨 꼭대기에서 굴리든, 밑바닥 근처에서 살살 굴리든, 목적지인 밑바닥에는 동일한 시간 후에 도착한다는 것이다! (공기 저항과 마찰을 무시하는 한)

이 말을 들으면 같이 떠오르는 게 있다. 바로 단진자의 운동.
학창 시절에 배운 기억으로는 단진자도 각도라든가 추의 무게와 무관하게 왕복 주기가 동일하다. 오로지 실의 길이 l과 중력가속도 g의 영향만 받아서 2π*sqrt(l/g)라는 공식이 나왔었다. (줄 길이의 제곱근에 비례)

그런데 공중에 대롱대롱 매달린 단진자가 그리는 궤적은 사이클로이드가 아니라 그냥 원호이다. 둘이 어떻게 동일할 수가 있지?
관련 자료를 다시 찾아보니, 단진자의 주기는 각도와 완전히 무관한 게 아니었다. 중등 수준에 맞게 난이도 보정을 하느라 θ가 충분히 작아서 sinθ가 θ와 얼추 같다고 “치고” 식을 유도한 것이었다. 원래 미분방정식은 θ'' + sinθ *g/l = 0 이렇게 나오는데, sinθ를 놔두면 미분방정식이 너무 어렵고 복잡해지기 때문이다.

언젠가 물리 시간에 뭔가 “sinθ가 θ와 얼추 같다고 치고” 식을 유도하는 걸 접한 기억은 있지만 그게 단진자 주기 공식이었다는 건 오랫동안 잊고 있었다.
sinθ를 살려서 문제를 해석학적으로 정확하게 풀려면 타원 적분이란 걸 동원해야 하며, 최종적인 답이 무한급수의 형태로 튀어나오게 된다. 이런 복잡한 항들이 2π*sqrt(l/g)의 뒤에 덕지덕지 추가적으로 곱해진다. 자세한 것은 타 사이트의 설명을 참고하시기 바란다. 단진자가 이 정도로 까다로운 물건이었다니..;;

그럼 60도짜리 진동과 15도짜리 진동이 일치할 정도가 되려면.. 진자도 원호 궤적이 아니라 사이클로이드 궤적을 따라 움직여야 한다. 사이클로이드는 원 궤적이 아니라는 것을 이미 앞서 설명한 바 있다. 진자의 궤적이 사이클로이드가 되게 바꾸는 방법은 다음과 같이 따로 있다..!

바로 이렇게.. 반지름 a/2짜리 원을 굴려서 얻어지는 사이클로이드 모양의 벽을 반반 쪼개서 양 옆에 배치하면 된다. 그러면 아무 지점에서 진동시켜도 사이클로이드 궤적을 그리면서 완벽한 등시 등주기가 실현된다. 얘는 각도와 무관하게 진동 주기에 π*sqrt(l/g) 다음으로 원운동 진자 같은 복잡한 항이 덕지덕지 붙지 않는다. 이것을 첫 발견자의 이름을 따서 호이겐스 진자라고 부른다.

사이클로이드 궤적은 어째서 등시 강하가 보장되며, 저렇게 진자를 진동을 보정하면 어떤 원리로 사이클로이드가 되는 걸까?

우린 임의의 궤적 함수가 주어졌을 때 구슬이 끝까지 다 굴러가는 데 걸리는 시간을 구하는 함수 T를 이미 구한 바 있다. 분자는 sqrt(1+f'(x)^2) dx이고, 분모는 sqrt(2*g*f(x))이던 그 식 말이다.
그걸 써먹으면 된다. 얘를 임의의 시작점을 집어넣어서 끝까지 적분을 하더라도 항상 동일한 값이 나온다는 것을 입증하면 등시 강하 특성을 입증할 수 있다.

그런데 사이클로이드는 매개변수 형태로 표현돼 있으므로..
x = r*(θ-sinθ), y = r*(1-cosθ) 로부터
dx/dθ = r*(1-cosθ), dy/dθ = r*sinθ 를 얻는다.
그러므로 dy/dx는 양변을 r로 나눈 sinθ / (1-cosθ)가 된다.
이 점을 염두에 두면 식의 분자 부분을 다음과 같이 깔끔하게 dx 대신 dθ로 치환할 수 있다.

다음으로 분모는? f(x), 아니 y에 해당하는 r*(1-cosθ)만 대입해 주면 되는데...
우리의 문제를 풀기 위해서는 약간의 변화를 줘야 된다.
일반적인 상황이라면 최고점인 y=0일 때 속도가 0이었지만, 우리는 y>0이고 더 낮은 임의의 지점에서 가속을 시작하여 구슬을 굴리더라도 끝까지 다 구르는 데 걸리는 시간이 동일하다는 걸 입증해야 한다.

그 임의의 지점을 단순 높이로 나타내건 각도로 나타내건 그건 중요하지 않다. 하지만 우리는 원을 0도에서 180도(pi!) 반원 각도만치 굴려서 사이클로이드 궤적을 굴리는 상황을 설정할 것이므로 각도 표기를 계속 사용할 것이다.
θ0이라는 각도에서 처음으로 공을 굴렸다면 속도가 어떻게 될까?
v = sqrt(2g*(y-y0)) 로부터 sqrt( 2g*r*( cos(θ0) - cos(θ) )) 이 나온다. 각도를 빼는 게 아니라 y축 관점에서 cos 결과값을 빼 주면 된다.

무슨 말이냐 하면, 사이클로이드의 y축 매개변수 r*(1-cosθ)에서 1이 cos(θ0)로 바뀐다는 것이다.

자, 그냥 제일 높은 곳에서 처음부터 구른 것이면 θ0는 0도이기 때문에 분자와 분모에서 2a*(1-cosθ)가 통째로 약분되어 없어진다. 적분식은 θ와 무관한 상수가 되며(= 극값), 사이클로이드가 고정된 최단 강하 곡선임이 증명된다.
다음으로 분모의 1을 cos(θ0)이라는 상수로 바꾸고, 얘를 θ0부터 pi까지 정적분한 값도 θ0와 무관하게 고정된 값이 나온다는 것을 입증하면.. 이건 구르는 소요 시간이 θ0의 값과 무관한 등시 곡선임을 추가로 증명할 수 있다.

얘도 다루기가 많이 까다로워 보이지만.. 그래도 끝자락이 삼각함수에서 비교적 취급하기 쉬운 상수인 pi이며, 삼각함수들을 반각/배각 공식을 이용해 제곱으로 치환하면 제곱근 근호를 걷어내고 적분을 그럭저럭 격파할 수 있다. 이렇게 말이다.

아.. 경이롭다. =_=;; 특히 마지막에 arcsin으로 치환되는 적분 부분이 정말 압권이다.
이를 통해 θ0은 싹 사라져 버리고, θ0가 있더라도 그냥 θ0=0일 때와 동일하게 상수 시간이 도출된다는 것을 알 수 있다.

이상이다. 사이클로이드를 제대로 마스터 하려면 호이겐스 진자의 진동 궤적이 왜 사이클로이드가 되는지도 해석학적으로 유도해야겠지만.. 그건 이 자리에서는 생략하겠다.
아무쪼록 (1) 사이클로이드가 원이나 타원의 호가 아닌 이유, (2) 저게 최단 강하 곡선인 이유, (3) 어느 위치에서건 등시 강하 곡선인 이유를 제대로 설명할 수 있다면 당신은 대학 나온 훌륭한 이과생 공돌이라고 불릴 자격이 될 것이다.

1600년대 중후반은 갈릴레오와 케플러에 이어 뉴턴, 라이프니츠, 베르누이 같은 사람이 등장하면서 미적분학이 태동하고 고전역학이 전성기를 맞이한 시기였다. 그 옛날 사람들이 사이클로이드라든가 단진자 운동의 특성을 갖고 저런 고민을 하고 있었던 것이다.

특히 저 때는 괘종시계가 현역이던 옛날이었기 때문에 단진자 운동을 정확하게 기술하고 등주기를 실현하는 것이 매우 중요했다. 그게 시계의 정확도와 직결됐기 때문이다.
그러다 1700년대에 가서는 오일러와 라그랑주가 나오면서 이 바닥은 더 세련된 계산법이 개발되고 발전을 이루게 되었다.;;

Posted by 사무엘

이전 글들에서는 현수선 얘기를 하면서 한창 쌍곡선함수 얘기를 늘어놨는데, 이 글에서는 분위기를 바꿔 정통(?) 삼각함수 얘기를 좀 하겠다.
수학에서 다루는 많고 많은 곡선 중에는 “원을 직선 위에서 굴릴 때 그 원에 놓인 정점이 그리는 궤적”이란 게 있다. 이걸 ‘사이클로이드’라고 한다.

반지름 r인 원이 (0, r) 위치에 놓여 있다. 그리고 궤적을 추적하고자 하는 점은 처음에 바닥의 원점 (0, 0)에 있다고 치자.
수학 좌표계는 x의 양의 방향이 오른쪽이고 y의 양의 방향은 위쪽이다. 이 상태로 원이 x의 양의 방향으로 굴러가려면 시계 방향으로 굴러가야 하며, 각도는 처음에 270도에서 시작했다가 줄어들어야 한다. 각도가 늘어나는 방향은 반시계 방향이기 때문이다.

즉, (cos(3π/2 - θ), sin(3π/2 - θ)+1) 정도 될 텐데..
cos(3π/2 - θ)는 sin(-θ) → -sin(θ)로 간소화시킬 수 있다.
sin(3π/2 - θ)도 -cos(θ)로 간소화시킬 수 있다.

이건 원이 제자리에서만 빙글빙글 도는 상황이다.
실제로 굴러간다는 조건을 추가하면 원은 아주 단순하게 x축의 양의 방향으로 r*θ만치만 이동하게 된다. 원의 전체 둘레는 2πr이고 θ는 각도 겸 원호의 둘레와 정확하게 대응하기 때문이다.

그러므로 사이클로이드의 궤적은 ( r*(θ-sinθ), r*(1-cosθ))라고 깔끔하게 나온다. 원이 헛돌지 않고 제대로 굴러가는 이상, 점이 아예 뒤로 후퇴하는 일은 어떤 경우에도 절대 없다는 걸 알 수 있다.

사이클로이드는 원이나 타원의 호와 비슷하게 생겼지만, 실제로는 이와 무관한 별개의 궤적이다. 아래의 그림 비교와 관련 설명을 참고하시라.

그리고 얘는 x와 y의 위치를 매개변수를 통해 제각기 나타낼 때는 식이 저렇게 아주 깔끔하게 구해지는 반면, y=f(x) 꼴의 단순한 양함수로 나타내는 건 굉장히 복잡하고 까다롭다. x축이 그냥 sin(t)가 아니라 t-sin(t)인데.. 이게 역함수를 구하기가 몹시 난감하기 때문이다.

그나마 x=g(y)라는 일종의 음함수에 가까운 형태로 나타내는 게 y=f(x)보다는 더 할 만하다. 원이 굴러간 거리에 대한 함수가 아니라 원의 범위 영역에 대한 함수 말이다.
사이클로이드는 이렇게 까다로운 면모가 있음에도 불구하고 한 사이클 궤적의 길이라든가 그 궤적 아래의 면적은 적분을 통해 깔끔하게 떨어지는 형태로 구할 수 있다. 이쪽으로 더 관심 있으신 분은 위키백과 등의 타 사이트를 참고하시라.

이 글에서 현수선 다음으로 사이클로이드를 소개한 이유는 얘도 현수선과 동급으로 역학적으로 아주 흥미로운 특성이 있는 물건이기 때문이다.

1. 최단 강하

먼저, 얘는 일명 ‘최단(최고속) 강하 곡선’이다. 높은 지점에서 낮은 지점으로 공을 굴리는데 무작정 최단 거리인 직선 경사로가 아니라, 사이클로이드를 상하 반전시킨 형태의 경사로를 만드는 게 좋다. 그러면 공이 중력 버프를 받아서 목적지에 더 빨리 도달하게 된다! (지면 마찰과 공기 저항 따위는 모두 무시)

그렇다면 비행기의 비상 탈출용 미끄럼틀도 직선이 아니라 이런 형태로 만드는 게 더 나을 것 같다. 물론 현실적으로는 무리이겠지만, 그래도 어린이용 미끄럼틀 중에는 이런 과학 고증(?)을 반영하여 사이클로이드 곡선을 흉내 낸 모양인 것도 있다.

그런데 의문이 든다. 이것도 까놓고 말해 그냥 포물선이라든가, 사분원이나 2차 베지어 곡선, 혹은 아예 역학적인 안정성이 검증돼 있는 현수선이 최단 강하일 수도 있을 텐데 왜 하필 사이클로이드가 당첨인 걸까? 직선보다 아래로 볼록한 모양이어야 하겠다는 건 수긍이 가지만 왜 하필 저 모양이 최적인지는 직관적으로 납득이 잘 안 된다.

먼저 우리는 구슬이 처음 놓여 있는 좌측 최상단 꼭대기가 원점 (0, 0)이라고 가정하자. 그리고 계산의 편의를 위해, y축은 중력이 향하는 방향(=아래)이 +로 증가하는 방향이라고 간주하도록 하겠다. 즉, 통상적인 수학 좌표가 아니라 컴퓨터 화면의 좌표계를 사용한다는 것이다.
그러면 구슬은 아래로 굴러가면서 x축과 y축 모두 값이 일관되게 증가한다고 생각할 수 있으며, 저 사이클로이드 공식도 y축의 부호를 뒤집을 필요 없이 그대로 사용할 수 있다.

이 문제를 푸는데 매우 중요하게 활용되는 단서는 (1) "역학적 에너지 보존의 법칙"이다. 구슬이 직선 경사로를 쫙 구르건 그 어떤 꼬불꼬불한 곡선을 타고 오르내렸건, y축 중력 가속도가 g인 상태에서 y라는 높이만치 내려가 있다면 그 당시에 물체의 속도는 mgy = mv^2 /2를 근거로 v=sqrt(2*g*y)가 된다. 위치 에너지가 그만치 운동 에너지로 바뀌었음을 뜻한다. 당연한 말이지만 x축과는 전혀 무관하다.

이 속도는 xyz 같은 축을 총체적으로 고려한 단위 시간당 이동량이다. 그 속도를 얻은 상태에서 수평 이동을 한다면 구슬은 공기 저항과 마찰이 없는 한, 영원히 등속 운동을 할 것이다. 상승하면 속도가 점차 느려질 것이고, 하강하면 속도가 붙어서 빨라질 것이다.

여기까지 준비가 완료됐으면 사이클로이드 문제는 크게 두 방법으로 풀 수 있다.
하나는 뉴턴이 터를 닦으시고 오일러-라그랑주가 끝장을 낸 (2) 변분법을 동원하는 것이다. 하드코어한 고전역학의 범주에서 끝을 보는 것이다.
다른 하나는 파동 분야에서 정립된 (3) 스넬의 법칙을 접목하는 것이다. 겉보기로 수식의 압박이 변분법보다는 '약간' 덜할지 모르지만, 그럼 스넬의 법칙은 왜 어째서 성립하는지를 일일이 따지고 든다면 난이도는 비슷하게 안드로메다로 치솟는다.

그럼 변분법 버전부터 먼저 살펴보도록 하겠다.
통상적인 미분이 함수의 극점을 구해서(= 도함수가 0) 특정 구간에서 함수의 최소/최대값을 구하는 데 쓰인다면, 변분법이란 범함수의 최소/최대값을 구하는 방법론을 말한다.
그럼 범함수란 무엇이냐? 프로그래밍에는 평범한 숫자나 객체가 아니라 함수 자체나 람다를 다른 함수의 인자로 넘겨주는 게 있다. 그것처럼 수학에도 함수를 받아들여서 스칼라 형태의 값을 되돌리는 함수가 있는데, 그걸 범함수(functional)라고 한다. 특정 함수가 특정 구간에서 미분 불가능한 지점의 개수.. 그런 것도 범함수의 일종이 될 수 있다.

변분법을 이용하면 두 점을 잇는 가장 짧은 경로가 직선인 이유처럼 공리 수준의 너무 당연해 보이는 것부터 시작해서..
둘레가 동일한 도형 중에 면적이 최대인 놈이 원의 형태가 되는 이유, 물방울· 비누 방울이나 잠수정이 전부 구형으로 만들어지는 이유(단위 면적당 압력 최소화) 같은 것도 다 수학적으로 엄밀하게 유도할 수 있다.

최단 강하 곡선 문제에서 우리가 구해야 하는 것은 길이나 면적 따위가 아니라 "도달 시간"의 최소값이다.
구슬이 위에서 아래로 굴러가는 궤적을 y = f(x)라는 함수로 나타낸다면, a라는 임의의 x축 지점에서 구슬의 진행 속도는 f(a)의 값으로부터 구할 수 있다. 속도는 앞서 공식을 구한 바와 같이 전적으로 y축에 의해서만 결정되고, y축의 값은 x축 대한 함수 f로부터 구할 수 있으니 말이다.

그렇다면 (2-1) 이 함수를 어째 잘 적분해서 f(x)라는 함수에 대해 구슬이 다 내려갈 때까지 걸리는 시간도 구할 수 있을 것이다.

이건 흔히 떠올릴 수 있는 상황은 아니다. 속도 함수를 적분해서 구하는 건 보통은 거리이기 때문이다. 시간에 대해 적분하면 끝이다. (예: 시속 100km로 2시간을 달린 거리는 200km)
그런데 반대로 소요 시간을 구하기 위해서는 반대로 함수 궤적의 거리가 주어져 있어야 한다. (예: 시속 100km로 300km를 달리는 데 걸리는 시간은 3시간) 그래도 이것도 적분이긴 하다.

앞서 sqrt(2*g*y)이라고 값을 구한 속도 v는 거리/시간, 즉 ds/dt라는 개념이다. 그런데 이 ds라는 건 역시 x 변화량 대비 y의 변화량을 거리화한 것이며, y 변화량은 곧 f(x)의 변화량과 같다. 그러므로 이것은 sqrt(1 + f'(x)^2 )로 나타낼 수 있다.
v = ds/dt 에서 sqrt(2*g*f(x)) = sqrt(1 + f'(x)^2 ) / dt 가 되고..
이 관계로부터 dt = sqrt(1 + f'(x)^2 ) / sqrt(2*g*f(x))이 된다.

오.. 그러므로 구슬이 구르는 궤적 함수가 f이고 중력 가속도는 g, 구슬이 다 구른 오른쪽 끝의 x축 지점이 a라 할 때, 구슬이 다 구르는 데 걸리는 시간 T는 다음과 같이 깔끔하게 구해진다.

여기까지 도달했으면 문제가 반 정도는 해결됐다고 볼 수 있다. 우리는 a와 g가 고정돼 있을 때 T를 최소화하는 f(x)를 구하면 된다.

1/sqrt(2*g)는 그냥 상수이므로 적분 기호 밖으로 옮겨도 아무 상관 없다. 그리고 범함수에서는 그 정의상 x라는 변수뿐만 아니라 그 x에 대한 함수 f(x)까지 범함수의 parameter가 될 수 있다. 그렇기 때문에 f(x) 대신 그냥 y로 표기하며, 도함수도 f'(x) 대신 dy/dx라고 표기하는 걸 더 일상적으로 보게 된다. 애초에 사이클로이드도 f(x) 대신 x, y축 따로 매개변수 형태로 표기하는 게 더 유리하기도 하고 말이다.

이 문제를 풀기 위해 쓰이는 도구는 바로 (2-2) 오일러-라그랑주 방정식이다.
x, y라는 변수가 있고 F가 y, y', x에 대한 범함수라고 하자. 우리가 푸는 문제에서는 x, y는 구슬의 궤적을 나타내며, F는 그 궤적으로 끝까지 구르는 데 걸리는 소요 시간을 구하는 구간 적분이 된다.
오일러-라그랑주 방정식에 대해 제일 간단하고 엉성하게 요점만 말하자면, F가 극값(최소 또는 최대)이 나올 때 이들은 다음 등식을 만족한다.

(왜 그런지 궁금하면 이변수함수의 연쇄법칙을 동원해서 저 식을 직접 유도하면 된다. 하지만 이 자리에서는 시간과 지면의 부족으로 증명을 생략하겠다. ㄲㄲ)

저 식에서 F에다가는.. 우리가 구한 처음 식 T에서 상수배 항과 dy 적분을 걷어낸 sqrt((1+x'^2)/y), 즉 순간 변화량을 대입한 뒤, 식을 x'에 대해 풀면 된다.
중요한 것은 우리의 경우, ∂F/∂x가 0이 보장된다는 것이다. 앞서 살펴보았듯이, 물체의 속도는 x축 위치와는 전혀 무관하고 y축 높이에 의해서만 결정되기 때문이다.

그러니 저 식에서 d/dy (∂F/∂x')도 같이 0이어야만 한다. F를 x'에 대해 편미분한 뒤, 양변에 sqrt(y)를 곱한 결과는 아래와 같이 된다. (x가 어떤 함수인지 정확한 정체를 모르기 때문에 x와 x'는 서로 독립변수로 간주됨) 이 식을 x'에 대해서 정리하면..

헥헥.. 이제 x'이 구해졌다. 그러면 함수 x는 y에 대해 적분을 해서 구할 수 있긴 하지만.. sqrt(x/(1-x))라는 함수는 부정적분을 호락호락 쉽게 구할 수 있는 놈이 아니다.
쟤를 수월하게 적분하는 방법은 x, y를 다른 변수 형태로 치환하는 것이다. 애초에 사이클로이드는 x, y 궤적을 매개변수 형태로 치환해야 표현하거나 취급하기 용이한 물건이다.

y를 sin(t)^2 / C^2 으로 치환하면 거추장스러운 C가 없어진다. 그리고 근호 안의 식이 sin(t)^2 / cos(t)^2로 바뀌고 근호도 없어져서 식 전체가 tan(t)로 바뀌는 마법이 펼쳐진다.
단, dy도 sin(t)^2의 도함수인 sin(t)*cos(t) dt로 바뀌기 때문에 저 tan(t)에다가 sin(t)와 cos(t)를 곱해 줘야 한다.
tan은 sin/cos이므로 cos는 서로 약분되어 없어지고.. 최종적으로 적분해야 하는 식은 sin(t)^2 dt가 된다.

그리고 sin(x)^2은 삼각함수 덧셈 정리로부터 유도된 반각 공식에 의거하여 (1-cos(2x))/2로 처분 가능하다. 이게 적분하기 훨씬 더 편하다.
사이클로이드를 기술하는데 각도가 t건 2t건 그건 중요하지 않으므로 x축의 궤적은 각도 θ에 대해 C*(θ-sin(θ))가 도출되며, y축은 이미 sin(t)^2의 간소화형인 C*(1-cos(θ))로 답이 나와 있다. 이것으로 유도 끝..

우리는 사이클로이드의 x, y축별 매개변수식에서 아주 중요한 특성을 하나 주목하게 된다. 바로 매개변수 t에 대해서 x축의 궤적 함수는 y축 궤적 함수의 부정적분이라는 것이다! 반대로 y축의 궤적 함수는 x축 궤적 함수의 도함수이다. 애초에 이런 관계였구나..

이것은 사이클로이드를 매개변수가 아닌 양/음함수로만 기술할 때는 간파할 수 없는 특성이다. 그리고 이런 특성이 존재하기 때문에 사이클로이드의 x, y 궤적은 앞서 제시되었던 미분 방정식을 만족하고 해당 범함수에 대해 오일러 방정식을 충족하고, 최단 강하 곡선 역할까지 한다고 볼 수 있다.

이상이다. 현수선하고 뭔가 비슷한 구석이 있는 것 같으면서도 현수선을 분석하는 것보다는 확실히 더 어려운 것 같다..;; 현수선이 역학적으로 자연스럽고 안정된 궤적이라면, 사이클로이드는 좀 더 인위적이고 최적화된 듯한 느낌이 드는 궤적이다. 다만, 둘 다 중력이 있기 때문에 존재한다는 공통점이 있다.

얘는 중력이 작용하는 지표면에서 공을 몇 도로 던졌을 때 제일 멀리 날아가냐 하는 문제와도 비슷한 느낌이 든다. 삼각함수의 반각 공식이 쓰였다는 공통점도 있으니 말이다. 하지만 사이클로이드는 그것보다도 당연히 훨씬 더 복잡하고 난해하다.
글이 이것만으로도 너무 길어졌으니, 최단 강하 증명의 다른 풀이법 등 나머지 얘기는 다음 시간에 계속하도록 하겠다. =_=;;

Posted by 사무엘