상미분 방정식 시스템 솔버

상미분 방정식 시스템 솔버 정보

상미분 방정식 시스템 솔버는 연립 선형 및 비선형 시스템을 위한 올인원 미분 방정식 도구 상자입니다. 2×2 또는 3×3 계수 행렬을 입력하면 도구가 전체 고윳값/고유벡터 분석을 수행하고, LaTeX 형식으로 폐쇄형 일반해 및 특수해를 작성하며, 원점의 평형점을 안장점, 마디, 나선점 또는 중심으로 분류하고, 애니메이션 궤적이 포함된 대화형 위상 초상화를 그립니다. 비선형 평면 시스템의 경우 임의의 우변 식 \(f(x,y)\) 및 \(g(x,y)\)를 입력하여 고정밀 RK4 위상 초상화를 생성할 수 있습니다.

상미분 방정식 시스템이란 무엇입니까?

상미분 방정식(ODE) 시스템은 여러 미지 함수의 도함수를 통해 하나의 독립 변수(보통 시간 \(t\))와 연결합니다. 가장 간결한 형태로는 다음과 같습니다.

상수 행렬 \(A\)에 대해 \(\mathbf{F}(t, \mathbf{x}) = A\mathbf{x}\)인 경우, 시스템은 선형 자율 시스템이 됩니다. 이 영역에서 이론은 가장 아름답게 펼쳐집니다. 전체적인 장기적 동작이 행렬 \(A\)의 고윳값에 의해 결정되기 때문입니다.

선형 시스템을 위한 고윳값 계산법

\(\mathbf{x}' = A\mathbf{x}\)에 대한 표준 방법은 다음과 같습니다.

1. 특성 다항식 \(\det(\lambda I - A) = 0\)을 계산합니다.
2. 고윳값 \(\lambda_1, \lambda_2, \dots\)을 구합니다.
3. 각 고윳값에 대해 \((A - \lambda I) v = 0\)을 풀어 고유벡터 \(v\)를 찾습니다.
4. 선형 결합으로 일반해를 구성합니다: \(\mathbf{x}(t) = c_1 v_1 e^{\lambda_1 t} + c_2 v_2 e^{\lambda_2 t} + \cdots\).
5. 초기 조건 \(\mathbf{x}(0)\)을 일반해에 대입하여 상수 \(c_i\)를 확정합니다.

2×2 시스템의 세 가지 사례

자취-행렬식 평면

자취 \(T = a_{11} + a_{22}\)와 행렬식 \(D = a_{11} a_{22} - a_{12} a_{21}\)을 갖는 2×2 행렬의 경우, 전체 분류는 하나의 다이어그램으로 요약됩니다.

이것이 결과 패널에서 \(T\), \(D\), \(\Delta = T^2 - 4D\)를 두드러지게 보여주는 이유입니다. 이 세 가지 숫자만으로 평형점의 이름을 정할 수 있습니다.

비선형 시스템과 위상 초상화

대부분의 실제 미분 방정식은 비선형이며 폐쇄형 해가 없습니다. 이 도구는 국부 절단 오차가 \(O(h^5)\)이며 매끄러운 벡터장에 대한 표준적인 방법인 4차 룬게-쿠타 (RK4) 방식으로 방정식을 수치 적분하여 이를 처리합니다.

위상 초상화 오버레이:

• 13×13 그리드에서 샘플링된 벡터장: 모든 지점에서의 흐름 방향을 보여줍니다.
• 초기 조건으로부터의 궤적: 빨간색으로 표시되며, 시간의 방향을 보여주는 주황색 애니메이션 러너가 포함됩니다.
• 여러 시작점 유선: 고리 모양의 시작점에서 뻗어나가는 선들을 통해 역학의 전역적인 그림을 제공합니다.
• 2×2 선형 시스템의 경우 고유벡터 축 (점선 청록색): 해가 지수적으로 미끄러지는 불변 방향입니다.

이 솔버 사용 방법

1. 폼 상단의 탭을 통해 2×2 선형, 3×3 선형 또는 비선형 2D 모드를 선택합니다.
2. 계수 또는 방정식을 입력합니다. 빠른 예제를 클릭하여 전형적인 시스템(안정된 마디, 중심, 안장점, 진자, 반 데르 폴 등)을 자동으로 채울 수 있습니다.
3. 초기 조건 \((x_0, y_0)\)과 시간 범위 \(T\)를 입력합니다. 일반적인 \(T\) 값은 진동자의 경우 6~20, 빠르게 감쇠하는 안정된 시스템의 경우 3~6입니다.
4. 계산 버튼을 클릭합니다. 분류, 고윳값, 고유벡터, 폐쇄형 해(선형 모드), 애니메이션 위상 초상화 및 시계열 플롯이 포함된 결과 페이지가 나타납니다.
5. 위상 초상화 아래의 버튼을 사용하여 궤적을 다시 재생하고 러너가 초기 조건 곡선을 횡단하는 모습을 다시 볼 수 있습니다.

작업 예제 — 감쇠 조화 진동자

감쇠 진동자 \(\ddot{x} + 2\zeta \omega \dot{x} + \omega^2 x = 0\)는 \(y = \dot{x}\)라고 둠으로써 2D 시스템으로 다시 쓸 수 있습니다.

\(\omega = 1\)이고 \(\zeta = 0.2\) (부족 감쇠)인 경우, 행렬은 \(A = \begin{pmatrix} 0 & 1 \\ -1 & -0.4 \end{pmatrix}\)입니다. 자취 \(T = -0.4\), 행렬식 \(D = 1\), 판별식 \(\Delta = 0.16 - 4 = -3.84 < 0\)이므로, 고윳값 \(-0.2 \pm 0.9798\,i\)를 갖는 안정된 나선점을 얻습니다. 궤적은 원점으로 소용돌이치며 들어가고, 시계열은 지수적으로 감쇠하는 사인파를 보여줍니다.

응용 분야

• 기계 시스템 — 결합된 용수철-질량 시스템, 진자, 자이로스코프.
• 전기 회로 — RLC 네트워크, 연산 증폭기 필터, 상태 공간 제어.
• 인구 역학 — 로트카-볼테라 포식자-피식자 모델, 경쟁 종, 역학 모델 (SIR, SIS).
• 화학 반응 속도론 — 반응 네트워크, 벨로우소프-자보틴스키 진동자.
• 신경과학 — 피츠휴-나구모 뉴런 모델, 호지킨-헉슬리 약축 모델.
• 제어 이론 — 선형화된 플랜트 모델, 관측기 설계, 안정성 여유.

팁 및 주의사항

• 궤적이 너무 빨리 발산한다면 시간 범위 T를 줄이세요. 불안정한 시스템은 단 몇 단위 시간 만에 뷰포트를 벗어날 수 있습니다.
• 중근 고윳값의 경우, 솔버가 \((A - \lambda I)w = v\)를 풀어 일반화된 고유벡터 \(w\)를 자동으로 찾으므로 수동 작업 없이 \(tv\) 항을 얻을 수 있습니다.
• 비선형 시스템의 경우 벡터장 화살표가 원점이 아닌 평형점도 청록색 점으로 보여줍니다. 화살표 크기가 0인 영역을 주시하세요.
• 3×3 시스템의 경우 위상 초상화는 제공되지 않지만(3D는 2D 페이지에 표시하기 어렵습니다), 폐쇄형 해와 안정성 판단은 여전히 적용됩니다.
• 초기 조건과 시간 범위는 분류와 무관합니다. 이를 변경해도 빨간색 궤적만 이동할 뿐 고윳값 판정 결과는 바뀌지 않습니다.

자주 묻는 질문

상미분 방정식 시스템이란 무엇입니까?

상미분 방정식(ODE) 시스템은 일반적으로 시간인 하나의 독립 변수에 대한 여러 미지 함수의 도함수를 연결하는 연립 방정식 세트입니다. 클래식한 형태는 \( \mathbf{x}'(t) = F(t, \mathbf{x}(t)) \)이며, 여기서 \( \mathbf{x} \)는 상태 벡터이고 \(F\)는 벡터장입니다. 선형 시스템은 \( \mathbf{x}' = A\mathbf{x} + \mathbf{b} \)로 간결하게 쓸 수 있으며, 그 동작은 거의 전적으로 계수 행렬 \(A\)의 고윳값에 의해 결정됩니다.

고윳값은 2×2 선형 시스템의 평형점을 어떻게 분류합니까?

2×2 시스템 \( \mathbf{x}' = A\mathbf{x} \)에서 원점은 \(A\)의 자취 \(T\)와 행렬식 \(D\)에 의해 분류됩니다. \(D < 0\)이면 안장점(불안정)입니다. \(D > 0\)이면서 \(T^2 > 4D\)이면 마디(\(T < 0\)이면 안정, \(T > 0\)이면 불안정)입니다. \(D > 0\)이면서 \(T^2 < 4D\)이면 나선점(\(T < 0\)이면 안정, \(T > 0\)이면 불안정, \(T = 0\)이면 순수 중심)입니다. 경계인 \(T^2 = 4D\)는 퇴화된 마디를 생성합니다.

고윳값이 복소수일 때 폐쇄형 해는 어떤 모양입니까?

\(A\)가 복소 고유벡터 \( v = p + iq \)와 함께 복소 공액 고윳값 \( \alpha \pm i\beta \)를 갖는 경우, 실수 일반해는 \( \mathbf{x}(t) = e^{\alpha t} \left[ c_1 (p \cos\beta t - q \sin\beta t) + c_2 (p \sin\beta t + q \cos\beta t) \right] \)입니다. 지수 \(e^{\alpha t}\)는 진폭(성장, 감쇠 또는 일정)을 제어하고 사인과 코사인은 회전을 담당합니다.

행렬에 중근 고윳값이 있으면 어떻게 됩니까?

행렬에 중근 고윳값 \(\lambda\)가 있지만 선형 독립인 고유벡터 \(v\)가 하나만 있는 경우, \( (A - \lambda I) w = v \)를 만족하는 일반화된 고유벡터 \(w\)가 추가로 필요합니다. 그러면 일반해는 \( \mathbf{x}(t) = c_1 v e^{\lambda t} + c_2 (tv + w) e^{\lambda t} \) 형식을 취합니다. 만약 고유 공간이 2차원이라면, 행렬은 해당 불변 부분 공간에서 단위 행렬의 스칼라 배수이며 해는 \( \mathbf{x}(t) = (c_1 v1 + c_2 v2) e^{\lambda t} \)로 단순화됩니다.

이 도구는 비선형 시스템을 기호적으로 풀 수 있습니까?

비선형 모드는 4차 룬게-쿠타(RK4) 적분기를 사용하여 시스템을 수치적으로 풀고 위상 초상화를 그립니다. 대부분의 비선형 시스템은 폐쇄형 해가 없으므로 이것이 표준적인 접근 방식입니다. 고정점에서의 자코비안을 계산하여 \(A\)로 입력하는 선형화를 통해 평형점 근처의 국부적 동작을 여전히 파악할 수 있으며, 이는 2×2 선형 모드에서 처리합니다.

위상 초상화란 무엇입니까?

위상 초상화는 \(x\)–\(y\) 평면에서 2D 시스템 해의 기하학적 그림입니다. 각 해는 궤적이라고 불리는 곡선을 그리며, 벡터장 화살표와 함께 궤적의 모음은 해가 안으로 소용돌이치는지, 안장점처럼 갈라지는지, 진동하는지 또는 평형점으로 수렴하는지와 같은 정성적 동작을 드러냅니다. 위상 초상화를 통해 시스템의 전역 구조를 한눈에 볼 수 있습니다.

더 읽어보기

• 미분 방정식 시스템 — Wikipedia
• 위상 초상화 — Wikipedia
• 고윳값과 고유벡터 — Wikipedia
• 룬게-쿠타 방법 — Wikipedia
• 반 데르 폴 진동자 — Wikipedia