공데셍의 전공 지식 저장소

수학/미분방정식 (Differential Equation)

10. 계수가 상수인 이계 선형 미분방정식 (1편) (Second-Order Linear ODE with Constant Coefficient)

Ball Dessin 2022. 8. 16. 03:34
반응형

일반적인 이계 선형 미분방정식의 풀이에 대해 설명하기 앞서서

계수가 상수인 경우에는 어떻게 푸는지 보일 것이다. 여기선 Homogeneous인 경우를 다룬다.

계수가 상수인 이계 선형 미분방정식은 다음의 형태를 갖는다.

ay+by+cy=0

 

 


 

 

예제1

다음 미분방정식을 풀어라.
yy=0 y(0)=2,y(0)=1

 

 

식을 이항하여 y=y 로 만들면 두 번 미분하여 자기 자신이 되는 y 를 찾는 문제라고 생각할 수 있다.

두 번 미분했을 때 형태가 잘 유지되는 함수가 무엇이 있는지 떠올려 보자.

 

y=sint 는 두 번 미분하면 y=sint 가 되므로 고려해봄직 하다. cost 도 마찬가지이다.

하지만 주어진 미분방정식은 두 번 미분하여 자기자신이 되는 것이지

자기자신에 음수배가 되는 것이 아니므로 사인 코사인은 후보에서 탈락시키자.

하지만 y=et 라면 두 번 미분해도 음수배가 되지 않는다.

이 뿐 아니라 y=et 인 경우에도 두 번 미분하면 마이너스가 두 번으로 상쇄되므로 해가 될 수 있다.

 

그런데 etet상수 c 를 곱해도 대입해보면 여전히 방정식을 만족함을 알 수 있다.

예를 들어 y=5et 인 경우 y=5et 이여서 이 둘은 같으므로 미분방정식의 해가 된다.

심지어 이들을 선형결합하여 만든 y=4et+3et 역시 두 번 미분해보면 같은 식이 나온다.

일반화하면  y=c1et+c2ety=y 의 해가 된다고 볼 수 있을 것이다.

y=ety=et 말고 다른 해가 존재하는지는 아직은 모르지만, 적어도 조금은 일반스러워 보이는 해를 얻었다.

그리고 사실 이게 일반해가 맞다는 것을 나중에 보일 것이다.

 

이제 초기값 문제를 풀자.

y=c1et+c2et 였으므로 y=c1etc2et 이다.

여기에 초기값들 y(0)=2, y(0)=1 을 대입하면

c1=12, c2=32 를 얻을 수 있다.

따라서 초기값 문제의 해는 다음과 같다.

y=12et+32et

 

 

 


 

 

 

위의 예제에서 추측할 수 있는 것은, 계수가 실수인 상수를 갖는 이계 선형 미분방정식

ay+by+cy=0

y=ert 꼴을 갖는 두가지 해의 선형결합으로 구성 할 수 있다는 것이다.

 

 

y=ert 이면 y=rert 이고 y=r2ert 이므로 이들을 미분방정식에 대입하면 다음과 같다.

(ar2+br+c)ert=0

ert0 이므로 양변에 나눠주면, 다음 이차방정식을 푸는 문제로 변한다.

ar2+br+c=0

이 식을 특성방정식(Characteristic Equation) 라고 부른다.

그리고 특성방정식은 이차방정식이므로 다음 세가지 경우로 나뉜다.

 

1. 두 개의 서로 다른 실근을 가짐

2. 하나의 중근을 가짐

3. 한 쌍의 켤레복소수 근을 가짐 (두 개의 허근)

 

 

 


 

 

 

이번 글에서는 서로 다른 실근을 가지는 경우에 대해 설명한다.

특성 방정식이 서로 다른 두 개의 근을 갖는다는 말은

r1r2r1, r2 가 존재하여 다음과 같은 해를 갖는다는 말이다.

(2)y=c1er1t+c2er2t

 

r1,r2 가 특성방정식의 근이라는 가정하에 위 식이 미분방정식의 해가 맞음을 증명하는 것은 어렵지 않다.

{y=c1er1t+c2er2ty=c1r1er1t+c2r2er2ty=c1r12er1t+c2r22er2t

이므로 이들을 ay+by+cy 에 대입하면 다음과 같이 정리된다.

ay+by+c=c1[ar12+br1+c]er1t+c2[ar22+br2+c]er2t

r1,r2특성방정식의 근이므로 ar12+br1+c=0 이고  ar22+br2+c=0 이다.

따라서 ay+by+c=0 이고 방정식을 잘 만족함을 볼 수 있다.

 

 

 

초기값이 다음과 같이 일반적인 형태로 주어진 경우에도 문제가 없는지 살펴보자.

y(t0)=y0,y(t0)=y0

이 값들을 식 (2) 에 대입하여 정리해보면 c1c2 는 다음과 같이 구해진다.

{c1=y0y0r2r1r2er1t0c2=y0r1y0r1r2er2t0

서로 다른 근이라고 하였으므로 r1r2 이므로 분모가 0 이 되지 않는다.

따라서 두 근의 값에 상관없이 c1c2 는 항상 잘 결정된다.

 

 

 


 

 

 

예제2

다음 미분방정식의 일반 해(General Solution)을 구하라.
y+5y+6y=0

 

y=ert 라고 가정하면 특성방정식은 다음과 같다.

r2+5r+6=(r+2)(r+3)=0

두 근은 2,3 이므로 일반해는 다음과 같다.

y=c1e2t+c2e3t

 


 

다음 글에서는 Homogeneous 2nd Order Linear ODE 에 관련된 정리들에 대해 다룰 것이다.

그 이후 특성방정식이 허근을 갖는 경우, 중근을 갖는 경우에 대해 다룬다.

반응형