공업수학 시리즈 6편: 고차 미분방정식과 초기값 문제

지금까지는 주로 1차와 2차 미분방정식을 봤습니다. 하지만 실제 시스템을 더 정교하게 모델링하면 3차 이상 고차 방정식도 자연스럽게 등장합니다. 핵심은 복잡해 보여도 기본 원리는 그대로라는 점입니다.

고차 미분방정식이란

최고차 도함수가 3차 이상인 식을 고차 미분방정식이라고 부릅니다. 예를 들면

$y''' - y'' - 2y' = 0$

같은 식입니다.

상수계수 선형 형태라면 여전히 지수함수 $e^{rx}$를 넣어 특성방정식을 만들 수 있습니다.

왜 초기조건 개수가 차수와 같을까

$n$차 미분방정식의 일반해에는 보통 독립적인 상수 $n$개가 들어갑니다. 따라서 유일한 해를 정하려면 조건도 $n$개가 필요합니다.

예를 들어 3차 방정식이면 보통

$y(0), \quad y'(0), \quad y''(0)$

처럼 세 개의 초기조건이 필요합니다.

이 관점은 매우 중요합니다. 미분방정식 풀이란 결국 상수들을 조건으로 결정해 하나의 실제 운동을 고르는 과정이기 때문입니다.

손으로 풀어보는 예제

다음 문제를 보겠습니다.

$y''' - y'' = 0$

특성방정식은

$r^3 - r^2 = 0$

즉

$r^2(r-1)=0$

이므로 근은 $r=0$이 중근, $r=1$이 단순근입니다. 따라서 일반해는

$y = C_1 + C_2 x + C_3 e^x$

입니다.

이제 조건을

$y(0)=1, \quad y'(0)=0, \quad y''(0)=2$

라고 두겠습니다.

미분하면

$y' = C_2 + C_3 e^x, \quad y'' = C_3 e^x$

입니다. $x=0$을 넣으면

$C_1 + C_3 = 1$

$C_2 + C_3 = 0$

$C_3 = 2$

이므로

$C_1 = -1, \quad C_2 = -2, \quad C_3 = 2$

따라서 해는

$y = -1 - 2x + 2e^x$

가 됩니다.

공학 응용

제어 시스템의 전달함수

제어공학에서는 3차, 4차 이상의 시스템이 흔합니다. 모터, 기어, 센서, 필터가 결합되면 동역학 차수가 올라갑니다.

보와 판의 변형

구조역학에서는 4차 미분방정식이 나오는 경우가 많습니다. 이는 곡률과 굽힘 모멘트가 연결되기 때문입니다.

신호 처리 필터

고차 필터는 입력에 대한 응답을 고차 미분방정식이나 동등한 상태공간 시스템으로 표현할 수 있습니다.

초기값 문제와 경계값 문제

입문 단계에서는 초기값 문제를 먼저 봅니다. 이는 한 시점의 상태가 모두 주어지는 경우입니다.

반면 경계값 문제는 서로 다른 위치나 시간 끝점에서 조건이 주어집니다. 예를 들어 막대 양끝의 온도나 보의 양끝 처짐 조건이 여기에 해당합니다. 뒤에서 편미분방정식을 배울 때 더 자주 나오게 됩니다.

자주 하는 실수

조건 수를 덜 준다

$n$차 식에 조건이 충분하지 않으면 해가 하나로 정해지지 않습니다.

중근 처리를 놓친다

중근이 있으면 $C_1 + C_2 x$처럼 $x$가 곱해진 항이 추가됩니다. 이 부분을 자주 빼먹습니다.

미분 계수를 잘못 정리한다

고차식에서는 미분 실수 하나가 전체 상수 계산을 망칩니다. 일반해를 쓰고, 차례대로 미분한 뒤, 마지막에 조건을 대입하는 흐름을 지키는 편이 좋습니다.

한 줄 요약

고차 미분방정식도 기본 원리는 같으며, 차수가 높아질수록 필요한 초기조건 수도 함께 늘어납니다.

다음 편 예고

다음 글에서는 여러 변수의 움직임을 한꺼번에 다루는 연립 미분방정식과 행렬 표현의 필요성을 살펴보겠습니다.

참고자료

• Erwin Kreyszig, Advanced Engineering Mathematics, 10th Edition
• Earl A. Coddington, An Introduction to Ordinary Differential Equations
• Gilbert Strang, Linear Algebra and Differential Equations 강의 자료