[논문 리뷰] Optimal linear stability condition for scalar differential equations with distributed delay

Samuel Bernard, Fabien Crauste|arXiv (Cornell University)|2014. 05. 28.

Stability and Controllability of Differential Equations참고 문헌 50인용 수 10

한 줄 요약

이 논문은 분포 지연을 가진 스칼라 선형 미분방정식에 대해, 해당 시스템이 단일 이산 지연을 가진 경우 안정할 경우, 동일한 평균을 가진 임의의 분포 지연에 대해서도 점점 안정성이 보장됨을 입증한다—이는 분포 형태에 영향을 받지 않는 강력한 안정성 기준을 제공한다. 핵심 기여는 이산 지연이 가장 불안정성을 유도하는 경우임을 보여주며, 따라서 이산 지연 하에서의 안정성이 동일한 평균을 가진 모든 분포 지연에 대해 안정성을 보장함을 의미한다.

ABSTRACT

Linear scalar differential equations with distributed delays appear in the study of the local stability of nonlinear differential equations with feedback, which are common in biology and physics. Negative feedback loops tend to promote oscillations around steady states, and their stability depends on the particular shape of the delay distribution. Since in applications the mean delay is often the only reliable information available about the distribution, it is desirable to find conditions for stability that are independent from the shape of the distribution. We show here that for a given mean delay, the linear equation with distributed delay is asymptotically stable if the associated differential equation with a discrete delay is asymptotically stable. We illustrate this criterion on a compartment model of hematopoietic cell dynamics to obtain sufficient conditions for stability.

연구 동기 및 목표

• 지연 분포의 형태에 관계없이 분포 지연을 가진 스칼라 선형 미분방정식의 안정성 조건을 규명하는 것.
• 생물학적 및 물리적 모델에서 일반적으로 지연의 평균만 알려져 있을 뿐 전체 분포가 알려져 있지 않은 점에 대응하는 것.
• 이러한 시스템에서 이산 지연이 가장 불안정성을 유도하는지 여부에 대한 열린 문제를 해결하는 것.
• 혈액 세포 생성 동역학과 같은 모델에 적용 가능한 실용적이고 계산 가능성이 높은 안정성 기준을 제공하는 것.
• 이산 지연 하에서의 안정성이 동일한 평균을 가진 임의의 분포 지연에 대해 안정성을 보장함을 보여주는 것.

제안 방법

• 분포 지연을 가진 스칼라 함수 미분방정식의 선형화된 형태를 분석: $ \dot{x} = -a x - b \int_0^\infty x(t - \tau) \, d\eta(\tau) $, 여기서 $ \eta $는 누적분포함수이다.
• 분포 지연 시스템의 안정성과 동일한 평균 $ \mu $를 가진 단일 이산 지연 $ \tau = \mu $을 가진 시스템의 안정성을 비교한다.
• 선형 지연된 함수 미분방정식 이론과 고유값 분석을 적용하여 점점 안정성 조건을 도출한다.
• 선형화된 시스템에서 유도된 특성방정식을 사용하여 이산 지연과 분포 지연의 경우에 고유값의 위치를 비교한다.
• 특정 경우(예: 지수분포 및 감마분포)에서 뉴퀴스트 기준과 소형 이득 정리를 적용하여 일반 결과를 검증한다.
• 혈액 줄기세포 동역학의 분할 모델을 사용하여 이론적 결과를 검증한다. 이 모델은 분포 지연과 이산 지연 모두를 포함한다.

실험 결과

연구 질문

• RQ1 분포 지연을 가진 선형 스칼라 미분방정식의 안정성은 지연 분포의 구체적 형태에 의존하는가?
• RQ2 전체 분포의 지식 없이도 평균 지연에만 의존하는 안정성 조건을 유도할 수 있는가?
• RQ3 이러한 시스템에서 이산 지연이 가장 불안정성을 유도하는 분포인가?
• RQ4 이산 지연 하에서의 안정성이 동일한 평균을 가진 모든 분포 지연에 대해 안정성을 보장하는가?
• RQ5 이 결과는 피드백 루프를 포함한 혈액 세포 생성 동역학과 같은 생물학적 모델에 어떻게 적용되는가?

주요 결과

• 주어진 평균 지연 $ \mu $에 대해, 단일 이산 지연 $ \tau = \mu $을 가진 선형 방정식이 점점 안정할 경우, 동일한 평균을 가진 임의의 분포 지연을 가진 방정식 역시 점점 안정하다.
• 이산 지연 분포는 가장 불안정성을 유도한다: 시스템이 이산 지연 하에서 불안정할 경우, 동일한 평균을 가진 모든 분포 지연에 대해서도 불안정하다.
• 중립성 또는 음성 협동성의 경우, 혈액 세포 모델에서 정상 상태는 조건 $ E_p := \sum_{i=1}^3 p_i \frac{q_i}{\beta_i} < \frac{k_0 \pi}{2 \alpha h(k_0 - \alpha)} $ 가 성립할 경우 선형적으로 점점 안정하다 (중립성의 경우 $ r=1 $).
• 음성 협동성의 경우 ($ 0 \leq r < 1 $), 정상 상태는 $ \bar{x}^{h+1-r} \leq \frac{k_0(1-r)}{\alpha h} $ 이거나, 이 부등식이 성립하지 않지만 조건 (31) 이 성립할 경우 안정하다.
• 수치 시뮬레이션 결과, 분포 지연 시스템은 이산 지연 시스템이 안정한 점들에서도 안정성을 유지하지만, 그 반대는 성립하지 않음—이는 이산 지연이 더 불안정성을 유도함을 확인한다.
• 결과적으로 실험적으로 관측된 주기적 진동을 재현하는 것은 분포 지연을 사용할 경우 이산 지연을 사용할 경우보다 더 어렵다는 것을 의미한다. 진동은 정상 상태의 불안정화를 필요로 하기 때문이다.