Physics World(물리,수학)

Euler-Cauchy equation(오일러-코시 방정식)Ch 2.5 본문

공학수학(미분방정식)/ODE(이변수미방)

Euler-Cauchy equation(오일러-코시 방정식)Ch 2.5

물리영이 2020. 10. 2. 18:00
728x90
반응형

-이전 글

physical-world.tistory.com/18

 

Homogeneous Linear ODE with constant coefficients(상수 계수 제차 2계 미방) Ch2.2

2계 미방부터는 해법이 간단하지 않다. 따라서, 다음과 같은 간단한 상수계수 제차2계미분방정식에 대해 알아본다. $$y''+ay'+by=r(x)$$ 제차 : \(r(x)=0\) $$\Rightarrow y''+ay'+by=0$$ 미분해서 상수만 튀어나..

physical-world.tistory.com

Euler-Cauchy equation(오일러-코시 방정식)Ch 2.5

다음과 같은 방정식을 Euler-Cauchy equation이라고 부른다.

$$x^2y''+axy'+by=0, \: a,b:상수$$

$$기법:y=x^m$$

$$\Rightarrow y'=mx^{m-1},y''=m(m-1)x^{m-2}$$

방정식에 대입하면,

$$m(m-1)x^{m-2}x^2+amx^{m-1}x+bx^m=0$$

$$ \Rightarrow m(m-1)^2x^m+amx^m+bx^m=0$$

$$x^m \neq0$$

$$\Rightarrow m(m-1)+am+b=m^2+(a-1)m+b=0$$

다음과 같은 방정식을 얻어낼 수 있다.

그리고 이 방정식의 해가 \(m_1,m_2\)라면, \(y_1=x^{m_1},y_2=x^{m_2}\)이다.

이차방정식이므로 당연히, 세가지 경우로 나눌 수 있다.

두 실근, 중근, 두 허근.

 

① 두 실근(\(m=m_1,m_2\))

중첩의 원리에 의해 일반해

$$y=c_1x^{m_1}+c_2x^{m_2}$$

 

②중근( \(m_1=m_2=\frac{1-a}{2}\) )

중근이므로 앞선 상수계수 미분방정식과 유사하게,

계수 감소법(Reduction of order)을 이용하여 \(y_2\)을 구한다.

참고:

https://physical-world.tistory.com/17

 

Linearly independent, Reduction of order(선형 독립, 계수 감소법)Ch2.1

2계 이상의 미분방정식은 해당하는 \(y_1,y_2\)처럼 해가 두 개이상 존재한다. 여기서 \(y_1,y_2\)의 선형 독립인지 선형 종속인지가 중요하다. 왜 중요한지는 글 하단 부분에서 알 수 있다. 정의: ①선

physical-world.tistory.com

$$y_1=x^{\frac{1-a}{2}},y_2=uy_1$$

오일러 코시방정식의 꼴을 바꿔보면

$$y''+\frac{a}{x}y'+\frac{b}{x^2}y=0$$

계수감소법에 의해, 

$$y_2=y_1 \int \frac{1}{y_1^2}e^{-\int p dx}dx$$

$$u=\int\frac{1}{x^{(\frac{1-a}{2})2}}e^{-\int\frac{a}{x} dx}dx=\int\frac{e^{-alnx}}{x^{1-a}}dx=lnx$$

일반해

$$y=(c_1+c_2lnx)x^{\frac{1-a}{2}}$$

여기서 \(m= \frac{1-a}{2} \)는 암기하려면 어려울 수 있으니, m에 대한 이차방정식에서 중근을 직접구하면 된다.

 

③두 허근(\(m=\frac{1-a}{2} \pm wi\))

m에 대한 이차방정식에서 두 허근을 근의공식으로 구하면 위와 같다.

일반해

$$y=c_1x^{\frac{1-a}{2}+wi}+c_2x^{\frac{1-a}{2}-wi}$$

위 일반해는 엄밀히 복소수 일반해이다.

 

여기서 선형결합으로 해를 조작해주면,

$$y_1,y_2=x^{\frac{1-a}{2}}x^{\pm wi}=x^{\frac{1-a}{2}}e^{\pm wi lnx}$$

\( e^{iwlnx}=cos(wlnx)+isin(wlnx)\)이므로

$$\frac{y_1+y_2}{2}=x^{\frac{1-a}{2}}cos(wlnx)$$

$$\frac{y_1+y_2}{2i}=x^{\frac{1-a}{2}}sin(wlnx)$$

이 두 식이 또한 독립이 되므로

실수 일반해

$$y=x^{\frac{1-a}{2}}[c_1cos(wlnx)+c_2sin(wlnx)]$$

이다.

 

2계 상수 미분방정식처럼,

복소수가 유용한지, 실수가 유용한지는 초기조건 또는 미분방정식마다 다르므로 적절히 선택해야한다.

참고: 2계 상수 미방

physical-world.tistory.com/18

 

Homogeneous Linear ODE with constant coefficients(상수 계수 제차 2계 미방) Ch2.2

2계 미방부터는 해법이 간단하지 않다. 따라서, 다음과 같은 간단한 상수계수 제차2계미분방정식에 대해 알아본다. $$y''+ay'+by=r(x)$$ 제차 : \(r(x)=0\) $$\Rightarrow y''+ay'+by=0$$ 미분해서 상수만 튀어나..

physical-world.tistory.com

 

※ \(x^2y''+axy'+by=0\) 꼴을 기억하는 법: 두 번 미분엔 두 번 \(x\)곱, 한 번 미분엔 한 번 \(x\)곱, 미분안한 항엔 상수만 있는 꼴이다. 이런 꼴이므로 \(x^m\)기법을 사용할 수 있는 것이다.

 

※ 양자역학 2D free space(V=0)에서 Near field(\(r \to 0\))일 때, 본 공식을 사용한다.

$$r \equiv l \xi $$

$$[\frac{\partial^2}{\partial \xi^2}+\frac{\partial}{\xi \partial \xi}-\frac{m^2}{\xi^2}-\xi^2+2\epsilon]\psi(\xi)=0$$

에서 

\(\xi \to 0\) 일 때, 근사하기를

$$\psi ''+\frac{\psi'}{\xi}-\frac{m^2}{\xi^2}\psi=0 \: for \: Near \: Field$$

$$\Rightarrow \psi(\xi) \thicksim \xi^{|m|}$$

 

-다음 글

physical-world.tistory.com/22

 

Existence and Uniqueness of solutions, Wronskian Ch2.6

이계 미분 방정식의 존재성과 유일성 정리에 대해 알아보자. 정리의 순서는 미분방정식의 책마다 다르겠지만, Kreysig 공업수학 책을 기준으로 정리한다. 우선 존재성과 유일성을 판단하기 이전

physical-world.tistory.com

 

반응형
Comments