점화식, 미분방정식, 선형대수학

수학노트
http://bomber0.myid.net/ (토론)님의 2010년 2월 12일 (금) 16:58 판 (피타고라스님이 이 페이지의 이름을 점화식 미분방정식 선형대수학로 바꾸었습니다.)
둘러보기로 가기 검색하러 가기

07 점화식

 

 

선형점화식

\(pa_{n+2} + qa_{n+1} + ra_n = 0\) 꼴의 점화식

점화식의 해는 2차원 벡터공간을 이루므로, 두 개의 선형독립인 해를 찾으면 된다.

 

특성방정식 \(px^2 + qx + r = 0 \) 가 서로 다른 두 근을 \(\alpha, \beta\) 를 갖는 경우.

수열 \(\alpha^{n}\)와 \(\beta^{n}\)는 선형독립인 두 해가 된다.

따라서 점화식의 일반해는 그 선형결합 \(a_n = A\alpha^{n} + B\beta^{n}\) 꼴로 주어진다.

 

특성방정식 \(px^2 + qx + r = 0 \) 가 중근 \(\alpha\) 를 가지는 경우

수열 \(\alpha^{n} \)와 \(n\alpha^{n} \)는 선형독립인 두 해가 된다.

따라서 점화식의 일반해는 \(a_n = A\alpha^{n} + Bn\alpha^{n}\) 꼴로 주어진다.

(증명)

수열 \(n\alpha^{n} \)이 점화식의 해가 되는지를 확인하면 된다.

\(p(n+2)\alpha^{n+2} + q(n+1)\alpha^{n+1} + rn\alpha^{n} =(n(p\alpha^2+q\alpha+r)+(2p+q)) \alpha^n=0\)

여기서 \(px^2 + qx + r = 0 \)가 중근 \(\alpha\)을 가지므로 \(p\alpha^2+q\alpha+r=0, 2p+q=0\)이다.

 

 

상수계수 이계 선형미분방정식

\(ay''+by'+cy=0\)

초기값을 가지는 이계 선형미분방정식의 해는 2차원 벡터공간을 이루므로, 두 개의 선형독립인 해를 찾으면 된다.

 

특성방정식 \(ax^2 + bx + c = 0 \) 가 서로 다른 두 근을 \(\alpha, \beta\) 를 갖는 경우.

함수 \(e^{\alpha t}\)와 \(e^{\beta t}\)는 선형독립인 두 해가 된다.

따라서 점화식의 일반해는 그 선형결합 \(y(t) = Ae^{\alpha t} + Be^{\beta t}\) 꼴로 주어진다.

 

특성방정식 \(ax^2 + bx + c = 0 \) 가 중근을 \(\alpha\) 를 갖는 경우.

 

함수 \(e^{\alpha t}\)와 \(te^{\beta t}\)는 선형독립인 두 해가 된다.

 

따라서 점화식의 일반해는 그 선형결합 \(y(t) = Ae^{\alpha t} + Bte^{\alpha t}\) 꼴로 주어진다.

 

 

(증명)

\(ax^2 + bx + c = 0 \)가 중근 \(\alpha\)을 가지므로 \(a\alpha^2+b\alpha+c=0, 2a+b=0\)이다.

\(y(t) = te^{\alpha t}\) 라 하자.

\(y'(t) = (\alpha t+1)e^{\alpha t}\)

\(y''(t) = (\alpha^2 t+2\alpha)e^{\alpha t}\)

미분방정식에 대입하면,

\(ay''(t)+by'(t)+cy =\{a(\alpha^2 t+2\alpha)+b(\alpha t+1)+ct\}e^{\alpha t}=\{(a\alpha^2 +b\alpha+c)t+(2a \alpha +b)\}e^{\alpha t}=0\) ■

 

 

 

 

 

벡터공간

\(f(0)=f(\pi)=0\) 을 만족시키는 \([0,\pi]\)에서 정의된 함수공간

 

내적

\((f,g)=\int_0^{\pi}f(x)g(x)\,dx\)

 

선형사상

\(L[y]=y''\)

 

(정리)

\(L[y]=y''\)은 Hermitian operator 이다.

즉 \((L[f],g)=(f'',g)=(f,g'')=(f,L[g])\)

(증명)

\((Lf,g)=\int_0^{\pi}f''(x)g(x)\,dx=[f'(x)g(x)]_0^{\pi}-\int_0^{\pi}f'(x)g'(x)\,dx=-[f(x)g'(x)]_0^{\pi}+\int_0^{\pi}f(x)g''(x)\,dx=(f,Lg)\)

 

http://www.wolframalpha.com/input/?i=+sin+3x+*+sin+4x+

http://www.wolframalpha.com/input/?i=integrate_0^(pi)+1/2+(cos(x)-cos(7+x))+dx