일 | 월 | 화 | 수 | 목 | 금 | 토 |
---|---|---|---|---|---|---|
1 | ||||||
2 | 3 | 4 | 5 | 6 | 7 | 8 |
9 | 10 | 11 | 12 | 13 | 14 | 15 |
16 | 17 | 18 | 19 | 20 | 21 | 22 |
23 | 24 | 25 | 26 | 27 | 28 |
- 공수 문제풀이
- 비제차 상미분 방정식
- homogeneous
- ODEs
- 미분방정식
- 문제풀이
- Problem set 2.7
- Conversation
- Nonhomogeneous ODEs
- 스피치 연습
- 코딩테스트
- speech
- 공수1
- Problem set 1.5
- 공수
- SW역량테스트
- 공학수학
- Advanced Engineering Mathematics
- English
- Problem Set 1.4
- kreyszig
- Homogeneous ODEs
- 공업수학
- vocabulary
- 백준
- 삼성SW역량테스트
- Python
- 미방
- Ode
- 영어회화
- Today
- Total
한걸음
Ch. 1 1st order ODEs : 1.5 Linear ODEs. 본문
Ch. 1 1st order ODEs : 1.5 Linear ODEs.
우당탕탕 할 수 있다!!! 2023. 11. 2. 12:23
드디어 1.5절! 정리시작해봅시다.
Intro.
아래와 같은 형태를 가질 때 First-order ODE는 선형(Linear)이라고 한다.
$$ {y}' + p(x)y = r(x) $$
y의 차수가 1차 이기 때문에 해당 식을 선형이라고 부른다.
y의 차수가 2차 이상이 되는 경우 비선형(NonLinear)이라고 한다.
이번 절에서는 Linear문제의 해를 구하는 것과, Nonlinear를 Linear로 바꾸어 해를 구하는 것을 공부한다.
1. Linear ODE
1.1 Homogeneous Linear ODE
r(x)가 0 인 경우 선형미분방정식은 Homogeneous Linear ODE라고 한다.
$$ {y}' + p(x)y = 0 $$
변수분리법(Sperating Variables)을 활용하여 적분하면
$$ \frac{dy}{y} = - p(x) dx , y(x) = ce^{-\int{p(x)dx}} $$
가 된다.
1.2 Nonhomogeneous Linear ODE
Nonhomogeneous 의 경우에는 Integrating Factor를 구해서 전개하면 된다.
임의의 Integrating Factor F(x)를 곱하면,
$$ F{y}' + pFy = rF $$
${y}' = \frac{dy}{dx}$ 로 표현하면,
$$ Fdy + (pFy - rF)dx = 0 $$
우리가 앞서 1.4절에서 공부했던 전미분 형태의 식이 된다. Exactness를 위해서는 다음과 같은 조건을 만족해야 한다.
$$ \frac{\partial (pFy - rF)}{\partial y} = \frac{\partial F}{\partial x} $$
$$ \frac {1}{F} = F' $$
$$ F = e^{\int{pdx}} $$
이 떄, $ h \equiv \int{pdx} $ 라고 하면,
$$ F = e^{h} $$
이고, 치환하면,
$$ e^h y' + h' e^h y = re^h $$
이다. 여기에서 잘 관찰하면 약간의 트릭을 발견할 수 있다. LHS 의 형태가 $ e^h y $ 를 x에 대해서 미분한 것과 같은 형태이므로,
$$ (e^h y)' = re^h $$
이다.
적분하면,
$$ e^h y = \int{e^h r dx} + c $$
가 된다. 정리하면,
$$ y(x) = e^{-h} (\int{e^h r dx} + c) , h = \int {p(x)dx} $$
책에서는 전체 출력을 초기 데이터와 입력 값에 대한 반응의 합으로 표현된다. (명심하자)
"Total output = Response to the Input r + Response to the Initial Data"
2. Nonlinear ODE
ODE의 y 차수가 2 이상인 경우 비선형이라고 하고 다음과 같은 형태로 나타낼 수 있다.
$$ y' + p(x)y = g(x) y^a $$
a = 0, a = 1 일 때는 1차 선형 미분 방정식이니 패스.
나머지의 경우에 대해서 풀어보자. 비선형을 선형방정식 형태로 바꾸는 것은 생각보다 간단하다. 다음과 같이 치환하고, x에 대해 미분하면 다음과 같다.
$$ u(x) = [y(x)]^{1-a} $$
$$ u' = (1-a)y^{-a} y' = (1 - a) y^{-a}(gy^a - py) $$
그리고, 식을 정리하면
$$ u' + (1 - a)pu = (1 - a) g $$
가 된다. 이제 Linear ODE형태가 되었으니 문제를 풀 수 있다!
Q. $ y^b $ 가 하나 더 생기면 어떻게 풀까?
Advanced Engineering Math. Problem set 1.5 : 2 ~ 4
Problem SET 1.5 : 2 ~ 4 Problem 2. Integration Constant Give a reason why you may choose the constant of integration in $ \int p dx $ to be zero. Nonhomogeneous 적분인자 h 를 구할 때 상수 C를 무시한 이유에 대해 설명하라는 문제이
aeromaster.tistory.com
Advanced Engineering Math. Problem set 1.5 : 5 ~ 7
Problem 5. $$ y' + ky = e^{-kx} $$ Let's say, $ p \equiv k , r(x) \equiv e^{-kx} $ $ h = \int{k dx} = - kx $ 이다. Substitute equation h in $ y(x) = e^{-h} (\int{e^h r dx} + C) $ $$ y = e^{ -kx } ( \int{e^{kx} e^{-kx} dx} + C ) $$ $$ y(x) = e^{-kx}(x + C
aeromaster.tistory.com
Advanced Engineering Math. Problem set 1.5 : 8 ~ 10
하루에 세 문제씩 풀고 있다. 이제는 적분인자 구하는 과정 및 최종 결과식이 자연스레 외워진다. Review : For Nonhomogeneous ODE, for instance, $ y' + p(x) y = r(x) $ $$ y(x) = e^{-h} (\int{e^h r(x) dx} + C) \ , \ h = \int{
aeromaster.tistory.com
'Engineering > Advanced Engineering Math.' 카테고리의 다른 글
Ch. 2 2nd order ODEs : 2.2 Homogeneous Linear ODEs with Constant Coefficients (1) | 2023.11.13 |
---|---|
Ch. 2 2nd order ODEs : 2.1 Homogeneous Linear ODEs for 2nd order (0) | 2023.11.10 |
Ch. 1 1st order ODEs : 1.4 Exact ODEs. Integrating Factors (2) | 2023.10.23 |
Ch. 1 1st order ODEs : 1.3 Separable ODEs. Modeling (0) | 2023.10.11 |
Ch. 1 1st order ODEs : Basic Concepts. Modeling (0) | 2023.10.11 |