1. Preliminaries
Cutting plane methods for optimization have a long history that goes back at least to a fundamental paper by Kelley [
1]. The theoretical approach of the analytic center cutting plane methods started from Gon and Vial [
2]. du Merle [
3] developed an implementation of the method of the prototype, which was successfully applied to solve several nontrivial convex optimization problems [
4,
5]. Some later developments of the analytic center cutting plane methods have been proposed for solving various variational inequalities, e.g., [
6,
7,
8,
9].
We present two proximal analytic center cutting plane algorithms for solving variational inequalities whose domains are normal regions.
Our proximal analytic center cutting plane algorithms are also constructive solutions to our Nash economic equilibrium problems.
This study contains a detailed description of computational schemes of algorithms and provides the theoretical proofs of their convergence to the desired solution.
Suppose
X is a non-empty subset of the
n-dimensional Euclidean space
, and
F:
X →
is a function. We call that a point
x* ∈
X is a solution of the variational inequality
VI [
F,
X] if
The point
x* ∈
X is a solution of the dual variational inequality
VID (
F,
X) if
Given
VI[
F,
X] (
VID[
F,
X]), the gap function is defined as
A point is said to be a ε-solution of the variational inequality (1) if , for given .
A function
F:
X →
is pseudomonotone on
X if ∀
x ∈
XIt is known that (see Auslender [
10]), if
F is continuous, then a solution
x* ∈
X of
VID (
F,
X) is a solution of
VI (
F,
X); and if
F is continuous pseudomonotone, then
x* ∈
X is a solution of
VI (
F,
X) if and only if it is a solution of
VID (
F,
X).
It is known that the following Lemma 1 holds.
Lemma 1. A point x* ∈ X is a solution of VI[F, X] (VID[F, X]) if and only if
().
The convex hull of a set
B ⊆
is the set
A polytope is a set P ⊆ which is the convex hull of a finite set.
A polyhedron is a set
where
b ∈
, and
A is an
m ×
n matrix.
A polytope is a polyhedron. The following intuitively clear but nontrivial to prove result is essentially due to Farkas [
11], Minkowski [
12], and Weyl [
13]:
Lemma 2. P is a polytope if and only if it is a bounded polyhedron.
In the sequel, we assume that a polytope always has a non-empty interior.
Definition 1. A subset X ⊆
is said to be a normal region if it is a closed bounded set and if there exists a sequence of polytopes {Cj}, which satisfies such that The proof of the following Theorem 1 is trivial.
Theorem 1. A closed, bounded, convex region X ⊆ is a normal region.
We denote the topological closure of X.
Theorem 2. A subset X ⊆ is a normal region if there exists a uniformly bounded sequence of polytopes {Cj}, Cj ⊆ X (j = 1, ⋯), such that Proof. Then, the sequence satisfies the conditions in Definition 1. □
Definition 2. A subset X ⊆ is a unbounded normal region if there exists a sequence of bounded normal regions (j = 1, 2, …), such that each Xj contains all boundary points of X, and Actually, a subset X ⊆ is a bounded or unbounded normal region if and only if it is a bounded or unbounded convex region.
2. Proximal Analytic Center Cutting Plane Algorithms
This section modifies the method in Shen and Pang [
6] and presents proximal analytic center cutting plane algorithms for solving variational inequality
VI[
F,
X], whose domains are normal regions.
From now on, we make the following assumptions: for each x, y ∈ X, given any , , where , we can always find and such that
- (i)
,
- (ii)
if , no matter the relationship between and ,
- (iii)
, where L is a constant.
Assume the auxiliary
be a mapping that is continuous in
x and
y, strong monotone with respect to
y with a constant
M > 0, i.e.,
Consider the auxiliary variational inequality associated with
, whose solution
satisfies
In view of the strong monotonicity of
with respect to y, this auxiliary variational inequality has a unique solution (Goffin, Luo, and Ye [
7]).
For any polytope
,
is associated with the potential function
. An approximate analytic center, introduced by Goffin, Luo, and Ye [
7], is the maximizer of the potential function
and the unique solution of the system
where
z is a dual vector, and
Z the diagonal matrix built upon
z.
Let
X ⊆
be a normal region, there exists a sequence of polytopes {
Cj} satisfying
Cj ⊆
Cj+1 (
j = 1, 2, ⋯) such that
Then, there is a sequence of variational inequalities VI[F, Cj] (j = 1, 2, ⋯) from the original variational inequality VI[F, X], where the polytope Cj is given by the linear inequalities , x, bj ∈ , and Aj is an m × n matrix, (j = 1, 2, ⋯).
The following Algorithm 1 extends Algorithm 3.1 in Shen and Pang [
6] from a feasible region of polytope to a normal region.
Algorithm 1 |
Assume X is a bounded normal region. Let and be two constant. Set k = 0, j = 1, Ak = Aj, bk = bj, εk = ε, and . Step 1 (Computation of the approximate analytic center). Find an approximate analytic center of . Step 2 (Stopping criterion). If () = 0, then stop. Else go to Step 3. Step 3 (Computation of a solution). Find of If then increase j by 1 and go to Step 1. Else go to Step 4. Step 4 (Solving the approximate auxiliary variational inequality problem). Find , such that , where satisfies , . Step 5 (construction of the approximate cutting plane). Let and , where lk,j is the smallest integer that satisfies , where satisfies
where Let , , . Increase k by 1 and go to Step 1. End of Algorithm 1. |
Theorem 3. Algorithm 1 either stops with a solution of the variational inequality VI(F, X) after a finite number of iterations, or there exists a sequence in X that converges to the solution of VI(F, X).
Proof. For any given
j, ∃
xj ∈
Cj, such that
, (
j = 1, 2, …), then
Because X is a bounded set, there exists a subsequence of {xj} and a point x* ∈ X such that the subsequence converges to . We may assume that .
Since
X is closed,
such that
=
. And noting that
,
such that
, and for given
jSo,
for
There is a subsequence of
, and without a loss of generosity, we may assume
itself, which satisfies that
. Therefore, for
On the other hand, the continuity of
implies that
Which concludes that x* is a solution of VI(F, X). The proof is complete. □
By use of our Algorithm 1, we are going to present the following Algorithm 2 to solve variational inequality VI(F, X) over an unbounded normal region X.
We can find a sequence of bounded normal regions
Xj (
j = 1, 2, …), such that each
Xj contains all boundary points of
X and
We note that each variational inequality
VI [
F,
Xj] has a unique solution
(
j = 1, 2, …). And, when
j is large enough, say
j >
N, the unique solution
of the
VI[
F,
X] satisfies that
.
Algorithm 2 |
Assume X is an unbunded normal region. Set j = 1. Step 1 (Computation of a sequence such that )) By use of Algorithm 1, computing a sequence such that and .), (and let , one has , although is still unknown yet). Step 2 (Stopping criterion). If () = 0 for any of , then stop. Else go to Step 1. End of Algorithm 2. |
Theorem 4. Algorithm 2 either stops with a solution of the variational inequality VI(F, X) after a finite number of iterations, or there exists a sequence in X that converges to the solution of VI(F, X).
Proof. Algorithm 2 stops with a solution of
VI(
F,
X), or we can find sequences
such that
(
j = 1, 2, …), and
if
j >
N, i.e.,
(although
is still unknown) if
j >
N. Therefore, from Step 1, one obtains
Then, take the sequence
(consisting of the first term of each sequence
,
j = 1, 2, …), and one has
Hence, is a Cauchy sequence and so is convergent. From (2), one obtains , which completes the proof. □
3. Nash Economic Equilibrium Application
Economic equilibrium is a condition or state in which economic forces are balanced. These economic variables remain unchanged from their equilibrium values in the absence of external influences. Economic equilibrium may also be defined as the point at which supply equals demand for a product, with the equilibrium price existing where the hypothetical supply and demand curves intersect. Economists can usually explain the past and sometimes predict the future, but not without help. One of the most important tools at their disposal is the Nash equilibrium, named after John Nash [
14], who won a Nobel Prize in 1994 for his discovery. There were plenty of discussions followed after Nash, e.g., some recent approaches can be seen in Fischer [
15], Faraci [
16], and Boilan [
17].
In this section, we explain that our proximal analytic center cutting plane algorithms can be used to solve some practical problems whose domains are normal regions, e.g., Nash economic equilibrium problems over normal regions.
Consider an oligopolistic economy in which a homogeneous product is supplied by
n firms. Let
p(
σ) denote the inverse demand function, which is the price at which consumers will purchase a quantity
σ. If each firm supplies
qi units of the product, then the total supply is
We denote by
the
i-th firm’s total cost of supplying
qi units of the product; the profit of the
i-th firm is given by
A vector
is a said to be a Nash equilibrium solution if it is an optimal solution to the problem
where K is the box-constrained set (see Konnov and Volotskaya [
18]),
aj are constants,
bj are either constants or +∞ (
j = 1, …,
n). And
(
i = 1, 2, …,
n).
The Clarke–Rokafellar generalized derivatives of
h at
p0 in the direction
d is defined by
where
,
means that
and
,
indicates that both
and
.
The Clarke–Rokafellar subdifferential of
h at
p0 is given by
From Golshtein and Tretyakov [
19], as well as Murphy, Sherali, and Soyster [
20], a vector is a Nash equilibrium solution in the oligopolistic economy if and only if it is a solution to the problem
where K =
.
Write
, then
Similarly, we can introduce the Nash equilibrium solution in the oligopolistic economy if the domain is a normal region X:
A Nash equilibrium solution in the oligopolistic economy can be defined as a solution of the variational inequality
VI[
F,
X]
if
X is a normal region.
Example 1. The following mathematical model is used to calculate the remaining loan balance of a fixed mortgage loan. The mortgage payment amount should be paid periodically for m periods on a mortgage amount L at a periodic interest rate of .
After periods for full amortization m periods (), the remaining balance of the loan is given by Assume one would like to consider a selling price
that satisfies
where
D is the down payment at the beginning of the mortgage.
If there are n = (l1, l2, …, ln) mortgage providers in the same mortgage rate, p = r1 = p(σ). Mortgage provider li provides a mortgage amount of qi (i = 1, 2, …, n). Then, qi are all functions of r1, r2, and r3.
Let
Y =
Y (
r1,
r2,
r3) be the convex region given by (5), where 0 <
q1 ≤ M < ∞ for a given M. Then,
which is an example of Nash equilibrium over a convex region.
Therefore, by Example 1, a domain of a Nash equilibrium problem may be a normal region, and so our proximal analytic center cutting plane methods can be applied to solve it.
A function
F:
X →
is said to be strictly pseudomonotone on
X if
Theorem 5. Let X ⊆ be a normal region and F: X → a continuous and strictly pseudomonotone function, then the Nash equilibrium (4) has one and only one solution.
Proof. ,
such that
and
. Without loss of generality, we suppose that
, then
. Due to the convexity of
, we have
Which means that
is convex. It is easy to see that the closure of any convex set is convex. Therefore,
is a convex and compact set in
. From Hartman and Stampacchia [
21], the Nash equilibrium (4) has solutions.
On the other hand, assume that
is a solution of the Nash equilibrium (4); then,
Due to the strict pseudomonotonicity of
F, we have
i.e.,
Which indicates that with is not a solution of (2). Therefore, the Nash equilibrium (2) has at most one solution. We complete the proof. □