적분 인자, Integrating factor 란 미분방정식이 완전(exact)하지 않을 때 완전미분방정식으로 만들어주기 위해서 곱하는 함수이다.
정의 - 적분인자
적분인자 구하는 법
PFdx + QFdy = 0 이 식에서 PF=M 이라 하고 QF=N 이라 하자. 이 미분방정식은 완전하기 때문에, 다음 식이 성립한다(Link).
M=PF, N=QF를 대입하자.
1. 우리가 구하는 적분인자 F가 x에 관한 함수라고 하자. |
두번째 항은 F를 y에 대해 미분했으므로 0이 된다. 따라서 다음과 같이 정리할 수 있다. 양변을 적분해주자. 그리고 편의상 P를 y로 편미분하면 Py로, Q를 x로 미분하면 Qx로 표시하겠다. 좌변을 적분하면, 양변에서 로그를 없애면, 이것이 x에 관한 적분인자 F의 답이 된다. |
2. 우리가 구하는 적분인자 F가 y에 관한 함수라고 하자. |
그러면 네번째 항은 F를 x에 대해 미분했으므로 0이 된다. 이후 위에서 했던 방식처럼 하면 F=F(y)를 구할 수 있는 식을 얻는다. |
적분인자를 사용하여 ODE를 풀어보자.
그리고 y(2)=1 이라는 조건이 주어졌다.
먼저 이 미분방정식이 exact 한지 알아보자. 먼저 P=x4+y2 이고, Q=-xy이다. 그래서 Py=2y 이고, Qx = -y가 된다. 두 함수가 같지 않으므로 exact 하지 않다.
그러면 적분인자를 구해보자. 적분인자가 x에 관한 함수인지, y에 관한 함수인지 모르므로 둘 다 구해본다.
먼저, F가 x에 관한 함수인 경우,
그리고 F가 y에 관한 함수인 경우,
적분인자 F가 x에 관함 함수일 경우 적분인자를 구할 수 있었으나, y에 관한 함수는 구할 수 없다. 따라서 적분인자F는 x에 관한 함수이고, 준식에 적분인자 F를 곱하자.
이제 준 식이 exact해졌다. 이제 완전 미분방정식 해법에 따라 풀어보자.
이제 u를 y로 미분한 것은 N이 되어야한다.
따라서 k'(y)=0 이다. 이제 일반해는 다음과 같다.
초기조건이 (2,1)로 주어졌다. c=2-1/8 = 15/8 이 된다. 즉 이 문제의 특수해는 다음과 같다.
'Mathematics > Ordinary Differential Equation' 카테고리의 다른 글
Bernoulli differential equation 베르누이 미분 방정식 (1) | 2012.01.07 |
---|---|
1계 선형 미분방정식, First Order Linear ODE (1) | 2012.01.06 |
완전 미분방정식, Exact ODEs (0) | 2010.05.02 |
변수 분리법, Separable ODEs (3) | 2010.03.31 |
미분 방정식 (Differential Equations) (0) | 2010.03.05 |