A Smoothing Neural Network Algorithm for Absolute Value Equations ()
1. Introduction
Consider the following absolute value problem [1] -[3] :
(1)
where, is absolute value of x, it is a subclass of absolute value equations which is proposed by Rohn [4] , and it is a NP-hard problem [1] .
The AVE has closed relation with some important problems, for example, the linear programming, Quadratic programming problem and the bimatrix game problem. The above problems can be transformed into the linear complementarity problem, and the linear complementarity problem can be transformed into the absolute value equations. Due to its simple and special structure and application value, the research on absolute value equation has drawn attention of many researchers. Mangasarian [5] pointed out the relationship between backpack feasibility problem and the AVE. The problem of AVE has been studied deeply by Yamashita and Fukushima [6] , and the results of the research on the problem of AVE are applied to the problem of location selection, good results are obtained. The numerical solution methods of AVE, such as Newton method, quasi-Newton method, are reachable in [7] -[12] .
In this paper, we present a smooth approximation function which is based on neural network method to solve the AVE. By using a smooth approximation function of, we turn it into a differentiable unconstrained optimization problem. Furthermore, we obtain the approximate solution of the original problem based on our established unconstrained optimization problem and the neural network model. Compared with the Newton method, the neural network model needs less requirement for the hardware of compute and the iterative process is real-time.
2. The Smoothing Reformulating of AVE
The absolute value Equation (1) is equivalent to the nonlinear equations:
(2)
where. Since it is a non smooth function, we construct a smooth function to approximate it.
Definition 1.1 Smoothing approximation function, given a function, smoothing function is called smoothing approximation function, if for any, there exists so that
where is not dependent on the x.
In this paper, we use the aggregate function [13] to give a smooth approximation of the absolute value equation:
Let, so every component of the absolute value function can be recorded as
For any, the definition of smoothing function is as follows
So the function of absolute value is obtained as follows
Thus the absolute value equation is transformed into the following smooth equations
(3)
where
We define the function as follows
where is the smoothing approximation of, is said as the energy function
of the neural network. Thus, the approximate solution of the absolute value equation is transformed to the global optimal solution of the optimization problem
3. Neural Network Model for Absolute Value Equation
Consider the following unconstrained optimization problem
(4)
the gradient can be calculated by the following formula:
where
now, we can give a neural network model for solving the absolute value equation, which is based on the steepest descent neural network model for (4).
(5)
where is a parameter represents that one can use a larger step size in the simulation, specific details can be referred to [14] -[16] . To simplify our analysis, we let throughout this paper. A block diagram (Figure 1) of the neural network is shown as follows.
4. Analysis of Stability and Existence
Next, we recall some materials about first order differential equations (ODE) [17] :
(6)
where is a mapping. We also introduce three kinds of stability that will be discussed later.
Definition 3.1 A point is called an equilibrium point or a steady state of the dynamic system (6)
if If the reisaneighborhood of such that a,
then is called an isolated equilibrium point.
Lemma 3.1 Assume that is a continuous mapping. Then, for any and, there exists a local solution for (6) with for some. If, in addition, H is locally Lipschitz continuous at, then the solution is unique.
Definition 3.2 (Asymptotic stability). An isolated equilibrium point is said to be asymptotically stable if
in addition to being Lyapunov stable, it has the property that as for all
Definition 3.3 (Lyapunov stability). Stability in the sense of Lyapunov Let be a solution for (6). An isolated equilibrium point is Lyapunov stable if for any and any there exists a
such that for all and
Definition 3.4 (Lyapunov function). Let be an open neighborhood of. A continuously differentiable function is said to be a Lyapunov function at the state over the set for Equation (6) if
Figure 1. The block diagram of neural network (5).
Lemma 3.2 a) An isolated equilibrium point is Lyapunov stable if there exists a Lyapunov function over some neighborhood of
b) An isolated equilibrium point is asymptotically stable if there is a Lyapunov function over some
neighborhood of such that for all
Lemma 3.3 [11] For any, , if then.
Theorem 3.1 is Lyapunov function over some neighborhood of
Proof. Let be the solution of the absolute value equation.
1) The function is obtained by our smooth approximation. So is continuous with respect to
. Obviously have continuous partial derivatives at all components of the.
2) Since, then
3) If then is always holds
So, by the Definition 3.4 we know that is Lyapunov function over some neighborhood of
Theorem 3.2 Each solution of the absolute value equation is the equilibrium point of the neural network (5).
Conversely, if, the equilibrium point of the neural network (5) is the solution of the absolute value equation.
Proof. Assume that is the solution of the absolute value equation, since, for any
we have if and only if is the solution of the absolute value equation.
Obviously, we got, so is the equilibrium point of the neural network (5). On the other
hand if and, then we get. So, the equilibrium point of the neural network (5) is the solution of the absolute value equation.
Next, we can prove that is not only Lyapunov stable and asymptotically stable.
Theorem 3.3. Let the be the isolated equilibrium of the neural network. is the Lyapunov stability and asymptotic stability for neural networks.
Proof. Since is the isolated equilibrium of the neural network, the solution of the absolute value eq-
uation is known by the Theorem 3.2. Therefore,. In addition, Since x is the isolated equilibrium
point, so and are hold over the neighborhood of. By
Theorem 3.1 we know that is Lyapunov function over some neighborhood of, so by Lemma 3.2 the isolated equilibrium is Lyapunov stable. Because is isolated, it is not difficult to compute:
Consequently, we have. By Lemma 3.2, is asymptotic stability.
5. Numerical Experiment
In this section we give some smooth of numerical tests of neural network algorithm, due to the complementarity problem can be transformed to absolute value equations, we consider the linear complementarity problem which is equivalent to the absolute value equations as test cases.
For a given matrix and vector, The linear complementarity problem is to find a vector to satisfy.
From the Theorem 2 in the literature [11] , if 1 is not the eigenvalues of the matrix, then is equivalent to the following absolute value equation:
where and is the solution of the absolute value equation.
Example 1 [11] . Consider the following linear complementary problem:
Since1 is not included in the eigenvalues of, then the linear complementary problem can be transformed into the following absolute value equation and they are equivalent:
where
We can find that is a solution of the absolute equation.
By using the neural network model, the initial point is generated by x0 = rand (n,1), and the program is performed under the environment of MATLAB7.11.0. The following two figures (Figure 2 and Figure 3) describe how the approximate solution of example 1 and the energy function varies with time.
Figure 2. Transient behavior of x(t) of example 1.
Figure 3. Transient behavior of energy function of example 1.
Since, then we can capture the solution of the linear complementary problem, the solution is.
Example 2 [11] . Consider the following linear complementary problem:
Through calculation, we can get one eigenvalue of is 1. By literature [11] , we can find that if 1 is the eigenvalue of matrix, then the and of the linear complementary problem need to be multiplied by a positive constant and makes 1 not the eigenvalue of (and the solution of the linear complementary pro- blem keeps invariant). Then we can transform linear complementary problem into absolute value equation by applied Theorem 2 and Theorem 3 in literature [11] .
Set, then we can find that 1 is not included in the eigenvalues of. And and have the common optimal solution, while z can be transformed into the absolute value equation by applying the Theorem 2. Then, we have
And the absolute value equation is, where:
Thus, we can get one solution of the absolute value equation whcih is, then the following two figures (Figure 4 and Figure 5) describe how the approximate solution of example 2 and the energy function varies with time
Since, then we can capture the solution of the linear complementary problem LCP(M,q), the solution is.
Example 3. Consider the following linear complementary problem:
Figure 4. Transient behavior of x(t) of example 2.
Figure 5. Transient behavior of energy function of example 2.
Through calculation we can get one eigenvalue of is 1. And the same as example 2, set, then, we can find that:
And the absolute value equation is:, where:
Thus, we can get the solution of the absolute value equation which is, then the following two figures (Figure 6 and Figure 7) describe how the approximate solution of example 3 and the energy function varies with time
Since, then we can capture the solution of the linear complementary problem LCP(M,q), the solution is z* = (2.5 2.5 0 2.5).
6. Conclusion
This paper adopted the aggregate function method to tackle the absolute value equation with smooth processing, and then turned the absolute value equation into a differentiable unconstrained optimization problem. In order to obtain the approximate solution of the original problem we use the proposed neural network model to solve the
Figure 6. Transient behavior of x(t) of example 3.
Figure 7. Transient behavior of energy function of example 3.
unconstrained optimization problem. At the same time, we propose one neural network which is based on different energy function. Through the transformation between linear complementary problem and absolute value equation, it can be used to solve the linear complementary problem, too. For the traditional energy function based on the NCP function, we can avoid a lot of matrix calculation. Numerical examples show that the algorithm is very effective for solving this kind of absolute value equation, and the accuracy of solution can be controlled by the parameters completely. In view of the fact that it is relatively difficult to solve the absolute value equation, the proposed method in this paper can be used to solve the absolute value problem effectively.
Acknowledgements
This work is supported by National Natural Science Foundation of China (No.11171221) and Innovation Program of Shanghai Municipal Education Commission (No.14YZ094).