1계 상미분방정식 해결사

1계 상미분방정식을 기호 및 수치적으로 해결합니다. 변수분리형, 선형, 완전 및 자율 형태를 자동으로 감지하여 적절한 기법을 적용하고, 해 곡선이 중첩된 대화형 방향장을 렌더링합니다.

1계 상미분방정식 해결사

1계 상미분방정식 해결사 정보

1계 상미분방정식 해결사는 dy/dx = f(x, y) 형태의 상미분방정식을 입력받아 구조(변수 분리형, 선형, 자율형, 완전형 또는 일반형)를 자동으로 분류하고, 가능한 경우 기호적 폐쇄형 해와 모든 경우에 대한 높은 정확도의 수치 해를 모두 생성합니다. 해 곡선이 중첩된 실시간 방향장 시각화는 방정식의 기하학적 의미를 즉각적으로 명확하게 해줍니다. 해는 모든 화살표에 접하는 곡선입니다.

1계 상미분방정식이란 무엇인가요?

1계 상미분방정식은 알 수 없는 함수 y(x)와 그 1계 도함수 y'(x)만을 포함합니다. 표준 명시적 형태는 다음과 같습니다:

dy/dx = f(x, y)

초기 조건 y(x₀) = y₀와 결합하면, 이는 초기값 문제(IVP)를 정의합니다. Picard-Lindelöf 정리는 f가 (x₀, y₀) 근처에서 y에 대해 립시츠 연속인 한 x₀의 일정 이웃 내에서 유일한 해를 보장합니다. 기하학적으로 IVP는 (x₀, y₀)를 통과하며 모든 지점의 기울기가 해당 지점의 f와 일치하는 유일한 곡선, 즉 방향장에 접하는 곡선을 찾는 것입니다.

해결사가 인식하는 6가지 클래스

클래스	형태	표준 풀이 기법	이 도구가 하는 일
순수 적분	dy/dx = f(x)	직접 적분: y = ∫f(x) dx + C	수치 적분 (RK4가 심슨 공식과 유사한 구적법으로 축소됨)
선형 (상수 계수)	dy/dx = a·y + b	적분 인자 또는 특성근을 통한 폐쇄형	전체 기호 해답 + 단계별 유도 과정
자율형	dy/dx = f(y)	변수 분리: ∫dy/f(y) = x + C	수치 해 + 방향장 시각화
변수 분리형	dy/dx = g(x)·h(y)	변수 분리: ∫dy/h(y) = ∫g(x) dx + C	교차 곱 테스트를 통한 형태 감지; 수치 해 표시
선형 (변수 계수)	dy/dx + P(x)·y = Q(x)	적분 인자 μ(x) = e^∫P(x) dx	유한 차분 선형성 테스트를 통한 형태 감지; 수치 해 표시
일반형	기타 모든 dy/dx = f(x, y)	수치 해석법 (RK4, RK45, BDF, …)	600개의 하위 단계를 갖는 고전적 Runge-Kutta

폐쇄형 메서드: 상수 계수를 갖는 선형 방정식

우변이 상수 a와 b를 사용하여 dy/dx = a·y + b로 단순화될 때, 적분 인자 μ(x) = e^(-a·x)가 정확한 해를 제공합니다. 일반해는 다음과 같습니다:

y(x) = -b/a + C · e^(a·x) (a ≠ 0) y(x) = b·x + C (a = 0)

초기 조건 y(x₀) = y₀를 적용하면 상수 C가 결정되어 고유한 특수해를 얻습니다. 이 단일 클래스는 수많은 교과서 문제를 다룹니다:

지수 성장 — dy/dx = k·y, 특수해 y(t) = y₀·e^(k·t).
지수 감쇠 — dy/dx = -k·y, 반감기 ln 2 / k.
뉴턴의 냉각 법칙 — dy/dx = -k·(y - T_amb), 체온이 주변 온도를 향해 지수적으로 완화됩니다.
RC 회로 충전 — dV/dt = (1/RC)·(V_in - V), 커패시터 전압이 전원 전압에 접근합니다.
약물 제거 — 제거율 k를 갖는 1계 약동학.

방향장 읽는 방법

모든 격자점 (x, y)에서 도구는 기울기가 f(x, y)와 같은 짧은 선분을 그립니다. 세 가지 유용한 관찰 사항:

평형점은 f(x, y) = 0인 지점이며 방향장이 수평입니다. 자율 방정식의 경우, 이는 f(y*) = 0을 만족하는 고정점 y*입니다. 근처의 궤적은 y*에 접근하거나(안정) 멀어집니다(불안정).
등경선(Isoclines)은 f(x, y)가 상수 c와 같은 곡선으로, 곡선을 따라 모든 화살표가 동일한 기울기 c를 가집니다.
해 곡선은 절대 교차하지 않음 (f가 립시츠인 경우) — 두 교차 곡선은 교차점에서 서로 다른 기울기를 가져야 하므로 시각적으로 분명합니다.

수치 해석법: 고전적 Runge-Kutta (RK4)

(x_n, y_n)이 주어지면 다음 값은 네 개의 기울기 추정치를 평균하여 계산됩니다:

k₁ = f(x_n, y_n) k₂ = f(x_n + h/2, y_n + h·k₁/2) k₃ = f(x_n + h/2, y_n + h·k₂/2) k₄ = f(x_n + h, y_n + h·k₃) y_{n+1} = y_n + (h/6) · (k₁ + 2k₂ + 2k₃ + k₄)

RK4는 국소 절단 오차 O(h⁵)와 전역 오차 O(h⁴)를 가지며, 비강성 방정식에 대해 기본 단계 수에서 약 6자리의 정확도를 제공합니다. 해석기는 초기 지점에서 두 x 방향으로 바깥쪽으로 적분하며, y의 크기가 10¹⁵을 초과하면 (dy/dx = y²과 같이 유한한 시간에 폭발하는 해의 전형적인 경우) 깔끔하게 중단됩니다.

이 계산기 사용 방법

우변 입력: dy/dx = ... 필드에 내용을 입력하세요. 변수로 x와 y를 사용하고, 곱셈은 *, 거듭제곱은 ^ 또는 **, 그리고 sin, cos, exp, log, sqrt와 같은 표준 함수를 사용합니다. 상수 pi와 e를 인식합니다.
초기 조건 지정 (x₀, y₀) — 유일한 해 곡선이 이 지점을 통과하게 됩니다.
x 범위 선택: 방향장과 해 곡선을 그릴 범위를 정하세요. y 범위는 적분된 해에서 자동으로 맞춰집니다.
풀이 및 시각화 클릭: 분류기가 먼저 실행됩니다. 방정식이 폐쇄형 패턴(상수 계수 선형)과 일치하면 기호 답안을 얻습니다. 방향장과 해 곡선은 항상 렌더링됩니다.
방향장 토글: 해 곡선에 집중하기 위해 방향장을 켜거나 끄고, 초기 지점에서 적분이 어떻게 진행되는지 보려면 곡선 그리기 애니메이션을 다시 재생하세요.

풀이 예시: 뉴턴의 냉각 법칙

80 °C의 커피 한 잔이 20 °C의 방에서 식습니다. 열 전달률은 온도 차이에 비례합니다:

dT/dt = -0.1 · (T - 20), T(0) = 80

이것은 상수 계수를 갖는 선형 방정식(a = -0.1, b = 2)입니다. 폐쇄형은 다음과 같습니다:

T(t) = 20 + 60 · e^(-0.1 t)

30분 후: T(30) = 20 + 60·e⁻³ ≈ 22.99 °C. 방향장 뷰는 제한 동작을 명확하게 보여줍니다. 시작 온도에 관계없이 모든 해 곡선은 수평선 T = 20으로 점근합니다.

일반적인 응용 분야

인구 동태 — 지수, 로지스틱, 알리 효과 모델.
약동학 — 약물 흡수 및 제거, 반감기 계산.
열 전달 — 뉴턴의 냉각 법칙, 집중 용량 모델.
RC 및 RL 회로 — 1계 선형 전기 과도 현상.
방사성 붕괴 — 단일 동위원소 붕괴 사슬.
혼합 탱크 — 유입/유출 하에서 용질의 농도.
항력을 받는 낙하 물체 — 종단 속도 분석 dv/dt = g - kv.

자주 묻는 질문

1계 상미분방정식이란 무엇인가요?

1계 상미분방정식(ODE)은 알 수 없는 함수 y(x)와 그 1계 도함수를 포함하는 dy/dx = f(x, y) 형태의 방정식입니다. ODE를 푼다는 것은 도함수가 우변과 일치하는 함수 y(x)를 찾는 것을 의미합니다. 초기 조건 y(x₀) = y₀가 주어지면, 가벼운 규칙성 가정(Picard-Lindelöf 정리) 하에 해는 유일합니다.

방향장(Slope field)이란 무엇인가요?

방향장(또는 기울기장)은 각 격자점 (x, y)에서 기울기가 f(x, y)와 같은 작은 선분을 그리는 그래프입니다. ODE의 해 곡선은 모든 지점에서 이러한 선분에 접하는 곡선입니다. 방향장은 방정식을 기호적으로 풀지 않고도 해의 전역적 동작에 대한 즉각적인 시각적 직관을 제공합니다.

이 도구는 어떤 종류의 1계 ODE를 풀 수 있나요?

이 도구는 방정식을 다음과 같이 자동 분류합니다: 적분 가능(x에만 의존, 직접 적분으로 풀이), 상수 계수 선형 y' = a·y + b(전체 폐쇄형 제공), 자율형(y에만 의존), 변수 분리형(g(x)·h(y)로 인수분해), 변수 계수 선형(P(x)·y + Q(x)) 또는 일반형. 모든 클래스에 대해 높은 정확도의 Runge-Kutta 수치 해와 방향장 시각화가 생성됩니다.

어떤 수치 해석법이 사용되나요?

고전적인 4차 Runge-Kutta 방법(RK4)이 초기 지점에서 각 방향으로 300개의 하위 단계를 사용하여 적용됩니다. RK4는 국소 절단 오차 O(h⁵)를 가지며 이 규모의 비강성 ODE에 대한 표준적인 방법입니다. 해석기는 발산(오버플로 또는 NaN)을 감지하고 그래프가 유효하게 유지되도록 적분을 깔끔하게 중단합니다.

선형 ODE의 적분 인자법이란 무엇인가요?

선형 1계 ODE y' + P(x)·y = Q(x)의 경우, 양변에 적분 인자 μ(x) = e^∫P(x) dx를 곱합니다. 왼쪽은 완전 도함수 d/dx[μ·y]가 되므로, y(x) = (1/μ(x)) · (∫μ(x)·Q(x) dx + C)가 됩니다. P와 Q가 상수일 때 이는 도구가 자동으로 반환하는 폐쇄형 y = -b/a + C·e^(a·x)로 축소됩니다.

이 도구가 강성 방정식이나 ODE 시스템을 처리할 수 있나요?

이 해석기는 비강성 1계 스칼라 ODE용으로 설계되었습니다. 매우 강성인 문제(해의 시간 척도가 여러 자릿수 차이 나는 경우)는 후방 오일러 또는 Rosenbrock과 같은 음함수 방법이 필요할 수 있으며, 결합 시스템은 벡터값 해석기가 필요합니다. 그러한 경우에는 SciPy의 solve_ivp 또는 전문적인 강성 ODE 해석기와 같은 전용 패키지를 사용하십시오.