[논문 리뷰] A Weiszfeld-like algorithm for a Weber location problem constrained to a closed and convex set

G. Torres|arXiv (Cornell University)|2012. 04. 04.

Computational Geometry and Mesh Generation참고 문헌 21인용 수 3

한 줄 요약

이 논문은 닫힘과 볼록 집합 내에서 제약 조건이 있는 웨이버 위치 문제를 위한 위즈펠드 유사 알고리즘을 제안하며, 반복값이 꼭짓점에 도달하더라도 수렴을 보장한다. 수직 투영과 수정된 위즈펠드 반복을 조합함으로써 타당성 있고 내림차순인 수열을 생성하며, KKT 조건을 만족하는 유일한 해로 수렴한다. 수치 실험에서 기존 비선형 소프트웨어보다 뛰어난 성능을 보였다.

ABSTRACT

The Weber problem consists of finding a point in $\mathbbm{R}^n$ that minimizes the weighted sum of distances from $m$ points in $\mathbbm{R}^n$ that are not collinear. An application that motivated this problem is the optimal location of facilities in the 2-dimensional case. A classical method to solve the Weber problem, proposed by Weiszfeld in 1937, is based on a fixed point iteration. In this work a Weber problem constrained to a closed and convex set is considered. A Weiszfeld-like algorithm, well defined even when an iterate is a vertex, is presented. The iteration function $Q$ that defines the proposed algorithm, is based mainly on an orthogonal projection over the feasible set, combined with the iteration function of a modified Weiszfeld algorithm presented by Vardi and Zhang in 2001. It can be seen that the proposed algorithm generates a sequence of feasible iterates that have descent properties. Under certain hypotheses, the limit of this sequence satisfies the KKT optimality conditions, is a fixed point of the iteration function that defines the algorithm, and is the solution of the constrained minimization problem. Numerical experiments confirmed the theoretical results.

연구 동기 및 목표

• 반복값이 꼭짓점에 도달할 경우 고전적 위즈펠드 알고리즘이 조기 수렴을 일으키는 문제를 해결하기 위해.
• 닫힘과 볼록 집합에 제약된 웨이버 문제를 위한 강건한 타당 반복 방법을 개발하기 위해.
• KKT 조건을 만족하는 최적 해로 내림차순 성질과 수렴을 보장하기 위해.
• 웨이버 함수의 미분 가능성 문제를 해결하고, 꼭짓점에서도 유효한 알고리즘을 제공하기 위해.

제안 방법

• 알고리즘은 내부 점에 대해 수정된 위즈펠드 반복 함수 T₀를 사용하며, 타당 집합 Ω에 수직 투영을 통해 조정된다.
• 반복값이 꼭짓점 aᵏ와 일치할 경우, T₀(aᵏ)와 aᵏ를 연결하는 직선을 따라 새로운 점을 정의하여 타당성을 확보한다.
• 반복 함수 Q는 투영과 꼭짓점 처리 논리를 통합하여 모든 반복값이 볼록 타당 집합 내에 있도록 보장한다.
• 알고리즘은 목적 함수 값 f(x⁽ˡ⁾)의 비증가 수열을 보장하며, x⁽ˡ⁾ ≠ Q(x⁽ˡ⁾)일 경우 엄격한 감소를 보인다.
• 알고리즘이 Q의 고정점로 수렴함을 증명하였으며, 이는 KKT 조건을 만족하는 유일한 해에 해당한다.
• 수치 검증은 MATLAB을 사용하여 m=50개의 꼭짓점과 임의의 가중치를 가진 1000개의 랜덤 인스턴스에서 fmincon과의 비교를 실시하였다.

실험 결과

연구 질문

• RQ1 제약 조건이 있는 웨이버 문제에서 반복값이 꼭짓점에 도달할 경우에도 강건한 위즈펠드 유사 알고리즘을 설계할 수 있는가?
• RQ2 닫힘과 볼록 집합에 제약된 조건에서 제안된 알고리즘이 내림차순 성질과 타당성을 유지하는가?
• RQ3 목적 함수가 꼭짓점에서 미분 가능하지 않을 경우에도 알고리즘이 진정한 해로 수렴할 수 있는가?
• RQ4 제약 조건이 있는 웨이버 문제에서 제안된 알고리즘의 성능은 fmincon과 같은 전통적 비선형 프로그래밍 소프트웨어보다 어떻게 비교되는가?
• RQ5 반복 함수 Q의 고정점은 주어진 조건 하에서 KKT 해와 동일한가?

주요 결과

• 제안된 알고리즘은 타당 반복값 수열을 생성하며, 내림차순 성질을 유지하여 최적 해로 수렴함을 보장한다.
• 반복값이 정확히 꼭짓점에 도달하더라도 알고리즘이 잘 정의되어 있으며, 고전적 위즈펠드 방법의 정체 문제를 피한다.
• 제약 조건이 있는 웨이버 문제의 유일한 해는 KKT 조건을 만족하며, 반복 함수 Q의 고정점이다.
• 1000개의 수치 실험에서 제안된 알고리즘이 fmincon보다 동일하거나 더 좋은 해를 모두 확보하였으며, 35개의 경우에서 차이가 0.01 초과였다.
• 알고리즘의 성능 우월성은 기울기 기반 방법에 의존하지 않고 꼭짓점에서의 비미분성을 효과적으로 다룰 수 있기 때문이다.