권찡's 공학이야기

1계 미분방정식(완전미분방정식형,Exact Differential Equation) 본문

mathematics/engineering math

1계 미분방정식(완전미분방정식형,Exact Differential Equation)

권찡 2018. 10. 28. 17:34

먼저 완전미분방정식이 어떤 것인가 정의해야합니다.



전미분 방정식(Total Differential Equation)에 대해서 알고 있다면 완전 미방을 이해하기 쉬울것입니다.



이변수 함수의 u(x,y)=c 의 전미분을 하면 아래의 식으로 정의됩니다.


du%5Cquad%20%3D%5Cfrac%20%7B%20%5Cpartial%20u%20%7D%7B%20%5Cpartial%20x%20%7Ddx%2B%5Cfrac%20%7B%20%5Cpartial%20u%20%7D%7B%20%5Cpartial%20y%20%7Ddy%20


이 식을 통해서 1계 미분방정식의 형태로 바꿀수가 있는데


%5Cfrac%20%7B%20dy%20%7D%7B%20dx%20%7D%3D-%5Cfrac%20%7B%20%5Cpartial%20u%2F%5Cpartial%20x%20%7D%7B%20%5Cpartial%20u%2F%5Cpartial%20y%20%7D%3Df(x%2Cy)%20



이때 아래와 같이 가정해봅시다.


M(x%2Cy)%3D%5Cfrac%20%7B%20%5Cpartial%20u%20%7D%7B%20%5Cpartial%20x%20%7D%5Cquad%20%2C%5Cquad%20N(x%2Cy)%3D%5Cfrac%20%7B%20%5Cpartial%20u%20%7D%7B%20%5Cpartial%20y%20%7D%20



이 식을 전미분방정식에 대입하면



du%3DM(x%2Cy)dx%2BN(x%2Cy)dy%20



이러한 전미분방정식으로 됩니다.



위 전미분 방정식이 M(x%2Cy)%3D%5Cfrac%20%7B%20%5Cpartial%20u%20%7D%7B%20%5Cpartial%20x%20%7D%5Cquad%20%2C%5Cquad%20N(x%2Cy)%3D%5Cfrac%20%7B%20%5Cpartial%20u%20%7D%7B%20%5Cpartial%20y%20%7D%20 를 만족하는 이변수함수 u가 존재하면 이를 두고 와전 미분방정식이라고 합니다.








우리가 마주하는 미방이 완전미방인지 아닌지 아는 방법은 간단합니다.




%5Cfrac%20%7B%20%5Cpartial%20M%20%7D%7B%20%5Cpartial%20y%20%7D%3D%5Cfrac%20%7B%20%5Cpartial%20N%20%7D%7B%20%5Cpartial%20x%20%7D%5C%5C%20%EB%8B%A4%EC%8B%9C%EB%A7%90%ED%95%B4%5Cquad%20%5Cfrac%20%7B%20%5Cpartial%20u%20%7D%7B%20%5Cpartial%20y%5Cpartial%20x%20%7D%3D%5Cfrac%20%7B%20%5Cpartial%20u%20%7D%7B%20%5Cpartial%20x%5Cpartial%20y%20%7D%20


위의 조건이 만족하면 완전 미분방정식이 성립하는 것이죠.






 예를 들어서 완전미방을 설명하도록 하겠습니다.


(%5Ccombi%20%5E%7B%203%20%7D%7B%20x%20%7D%2B3x%5Ccombi%20%5E%7B%202%20%7D%7B%20y%20%7D)dx%5Cquad%20%2B%5Cquad%20(3%5Ccombi%20%5E%7B%202%20%7D%7B%20x%20%7Dy%5Cquad%20%2B%5Ccombi%20%5E%7B%203%20%7D%7B%20y%20%7D)dy%5Cquad%20%3D0%20


위와 같은 미분방정식이 있을때 완전 미방임을 확인해 봅시다.


dx 앞에 식을 M(x,y) 로 보고 dy앞에 식을 N(x,y)라고 본다면


M을 y로 편미분, N을 x로 편미분해서 확인합니다.

%5C%5C%20N%EC%9D%84%5Cquad%20x%EB%A1%9C%5Cquad%20%ED%8E%B8%EB%AF%B8%EB%B6%84%ED%95%9C%5Cquad%20%EA%B0%92%5Cquad%20%3A%5Cquad%206xy%5C%5C%20M%EC%9D%84%5Cquad%20y%EB%A1%9C%5Cquad%20%ED%8E%B8%EB%AF%B8%EB%B6%84%ED%95%9C%5Cquad%20%EA%B0%92%5Cquad%20%3A%5Cquad%206xy%20


값은 같은 것을 확인했습니다. 이럴 경우 완전미분방정식의 해법으로 푸는 것입니다.


일반해의 형태부터 보겠습니다.


%5Cint%20%7B%20M(x%2Cy)%20%7Ddx%2B%5Cint%20%7B%20%5C%7BN(x%2Cy)-%5Cfrac%20%7B%20%5Cpartial%20%20%7D%7B%20%5Cpartial%20y%20%7D%5Cint%20%7B%20M(x%2Cy)dx%20%7D%20%7D%5C%7D%5Cquad%20dy%5Cquad%20%3D%5Cquad%20c(%EC%A0%81%EB%B6%84%EC%83%81%EC%88%98)%20



복잡해보이는 식입니다.



그러나 위 식을 자세하게 살펴보면 위식은 원래의 이변수 함수 u(x,y)를 구해내는 식임을 알수 있습니다.  이때 u함수를 potential function이라고 부르며


영어 의미 그대로 원래의 본함수라는 의미입니다. 즉 완전미방의 해는 potential function이다 라는 것이죠.




쉬운것 하나 예를 들어보겠습니다.


(%5Ccombi%20%5E%7B%203%20%7D%7B%20x%20%7D%2B3x%5Ccombi%20%5E%7B%202%20%7D%7B%20y%20%7D)dx%5Cquad%20%2B%5Cquad%20(3%5Ccombi%20%5E%7B%202%20%7D%7B%20x%20%7Dy%5Cquad%20%2B%5Ccombi%20%5E%7B%203%20%7D%7B%20y%20%7D)dy%5Cquad%20%3D0%20



이 식이 완전미방인지 확인했으니 potential function를 구해봅시다.



M(x,y)가 potential function의 x편미분식이고, N(x,y)가 potential function의 y편미분식입니다.



반대로 연산한다고 생각해보죠. 각각 변수에 맞게 적분 한다는 것이죠.



M%EC%9D%84%5Cquad%20x%EB%A1%9C%5Cquad%20%EC%A0%81%EB%B6%84%5Cquad%20%3A%5Cquad%20%5Cfrac%20%7B%201%20%7D%7B%204%20%7D%5Ccombi%20%5E%7B%204%20%7D%7B%20x%20%7D%2B%5Cfrac%20%7B%203%20%7D%7B%202%20%7D%5Ccombi%20%5E%7B%202%20%7D%7B%20x%20%7D%5Ccombi%20%5E%7B%202%20%7D%7B%20y%20%7D%5C%5C%20N%EC%9D%84%5Cquad%20y%EB%A1%9C%5Cquad%20%EC%A0%81%EB%B6%84%5Cquad%20%3A%5Cquad%20%5Cfrac%20%7B%203%20%7D%7B%202%20%7D%5Ccombi%20%5E%7B%202%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%204%20%7D%5Ccombi%20%5E%7B%204%20%7D%7B%20y%20%7D%20





공통된 부분은 중복된다 보고 따로 있는것은 그냥 씁니다.



%5Cfrac%20%7B%201%20%7D%7B%204%20%7D%5Ccombi%20%5E%7B%204%20%7D%7B%20x%20%7D%2B%5Cfrac%20%7B%203%20%7D%7B%202%20%7D%5Ccombi%20%5E%7B%202%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%204%20%7D%5Ccombi%20%5E%7B%204%20%7D%7B%20y%20%7D%20




이것이 potential function 함수가 되는 것이죠.


즉 위 미방의 해는


u(x,y)=c(적분상수)





조금 복잡해보이는 미분방정식을 보여드리면


(%5Ccombi%20%5E%7B%202%20%7D%7B%20y%20%7Dcosx-%5Cquad%203%5Ccombi%20%5E%7B%202%20%7D%7B%20x%20%7Dy-2x)dx%2B(2ysinx-%5Ccombi%20%5E%7B%203%20%7D%7B%20x%20%7D%2Bln%5Ccombi%20%7B%20y%20%7D)dy%5Cquad%20%3D0%20



이 문제의 해를 구해봅시다.



앞서 배운 변수분리법으로 풀리지 않을것 같습니다. 그렇다면 완전미방인지 확인해보겠습니다.



(2ysinx-%5Ccombi%20%5E%7B%203%20%7D%7B%20x%20%7D%2Bln%5Ccombi%20%7B%20y%20%7D)%EC%9D%84%5Cquad%20x%EB%A1%9C%5Cquad%20%EB%AF%B8%EB%B6%84%5Cquad%20%5Cto%20%5Cquad%202ycosx-3%5Ccombi%20%5E%7B%202%20%7D%7B%20x%20%7D%5C%5C%20(%5Ccombi%20%5E%7B%202%20%7D%7B%20y%20%7Dcosx-%5Cquad%203%5Ccombi%20%5E%7B%202%20%7D%7B%20x%20%7Dy-2x)%EC%9D%84%5Cquad%20y%EB%A1%9C%5Cquad%20%EB%AF%B8%EB%B6%84%5Cquad%20%5Cto%20%5Cquad%202ycosx%5Cquad%20-3%5Ccombi%20%5E%7B%202%20%7D%7B%20x%20%7D%20



완전미분방정식이네요.



그렇다면 해를 구하는 것은 간단합니다. potential function를 구하면 됩니다.



M%EC%9D%84%5Cquad%20x%EB%A1%9C%5Cquad%20%EC%A0%81%EB%B6%84%3A%5Cquad%20%5Ccombi%20%5E%7B%202%20%7D%7B%20y%20%7Dsinx-%5Ccombi%20%5E%7B%203%20%7D%7B%20x%20%7Dy-%5Ccombi%20%5E%7B%202%20%7D%7B%20x%20%7D%5C%5C%20N%EC%9D%84%5Cquad%20y%EB%A1%9C%5Cquad%20%EC%A0%81%EB%B6%84%3A%5Cquad%20%5Ccombi%20%5E%7B%202%20%7D%7B%20y%20%7Dsinx-%5Ccombi%20%5E%7B%203%20%7D%7B%20x%20%7Dy%2Byln%5Ccombi%20%7B%20y%20%7D-y%5C%5C%20%5C%5C%20%5C%5C%20%5C%5C%20%5Ccombi%20%5E%7B%202%20%7D%7B%20y%20%7Dsinx-%5Ccombi%20%5E%7B%203%20%7D%7B%20x%20%7Dy-%5Ccombi%20%5E%7B%202%20%7D%7B%20x%20%7D%2Byln%5Ccombi%20%7B%20y%20%7D-y%5Cquad%20%3D%5Cquad%20%ED%8F%AC%ED%85%90%EC%85%9C%5Cquad%20%ED%95%A8%EC%88%98%5Cquad%20U%5C%5C%20%5C%5C%20%EC%9D%B4%EA%B2%83%EC%9D%B4%5Cquad%20%EC%9D%BC%EB%B0%98%ED%95%B4%EC%9D%B4%EA%B8%B0%EB%95%8C%EB%AC%B8%EC%97%90%5Cquad%20%5C%5C%20%5C%5C%20%5Ccombi%20%5E%7B%202%20%7D%7B%20y%20%7Dsinx-%5Ccombi%20%5E%7B%203%20%7D%7B%20x%20%7Dy-%5Ccombi%20%5E%7B%202%20%7D%7B%20x%20%7D%2Byln%5Ccombi%20%7B%20y%20%7D-y%5Cquad%20%3DU(x%2Cy)%5Cquad%20%3DC%20




완전 간단한 미분방정식 풀이네요.