베르누이 미분방정식 계산기

베르누이 미분방정식 계산기 정보

베르누이 미분방정식 계산기는 가장 유명한 1계 비선형 미분방정식 중 하나인 베르누이 방정식 y' + P(x)y = Q(x)yⁿ을 다룹니다. 교과서에 나오는 전형적인 유도 과정을 대화형 단계별 설명으로 변환해 줍니다. v = y¹⁻ⁿ 치환을 통해 방정식을 선형화하고, 적분 인자 μ(x)를 생성하며, 그 결과로 나온 폐쇄형 곡선을 RK4 수치해 및 방향장 위에 겹쳐서 모든 세부 사항을 한눈에 볼 수 있게 합니다.

베르누이 미분방정식이란 무엇인가요?

1695년 야코프 베르누이가 소개한 베르누이 방정식은 다음과 같은 형태의 1계 상미분 방정식입니다.

n = 0일 때 방정식은 이미 선형이며, n = 1일 때는 변수분리형입니다. 그 외의 모든 실수 n에 대해 방정식은 비선형이지만, 전형적인 v = y¹⁻ⁿ 치환을 통해 v에 대한 선형 상미분 방정식으로 변환할 수 있으며, 이는 표준적인 적분 인자 기법으로 풀 수 있습니다.

6단계 베르누이 방법

y' + P(x)y = Q(x)yⁿ에서 시작합니다:

1. yⁿ으로 나누기: \( y^{-n}y' + P(x)\,y^{1-n} = Q(x) \).
2. v = y¹⁻ⁿ 치환: \( v' = (1-n)y^{-n}y' \)이므로, \( y^{-n}y' = v'/(1-n) \)이 됩니다.
3. 선형화: \( v' + (1-n)P(x)\,v = (1-n)Q(x) \) — v에 대한 1계 선형 상미분 방정식입니다.
4. 적분 인자: \( \mu(x) = \exp\!\left(\int (1-n)P(x)\,dx\right) \)이므로, \( (\mu v)' = \mu(1-n)Q(x) \)가 됩니다.
5. v(x)에 대해 풀기: \( v(x) = \frac{1}{\mu(x)}\left[\mu(x_0)v_0 + \int_{x_0}^{x}\mu(t)(1-n)Q(t)\,dt\right] \).
6. 역치환: \( y(x) = v(x)^{1/(1-n)} \).

관련 적분이 초등 함수인 경우 깔끔한 폐쇄형 솔루션을 얻을 수 있습니다. 그렇지 않은 경우, 계산기는 심슨 공식을 사용하여 수치적으로 계산하여 솔루션 곡선을 그립니다.

자동으로 처리되는 특수 사례

풀이 예제 — n = 2, 로지스틱 스타일

초기 조건 y(1) = 1인 y' + y/x = x·y² 방정식을 고려해 봅시다. 여기서 P(x) = 1/x, Q(x) = x, n = 2이므로 1 − n = −1입니다.

1. v = y⁻¹ = 1/y로 치환합니다. 그러면 v' = −y⁻²y'이 되고 방정식은 v' − (1/x)v = −x가 됩니다.
2. 적분 인자: μ(x) = exp(∫−1/x dx) = 1/x.
3. (μ·v)' = μ·(−x) = −1. 적분하면: (1/x)·v = −x + C, 즉 v = −x² + Cx가 됩니다.
4. 초기 조건 적용: x = 1일 때 v = 1/1 = 1이므로, 1 = −1 + C ⇒ C = 2입니다. 따라서 v(x) = −x² + 2x입니다.
5. 역치환: y(x) = 1/v(x) = 1/(2x − x²) = 1/(x(2 − x))입니다.

폐쇄형 솔루션 y = 1/(x(2−x))는 x = 0과 x = 2에서 수직 점근선을 가집니다. 이는 방향장을 통해 한눈에 명확하게 확인할 수 있습니다.

이 계산기 사용 방법

1. 방정식 빌더를 채웁니다. 파란색 슬롯에 P(x)와 Q(x)를 입력하고, 작은 위첨자 상자에 지수 n을 입력합니다. 레이아웃은 표준 형태인 y' + P(x)y = Q(x)yⁿ을 따릅니다.
2. 초기 조건을 설정합니다. (x₀, y₀)와 그래프 범위 [x 최소값, x 최대값]을 입력합니다. 범위는 x₀를 포함해야 합니다.
3. 계산하기를 클릭합니다. 계산기는 특수 사례(n = 0 또는 n = 1)인지 감지하여 그에 맞는 유도 과정을 보여줍니다. 그렇지 않은 경우 MathJax로 렌더링된 수식과 함께 전체 6단계 베르누이 치환 과정을 실행합니다.
4. 그래프를 읽습니다. 주황색 곡선은 RK4 수치해입니다. 파란색 점선은 적분 인자를 통해 계산된 폐쇄형 솔루션입니다. 화살표 장은 모든 곳에서의 y'을 보여주므로 다른 솔루션의 흐름도 가늠해 볼 수 있습니다.
5. 샘플 포인트의 CSV를 복사합니다. 궤적 데이터를 다른 프로그램으로 가져오고 싶은 경우 사용하세요.

팁, 주의사항 및 예외 사례

• 정수가 아닌 n은 y > 0을 요구합니다. 솔루버는 n = 1/2이면서 y₀ ≤ 0인 경우와 같이 yⁿ이 복소수가 되는 조합을 감지하여 알립니다.
• y₀ = 0은 종종 특이점입니다. Q ≠ 0이고 n > 0인 모든 베르누이 방정식은 자명한 해 y ≡ 0을 갖지만, 이는 일반적으로 사용자가 원하는 분기가 아닙니다.
• x₀ 근처에서 P(x)의 발산을 피하세요. 1/x와 같은 식은 x₀ ≠ 0이어야 합니다. 솔루버는 실행 전 이를 검증합니다.
• 큰 지수 (|n| > 20)는 오버플로 방지를 위해 거부됩니다. 실제로 n이 이 정도로 큰 베르누이 방정식은 실무에서 거의 나타나지 않습니다.
• 수직 점근선. RK4가 발산하면 솔루션이 유한하게 유지되는 x₀ 쪽으로 x 범위를 좁혀 보세요.

베르누이 방정식이 나타나는 분야

• 인구 동태학 — 로지스틱 방정식 y' = ry(1 − y/K)는 정리하면 n = 2인 베르누이 방정식입니다.
• 화학 반응 속도론 — 자기 촉매 반응은 종종 y' ∝ y − y²을 따릅니다.
• 전기 회로 — 특정 비선형 저항을 가진 RL 회로는 베르누이 형태를 띱니다.
• 유체 역학 — 유사성 감소 후의 경계층 방정식.
• 전염병 모델 — SIR 모델의 취약계층 분율은 베르누이 형태로 환원될 수 있습니다.
• 경제 성장 — 저축률이 일정한 솔로우-스완 모델은 n = α인 베르누이 형태입니다.

자주 묻는 질문 (FAQ)

베르누이 미분방정식이란 무엇인가요?

베르누이 방정식은 y' + P(x)y = Q(x)yⁿ 형태의 1계 상미분 방정식으로, 여기서 P와 Q는 연속 함수이고 n은 임의의 실수입니다. v = y¹⁻ⁿ 치환을 통해 선형 방정식으로 변환할 수 있는 비선형 상미분 방정식의 전형적인 예입니다.

v = y¹⁻ⁿ 치환은 어떻게 작동하나요?

원래 방정식에 y⁻ⁿ을 곱하여 모든 y 항이 y¹⁻ⁿ 또는 y⁻ⁿy'이 되도록 합니다. v = y¹⁻ⁿ으로 설정하면 v' = (1−n)y⁻ⁿy'이 됩니다. 이를 대입하면 베르누이 방정식은 v' + (1−n)P(x)v = (1−n)Q(x)로 변환되며, 이는 v에 대한 선형 방정식이 되어 적분 인자로 풀 수 있습니다.

n = 0 또는 n = 1일 때 어떻게 되나요?

n = 0이면 방정식은 이미 1계 선형 방정식이므로 치환이 필요하지 않습니다. n = 1이면 베르누이 공식에서 1 − n = 0으로 나누게 되므로 이를 별도로 처리합니다. 방정식은 y' = (Q(x) − P(x))·y로 축소되며, 이는 y = y₀·exp(∫(Q−P) dx) 형태의 폐쇄형 솔루션을 갖는 변수분리형 방정식입니다.

베르누이 방정식은 항상 폐쇄형으로 풀 수 있나요?

원칙적으로는 그렇지만, 적분 인자를 포함하는 결과 적분이 초등 함수로 표현되지 않을 수 있습니다. 이 경우 계산기는 심슨 공식을 사용하여 수치적으로 계산하고 솔루션 곡선을 그립니다. 방법 자체는 항상 베르누이 상미분 방정식을 구적법으로 환원합니다.

음수 y와 정수가 아닌 n이 문제가 되는 이유는 무엇인가요?

n이 정수가 아니면 yⁿ은 exp(n·ln y)로 정의되며 y > 0일 때만 실수입니다. 음수 y를 입력하면 복소수가 생성됩니다. 솔루버는 이 상황을 감지하여 솔루션이 실수 값을 유지할 수 있도록 y₀ > 0 또는 정수 지수를 요구합니다.

방향장은 무엇을 보여주나요?

방향장은 해당 (x, y) 점에서의 y'과 각도가 같은 아주 작은 접선 조각들의 격자입니다. 모든 솔루션 곡선은 이 접선들을 따라가야 하므로, 방향장을 통해 모든 솔루션의 질적인 모양을 한눈에 볼 수 있으며 초기 조건에 따라 특정 곡선이 결정됩니다.

더 읽어보기

• 베르누이 미분 방정식 — 위키백과
• 적분 인자 — 위키백과
• 로지스틱 함수 — 위키백과
• Slope field — Wikipedia (영문)