권찡's 공학이야기

특수한 미분방정식(리카티 미분방정식) 본문

mathematics/engineering math

특수한 미분방정식(리카티 미분방정식)

권찡 2018. 10. 30. 22:16

리카티(Riccati) 미분방정식은 아래와 같은 식입니다.


y'%5Cquad%20%3D%5Cquad%20p(x)%5Ccombi%20%5E%7B%202%20%7D%7B%20y%20%7D%2Bq(x)y%2BR(x)%20




위와 같은 형태의 미분 방정식을 보고 리카티 미분방정식이라고 합니다. 지금까지 정리한 내용으로 이런 미분방정식을 푸는 법이 떠오르지가 않습니다.



상당히 어려운 미분방정식입니다. 왜냐?



일단 풀이가 이미 알고 있는 particular Solution이 있어야 됩니다.




??????




무슨 말이냐면, 하나의 특수해를 알고 있어야 다음의 풀이가 가능하다는 것인데, 이 특수해를 얻는 방법이 따로 있는 것이 아니며, 직관적으로 혹은 여러번의 시도로 얻어야 합니다.



가령 y=2x 넣으면 성립할 것 같은데?? 혹은 y=1 넣으면 성립하는데??



와 같은 방식으로 찾아야 하기 때문입니다.



리카티 미분방정식의 예를 하나 들어볼까요??


y'%5Cquad%20%3D%5Cquad%20%5Cfrac%20%7B%201%20%7D%7B%20x%20%7D%5Ccombi%20%5E%7B%202%20%7D%7B%20y%20%7D%2B%5Cfrac%20%7B%201%20%7D%7B%20x%20%7Dy-%5Cfrac%20%7B%202%20%7D%7B%20x%20%7D%5Cquad%20%EA%B2%BD%EC%9A%B0%5Cquad%20y%3D1%5C%5C%20y'%5Cquad%20%3D%5Cquad%20%5Cfrac%20%7B%201%20%7D%7B%20%5Ccombi%20%5E%7B%202%20%7D%7B%20x%20%7D%20%7D%5Ccombi%20%5E%7B%202%20%7D%7B%20y%20%7D-%5Cfrac%20%7B%201%20%7D%7B%20x%20%7Dy%2B1%5Cquad%20%EA%B2%BD%EC%9A%B0%5Cquad%20y%3Dx%20



이런식으로 하나의 해를 때려 맞춰야 됩니다.






자, 만약에 위와 같이 때려 맞춰서 하나의 해를 얻었다고 해봅시다. 이 해는 y1 이라고 칭하겠습니다.


그리고 일반해를 y=y1+u(x) 라 가정합니다.



이 일반해를 리카티 미분방정식에 대입합니다.




y%3D%5Ccombi%20_%7B%201%20%7D%7B%20y%20%7D%2Bu%5C%5C%20y'%3Dp(x)%5Ccombi%20%5E%7B%202%20%7D%7B%20y%20%7D%2Bq(x)y%2BR(x)%5C%5C%20%5Cdownarrow%20%5C%5C%20%5Ccombi%20_%7B%201%20%7D%7B%20y%20%7D'%2Bu'%5Cquad%20%3Dp(x)(%5Ccombi%20%5E%7B%202%20%7D%7B%20%5Ccombi%20_%7B%201%20%7D%7B%20y%20%7D%20%7D%2B2u%5Ccombi%20_%7B%201%20%7D%7B%20y%20%7D%2B%5Ccombi%20%5E%7B%202%20%7D%7B%20u%20%7D)%2Bq(x)(%5Ccombi%20_%7B%201%20%7D%7B%20y%20%7D%2Bu)%2BR(x)%5C%5C%20%5Cdownarrow%20%EC%9A%B0%EB%B3%80%EC%9D%84%5Cquad%20%EC%A0%95%EB%A6%AC%5C%5C%20%5Ccombi%20_%7B%201%20%7D%7B%20y%20%7D'%2Bu'%5Cquad%20%3D%5Bp(x)%5Ccombi%20%5E%7B%202%20%7D%7B%20%5Ccombi%20_%7B%201%20%7D%7B%20y%20%7D%20%7D%2Bq(x)%5Ccombi%20_%7B%201%20%7D%7B%20y%20%7D%2BR(x)%5D%2B%5B2p(x)%5Ccombi%20_%7B%201%20%7D%7B%20y%20%7D%2Bq(x)%5Du%2B%5Bp(x)%5D%5Ccombi%20%5E%7B%202%20%7D%7B%20u%20%7D%20





이때 y1이 주어진 미분방정식의 해라는 것은 이미 알고 있으므로


%5Ccombi%20_%7B%201%20%7D%7B%20y%20%7D'%3Dp(x)%5Ccombi%20%5E%7B%202%20%7D%7B%20%5Ccombi%20_%7B%201%20%7D%7B%20y%20%7D%20%7D%2Bq(x)%5Ccombi%20_%7B%201%20%7D%7B%20y%20%7D%2BR(x)%5Cquad%20%EA%B0%80%5Cquad%20%EC%84%B1%EB%A6%BD%20



그러면 위의 식에서 u가 없는 항이 정리가 됩니다.


u'%3D%5Cquad%20%5B2p(x)%5Ccombi%20_%7B%201%20%7D%7B%20y%20%7D%2Bq(x)%5Du%2B%5Bp(x)%5D%5Ccombi%20%5E%7B%202%20%7D%7B%20u%20%7D%20



정리가 되면 이 식은 베르누이 미분방정식과 같은 형태가 됩니다. 이것을 베르누이 미분방정식 같이 풀어내는 것입니다.



예를 들어서 하나 문제를 풀어보죠



y'%3D2-2xy%2B%5Ccombi%20%5E%7B%202%20%7D%7B%20y%20%7D%20 



무엇을 대입해야 위식이 만족할까요???



답은 2x입니다.



그렇다면 일반해는


y%3D2x%2Bu(x)%3D2x%2Bu%5C%5C%20y'%3D2%2Bu'%20 


위 식을 미분방정식에 대입해봅시다.


y'%3D2-2xy%2B%5Ccombi%20%5E%7B%202%20%7D%7B%20y%20%7D%5C%5C%202%2Bu'%3D2-2x(2x%2Bu)%2B%5Ccombi%20%5E%7B%202%20%7D%7B%20(2x%2Bu)%20%7D%5C%5C%20%5Cto%20u'%3D2xu%2B%5Ccombi%20%5E%7B%202%20%7D%7B%20u%20%7D%5C%5C%20%5Cto%20u'-2xu%3D%5Ccombi%20%5E%7B%202%20%7D%7B%20u%20%7D%5Cquad%20%2F%2F%5Cquad%20u%EC%97%90%5Cquad%20%EA%B4%80%ED%95%9C%5Cquad%20%EB%B2%A0%EB%A5%B4%EB%88%84%EC%9D%B4%5Cquad%20%EB%AF%B8%EB%B0%A9%5C%5C%20%5Cto%20%5Ccombi%20%5E%7B%20-2%20%7D%7B%20u%20%7Du'-2x%5Ccombi%20%5E%7B%20-1%20%7D%7B%20u%20%7D%3D1%5Cquad%20%5Cquad%20%2F%2F%5Cquad%20%5Ccombi%20%5E%7B%20-1%20%7D%7B%20u%20%7D%3Dt%EB%A1%9C%5Cquad%20%EC%B9%98%ED%99%98%5Cquad%20%2Ct'%3D%5Cfrac%20%7B%20dt%20%7D%7B%20dx%20%7D%3D-%5Ccombi%20%5E%7B%20-2%20%7D%7B%20u%20%7D%5Cfrac%20%7B%20du%20%7D%7B%20dx%20%7D%3D-%5Ccombi%20%5E%7B%20-2%20%7D%7B%20u%20%7Du'%5C%5C%20%5Cto%20%5Cquad%20-t'-2xt%3D1%5Cquad%20%5C%5C%20%5Cto%20%5Cquad%20t'%2B2xt%3D-1%5C%5C%20%20 



이제 1계 선형 미분방정식처럼 풀면 되겠죠


t'%2B2xt%3D-1%5C%5C%20t(x)%3D%5Ccombi%20%5E%7B%20-%5Cint%20%7B%202x%5Cquad%20dx%20%7D%20%7D%7B%20e%20%7D%5B%5Cint%20%7B%20-%5Ccombi%20%5E%7B%20%5Cint%20%7B%202x%5Cquad%20dx%20%7D%20%7D%7B%20e%20%7Ddx%2BC%20%7D%5D%5C%5C%20t(x)%3D%5Ccombi%20%5E%7B%20-%5Ccombi%20%5E%7B%202%20%7D%7B%20x%20%7D%20%7D%7B%20e%20%7D%5B%5Cint%20%7B%20-%5Ccombi%20%5E%7B%20%5Ccombi%20%5E%7B%202%20%7D%7B%20x%20%7D%20%7D%7B%20e%20%7Ddx%2BC%20%7D%5D%5C%5C%20u%3D%5Cfrac%20%7B%201%20%7D%7B%20t%20%7D%3D%5Cfrac%20%7B%20%5Ccombi%20%5E%7B%20%5Ccombi%20%5E%7B%202%20%7D%7B%20x%20%7D%20%7D%7B%20e%20%7D%20%7D%7B%20C-%5Cint%20%7B%20%5Ccombi%20%5E%7B%20%5Ccombi%20%5E%7B%202%20%7D%7B%20x%20%7D%20%7D%7B%20e%20%7Ddx%20%7D%20%7D%5C%5C%20%20 


따라서 최종적인 일반해는 


y%3D%5Ccombi%20_%7B%201%20%7D%7B%20y%20%7D%2Bu(x)%3D2x%2B%5Cfrac%20%7B%20%5Ccombi%20%5E%7B%20%5Ccombi%20%5E%7B%202%20%7D%7B%20x%20%7D%20%7D%7B%20e%20%7D%20%7D%7B%20C-%5Cint%20%7B%20%5Ccombi%20%5E%7B%20%5Ccombi%20%5E%7B%202%20%7D%7B%20x%20%7D%20%7D%7B%20e%20%7Ddx%20%7D%20%7D%20 



이런 식이 됩니다.