Processing math: 100%

반응형

12. 고계 미분방정식 요약

 

y(n)+pn1(x)y(n1)++p1(x)y+p0(x)y=r(x): n계 비동차 선형미분방정식 (r(x)=0이면 동차)

y(n)+an1y(n1)++a1y+a0y=r(x): 상수계수 n계 비동차 선형미분방정식 (r(x)=0이면 동차)

xny(n)+an1xn1y(n1)++a1xy+a0y=r(x): 오일러-코시 방정식  

 

2계 미분방정식 y+f(x)y+g(x)y=0의 알려진 해가 y1이면 계수축소법을 이용하여 다른 해 y2를 구한다.

y2=u(x)y1(x)이고 미분방정식에 직접 대입하여 구하면 y2=y11{y1(x)}2ef(x)dxdx  

 

상수계수 동차미분방정식의 경우 y=eλx를 대입하여 특성방정식을 얻고 오일러-코시 방정식의 경우 y=xm(x>0)를 대입하여 특성방정식을 얻는다.  

 

특성방정식의 해에 따른 동차해의 형태 

 

상수계수 동차 미분방정식

오일러-코시 방정식

서로다른 실근

y=c1eλ1x++cneλnx

y=c1xm1++cnxmn

m중근

(계수축소법을 이용하여 얻음)

y=(c0+c1x++cmxm)ys+

(ysm번 중복되는 해)

y=(c0+c1lnx++cm(lnx)m)ys+

(ysm번 중복되는 해)

켤레복소수근

(오일러공식을 이용하여 얻음)

y=eαx(c1cosβx+c2sinβx)+

y=xαx{c1cos(βlnx)+c2sin(βlnx)}+

 

n계 비동차 선형미분방정식에서 비동차해를 구하는 방법

r(x)가 다항함수, 지수함수, 삼각함수인 경우, 미정계수법을 이용하여 비동차해 yp를 구한다.

r(x)의 항

선택할 yp

kepx

Cepx

kxn(n=0,1,...)

Knxn+Kn1xn1++K1x+K0

ksinqx또는 kcosqx

Kcosqx+Msinqx

 

비동차해를 미정계수법으로 구할 수 있고 r(x)가 2개 이상의 함수들의 합으로 나타나는 경우는 다음의 중첩원리를 이용한다.

중첩원리

y(n)+pn1y(n1)++p1y+p0y=r1의 비동차해를 yp1, y(n)+pn1y(n1)++p1y+p0y=r2의 비동차해를 yp2라고 하자. 그러면 y(n)+pn1y(n1)++p1y+p0y=αr1+βr2의 비동차해는 y=αyp1+βyp2이다.

 

미정계수법으로 대응되는 함수가 동차해의 함수와 중복되는 경우, 중복도에 따라 x 또는 x2을 곱한다.

(오일러-코시 방정식의 경우는 lnx 또는 (lnx)2)

 

미정계수법으로 비동차해를 구할 수 없는 경우, 동차해를 y1,...,yn이라고 하면 비동차해는yp=y1W1Wr(x)dx++ynWnWr(x)dx이고 여기서 W=W[y1,...,yn]y1,...,yn의 론스키안이고 Wi(i=1,...,n)Wi번째 열을 열벡터 (01)로 바꾼 행렬식이다. 오일러-코시 방정식의 경우는 y(n)+an1xy(n1)++a1xn1y+a0xny=r(x)xn형태로 고치고 구한다.

 

 

 

 

 

반응형
Posted by skywalker222