미분방정식에서 구조적 인과 모델로 결정적 경우

초록

본 논문은 1차 결정론적 미분방정식(ODE)의 평형 상태를 구조적 인과 모델(SCM)로 변환하는 방법을 제시한다. 안정성 가정 하에 평형 방정식을 라벨링하여 완전한 ‘do’ 개입을 정확히 기술하고, 이를 통해 순환 인과 구조도 일관된 인과 의미를 갖는다고 설명한다.

상세 분석

이 논문은 인과 그래프 이론이 주로 비순환(acyclic) 시스템에 적용되어 왔던 전통을 넘어, 순환(cyclic) 구조를 포함한 동적 시스템을 인과 모델로 해석하는 새로운 프레임워크를 제시한다. 핵심 아이디어는 1차 ODE 시스템이 충분히 매끄러운 함수 f_i 로 정의되고, 모든 변수 X_i 가 시간에 따라 변화하는 미분 방정식 \dot X_i = f_i(X_{pa(i)}) 로 표현될 때, 시스템이 안정(stable)하면 t→∞ 에서 고유한 평형점 X* 로 수렴한다는 점이다. 평형 조건 0 = f_i(X_{pa(i)}) 은 변수들 간의 정적 관계를 나타내며, 이를 라벨링된 평형 방정식 집합 E = {(i, E_i)} 로 정리한다. 라벨링은 각 방정식이 어떤 변수에 대응하는지를 명시함으로써, ‘do(X_I = ξ_I)’와 같은 완전 개입이 발생했을 때 어떤 방정식이 교체되어야 하는지를 정확히 지정한다. 개입 후의 평형 방정식은 교체된 변수에 대해 0 = X_i – ξ_i 형태가 되고, 나머지 변수는 원래 방정식을 유지한다. 따라서, 개입 전후의 평형 해는 라벨링된 방정식 집합만으로도 완전히 기술될 수 있다.

논문은 두 가지 구체적 예시를 통해 이론을 검증한다. 첫 번째는 포식자‑피식자 관계를 모델링하는 Lotka‑Volterra 방정식이다. 이 시스템은 비안정적이며 두 개의 평형점(전부 0, 그리고 양의 고정점)을 가진다. ‘do(X_2 = ξ_2)’ 개입을 가하면 시스템은 일반적으로 안정적인 평형(0, ξ_2) 로 수렴한다. 두 번째 예시는 마찰이 포함된 다중 질량‑스프링 시스템이다. 각 질량‑위치‑운동량 쌍을 변수로 두고, 마찰에 의해 모든 초기 조건이 하나의 고유 평형으로 수렴한다. 특정 질량을 고정하는 개입(do(Q_i = ξ_i, P_i = 0))을 적용해도 마찬가지로 고유 평형이 존재한다. 이러한 예시들은 동적 시스템이 서로 다른 미분식이라도 동일한 라벨링된 평형 방정식, 즉 동일한 SCM을 생성할 수 있음을 보여준다.

핵심 통찰은 다음과 같다. (1) 평형 방정식 자체는 인과 구조를 완전히 드러내지 못하지만, 라벨링을 통해 ‘어떤 변수에 대한 방정식이 바뀌는가’를 명시하면 인과적 개입을 정확히 모델링할 수 있다. (2) 순환 인과 그래프는 동적 시스템의 평형에서 자연스럽게 등장하며, 이때 인과 의미는 시간 흐름이 사라진 정적 관계로 보존된다. (3) 안정성 가정은 필수적이며, 이는 모든 가능한 개입에 대해 고유 평형이 존재하고 수렴함을 보장한다. (4) 이 접근법은 기존의 ‘시간‑시점 고정점’ 해석과 차별화되며, 연속시간 동역학을 직접 모델링하지 않고도 인과적 효과를 예측한다. 따라서, 동적 시스템을 연구하는 물리·생물·공학 분야에서 인과 추론을 수행할 수 있는 강력한 이론적 기반을 제공한다.