2계 상미분방정식 해결사

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

2계 상미분방정식 해결사는 상수 실수 계수를 갖는 a·y″ + b·y′ + c·y = g(x) 형태의 선형 상미분방정식을 입력받아 특성 방정식을 자동으로 유도하고, 댐핑 체계(과감쇠, 임계 감쇠, 미달 감쇠, 무감쇠 또는 불안정)를 분류하며, 기호적 폐쇄형 해와 고정밀 수치 해를 모두 생성합니다. 대화형 출력은 y(x) 및 y′(x)의 이중 곡선 시간 그래프와 (y, y′)의 위상 평면 궤적을 함께 제공합니다. 이 뷰를 통해 미달 감쇠의 경우 안쪽으로의 나선형, 과감쇠의 경우 마디, 무감쇠의 경우 폐쇄 루프, 불안정의 경우 바깥쪽으로의 나선형 등 체계를 한눈에 파악할 수 있습니다.

상수 계수를 갖는 2계 선형 ODE란 무엇입니까?

실수 상수 계수를 갖는 2계 선형 상미분방정식(Ordinary Differential Equation)은 다음과 같은 형태의 방정식입니다.

여기서 a ≠ 0이고 b, c는 실수 상수이며 g(x)는 강제항입니다. 두 개의 초기 조건 y(x₀) = y₀ 및 y′(x₀) = y′₀를 추가하면 x₀ 근처에서 유일한 해를 갖는 초기값 문제(Initial-Value Problem)가 됩니다. 이는 동등한 1계 시스템에 적용된 피카르-린델뢰프 정리에 의해 보장됩니다.

만약 g(x) = 0이면 이 방정식은 제차(homogeneous) 방정식입니다. 그렇지 않으면 비제차(non-homogeneous)이며, 전체 해는 다음과 같이 분해됩니다.

여기서 y_h는 연관된 제차 방정식의 일반해(두 개의 임의 상수 포함)이고, y_p는 전체 방정식의 임의의 특수해입니다. 두 초기 조건을 적용하면 두 임의 상수가 고정됩니다.

특성 방정식

제차 방정식에서 y = e^(r·x)라고 가정하면 특성(또는 보조) 방정식을 얻게 됩니다.

판별식 Δ = b² − 4ac를 갖는 이 2차 방정식이 전체적인 정성적 거동을 제어합니다.

세 가지 근의 경우 및 댐핑 체계

미정계수법 (비제차의 경우)

g(x)가 다음과 같은 단순한 형태를 가질 때, 미정계수법은 미지의 계수를 포함한 동일한 형태의 시도해를 가정하고 이를 풀어 특수해를 제공합니다.

• 상수 g(x) = k. 시도해: y_p = K. 만약 c = 0이면 x를 곱하고, b = 0까지 만족하면 x²을 곱합니다.
• n차 다항식. 시도해: 일반적인 n차 다항식. 상수항이나 선형항이 공진하는 경우 x 또는 x²을 곱합니다.
• 지수함수 g(x) = A·e^(k·x). 시도해: y_p = K·e^(k·x). k가 특성근과 일치하면 x(단일근) 또는 x²(중근)을 곱하며, 이를 공진(resonance)이라고 합니다.
• 정현파 g(x) = A·cos(ω·x) + B·sin(ω·x). 시도해: y_p = K₁·cos(ω·x) + K₂·sin(ω·x). 만약 iω가 근이면 x를 곱합니다(순수 주파수 공진).
• 곱과 합은 선형성과 곱셈 법칙에 따릅니다.

위상 평면 해석하기

동등한 1계 시스템은 u = y, v = y′이며 u′ = v 및 v′ = (g(x) − b·v − c·u)/a입니다. x에 따라 매개변수적으로 u에 대한 v를 플롯하면 위상 평면 궤적이 됩니다. 제차 자율 시스템(g에 x가 없음)의 경우, 궤적은 시작점(y₀, y′₀)에 의해 유일하게 결정되며 체계를 한눈에 보여줍니다.

• 미달 감쇠: 궤적이 원점을 향해 안쪽으로 나선을 그리며 수렴합니다.
• 과감쇠: 궤적이 불변선(느린 고유벡터)을 따라 원점으로 접근합니다.
• 임계 감쇠: 퇴화된 마디 형태로, 궤적이 단일 고유벡터에 접하며 수렴합니다.
• 무감쇠: 원점을 둘러싼 닫힌 타원 — 영구적인 진동을 의미합니다.
• 불안정: 궤적이 나선을 그리거나 바깥쪽 무한대로 발산합니다.

풀이 예제: 강제 감쇠 조화 진동자

y(0) = 0, y′(0) = 0인 방정식 y″ + 2·y′ + 5·y = 10을 고려해 봅시다. 이는 강제력이 있는 미달 감쇠 시스템입니다.

1. 특성 방정식: r² + 2r + 5 = 0 → Δ = 4 − 20 = −16 → r = −1 ± 2i.
2. 제차해: y_h = e^(−x)·(C₁·cos 2x + C₂·sin 2x).
3. 특수해: 상수 강제력 g = 10에 대해 y_p = K를 대입하면 5K = 10이므로 y_p = 2입니다.
4. 초기 조건 적용: y(0) = 0 → C₁ + 2 = 0 → C₁ = −2. y′(0) = 0 → −C₁ + 2C₂ = 0 → C₂ = −1.
5. 최종 답안: y(x) = 2 − e^(−x)·(2·cos 2x + sin 2x) — 감쇠 봉투선과 함께 진동하며 y → 2로 수렴합니다.

이 계산기 사용 방법

1. 계수 입력: 첫 번째 줄에 a, b, c를 입력합니다. a는 0이 아니어야 합니다(그렇지 않으면 1계 방정식이 됩니다).
2. 강제항 g(x) 입력: 제차 문제의 경우 0으로 두거나 비워둡니다. 상수, 2차 이하의 다항식, 공진을 포함한 단일 지수 A·e^(k·x)에 대해 폐쇄형 특수해를 유도합니다.
3. 초기 조건 제공: (x₀, y₀, y′₀)를 입력합니다. 2계 방정식이므로 x₀에서의 y와 y′이 모두 지정되어야 합니다.
4. x 범위 선택: 그래프를 그릴 범위를 정합니다. 해석기는 RK4를 사용하여 x₀에서 양방향으로 통합을 수행합니다.
5. 풀이 및 시각화 클릭: 복소 평면상의 근과 함께 특성 방정식, 댐핑 체계 분류, 제차 및 특수 폐쇄형 해, y와 y′의 이중 곡선 시간 그래프 및 위상 평면 궤적을 얻을 수 있습니다.

일반적인 응용 분야

• 기계적 스프링-질량-댐퍼 시스템: m·x″ + c·x′ + k·x = F(t). 과감쇠, 임계 감쇠, 미달 감쇠는 감쇠비 ζ = c/(2·√(m·k))에 대응됩니다.
• RLC 전기 회로: 직렬 RLC 회로는 L·Q″ + R·Q′ + Q/C = V(t)를 따르며, 이는 기계 시스템과 동일한 구조를 가집니다.
• 진자 (작은 각도): θ″ + (g/L)·θ = 0은 단순 조화 운동을 나타내며, 공기 저항을 추가하면 감쇠 진동이 됩니다.
• 지진에 대한 건물 반응: 지반 가속도를 강제항으로 하는 단자유도 구조물 해석.
• PID 제어 서보 시스템: 폐쇄 루프 오차 역학은 댐핑 비율에 따라 오버슈트가 결정되는 2계 ODE로 축소됩니다.
• 관성이 있는 인구 모델: 자본 축적 시차가 있는 경제 성장 또는 지연 반응이 있는 생태계 모델.

수치적 방법 — 2D 시스템에 대한 고전적 룬게-쿠타(RK4)

이 도구는 a·y″ + b·y′ + c·y = g(x)를 다음과 같은 1계 시스템으로 축소합니다.

여기서 u(x₀) = y₀, v(x₀) = y′₀입니다. 그런 다음 벡터 상태 (u, v)에 4단계 룬게-쿠타를 적용합니다. RK4는 O(h⁵)의 국부 절단 오차와 O(h⁴)의 전역 오차를 가지며, 각 방향으로 기본 400개 하위 단계는 비강성 문제에 대해 약 6자리의 정확도를 제공합니다.

자주 묻는 질문

상수 계수를 갖는 2계 선형 ODE란 무엇입니까?

상수 계수를 갖는 2계 선형 ODE는 a·y″ + b·y′ + c·y = g(x) 형태를 가지며, 여기서 a, b, c는 실수 상수이고 g(x)는 강제(비제차) 항입니다. 두 개의 초기 조건 y(x₀) = y₀ 및 y′(x₀) = y′₀가 주어지면 해는 유일합니다. 제차형인 g(x) = 0인 경우에는 항상 특성 방정식 a·r² + b·r + c = 0을 통해 폐쇄형 해를 구할 수 있으며, 비제차형의 경우 y(x) = y_h(x) + y_p(x)로 풀이됩니다.

특성 방정식이란 무엇입니까?

a·y″ + b·y′ + c·y = 0에 대해 y = e^(r·x)라는 가설을 대입하면 특성 방정식 또는 보조 방정식인 a·r² + b·r + c = 0을 얻게 됩니다. 이 방정식의 근이 제차해의 형태를 결정합니다. 서로 다른 두 실근은 y_h = C₁·e^(r₁·x) + C₂·e^(r₂·x); 중근 r은 y_h = (C₁ + C₂·x)·e^(r·x); 켤레 복소수근 α ± β·i는 y_h = e^(α·x)·(C₁·cos(β·x) + C₂·sin(β·x)) 형태를 가집니다.

미달 감쇠, 임계 감쇠, 과감쇠는 무엇을 의미합니까?

이 용어들은 스프링-질량-댐퍼 모델 m·x″ + c·x′ + k·x = 0에서 유래되었습니다. 과감쇠(판별식 > 0, 두 실근)는 시스템이 진동 없이 천천히 평형 상태로 돌아감을 의미합니다. 임계 감쇠(판별식 = 0, 중근)는 오버슈트 없이 가장 빠르게 복귀하는 상태입니다. 미달 감쇠(판별식 < 0, 복소수근)는 감쇠 진동을 일으킵니다. 무감쇠(b = 0, c/a > 0)는 영원히 순수 정현파 진동을 유지합니다.

미정계수법이란 무엇입니까?

상수, 다항식, 지수, 사인, 코사인 및 그 곱과 같은 단순 강제력 g(x)의 경우, 특수해 y_p는 미지의 계수를 가진 g와 동일한 형태를 갖는다고 가정하고, 이를 ODE에 대입하여 계수를 결정합니다. g(x)가 특성근과 공진하는 경우 시도해에 x(또는 중근의 경우 x²)를 곱해야 합니다.

위상 평면이란 무엇입니까?

2D 시스템(y, y′)으로 축소된 2계 방정식의 경우, 위상 평면은 x가 진행됨에 따라 y에 대한 y′을 플롯합니다. 위상 평면의 솔루션 곡선은 체계를 한눈에 보여줍니다. 미달 감쇠의 경우 감쇠 나선형, 과감쇠의 경우 안쪽 마디, 무감쇠 조화 운동의 경우 닫힌 타원, 불안정 진동의 경우 바깥쪽 나선형을 나타냅니다. 이는 특성 방정식 근 선도의 기하학적 대응물입니다.

이 도구는 어떤 수치적 방법을 사용합니까?

고전적인 4계 룬게-쿠타(RK4) 방법이 u = y, v = y′인 동등한 1계 시스템에 적용됩니다(u′ = v 및 v′ = (g(x) − b·v − c·u)/a). RK4는 O(h⁵)의 국부 절단 오차를 가지며, 방향당 기본 400개 하위 단계는 선택한 창에서 비강성 방정식에 대해 대략 6자리의 정확도를 제공합니다.

더 읽어보기

• 선형 미분 방정식 — Wikipedia
• 특성 방정식 — Wikipedia
• 미정계수법 — Wikipedia
• 조화 진동자 — Wikipedia
• 위상 평면 — Wikipedia
• 룬게-쿠타 방법 — Wikipedia