Euler 방법의 Lyapunov 안정성 연구

초록

본 논문은 미분 포함 시스템의 Euler 이산화에 대해 Lyapunov 안정성 기준을 확장한다. 연속시간 Lyapunov 함수와 이산시간 Lyapunov 함수를 쌍으로 사용해 비고립적인 지역 최소점의 d‑안정성을 충분조건으로 제시하고, Bouligand 서브그라디언트 방법에 적용한다.

상세 분석

논문은 먼저 전통적인 Lyapunov 안정성 개념을 미분 포함 (F:\mathbb R^{n}\rightrightarrows\mathbb R^{n}) 에 정의하고, 이를 Euler 이산화 (x_{k+1}\in x_{k}+ \alpha_{k}F(x_{k})) 에 그대로 옮기려는 시도를 비판한다. 기존의 d‑stability 정의는 단계 크기 ({\alpha_{k}}) 가 임의로 작아도 모든 (\varepsilon) 에 대해 초기값을 충분히 작은 (\delta) 로 잡을 수 있음을 요구하지만, 실제로는 (F) 가 비연속이거나 비단조적인 경우 이 조건을 만족하지 못한다. 이를 해결하기 위해 저자는 두 개의 Lyapunov 함수 (f) (연속시간)와 (g) (이산시간)를 도입한다. (f) 는 보존적 필드 (D) 에 대한 potential이며, (\sup_{u\in F(x)}\min_{v\in D(x)}\langle u,v\rangle\le 0) 를 만족한다. (g) 는 “d‑Lyapunov” 함수로, 모든 허용된 단계 크기 (\alpha\le\bar\alpha) 에 대해 (g(x_{k+1})\le g(x_{k})) 가 성립한다.

주요 정리 Theorem 1.2는 다음 네 가지 가정을 전제로 한다. (i) (F) 는 상반연속이며 비어 있지 않은 콤팩트 값을 가진다. (ii) (f) 는 보존적 필드 (D) 의 potential이며, 지역적 강한 감소 조건 (\max_{u\in\operatorname{co}F(x)}\min_{v\in D(x)}\langle u,v\rangle\le -\min_{w\in D(x)}\kappa(w)) 를 만족하고, (f) 의 최소값은 격리된 (D)‑임계값이다. (iii) (g) 는 연속이고 (x) 근처에서 d‑Lyapunov이다. (iv) (x) 가 (f+g) 의 엄격한 지역 최소점이면 (x) 는 d‑stable한다.

이 정리는 기존의 반정밀(semialgebraic) 가정이나 일정한 단계 크기 제한을 완전히 없앤다. 특히, 변수 단계 크기 (\alpha_{k}=\beta(k+1)^{-\frac1\gamma}) ( (\gamma\ge p) )를 허용함으로써 실제 최적화 알고리즘에서 흔히 쓰이는 감소 스케줄을 이론적으로 정당화한다. 저자는 (f) 와 (g) 가 서로 보완적인 역할을 한다고 강조한다. (f) 는 값이 유계하게 유지하면서 (g) 는 단조 감소한다. 따라서 (f+g) 가 최소값을 갖는 영역 안에 모든 궤적이 머무르게 되고, 이는 d‑stability를 보장한다.

또한, 논문은 (g) 의 감소 속도를 (-\omega\alpha^{p}) 형태로 강제하면 (\sum\alpha^{p}=+\infty) 조건 하에 (f(x_{k})\to f(x)) 와 (g(x_{k})\to g(x)) 가 동시에 성립함을 보인다. 이는 전통적인 Lyapunov 이론에서 요구되는 단조성(monotonicity) 없이도 수렴을 확보할 수 있음을 의미한다.

기술적 측면에서 저자는 (F) 가 비연속이거나 비보존적일 때도 적용 가능한 충분조건을 제시한다. 예를 들어 (F(x)=-\operatorname{sign}(x)|x|^{p}) 와 같이 (F) 가 (0) 에서 점프하는 경우에도, 적절히 작은 (\alpha_{k}) 와 (g) 의 설계만으로 d‑stability를 확보한다는 점을 시연한다.

마지막으로, Bouligand 서브그라디언트 (\partial_{B}f) 를 (F) 로 잡은 경우, (f) 가 비광택(non‑smooth)이라도 위의 조건을 만족하면 해당 비정규 최적화 문제의 Euler 이산화가 지역 최소점 근처에서 안정적으로 수렴한다는 실용적인 결론을 도출한다.