오일러 미분 방정식의 풀이
📂상미분방정식오일러 미분 방정식의 풀이
정의
다음과 같은 꼴의 미분 방정식을 오일러 미분 방정식 혹은 오일러-코시 방정식이라 한다.
a2x2dx2d2y+a1xdxdy+a0y=0
설명
우변이 0이 아닌 비동차 방정일 경우에는 x=ez로 치환해서 풀면 된다.
풀이
계산의 편의를 위해 (1)의 양 변을 a2로 나누고 나머지 두 항의 계수를 다시 a1, a0라고 하자. 그러면
x2dx2d2y+a1xdxdy+a0y=0
미분 방정식을 잘 보면 두번미분하고 2차항을 곱한 항과 1번 미분하고 1차항을 곱한 항과 원래 함수를 더해서 0이된다. 따라서 해를 다음과 같다고 둘 수 있다.
y=xr
미분 방정식에 대입하면
r(r−1)xr+a1rxr+a0xr=0[r(r−1)+a1r+a0]xr=0[r2−(a1−1)r+a0]xr=0
xr=0이므로 r2−(a1−1)r+a0=0이다. 이는 간단한 2차 방정식으로 그 해는
r=2−(a1−1)±(a1−1)2−4a0
두 해를 각각 r1, r2라고 하자. 두 근의 상태에 따라서 미분 방정식의 해가 달라진다.
Case 1. r1, r2가 서로 다른 실수
방정식의 두 해는 y1=xr1과 y2=xr2이다. 론스키안을 확인해보면
W[y1,y2]=(r2−r1)xr1+r2−1
r1=r2이므로 x>0일 때 반드시 W[x1,r2]=0임을 알 수 있다. 따라서 두 해는 기본해집합을 이루므로 일반해는
y=c1xr1+c2xr2,x>0
■
Case 2. r1, r2가 서로 같은 실수
위 경우에서 y1, y2이므로 세컨드 솔루션을 찾아야한다. 첫번째 해를 y1=xr1이라 하고 미분 연산자 L을 아래와 같다고 하자.
L[y]=x2y′′+xy′+y
그러면
L[xr]=[r2−(a1−1)r+a0]xr=0
이때 r에 대한 2차 방정식이 중근을 가지는 경우이므로 아래와 같이 완전제곱꼴로 표현할 수 있다.
L[xr]=(r−r1)2xr=0
0을 미분하면 0이므로 양변을 r로 미분하면
∂r∂L[xr]=0
이고 x와 r의 미분 순서를 바꿔도 상관 없으므로
∂r∂L[xr]=L[∂r∂xr]=L[xrlnx]=0
따라서 y2=xr1lnx가 두번째 해이다. 론스키안을 계산해보면 W[xr1,x1lnx]=x2r1−1=0이므로 두 해는 기본해집합을 이룬다. 따라서 일반해는
y=c1xr1+c2xr1lnx,x>0
■
Case 3. r1, r2가 서로 다른 복소수
r1=λ+iμ, r2=λ−iμ라고 하자. 그러면 두 해는
y1=xλ+iμ,y2=xλ−iμ
따라서 기본해는
y=c1xλ+iμ+c2xλ−iμ
하지만 복소함수의 경우 삼각함수로 표현하는 것이 일반적이다. 오일러 공식에 의해서 아래의 식이 성립한다.
xλ+iμ=xλxiμ=xλelnxiμ== xλeiμlnx xλ[cos(μlnx)+isin(μlnx)],x>0
따라서 일반해는 복소수 상수 c1, c2에 대해서 다음과 같이 표현된다.
y=c1xλcos(μlnx)+c2xλsin(μlnx),x>0
cos과 sin은 독립적이므로 론스키안은 계산해보지 않아도 반드시 0이 되지 않음을 알 수 있다. 굳이 계산해보면
W[xλcos(μlnx),xλsin(μlnx)]=μx2λ−1=0
■